anchor.diff:          ! save the old U and V values so we can restore them
           allocate(savedU(nxu,nz,ny))
           allocate(savedV(nx,nz,nyv))
  +        
           savedU=domain%u
           savedV=domain%v
           
  
positive.code:         endif
          if (reverse) then
              write(*,*) "ERROR: reversing linear winds not set up for parallel fftw computation yet"
              stop
          endif
  
  
negative.code:     real, parameter :: nsqmax=log(6e-4)
      real, parameter :: nsqmin=log(1e-7)
      
      integer, parameter :: n_dir_values=25
      integer, parameter :: n_nsq_values=6
      integer, parameter :: n_spd_values=6
      
      ! real,parameter::pi=3.1415927
      complex,parameter :: j= (0,1)
  