anchor.diff:      v(:, :, k) = v(:, :, k) * hfacS * wetmask(:, :)
     end do
   
  -  ! ! If the winds are static, then set wind_ = base_wind_
  -  ! if (.not. UseSinusoidWind .and. .not. UseStochWind)  then
  -  !   wind_x = base_wind_x
  -  !   wind_y = base_wind_y
  -  ! end if
  -
  -  ! ! Initialise random numbers for stochastic winds
  -  ! if (UseStochWind) then
  -  !   ! This ensures a different series of random numbers each time the
  -  !   ! model is run.
  -  !   call ranseed()
  -  !   ! Number of timesteps between updating the perturbed wind field.
  -  !   n_stoch_wind = int(wind_period/dt)
  -  ! end if
   
     if (.not. RedGrav) then
       ! Initialise arrays for pressure solver
  
positive.code:   namelist /INITIAL_CONDITONS/ initUfile, initVfile, initHfile, initEtaFile
  
    namelist /EXTERNAL_FORCING/ zonalWindFile, meridionalWindFile, &
        ! UseSinusoidWind, UseStochWind, wind_alpha, wind_beta, &
        ! wind_period, wind_t_offset, 
        DumpWind, wind_mag_time_series_file
  
    open(unit=8, file="parameters.in", status='OLD', recl=80)
  
negative.code:   ! Zero vector - for internal use only
    zeros = 0d0
  
  
    ! Read in arrays from the input files
    call read_input_fileU(initUfile, u, 0.d0, nx, ny, layers)
  