anchor.diff:  subroutine read_input_fileH(name, array, default, nx, ny, layers)
     implicit none
   
  -  character(30), intent(in) :: name
  +  character(60), intent(in) :: name
     double precision, intent(out) :: array(0:nx+1, 0:ny+1, layers)
     double precision, intent(in) :: default(layers)
     integer, intent(in) :: nx, ny, layers
  
positive.code:   double precision, dimension(:,:),   allocatable :: fu
    double precision, dimension(:,:),   allocatable :: fv
    ! File names to read them from
    character(30) :: fUfile, fVfile
    character(30) :: wetMaskFile
    ! Numerics
    double precision :: dt
    double precision :: au, ar, botDrag
  
negative.code:   double precision, dimension(:,:,:), allocatable :: spongeH
    double precision, dimension(:,:,:), allocatable :: spongeU
    double precision, dimension(:,:,:), allocatable :: spongeV
    character(60) :: spongeHTimeScaleFile
    character(60) :: spongeUTimeScaleFile
    character(60) :: spongeVTimeScaleFile
    character(60) :: spongeHfile
    character(60) :: spongeUfile
    character(60) :: spongeVfile
    ! Main input files
    character(60) :: initUfile, initVfile, initHfile, initEtaFile
    character(60) :: zonalWindFile, meridionalWindFile
  
    ! Set default values here
  
  