anchor.diff:      public :: surface_physics_class
       public :: surface_physics_par_load, surface_alloc, surface_dealloc
       public :: surface_boundary_define, surface_physics_average
  +    public :: print_param, print_boundary_opt
   
   contains
   
  
positive.code:     end subroutine surface_boundary_define
  
  
  end module surface_physics
  
negative.code:         integer, intent(in) :: day, year
          
          ! physical parameters used in the different parameterizations
          double precision :: tmin, albr, albl, alb_smax, alb_smin, hcrit, csh, &
              ceff, tstic, amp, clh, &
              pdd_sigma, pdd_b, &
              itm_cc, itm_t, &
  