/[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 225 by guez, Mon Oct 16 12:35:41 2017 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, intent(in):: u1lay(:), v1lay(:) ! (knon)
21      ! u1lay----input-R- vent u de la premiere couche (m/s)      ! vent de la premiere couche (m/s)
     ! v1lay----input-R- vent v de la premiere couche (m/s)  
22    
23      REAL, intent(in):: coef(:, :) ! (knon, klev)      REAL, intent(in):: coef(:, :) ! (knon, klev)
24      ! Coefficient d'echange (m**2/s) multiplié par le cisaillement du      ! Coefficient d'echange (m**2/s) multiplié par le cisaillement du
# Line 36  contains Line 34  contains
34      ! delp-----input-R- epaisseur de couche (Pa)      ! delp-----input-R- epaisseur de couche (Pa)
35      REAL d_ven(klon, klev)      REAL d_ven(klon, klev)
36      ! d_ven----output-R- le changement de "ven"      ! d_ven----output-R- le changement de "ven"
37      REAL flux_v(klon, klev)  
38      ! flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s)      REAL, intent(out):: flux_v(:) ! (knon)
39        ! (diagnostic) flux du vent à la surface, en (kg m/s)/(m**2 s)
40        ! flux_v est le flux de moment angulaire (positif vers bas)
41    
42      ! Local:      ! Local:
43      INTEGER i, k      INTEGER i, k
# Line 99  contains Line 99  contains
99         ENDDO         ENDDO
100      ENDDO      ENDDO
101    
     ! flux_v est le flux de moment angulaire (positif vers bas) dont  
     ! l'unite est: (kg m/s)/(m**2 s)  
102      DO i = 1, knon      DO i = 1, knon
103         flux_v(i, 1) = zx_coef(i, 1)/(RG*dtime) &         flux_v(i) = zx_coef(i, 1)/(RG*dtime) &
104              *(local_ven(i, 1)*zx_alf1(i) &              *(local_ven(i, 1)*zx_alf1(i) &
105              +local_ven(i, 2)*zx_alf2(i))              +local_ven(i, 2)*zx_alf2(i))
106      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  
107    
108      DO k = 1, klev      DO k = 1, klev
109         DO i = 1, knon         DO i = 1, knon
Legend:
Removed from v.82  
changed lines
  Added in v.225
  ViewVC Help
Powered by ViewVC 1.1.21