/[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 229 by guez, Mon Nov 6 17:20:45 2017 UTC revision 231 by guez, Mon Nov 6 18:00:54 2017 UTC
# Line 4  module clvent_m Line 4  module clvent_m
4    
5  contains  contains
6    
7    SUBROUTINE clvent(dtime, u1lay, v1lay, coef, t, ven, paprs, pplay, delp, &    SUBROUTINE clvent(dtime, u1lay, v1lay, coef, cdrag, t, ven, paprs, pplay, &
8         d_ven, flux_v)         delp, d_ven, flux_v)
9    
10      ! Author: Z. X. Li (LMD/CNRS)      ! Author: Z. X. Li (LMD/CNRS)
11      ! Date: 1993/08/18      ! Date: 1993/08/18
# Line 20  contains Line 20  contains
20      REAL, intent(in):: u1lay(:), v1lay(:) ! (knon)      REAL, intent(in):: u1lay(:), v1lay(:) ! (knon)
21      ! vent de la premiere couche (m/s)      ! vent de la premiere couche (m/s)
22    
23      REAL, intent(in):: coef(:, :) ! (knon, klev)      REAL, intent(in):: coef(:, 2:) ! (knon, 2: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)
26      ! unité).  
27        REAL, intent(in):: cdrag(:) ! (knon) sans unité
28    
29      REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K)      REAL, intent(in):: t(:, :) ! (knon, klev) ! temperature (K)
30      REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m/s)      REAL, intent(in):: ven(:, :) ! (knon, klev) vitesse horizontale (m/s)
# Line 53  contains Line 54  contains
54      local_ven = ven      local_ven = ven
55    
56      DO i = 1, knon      DO i = 1, knon
57         zx_coef(i, 1) = coef(i, 1) * (1. + SQRT(u1lay(i)**2 + v1lay(i)**2)) &         zx_coef(i, 1) = cdrag(i) * (1. + SQRT(u1lay(i)**2 + v1lay(i)**2)) &
58              * pplay(i, 1) / (RD * t(i, 1)) * dtime * RG              * pplay(i, 1) / (RD * t(i, 1)) * dtime * RG
59      ENDDO      ENDDO
60    
Legend:
Removed from v.229  
changed lines
  Added in v.231
  ViewVC Help
Powered by ViewVC 1.1.21