118 |
c |
c |
119 |
REAL dtime |
REAL dtime |
120 |
real date0 |
real date0 |
121 |
integer itap |
integer, intent(in):: itap |
122 |
REAL t(klon,klev), q(klon,klev) |
REAL t(klon,klev), q(klon,klev) |
123 |
REAL u(klon,klev), v(klon,klev) |
REAL u(klon,klev), v(klon,klev) |
124 |
cIM 230604 BAD REAL radsol(klon) ??? |
cIM 230604 BAD REAL radsol(klon) ??? |
125 |
REAL, intent(in):: paprs(klon,klev+1) |
REAL, intent(in):: paprs(klon,klev+1) |
126 |
real pplay(klon,klev) |
real, intent(in):: pplay(klon,klev) |
127 |
REAL, intent(in):: rlon(klon), rlat(klon) |
REAL, intent(in):: rlon(klon), rlat(klon) |
128 |
real cufi(klon), cvfi(klon) |
real cufi(klon), cvfi(klon) |
129 |
REAL d_t(klon, klev), d_q(klon, klev) |
REAL d_t(klon, klev), d_q(klon, klev) |
356 |
ENDDO |
ENDDO |
357 |
ENDDO |
ENDDO |
358 |
DO i = 1, klon ! vent de la premiere couche |
DO i = 1, klon ! vent de la premiere couche |
|
ccc zx_alf1 = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2)) |
|
359 |
zx_alf1 = 1.0 |
zx_alf1 = 1.0 |
360 |
zx_alf2 = 1.0 - zx_alf1 |
zx_alf2 = 1.0 - zx_alf1 |
361 |
u1lay(i) = u(i,1)*zx_alf1 + u(i,2)*zx_alf2 |
u1lay(i) = u(i,1)*zx_alf1 + u(i,2)*zx_alf2 |