anchor.diff:  ! hard-coded numbers
         surface_physics_nml = 'semic.nml'
         forc_file = 'forcing.nc'
  +      mask_file = "glac.nc"
         nx = 192
         ny = 96
         nt = 365
  
positive.code: 
        end subroutine update
  
        subroutine read_forcing(fnm_in,nx,ny,nt,wind,tt,sf,rf,qq,rhoa,sp,lwd,swd)
            
            integer, intent(in) :: nx !< x dimension
  
negative.code:       subroutine init(dom,domname)
        implicit none
        character(len=*), intent(in) :: domname !< domain name
        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
  