double precision :: area_shelf_streamice_ad(1-olx:snx+olx,1-oly:

      double precision :: b_glen_ad(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)

      double precision :: bdot_streamice_ad(1-olx:snx+olx,1-oly:sny+oly,

      double precision :: c_basal_friction_ad(1-olx:snx+olx,1-oly:sny+

      double precision :: float_frac_streamice_ad(1-olx:snx+olx,1-oly:

      double precision :: h_streamice_ad(1-olx:snx+olx,1-oly:sny+oly,

      double precision :: surf_el_streamice_ad(1-olx:snx+olx,1-oly:sny+

      double precision :: tau_beta_eff_streamice_ad(1-olx:snx+olx,1-oly:

      double precision :: u_streamice_ad(1-olx:snx+olx,1-oly:sny+oly,

      double precision :: v_streamice_ad(1-olx:snx+olx,1-oly:sny+oly,

      double precision :: visc_streamice_ad(1-olx:snx+olx,1-oly:sny+oly,

      double precision :: r_low_si_ad(1-olx:snx+olx,1-oly:sny+oly,nsx,

      double precision, dimension(:), allocatable ::  tmp

        double precision npart,xpart,ypart,kpart,kfloat,iup,itop

        double precision depth(Nx,Ny),xc(Nx),yc(Ny),degX,degY

      double precision fac

      double precision ferr, berr

      double precision phtmpin(maxn), phtmpout(maxn)

      double precision phxout(maxn)

      double precision tmpvec1(maxn), tmpvec2(maxn)

      double precision phrwork(3*maxn)

      double precision sign, den, anum

      double precision z1, b, discri

      double precision    t, f, fp, ta, fa, fpa, tlower, tupper

      double precision phtmpadmtlm(maxn)

      double precision    gold(nn), xdiff(nn), diag(nn)

      double precision    r1, rmin, den

      double precision    x(n)

      double precision dd(nn), alpha(100), xdiff(nn)

      double precision     DDOT

      double precision    r

      double precision    fc, gnorm0

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

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

      double precision xpara1, xpara2

      double precision    factor

      double precision    barmax

      double precision    barmul

      double precision    barmin

      double precision    tg, fg, td, ta

      double precision    fa, fpa, fp

      double precision    fnew, fdiff

      double precision    z, test, barr

      double precision    left, right, told

      double precision     DDOT

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

      double precision  dd(nn), gold(nn), xdiff(nn)

      double precision phff

      double precision   fupd

      double precision   r1, tmin, tmax, tact, gnorm, gnorm0, eps1

      double precision   fold, ys

      double precision   dotdg

      double precision     DDOT, DNRM2

      double precision rmin

      double precision    ff, fmin, fold, gnorm0, dotdg

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

      double precision    tmin, tmax, tact, epsx

      double precision    fdiff, preco

      double precision     DDOT

      double precision    fc, gnorm0