xx=trant(i,j)*dlayer(i,j,k2)

         flcu(i,j,k1) =flcu(i,j,k1)+xx

         flxu(i,j,k1) =flxu(i,j,k1)+xx*fclr(i,j)

         xx=trant(i,j)*dlayer(i,j,k1)

         flcd(i,j,k2) =flcd(i,j,k2)+xx

         flxd(i,j,k2) =flxd(i,j,k2)+xx*fclr(i,j)

      _RL xx,xx1,xx2,xx3

           xx=dh2o(i,j,k)*(pa(i,j,k)/500.0)

           h2oexp(i,j,k,1)=exp(-xx*0.10624)

           xx=h2oexp(i,j,k,1)*h2oexp(i,j,k,1)

           xx=xx*xx

           h2oexp(i,j,k,2)=xx*xx

           xx=h2oexp(i,j,k,2)*h2oexp(i,j,k,2)

           xx=xx*xx

           h2oexp(i,j,k,3)=xx*xx

           xx=h2oexp(i,j,k,3)*h2oexp(i,j,k,3)

           xx=xx*xx

           h2oexp(i,j,k,4)=xx*xx

           xx=h2oexp(i,j,k,4)*h2oexp(i,j,k,4)

           xx=xx*xx

           h2oexp(i,j,k,5)=xx*xx

           xx=h2oexp(i,j,k,5)*h2oexp(i,j,k,5)

           xx=xx*xx

           h2oexp(i,j,k,6)=xx*xx

           xx=dco2(i,j,k)*(pa(i,j,k)/300.0)**0.5

           co2exp(i,j,k,1,1)=exp(-xx*2.656e-5)

           xx=co2exp(i,j,k,1,1)*co2exp(i,j,k,1,1)

           xx=xx*xx

           co2exp(i,j,k,2,1)=xx*xx

           xx=co2exp(i,j,k,2,1)*co2exp(i,j,k,2,1)

           xx=xx*xx

           co2exp(i,j,k,3,1)=xx*xx

           xx=co2exp(i,j,k,3,1)*co2exp(i,j,k,3,1)

           xx=xx*xx

           co2exp(i,j,k,4,1)=xx*xx

           xx=co2exp(i,j,k,4,1)*co2exp(i,j,k,4,1)

           xx=xx*xx

           co2exp(i,j,k,5,1)=xx*xx

           xx=co2exp(i,j,k,5,1)*co2exp(i,j,k,5,1)

           xx=xx*xx

           co2exp(i,j,k,6,1)=xx*xx

           xx=dn2o(i,j,k)*(1.+(1.4476e-3+3.6656e-6*dt(i,j,k))*dt(i,j,k))

           n2oexp(i,j,k,1)=exp(-xx*0.25238)

           xx=n2oexp(i,j,k,1)*n2oexp(i,j,k,1)

           xx1=xx*xx

           n2oexp(i,j,k,2)=xx*xx1*xx2*xx3

      _RL xx

           xx=   0.3153*th2o(i,j,1)

           xx=xx+0.4604*th2o(i,j,2)

           xx=xx+0.1326*th2o(i,j,3)

           xx=xx+0.0798*th2o(i,j,4)

           xx=xx+0.0119*th2o(i,j,5)

           tran(i,j)=tran(i,j)*xx

           xx=    0.2673*tco2(i,j,1,1)

           xx=xx+ 0.2201*tco2(i,j,2,1)

           xx=xx+ 0.2106*tco2(i,j,3,1)

           xx=xx+ 0.2409*tco2(i,j,4,1)

           xx=xx+ 0.0196*tco2(i,j,5,1)

           xx=xx+ 0.0415*tco2(i,j,6,1)

           tran(i,j)=tran(i,j)*xx

           xx=   0.970831*tn2o(i,j,1)

           xx=xx+0.029169*tn2o(i,j,2)

           tran(i,j)=tran(i,j)*(xx-1.0)

      _RL dp,xx,p1,dwe,dpe,a1,b1,fk1,a2,b2,fk2

         xx=pa(i,j,k)*0.001618*wa(i,j,k)*wa(i,j,k)*dp

         dcont(i,j,k) = xx*exp(1800./ta(i,j,k)-6.081)+1.e-10

         print *, 'pathei-lsopt xx(1), gg(1) ', xx(1), gg(1)

         call SIMUL( indic, nn, xx, ff, gg )

         print *, 'pathei: nach simul:   xx(1), gg(1) = ', xx(1), gg(1)

         call DOSTORE( nn, xx, .true., 1 )

         call LSWRI( isize, iiter, nn, xx, gg, lphprint )

      subroutine LSOPT_TOP( nn, xx, ff, gg, simul, optline

         call DOSTORE( nn, xx, .false., 1 )

     $           , DNRM2( nn, xx, 1 ), 0.

     &     , gg, dd, xx, xdiff

     &     , dd, gg, xx, xdiff

      call DOSTORE( nn, xx, .true., 1 )

      call LSWRI( isize, iiter, nn, xx, gg, lphprint )

     $     , DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 )

     &        , gg, dd, xx, xdiff

      double precision  xx(nn), ff, gg(nn), epsx, fmin, epsg

     $        , '  norm of x...................', DNRM2( nn, xx, 1 )

      real*4 xx(NN), gg(NN)

      subroutine LSWRI( isize, phniter, NN, xx, gg, lphprint )

      write(95,REC=phniter+1) xx

         xdiff(i) = xx(i) + tact*dd(i)

      double precision    gg(nn), dd(nn), xx(nn), xdiff(nn)

     &     , gg, dd, xx, xdiff