/[lmdze]/trunk/Sources/phylmd/clvent.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/clvent.f

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

revision 202 by guez, Wed Jun 8 12:23:41 2016 UTC revision 206 by guez, Tue Aug 30 12:52:46 2016 UTC
# Line 35  contains Line 35  contains
35      ! delp-----input-R- epaisseur de couche (Pa)      ! delp-----input-R- epaisseur de couche (Pa)
36      REAL d_ven(klon, klev)      REAL d_ven(klon, klev)
37      ! d_ven----output-R- le changement de "ven"      ! d_ven----output-R- le changement de "ven"
38      REAL flux_v(klon, klev)  
39      ! flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)      REAL, intent(out):: flux_v(:) ! (knon)
40        ! (diagnostic) flux du vent à la surface, en (kg m/s)/(m**2 s)
41        ! flux_v est le flux de moment angulaire (positif vers bas)
42    
43      ! Local:      ! Local:
44      INTEGER i, k      INTEGER i, k
# Line 98  contains Line 100  contains
100         ENDDO         ENDDO
101      ENDDO      ENDDO
102    
     ! flux_v est le flux de moment angulaire (positif vers bas) dont  
     ! l'unite est: (kg m/s)/(m**2 s)  
103      DO i = 1, knon      DO i = 1, knon
104         flux_v(i, 1) = zx_coef(i, 1)/(RG*dtime) &         flux_v(i) = zx_coef(i, 1)/(RG*dtime) &
105              *(local_ven(i, 1)*zx_alf1(i) &              *(local_ven(i, 1)*zx_alf1(i) &
106              +local_ven(i, 2)*zx_alf2(i))              +local_ven(i, 2)*zx_alf2(i))
107      ENDDO      ENDDO
     DO k = 2, klev  
        DO i = 1, knon  
           flux_v(i, k) = zx_coef(i, k)/(RG*dtime) &  
                * (local_ven(i, k)-local_ven(i, k-1))  
        ENDDO  
     ENDDO  
108    
109      DO k = 1, klev      DO k = 1, klev
110         DO i = 1, knon         DO i = 1, knon
Legend:
Removed from v.202  
changed lines
  Added in v.206
  ViewVC Help
Powered by ViewVC 1.1.21