anchor.diff:    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)
  
positive.code:   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
    do n = 1, nTimeSteps
  
      ! ! Time varying winds
      ! if (UseSinusoidWind .eqv. .true.) then
      !   wind_x = base_wind_x*(wind_alpha +  &
      !       wind_beta*sin(((2d0*pi*n*dt)/wind_period) - wind_t_offset))
      !   wind_y = base_wind_y*(wind_alpha +  &
      !       wind_beta*sin(((2d0*pi*n*dt)/wind_period) - wind_t_offset))
      ! else if (UseStochWind .eqv. .true.) then
      !   if (mod(n-1, n_stoch_wind).eq.0) then
      !     ! Gives a pseudorandom number in range 0 <= x < 1
      !     call random_number(stoch_wind_mag)
      !     ! Convert to -1 <= x < 1
      !     stoch_wind_mag = (stoch_wind_mag - 0.5d0)*2d0
      !     wind_x = base_wind_x*(wind_alpha +  &
      !         wind_beta*stoch_wind_mag)
      !     wind_y = base_wind_y*(wind_alpha +  &
      !         wind_beta*stoch_wind_mag)
      !   end if
      ! end if
  
      wind_x = base_wind_x*wind_mag_time_series(n)
      wind_y = base_wind_y*wind_mag_time_series(n)
  
negative.code:     v(:, :, k) = v(:, :, k) * hfacS * wetmask(:, :)
    end do
  
  
    if (.not. RedGrav) then
      ! Initialise arrays for pressure solver
  