anchor.diff:  
     call read_input_fileU(zonalWindFile, base_wind_x, 0.d0, nx, ny, 1)
     call read_input_fileV(meridionalWindFile, base_wind_y, 0.d0, nx, ny, 1)
  -  
  +
     allocate(wind_mag_time_series(nTimeSteps))
     call read_input_file_time_series(wind_mag_time_series_file, &
  -   wind_mag_time_series, 1.d0, nTimeSteps)
  +      wind_mag_time_series, 1.d0, nTimeSteps)
   
     wind_x = base_wind_x*wind_mag_time_series(1)
     wind_y = base_wind_y*wind_mag_time_series(1)
  
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)
    read(unit=8, nml=NUMERICS)
  
negative.code:   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
  
  