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) |
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) |