anchor.diff:          integer :: buffer                   ! number of grid cells to buffer around the domain MUST be >=1
           integer :: stability_window_size    ! window to average nsq over
           real :: max_stability               ! limits on the calculated Brunt Vaisala Frequency
  -        real :: min_stability               ! these may need to be a little narrower. 
  +        real :: min_stability               ! these may need to be a little narrower.
           real :: linear_contribution         ! multiplier on uhat,vhat before adding to u,v
           real :: linear_update_fraction      ! controls the rate at which the linearfield updates (should be calculated as f(in_dt))
  -    
  +
           real :: N_squared                   ! static Brunt Vaisala Frequency (N^2) to use
           logical :: remove_lowres_linear     ! attempt to remove the linear mountain wave from the forcing low res model
           real :: rm_N_squared                ! static Brunt Vaisala Frequency (N^2) to use in removing linear wind field
           real :: rm_linear_contribution      ! fractional contribution of linear perturbation to wind field to remove from the low-res field
  -    
  +
           logical :: spatial_linear_fields    ! use a spatially varying linear wind perturbation
           logical :: linear_mask              ! use a spatial mask for the linear wind field
           logical :: nsq_calibration          ! use a spatial mask to calibrate the nsquared (brunt vaisala frequency) field
  -    
  +
           ! Look up table generation parameters
           real :: dirmax, dirmin
           real :: spdmax, spdmin
  
positive.code:                 write(*,*) "WARNING WARNING WARNING"
              endif
          endif
      
      end subroutine options_check
      
      subroutine init_restart_options(filename, options)
          ! initialize the restart specifications
          ! read in the namelist, and calculate the restart_step if appropriate
          character(len=*),  intent(in)    :: filename    ! name of the file containing the restart namelist
          type(options_type),intent(inout) :: options     ! options data structure to store output
          
          ! Local variables
          character(len=MAXFILELENGTH) :: restart_file    ! file name to read restart data from
          integer :: restart_step                         ! time step relative to the start of the restart file
  
negative.code:                 warning_level=5
              endif
          endif
  
  
          if (smooth_wind_distance.eq.(-9999)) then
              smooth_wind_distance=dxlow*2
              write(*,*) " Default smoothing distance = lowdx*2 = ", smooth_wind_distance
          endif
  
          options%t_offset=t_offset
          if (smooth_wind_distance<0) then
              write(*,*) " Wind smoothing must be a positive number"
              write(*,*) " smooth_wind_distance = ",smooth_wind_distance
              if (warning_level>4) then
  