anchor.diff:              end if
   
               ! refreezing as fraction of melt (increases with snow height)
  -            f_rz = hsnow/(hsnow+hcrit)
  +            f_rz = hsnow/(hsnow+rcrit)
               if (.not. bnd%refr) then
                   !f_rz = (1._dp-dexp(-hsnow))
                   !refr = f_rz*dmin1(qcold/(rhow*clm),melt+rf)
  
positive.code:         alb_smax  = dom%par%alb_smax
          alb_smin  = dom%par%alb_smin
          hcrit     = dom%par%hcrit
          csh       = dom%par%csh
          tstic     = dom%par%tstic
          ceff      = dom%par%ceff
  
negative.code: 
          ! auxillary variables
          double precision, dimension(dom%par%nx) :: esat_sur, lwu, qsb, &
              qmelt, acc, runoff, qcold, below, above, shum_sat, f_rz, &
              melted_snow, melted_ice, new_ice, refreezing_max, refrozen_rain, &
              refrozen_snow, refrozen_snow_max, runoff_rain, runoff_snow, &
              shf_enhnc
  