anchor.diff:    ! Main input files
     character(60) :: initUfile, initVfile, initHfile, initEtaFile
     character(60) :: zonalWindFile, meridionalWindFile
  -  logical :: RelativeWind
  +  logical       :: RelativeWind
     double precision :: Cd
   
     ! External pressure solver variables
  
positive.code:   character(60) :: fUfile, fVfile
    character(60) :: wetMaskFile
    ! Numerics
    double precision :: dt
    double precision :: au, ar, botDrag
    double precision :: kh(layerwise_input_length), kv
    double precision :: slip, hmin
    integer          :: niter0, nTimeSteps
    double precision :: dumpFreq, avFreq, checkpointFreq, diagFreq
    double precision, dimension(:),     allocatable :: zeros
    integer maxits
    double precision :: eps, freesurfFac, thickness_error
    integer          :: debug_level
    integer          :: hAdvecScheme
    integer          :: TS_algorithm
    integer          :: AB_order
  
    ! Model
    double precision :: hmean(layerwise_input_length)
    ! 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
    character(60) :: wind_mag_time_series_file
    double precision, dimension(:),     allocatable :: wind_mag_time_series
    ! Sponge regions
  
negative.code: 
    integer, parameter :: layerwise_input_length = 10000
    ! Resolution
    integer :: nx ! number of x grid points
    integer :: ny ! number of y grid points
    integer :: layers ! number of active layers in the model
    ! Layer thickness (h)
    double precision, dimension(:,:,:), allocatable :: h
    ! Velocity component (u)
  