anchor.diff:        subroutine init(dom,domname)
         implicit none
         character(len=*), intent(in) :: domname !< domain name
  -      character(len=10) :: forc_file !< name of file is allocated, fuck you
  +      character(len=10) :: forc_file !< name of file is allocated
  +      character(len=7) :: mask_file !< name the mask file
         type(smb_class), intent(in out) :: dom  !< domain object
         integer :: nx, ny, nt, nloop
         character(len=256) :: surface_physics_nml
  
positive.code: ! hard-coded numbers
        surface_physics_nml = 'semic.nml'
        forc_file = 'forcing.nc'
        nx = 192
        ny = 96
        nt = 365
  
negative.code:       dom%surface%now%tsurf = 260.0
        dom%surface%now%alb_snow = 0.8
  
        call read_mask(trim(mask_file), nx, ny, dom%surface%now%mask)
        !dom%surface%now%mask = mask ! all points read from glac mask. 
  
        ! read model parameters from namelist file
        dom%surface_physics_nml = surface_physics_nml
  