anchor.diff:          alb_smax  = dom%par%alb_smax
           alb_smin  = dom%par%alb_smin
           hcrit     = dom%par%hcrit
  +        rcrit     = dom%par%rcrit
           csh       = dom%par%csh
           tstic     = dom%par%tstic
           ceff      = dom%par%ceff
  
positive.code:         character (len=256) :: name, boundary(30), alb_scheme
          character (len=3)   :: method
          integer             :: nx
          double precision    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,amp,csh,tmin,tstic,clh
          double precision    :: pdd_sigma, pdd_b
          double precision    :: itm_cc, itm_t
          double precision    :: tau_a, tau_f, w_crit
  
negative.code:         integer, intent(in) :: day, year
          
          ! physical parameters used in the different parameterizations
          double precision :: tmin, albr, albl, alb_smax, alb_smin, hcrit, rcrit, &
                              ceff, tstic, amp, clh, csh, &
                              pdd_sigma, pdd_b, &
                              itm_cc, itm_t, &
                              tau_a, tau_f, w_crit, &
                              mmfac, Pmaxfrac
  
          ! ocean/land/ice mask
          double precision, dimension(dom%par%nx) :: land_ice_ocean
  