anchor.diff:                                 itm_cc, itm_t
   
           namelist /surface_physics/ boundary,tstic,ceff,csh,clh,alb_smax,alb_smin,&
  -                                   albr,albl,tmin,hcrit,amp,&
  +                                   albr,albl,tmin,hcrit,rcrit,amp,&
                                      tau_a, tau_f, w_crit, &
                                      pdd_sigma,pdd_b,&
                                      itm_cc, itm_t, &
  
positive.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)') 'amp        ', par%amp
          write(*,'(a,g13.6)') 'csh        ', par%csh
          write(*,'(a,g13.6)') 'tmin       ', par%tmin
  
negative.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) for surface albedo
          par%rcrit      = rcrit          ! critical snow height (m) for refreezing fraction
          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)
  