171 |
REAL ypaprs(klon, klev + 1), ypplay(klon, klev), ydelp(klon, klev) |
REAL ypaprs(klon, klev + 1), ypplay(klon, klev), ydelp(klon, klev) |
172 |
REAL ycoefm0(klon, 2:klev), ycoefh0(klon, 2:klev) |
REAL ycoefm0(klon, 2:klev), ycoefh0(klon, 2:klev) |
173 |
REAL yzlay(klon, klev), zlev(klon, klev + 1), yteta(klon, klev) |
REAL yzlay(klon, klev), zlev(klon, klev + 1), yteta(klon, klev) |
|
REAL ykmm(klon, klev + 1), ykmn(klon, klev + 1) |
|
174 |
REAL yq2(klon, klev + 1) |
REAL yq2(klon, klev + 1) |
175 |
REAL delp(klon, klev) |
REAL delp(klon, klev) |
176 |
INTEGER i, k, nsrf |
INTEGER i, k, nsrf |
374 |
ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), ycdragm(:knon)) |
ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), ycdragm(:knon)) |
375 |
CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), & |
CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), & |
376 |
yu(:knon, :), yv(:knon, :), yteta(:knon, :), yq2(:knon, :), & |
yu(:knon, :), yv(:knon, :), yteta(:knon, :), yq2(:knon, :), & |
377 |
ykmm(:knon, :), ykmn(:knon, :), ustar(:knon)) |
ycoefm(:knon, :), ycoefh(:knon, :), ustar(:knon)) |
|
ycoefm(:knon, :) = ykmm(:knon, 2:klev) |
|
|
ycoefh(:knon, :) = ykmn(:knon, 2:klev) |
|
378 |
END IF |
END IF |
379 |
|
|
380 |
CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), ycoefm(:knon, :), & |
CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), ycoefm(:knon, :), & |
501 |
|
|
502 |
CALL stdlevvar(klon, knon, nsrf, u1(:knon), v1(:knon), tair1(:knon), & |
CALL stdlevvar(klon, knon, nsrf, u1(:knon), v1(:knon), tair1(:knon), & |
503 |
qair1, zgeo1, tairsol, qairsol, rugo1, psfce, patm, yt2m, & |
qair1, zgeo1, tairsol, qairsol, rugo1, psfce, patm, yt2m, & |
504 |
yq2m, yt10m, yq10m, wind10m(:knon), ustar) |
yq2m, yt10m, yq10m, wind10m(:knon), ustar(:knon)) |
505 |
|
|
506 |
DO j = 1, knon |
DO j = 1, knon |
507 |
i = ni(j) |
i = ni(j) |