trunk/dyn3d/covnat.f

SUBROUTINE covnat(klevel,ucov, vcov, unat, vnat )

  ! From LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07

  use dimens_m
  use paramet_m
  use comgeom

  IMPLICIT NONE

  !   Auteur:  F Hourdin Phu LeVan
  !   Objet:
  !    calcul des compos. naturelles a partir des comp.covariantes

  INTEGER klevel
  REAL ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
  REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
  INTEGER   l,ij

  !------------------------------------------------------------------

  DO l = 1,klevel
     DO ij = 1, iip1
        unat (ij,l) =0.
     END DO

     DO ij = iip2, ip1jm
        unat( ij,l ) = ucov( ij,l ) / cu(ij)
     ENDDO

     DO ij = ip1jm+1, ip1jmp1  
        unat (ij,l) =0.
     END DO

     DO ij = 1,ip1jm
        vnat( ij,l ) = vcov( ij,l ) / cv(ij)
     ENDDO
  ENDDO

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