--- trunk/Sources/phylmd/CV3_routines/cv3_yield.f 2015/04/29 15:47:56 134 +++ trunk/Sources/phylmd/CV3_routines/cv3_yield.f 2016/03/11 18:47:26 178 @@ -54,7 +54,7 @@ real rat, awat, delti 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)