anchor.diff:          debug=.False.
           
       end subroutine linear_perturb
  -end module linear_theory_winds+end module linear_theory_winds
  
positive.code:         ! 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
          
  
negative.code:         endif
          if (reverse) then
              write(*,*) "ERROR: reversing linear winds not set up for parallel fftw computation yet"
              write(*,*) "also not set up for the new spatial wind LUTs"
              stop
          endif
  
  