anchor.diff:    call read_input_fileV(initVfile, v, 0.d0, nx, ny, layers)
     call read_input_fileH(initHfile, h, hmean, nx, ny, layers)
   
  -  if (.not. RedGrav) then
  -    call read_input_fileH_2D(depthFile, depth, H0, nx, ny)
  -    call read_input_fileH_2D(initEtaFile, eta, 0.d0, nx, ny)
  -    ! Check that depth is positive - it must be greater than zero
  -    if (minval(depth) .lt. 0) then
  -      write(17, "(A)") "Depths must be positive."
  -      call clean_stop(0, .FALSE.)
  -    end if
  -  end if
  -
     call read_input_fileU(fUfile, fu, 0.d0, nx, ny, 1)
     call read_input_fileV(fVfile, fv, 0.d0, nx, ny, 1)
   
  
positive.code:   double precision array_small(nx, ny, layers)
    integer k
  
    if (name.ne.'') then
      open(unit=10, form='unformatted', file=name)
      read(10) array_small
  
negative.code:         spongeVTimeScale, spongeV, &
          nx, ny, layers, n, debug_level)
  
  
      ! Use dh/dt, du/dt and dv/dt to step h, u and v forward in time with
      ! the Adams-Bashforth third order linear multistep method
  
  