My Project
 All Classes Files Functions Variables Macros
covnat_p.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE covnat_p(klevel,ucov, vcov, unat, vnat )
5  USE parallel
6  IMPLICIT NONE
7 
8 c=======================================================================
9 c
10 c Auteur: F Hourdin Phu LeVan
11 c -------
12 c
13 c Objet:
14 c ------
15 c
16 c *********************************************************************
17 c calcul des compos. naturelles a partir des comp.covariantes
18 c ********************************************************************
19 c
20 c=======================================================================
21 
22 #include "dimensions.h"
23 #include "paramet.h"
24 #include "comgeom.h"
25 
26  INTEGER klevel
27  REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel )
28  REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
29  INTEGER l,ij
30  INTEGER :: ijb,ije
31 
32 
33  ijb=ij_begin
34  ije=ij_end
35 
36  if (pole_nord) then
37  DO l = 1,klevel
38  DO ij = 1, iip1
39  unat(ij,l) =0.
40  END DO
41  ENDDO
42  endif
43 
44  if (pole_sud) then
45  DO l = 1,klevel
46  DO ij = ip1jm+1, ip1jmp1
47  unat(ij,l) =0.
48  END DO
49  ENDDO
50  endif
51 
52  ijb=ij_begin
53  ije=ij_end
54  if (pole_nord) ijb=ij_begin+iip1
55  if (pole_sud) ije=ij_end-iip1
56 
57  DO l = 1,klevel
58  DO ij = ijb, ije
59  unat( ij,l ) = ucov( ij,l ) / cu(ij)
60  ENDDO
61  END DO
62 
63  ijb=ij_begin-iip1
64  ije=ij_end
65  if (pole_nord) ijb=ij_begin
66  if (pole_sud) ije=ij_end-iip1
67 
68  DO l = 1,klevel
69  DO ij = ijb,ije
70  vnat( ij,l ) = vcov( ij,l ) / cv(ij)
71  ENDDO
72 
73  ENDDO
74 
75  RETURN
76  END