source:
trunk/mars/libf/dyn3d/covnat.F
@
80
Last change on this file since 80 was 38, checked in by , 14 years ago | |
---|---|
File size: 1.1 KB |
Rev | Line | |
---|---|---|
[38] | 1 | SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat ) |
2 | IMPLICIT NONE | |
3 | ||
4 | c======================================================================= | |
5 | c | |
6 | c Auteur: F Hourdin Phu LeVan | |
7 | c ------- | |
8 | c | |
9 | c Objet: | |
10 | c ------ | |
11 | c | |
12 | c ********************************************************************* | |
13 | c calcul des compos. naturelles a partir des comp.covariantes | |
14 | c ******************************************************************** | |
15 | c | |
16 | c======================================================================= | |
17 | ||
18 | #include "dimensions.h" | |
19 | #include "paramet.h" | |
20 | #include "comgeom.h" | |
21 | ||
22 | INTEGER klevel | |
23 | REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel ) | |
24 | REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel ) | |
25 | INTEGER l,ij | |
26 | ||
27 | ||
28 | DO l = 1,klevel | |
29 | DO ij = 1, iip1 | |
30 | unat (ij,l) =0. | |
31 | END DO | |
32 | ||
33 | DO ij = iip2, ip1jm | |
34 | unat( ij,l ) = ucov( ij,l ) / cu(ij) | |
35 | ENDDO | |
36 | DO ij = ip1jm+1, ip1jmp1 | |
37 | unat (ij,l) =0. | |
38 | END DO | |
39 | ||
40 | DO ij = 1,ip1jm | |
41 | vnat( ij,l ) = vcov( ij,l ) / cv(ij) | |
42 | ENDDO | |
43 | ||
44 | ENDDO | |
45 | RETURN | |
46 | END |
Note: See TracBrowser
for help on using the repository browser.