--- trunk/phylmd/CV3_routines/cv3_yield.f 2014/04/25 14:58:31 97 +++ trunk/Sources/phylmd/CV30_routines/cv3_yield.f 2016/03/16 14:50:46 184 @@ -52,9 +52,9 @@ ! local variables: integer i,k,il,n,j,num1 real rat, awat, delti - real ax, bx, cx, dx, ex + real ax, bx, cx, dx real cpinv, rdcp, dpinv - real lvcp(nloc,na), mke(nloc,na) + real lvcp(nloc,na) real am(nloc), work(nloc), ad(nloc), amp1(nloc) !!! real up1(nloc), dn1(nloc) real up1(nloc,nd,nd), dn1(nloc,nd,nd) @@ -627,12 +627,6 @@ !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc do i=1,nd - do il=1,ncum - mke(il,i)=upwd(il,i)+dnwd(il,i) - enddo - enddo - - do i=1,nd DO il=1,ncum rdcp=(rrd*(1.-rr(il,i))-rr(il,i)*rrv) & /(cpd*(1.-rr(il,i))+rr(il,i)*cpv)