305 |
* (ypplay(:knon, k - 1) - ypplay(:knon, k)) |
* (ypplay(:knon, k - 1) - ypplay(:knon, k)) |
306 |
ENDDO |
ENDDO |
307 |
|
|
308 |
CALL clcdrag(nsrf, yu(:knon, 1), yv(:knon, 1), yt(:knon, 1), & |
CALL clcdrag(nsrf, sqrt(yu(:knon, 1)**2 + yv(:knon, 1)**2), & |
309 |
yq(:knon, 1), zgeop(:knon, 1), yts(:knon), yqsurf(:knon), & |
yt(:knon, 1), yq(:knon, 1), zgeop(:knon, 1), yts(:knon), & |
310 |
yrugos(:knon), ycdragm(:knon), ycdragh(:knon)) |
yqsurf(:knon), yrugos(:knon), ycdragm(:knon), ycdragh(:knon)) |
311 |
|
|
312 |
IF (iflag_pbl == 1) THEN |
IF (iflag_pbl == 1) THEN |
313 |
ycdragm(:knon) = max(ycdragm(:knon), 0.) |
ycdragm(:knon) = max(ycdragm(:knon), 0.) |