anchor.diff:    double precision :: slip, hmin
     integer nTimeSteps
     double precision :: dumpFreq, avFreq
  -  double precision, dimension(:), allocatable :: zeros
  +  double precision, dimension(:),     allocatable :: zeros
     integer maxits
     double precision :: eps, freesurfFac, thickness_error
     ! Model
  
positive.code:   double precision :: g_vec(layerwise_input_length)
    double precision :: rho0
    ! Wind
    double precision, dimension(:,:), allocatable :: base_wind_x
    double precision, dimension(:,:), allocatable :: base_wind_y
    logical :: DumpWind
    character(30) :: wind_mag_time_series_file
    double precision, dimension(:), allocatable :: wind_mag_time_series
    ! Sponge regions
    double precision, dimension(:,:,:), allocatable :: spongeHTimeScale
    double precision, dimension(:,:,:), allocatable :: spongeUTimeScale
  
negative.code:   ! Velocity component (v)
    double precision, dimension(:,:,:), allocatable :: v
    ! Free surface (eta)
    double precision, dimension(:,:),   allocatable :: eta
    ! Bathymetry
    character(30) :: depthFile
    double precision, dimension(:,:),   allocatable :: depth
    double precision :: H0 ! default depth in no file specified
    ! Grid
    double precision :: dx, dy
    double precision, dimension(:,:),   allocatable :: wetmask
    ! Coriolis parameter at u and v grid-points respectively
    double precision, dimension(:,:),   allocatable :: fu
    double precision, dimension(:,:),   allocatable :: fv
    ! File names to read them from
    character(30) :: fUfile, fVfile
    character(30) :: wetMaskFile
  