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

Diff of /trunk/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 227 by guez, Thu Nov 2 15:47:03 2017 UTC
# Line 17  contains Line 17  contains
17      INTEGER knon      INTEGER knon
18      REAL, intent(in):: dtime ! intervalle de temps (en s)      REAL, intent(in):: dtime ! intervalle de temps (en s)
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
25      ! vent (dV/dz). La première valeur indique la valeur de Cdrag (sans      ! vent (dV/dz). La première valeur indique la valeur de Cdrag (sans
26      ! unité).      ! unité).
27    
28      REAL t(klon, klev), ven(klon, klev)      REAL, intent(in):: t(klon, klev), ven(klon, klev)
29      ! t--------input-R- temperature (K)      ! t--------input-R- temperature (K)
30      ! ven------input-R- vitesse horizontale (m/s)      ! ven------input-R- vitesse horizontale (m/s)
31      REAL paprs(klon, klev+1), pplay(klon, klev), delp(klon, klev)      REAL paprs(klon, klev+1), pplay(klon, klev), delp(klon, klev)
# Line 35  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 98  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.202  
changed lines
  Added in v.227
  ViewVC Help
Powered by ViewVC 1.1.21