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:         alb_smax  = dom%par%alb_smax
          alb_smin  = dom%par%alb_smin
          hcrit     = dom%par%hcrit
          heff      = dom%par%heff
          csh       = dom%par%csh
          tstic     = dom%par%tstic
          ceff      = dom%par%ceff
  
negative.code:         character (len=256) :: name, boundary(30), alb_scheme
          character (len=3)   :: method
          integer             :: nx
          double precision    :: ceff,albr,albl,alb_smax,alb_smin,hcrit,amp,csh,tmin,tstic,clh
          double precision    :: pdd_sigma, pdd_b
          double precision    :: itm_cc, itm_t
          double precision    :: tau_a, tau_f, w_crit
  