fg     = ff

      fa     = ff

         fdiff = fnew - ff

         ff = fnew

            ff = fg

      ff    = fg

      double precision ff, dotdg, tmin, tmax, tact, epsx

     &     , ff, dotdg

      _RL  zth,ff,xx,taup,sscp,gp,gm1,gm2,gm3,akk,alf1,alf2,

                ff  = g0*g0

                xx  = one-ff*ssc

                sscp= ssc*(one-ff)/xx

           ff=0.5+(0.3739+(0.0076+0.1185*gg)*gg)*gg

           taux=taux*(1.-ww*ff)

         ff=0.5+(0.3739+(0.0076+0.1185*asyal(i,j,k2-1,ib))

         taux=taual(i,j,k2-1,ib)*(1.-ssaal(i,j,k2-1,ib)*ff)

      _RL w1,w2,w3,g1,g2,g3,ww,gg,ff,taux,reff1,reff2

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

         print *, 'pathei: nach simul:   nn, ff = ', nn, ff

         phff = ff

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

         ff = fupd

         phff = ff

     &        ff

            print iform, iiter, tact, ifunc, ff, gnorm0

      fold = ff

      phff = ff

     &     , ff, fmin, fold, gnorm0, dotdg

     &     , ff, dotdg

      print iform, iiter, tact, ifunc, ff, gnorm

      call OUTSTORE( nn, ff, gnorm0, nupdate, jmin, jmax )

     &        , ff, fmin, fold, gnorm0, dotdg

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

     $        , '  cost function...............', ff

      double precision    ff, fmin, fold, gnorm0, dotdg

     &     , ff, fmin, fold, gnorm0, dotdg

         fdiff = fold - ff