anchor.diff:  
     call read_input_file_time_series(wind_mag_time_series_file, &
         wind_mag_time_series, 1.d0, nTimeSteps)
  -
  -  wind_x = base_wind_x*wind_mag_time_series(1)
  -  wind_y = base_wind_y*wind_mag_time_series(1)
   
     call read_input_fileH(spongeHTimeScaleFile, spongeHTimeScale, &
         zeros, nx, ny, layers)
  
positive.code:   character(30) :: wetMaskFile
  
    ! Numerics
    double precision :: pi, dt
    double precision :: au, ar, botDrag
    double precision :: ah(layerwise_input_length)
    double precision :: slip
    double precision :: hmin
    integer nTimeSteps
    double precision :: dumpFreq, avFreq
    integer nwrite, avwrite
    double precision, dimension(:), allocatable :: zeros
    integer maxits
    double precision :: eps, rjac
    double precision :: freesurfFac
    double precision :: thickness_error
  
  
negative.code:   ! Switch for using n + 1/2 layer physics, or using n layer physics
    logical :: RedGrav
  
  
    ! Physics
    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
  