anchor.diff:    double precision :: wind_y(0:nx+1, 0:ny+1)
     double precision :: base_wind_x(0:nx+1, 0:ny+1)
     double precision :: base_wind_y(0:nx+1, 0:ny+1)
  -  ! logical :: UseSinusoidWind
  -  ! logical :: UseStochWind
     logical :: DumpWind
  -  ! double precision :: wind_alpha, wind_beta, wind_period, wind_t_offset
  -  ! integer :: n_stoch_wind
  -  ! double precision :: stoch_wind_mag
     character(30) :: wind_mag_time_series_file
     double precision, dimension(:), allocatable :: wind_mag_time_series
   
  
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:   ! 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)
  