anchor.diff:  
     return
   end subroutine wrap_fields_2D
  -
  -!-----------------------------------------------------------------
  -!> Robustly produce a different random seed for each run.
  -!! Code adapted from
  -!! http://web.ph.surrey.ac.uk/fortweb/glossary/random_seed.html
  -
  -subroutine ranseed()
  -  implicit none
  -
  -  ! ----- variables for portable seed setting -----
  -  integer :: i_seed
  -  integer, dimension(:), allocatable :: a_seed
  -  integer, dimension(1:8) :: dt_seed
  -  ! ----- end of variables for seed setting -----
  -
  -  ! ----- Set up random seed portably -----
  -  call random_seed(size=i_seed)
  -  allocate(a_seed(1:i_seed))
  -  call random_seed(get=a_seed)
  -  call date_and_time(values=dt_seed)
  -  a_seed(i_seed) = dt_seed(8); a_seed(1) = dt_seed(8)*dt_seed(7)*dt_seed(6)
  -  call random_seed(put=a_seed)
  -  return
  -end subroutine ranseed
  
positive.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
  
  
negative.code:   namelist /INITIAL_CONDITONS/ initUfile, initVfile, initHfile, initEtaFile
  
    namelist /EXTERNAL_FORCING/ zonalWindFile, meridionalWindFile, &
        DumpWind, wind_mag_time_series_file
  
    open(unit=8, file="parameters.in", status='OLD', recl=80)
  