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

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

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

revision 225 by guez, Mon Oct 16 12:35:41 2017 UTC revision 232 by guez, Tue Nov 7 10:23:25 2017 UTC
# Line 4  module cltracrn_m Line 4  module cltracrn_m
4    
5  contains  contains
6    
7    SUBROUTINE cltracrn(itr, dtime, u1lay, v1lay, coef, t, ftsol, pctsrf, tr, &    SUBROUTINE cltracrn(itr, dtime, u1lay, v1lay, coef, cdragh, t, ftsol, &
8         trs, paprs, pplay, delp, masktr, fshtr, hsoltr, tautr, vdeptr, lat, &         pctsrf, tr, trs, paprs, pplay, delp, masktr, fshtr, hsoltr, tautr, &
9         d_tr, d_trs)         vdeptr, lat, d_tr, d_trs)
10    
11      ! From phylmd/cltracrn.F, version 1.2 2005/05/25 13:10:09      ! From phylmd/cltracrn.F, version 1.2 2005/05/25 13:10:09
12    
# Line 32  contains Line 32  contains
32      ! dtime----input-R- intervalle de temps (en second)      ! dtime----input-R- intervalle de temps (en second)
33      REAL, intent(in):: u1lay(klon), v1lay(klon) ! vent de la premiere      REAL, intent(in):: u1lay(klon), v1lay(klon) ! vent de la premiere
34                                                  ! couche (m/s)                                                  ! couche (m/s)
35      REAL coef(klon, klev)      REAL coef(:, 2:) ! (klon, 2:klev)
36      ! coef-----input-R- le coefficient d'echange (m**2/s) l>1      ! coef-----input-R- le coefficient d'echange (m**2/s) l>1
37        real cdragh(:) ! klon
38      REAL, intent(in):: t(klon, klev) ! temperature (K)      REAL, intent(in):: t(klon, klev) ! temperature (K)
39      real, intent(in):: ftsol(klon, nbsrf), pctsrf(klon, nbsrf)      real, intent(in):: ftsol(klon, nbsrf), pctsrf(klon, nbsrf)
40      ! ftsol----input-R- temperature du sol (en Kelvin)      ! ftsol----input-R- temperature du sol (en Kelvin)
# Line 115  contains Line 116  contains
116      ENDDO      ENDDO
117    
118      DO i = 1, klon      DO i = 1, klon
119         zx_coef(i, 1) = coef(i, 1) &         zx_coef(i, 1) = cdragh(i) &
120              * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) &              * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) &
121              * pplay(i, 1)/(RD*t(i, 1))              * pplay(i, 1)/(RD*t(i, 1))
122         zx_coef(i, 1) = zx_coef(i, 1) * dtime*RG         zx_coef(i, 1) = zx_coef(i, 1) * dtime*RG
Legend:
Removed from v.225  
changed lines
  Added in v.232
  ViewVC Help
Powered by ViewVC 1.1.21