/[lmdze]/trunk/dyn3d/covnat.f
ViewVC logotype

Diff of /trunk/dyn3d/covnat.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/covnat.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/dyn3d/covnat.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC
# Line 1  Line 1 
1  !  module covnat_m
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $  
3  !    IMPLICIT NONE
4        SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )  
5        use dimens_m  contains
6        use paramet_m  
7        use comgeom    SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
8        IMPLICIT NONE  
9        ! From LMDZ4/libf/dyn3d/covnat.F, version 1.1.1.1 2004/05/19 12:53:07
10  c=======================================================================  
11  c      USE comgeom, ONLY: cu, cv
12  c   Auteur:  F Hourdin Phu LeVan      USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
13  c   -------  
14  c      ! Authors: F. Hourdin, Phu Le Van
15  c   Objet:  
16  c   ------      ! Objet : calcul des composantes naturelles à partir des
17  c      ! composantes covariantes.
18  c  *********************************************************************  
19  c    calcul des compos. naturelles a partir des comp.covariantes      INTEGER, intent(in):: klevel
20  c  ********************************************************************      REAL, intent(in):: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
21  c      REAL, intent(out):: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)
22  c=======================================================================  
23        ! Local:
24        INTEGER l, ij
25        INTEGER klevel  
26        REAL ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )      !------------------------------------------------------------------
27        REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )  
28        INTEGER   l,ij      DO l = 1, klevel
29           DO ij = 1, iip1
30              unat(ij, l) =0.
31        DO l = 1,klevel         END DO
32           DO ij = 1, iip1  
33              unat (ij,l) =0.         DO ij = iip2, ip1jm
34           END DO            unat(ij, l) = ucov(ij, l) / cu(ij)
35           ENDDO
36           DO ij = iip2, ip1jm  
37              unat( ij,l ) = ucov( ij,l ) / cu(ij)         DO ij = ip1jm+1, ip1jmp1
38           ENDDO            unat(ij, l) =0.
39           DO ij = ip1jm+1, ip1jmp1           END DO
40              unat (ij,l) =0.  
41           END DO         DO ij = 1, ip1jm
42              vnat(ij, l) = vcov(ij, l) / cv(ij)
43           DO ij = 1,ip1jm         ENDDO
44              vnat( ij,l ) = vcov( ij,l ) / cv(ij)      ENDDO
45           ENDDO  
46      END SUBROUTINE covnat
47        ENDDO  
48        RETURN  end module covnat_m
       END  
Legend:
Removed from v.3  
changed lines
  Added in v.76
  ViewVC Help
Powered by ViewVC 1.1.21