/[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/libf/dyn3d/covnat.f90 revision 45 by guez, Wed Apr 27 13:00:12 2011 UTC
# Line 1  Line 1 
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  
Legend:
Removed from v.3  
changed lines
  Added in v.45
  ViewVC Help
Powered by ViewVC 1.1.21