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 |
|
|
30 |
|
|
31 |
REAL, intent(in):: dtime |
REAL, intent(in):: dtime |
32 |
! dtime----input-R- intervalle de temps (en second) |
! dtime----input-R- intervalle de temps (en second) |
33 |
REAL u1lay(klon), v1lay(klon) |
REAL, intent(in):: u1lay(klon), v1lay(klon) ! vent de la premiere |
34 |
! u1lay----input-R- vent u de la premiere couche (m/s) |
! couche (m/s) |
35 |
! v1lay----input-R- vent v de la premiere couche (m/s) |
REAL coef(:, 2:) ! (klon, 2:klev) |
|
REAL coef(klon, 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) |
107 |
|
|
108 |
local_trs = trs |
local_trs = trs |
109 |
|
|
110 |
! Attention si dans clmain zx_alf1(i) = 1. |
! Attention si dans pbl_surface zx_alf1(i) = 1. |
111 |
! Il doit y avoir coherence (donc la meme chose ici) |
! Il doit y avoir coherence (donc la meme chose ici) |
112 |
|
|
113 |
DO i = 1, klon |
DO i = 1, klon |
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 |