anchor.diff:          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    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,rcrit,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
  
positive.code:         par%albl       = albl           ! background bare land (bare land)
          par%alb_smax   = alb_smax       ! maximum snow albedo (fresh snow)
          par%alb_smin   = alb_smin       ! minimum snow albedo (wet/old snow)
          par%hcrit      = hcrit          ! critical snow height (m)
          par%amp        = amp            ! diurnal cycle amplitude (K)
          par%csh        = csh            ! turbulent heat exchange coeffcicient (sensible and latent heat)
          par%clh        = clh            ! turbulent heat exchange coeffcicient (sensible and latent heat)
  
negative.code:         write(*,'(a,g13.6)') 'alb_smax   ', par%alb_smax
          write(*,'(a,g13.6)') 'alb_smin   ', par%alb_smin
          write(*,'(a,g13.6)') 'hcrit      ', par%hcrit
          write(*,'(a,g13.6)') 'rcrit      ', par%rcrit
          write(*,'(a,g13.6)') 'amp        ', par%amp
          write(*,'(a,g13.6)') 'csh        ', par%csh
          write(*,'(a,g13.6)') 'tmin       ', par%tmin
  