/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_yield.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_yield.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 54  contains Line 54  contains
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)
# Line 627  contains Line 627  contains
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)
Legend:
Removed from v.134  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21