1 |
! |
SUBROUTINE covnat(klevel,ucov, vcov, unat, vnat ) |
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $ |
|
3 |
! |
! From LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07 |
4 |
SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat ) |
|
5 |
use dimens_m |
use dimens_m |
6 |
use paramet_m |
use paramet_m |
7 |
use comgeom |
use comgeom |
8 |
IMPLICIT NONE |
|
9 |
|
IMPLICIT NONE |
10 |
c======================================================================= |
|
11 |
c |
! Auteur: F Hourdin Phu LeVan |
12 |
c Auteur: F Hourdin Phu LeVan |
! Objet: |
13 |
c ------- |
! calcul des compos. naturelles a partir des comp.covariantes |
14 |
c |
|
15 |
c Objet: |
INTEGER klevel |
16 |
c ------ |
REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel ) |
17 |
c |
REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel ) |
18 |
c ********************************************************************* |
INTEGER l,ij |
19 |
c calcul des compos. naturelles a partir des comp.covariantes |
|
20 |
c ******************************************************************** |
!------------------------------------------------------------------ |
21 |
c |
|
22 |
c======================================================================= |
DO l = 1,klevel |
23 |
|
DO ij = 1, iip1 |
24 |
|
unat (ij,l) =0. |
25 |
INTEGER klevel |
END DO |
26 |
REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel ) |
|
27 |
REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel ) |
DO ij = iip2, ip1jm |
28 |
INTEGER l,ij |
unat( ij,l ) = ucov( ij,l ) / cu(ij) |
29 |
|
ENDDO |
30 |
|
|
31 |
DO l = 1,klevel |
DO ij = ip1jm+1, ip1jmp1 |
32 |
DO ij = 1, iip1 |
unat (ij,l) =0. |
33 |
unat (ij,l) =0. |
END DO |
34 |
END DO |
|
35 |
|
DO ij = 1,ip1jm |
36 |
DO ij = iip2, ip1jm |
vnat( ij,l ) = vcov( ij,l ) / cv(ij) |
37 |
unat( ij,l ) = ucov( ij,l ) / cu(ij) |
ENDDO |
38 |
ENDDO |
ENDDO |
39 |
DO ij = ip1jm+1, ip1jmp1 |
|
40 |
unat (ij,l) =0. |
END SUBROUTINE covnat |
|
END DO |
|
|
|
|
|
DO ij = 1,ip1jm |
|
|
vnat( ij,l ) = vcov( ij,l ) / cv(ij) |
|
|
ENDDO |
|
|
|
|
|
ENDDO |
|
|
RETURN |
|
|
END |
|