anchor.diff:  
           ! Declaration of namelist parameters
           double precision    :: ceff, albi, albl, alb_smax, alb_smin, hcrit, rcrit, &
  -                               amp, csh, clh, shf_enh, lhf_enh, tmin, tmax, &
  +                               amp, csh, clh, tmin, tmax, &
                                  tstic, &
                                  afac, tmid, &
                                  tau_a, tau_f, w_crit, mcrit
           integer :: n_ksub
   
  -        namelist /surface_physics/ boundary, tstic, ceff, csh, clh, shf_enh, lhf_enh, &
  +        namelist /surface_physics/ boundary, tstic, ceff, csh, clh, &
                                      alb_smax, alb_smin, albi, albl, &
                                      tmin, tmax, hcrit, rcrit, amp, &
                                      tau_a, tau_f, w_crit, mcrit, &
  
positive.code:         tmid          = par%tmid 
  
          n_ksub        = par%n_ksub
          shf_enh       = par%shf_enh
          lhf_enh       = par%lhf_enh
  
          ! Read parameters from input namelist file
          open(7,file=trim(filename))
  
negative.code:         par%tmid       = tmid
  
          par%n_ksub     = n_ksub
  
          ! initialize sub-daily time step tsticsub
          par%tsticsub = par%tstic / dble(par%n_ksub)
  