14 |
|
|
15 |
! L. Fairhead, April 2000 |
! L. Fairhead, April 2000 |
16 |
|
|
17 |
USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats |
USE fcttre, ONLY: foede, foeew |
18 |
use nr_util, only: assert_eq |
use nr_util, only: assert_eq |
19 |
USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt |
USE suphec_m, ONLY: rcpd, rd, retv, rlstt, rlvtt, rtt |
20 |
USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2 |
USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2 |
90 |
sl = merge(RLSTT, RLVTT, tsurf < RTT) |
sl = merge(RLSTT, RLVTT, tsurf < RTT) |
91 |
|
|
92 |
! Q |
! Q |
93 |
oq = 1. - (beta * coef * peqBcoef * dtime) |
oq = 1. - beta * coef * peqBcoef * dtime |
94 |
mq = beta * coef * (peqAcoef - qsat + dq_s_dt * tsurf) / oq |
mq = beta * coef * (peqAcoef - qsat + dq_s_dt * tsurf) / oq |
95 |
nq = beta * coef * (- 1. * dq_s_dt) / oq |
nq = beta * coef * (- 1. * dq_s_dt) / oq |
96 |
|
|