anchor.diff:  
           ! Declaration of namelist parameters
           double precision    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,rcrit,  &
  -                               amp,csh,clh,tmin,tmax,&
  +                               amp,csh,clh,shf_enh,lhf_enh,tmin,tmax,&
                                  tstic, Ts_wt(2), &
                                  afac,tmid, &
                                  tau_a, tau_f, w_crit, mcrit, &
  
positive.code:                     subl = 0.0_dp
                      evap = 0.0_dp
                      lhf  = 0.0_dp
                      call latent_heat_flux(tsurf_new, usurf, shum, psurf, rhoatm, clh, eps, cls, clv, lhf, subl, evap)
                  end if
  
                  ! outgoing long-wave flux following Stefan-Boltzmann's law (W/m^2)
  
negative.code:         
          ! physical parameters used in the different parameterizations
          double precision :: tmin, tmax, albr, albl, alb_smax, alb_smin, hcrit, rcrit, &
                              ceff, tstic, amp, clh, csh, shf_enh, lhf_enh, &
                              pdd_sigma, pdd_b, &
                              itm_cc, itm_t, &
                              tau_a, tau_f, w_crit, mcrit, &
                              mmfac, Pmaxfrac, &
                              Ts_wt(2), afac,tmid 
          double precision :: tsticsub
          integer :: ksub, n_ksub
  
  