1 |
module cv3_yield_m |
module cv30_yield_m |
2 |
|
|
3 |
implicit none |
implicit none |
4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE cv3_yield(nloc,ncum,nd,na & |
SUBROUTINE cv30_yield(nloc,ncum,nd,na & |
8 |
,icb,inb,delt & |
,icb,inb,delt & |
9 |
,t,rr,u,v,gz,p,ph,h,hp,lv,cpn,th & |
,t,rr,u,v,gz,p,ph,h,hp,lv,cpn,th & |
10 |
,ep,clw,m,tp,mp,rp,up,vp & |
,ep,clw,m,tp,mp,rp,up,vp & |
14 |
,iflag,precip,VPrecip,ft,fr,fu,fv & |
,iflag,precip,VPrecip,ft,fr,fu,fv & |
15 |
,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd) |
,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd) |
16 |
use conema3_m |
use conema3_m |
17 |
use cv3_param_m |
use cv30_param_m |
18 |
use cvthermo |
use cvthermo |
19 |
use cvflag |
use cvflag |
20 |
|
|
54 |
real rat, awat, delti |
real rat, awat, delti |
55 |
real ax, bx, cx, dx |
real ax, bx, cx, dx |
56 |
real cpinv, rdcp, dpinv |
real cpinv, rdcp, dpinv |
57 |
real lvcp(nloc,na), mke(nloc,na) |
real lvcp(nloc,na) |
58 |
real am(nloc), work(nloc), ad(nloc), amp1(nloc) |
real am(nloc), work(nloc), ad(nloc), amp1(nloc) |
59 |
!!! real up1(nloc), dn1(nloc) |
!!! real up1(nloc), dn1(nloc) |
60 |
real up1(nloc,nd,nd), dn1(nloc,nd,nd) |
real up1(nloc,nd,nd), dn1(nloc,nd,nd) |
230 |
enddo |
enddo |
231 |
if(num1.le.0) cycle |
if(num1.le.0) cycle |
232 |
|
|
233 |
call zilch(amp1,ncum) |
amp1(:ncum) = 0. |
234 |
call zilch(ad,ncum) |
ad(:ncum) = 0. |
235 |
|
|
236 |
do k=i+1,nl+1 |
do k=i+1,nl+1 |
237 |
do il=1,ncum |
do il=1,ncum |
627 |
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
628 |
|
|
629 |
do i=1,nd |
do i=1,nd |
|
do il=1,ncum |
|
|
mke(il,i)=upwd(il,i)+dnwd(il,i) |
|
|
enddo |
|
|
enddo |
|
|
|
|
|
do i=1,nd |
|
630 |
DO il=1,ncum |
DO il=1,ncum |
631 |
rdcp=(rrd*(1.-rr(il,i))-rr(il,i)*rrv) & |
rdcp=(rrd*(1.-rr(il,i))-rr(il,i)*rrv) & |
632 |
/(cpd*(1.-rr(il,i))+rr(il,i)*cpv) |
/(cpd*(1.-rr(il,i))+rr(il,i)*cpv) |
697 |
enddo ! cld |
enddo ! cld |
698 |
enddo ! cld |
enddo ! cld |
699 |
|
|
700 |
end SUBROUTINE cv3_yield |
end SUBROUTINE cv30_yield |
701 |
|
|
702 |
end module cv3_yield_m |
end module cv30_yield_m |