20 |
|
|
21 |
|
|
22 |
! inputs: |
! inputs: |
23 |
integer len, nd |
integer, intent(in):: len, nd |
24 |
integer nk(len), icb(len) |
integer nk(len), icb(len) |
25 |
real, intent(in):: t(len,nd) |
real, intent(in):: t(len,nd) |
26 |
real q(len,nd), qs(len,nd), gz(len,nd) |
real, intent(in):: q(len,nd), qs(len,nd), gz(len,nd) |
27 |
real p(len,nd) |
real, intent(in):: p(len,nd) |
28 |
real plcl(len) ! convect3 |
real plcl(len) ! convect3 |
29 |
|
|
30 |
! outputs: |
! outputs: |
33 |
! local variables: |
! local variables: |
34 |
integer i, k |
integer i, k |
35 |
integer icb1(len), icbs(len), icbsmax2 ! convect3 |
integer icb1(len), icbs(len), icbsmax2 ! convect3 |
36 |
real tg, qg, alv, s, ahg, tc, denom, es, rg |
real tg, qg, alv, s, ahg, tc, denom, es |
37 |
real ah0(len), cpp(len) |
real ah0(len), cpp(len) |
38 |
real tnk(len), qnk(len), gznk(len), ticb(len), gzicb(len) |
real tnk(len), qnk(len), gznk(len), ticb(len), gzicb(len) |
39 |
real qsicb(len) ! convect3 |
real qsicb(len) ! convect3 |
176 |
clw(i,icbs(i))=qnk(i)-qg |
clw(i,icbs(i))=qnk(i)-qg |
177 |
clw(i,icbs(i))=max(0.0,clw(i,icbs(i))) |
clw(i,icbs(i))=max(0.0,clw(i,icbs(i))) |
178 |
|
|
|
rg=qg/(1.-qnk(i)) |
|
|
! ori tvp(i,icb(i))=tp(i,icb(i))*(1.+rg*epsi) |
|
179 |
! convect3: (qg utilise au lieu du vrai mixing ratio rg) |
! convect3: (qg utilise au lieu du vrai mixing ratio rg) |
180 |
tvp(i,icbs(i))=tp(i,icbs(i))*(1.+qg/eps-qnk(i)) !whole thing |
tvp(i,icbs(i))=tp(i,icbs(i))*(1.+qg/eps-qnk(i)) !whole thing |
181 |
|
|
273 |
clw(i,icb(i)+1)=qnk(i)-qg |
clw(i,icb(i)+1)=qnk(i)-qg |
274 |
clw(i,icb(i)+1)=max(0.0,clw(i,icb(i)+1)) |
clw(i,icb(i)+1)=max(0.0,clw(i,icb(i)+1)) |
275 |
|
|
|
rg=qg/(1.-qnk(i)) |
|
|
! ori tvp(i,icb(i))=tp(i,icb(i))*(1.+rg*epsi) |
|
276 |
! convect3: (qg utilise au lieu du vrai mixing ratio rg) |
! convect3: (qg utilise au lieu du vrai mixing ratio rg) |
277 |
tvp(i,icb(i)+1)=tp(i,icb(i)+1)*(1.+qg/eps-qnk(i)) !whole thing |
tvp(i,icb(i)+1)=tp(i,icb(i)+1)*(1.+qg/eps-qnk(i)) !whole thing |
278 |
|
|