anchor.diff:    ! Main input files
     character(60) :: initUfile, initVfile, initHfile, initEtaFile
     character(60) :: zonalWindFile, meridionalWindFile
  +  logical :: RelativeWind
  +  double precision :: Cd
   
     ! External pressure solver variables
     integer :: nProcX, nProcY
  
positive.code:   logical,          intent(in) :: RedGrav
    ! Whether to write computed wind in the output
    logical,          intent(in) :: DumpWind
  
    double precision :: dhdt(0:nx+1, 0:ny+1, layers)
    double precision :: dhdtold(0:nx+1, 0:ny+1, layers)
  
negative.code:     dx, dy, wetmask, hfacW, hfacE, hfacN, hfacS, fu, fv, &
      au, ar, botDrag, ah, slip, &
      RedGrav, g_vec, rho0, wind_x, wind_y, &
      RelativeWind, Cd, &
      spongeHTimeScale, spongeH, &
      spongeUTimeScale, spongeU, &
      spongeVTimeScale, spongeV, &
  