libf/dyn3d/covnat.f90

module covnat_m

  IMPLICIT NONE

contains

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

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

    USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
    USE comgeom, ONLY: cu, cv

    ! Authors: F. Hourdin, Phu Le Van

    ! Objet : calcul des composantes naturelles à partir des
    ! composantes covariantes.

    INTEGER, intent(in):: klevel
    REAL, intent(in):: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
    REAL, intent(out):: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)

    ! Local:
    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

end module covnat_m
1	guez	61	module covnat_m
2	guez	3
3	guez	61	IMPLICIT NONE
4	guez	3
5	guez	61	contains
6	guez	3
7	guez	61	SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
8	guez	3
9	guez	61	! From LMDZ4/libf/dyn3d/covnat.F, version 1.1.1.1 2004/05/19 12:53:07
10	guez	3
11	guez	61	USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
12			USE comgeom, ONLY: cu, cv
13	guez	3
14	guez	61	! Authors: F. Hourdin, Phu Le Van
15	guez	3
16	guez	61	! Objet : calcul des composantes naturelles à partir des
17			! composantes covariantes.
18	guez	3
19	guez	61	INTEGER, intent(in):: klevel
20			REAL, intent(in):: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
21			REAL, intent(out):: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)
22	guez	45
23	guez	61	! Local:
24			INTEGER l, ij
25	guez	45
26	guez	61	!------------------------------------------------------------------
27	guez	45
28	guez	61	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
37			DO ij = ip1jm+1, ip1jmp1
38			unat(ij, l) =0.
39			END DO
40
41			DO ij = 1, ip1jm
42			vnat(ij, l) = vcov(ij, l) / cv(ij)
43			ENDDO
44			ENDDO
45
46			END SUBROUTINE covnat
47
48			end module covnat_m