/[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

trunk/phylmd/clvent.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/phylmd/clvent.f revision 206 by guez, Tue Aug 30 12:52:46 2016 UTC
# Line 15  contains Line 15  contains
15      USE suphec_m, ONLY: rd, rg      USE suphec_m, ONLY: rd, rg
16    
17      INTEGER knon      INTEGER knon
18      REAL, intent(in):: dtime      REAL, intent(in):: dtime ! intervalle de temps (en s)
     ! dtime----input-R- intervalle du temps (en second)  
19    
20      REAL u1lay(klon), v1lay(klon)      REAL u1lay(klon), v1lay(klon)
21      ! u1lay----input-R- vent u de la premiere couche (m/s)      ! u1lay----input-R- vent u de la premiere couche (m/s)
# Line 36  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 99  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.82  
changed lines
  Added in v.206
  ViewVC Help
Powered by ViewVC 1.1.21