anchor.diff:          ! default assumes no errors in reading the LUT
           error = 0
   
  -        print*, "WARNING: minimum_layer_size hard coded at 100 m"
  -        ! minimum_layer_size = optione%minimum_layer_size
  +        ! minimum_layer_size = options%minimum_layer_size
           minimum_layer_size = 100
  +        print*, "WARNING: minimum_layer_size hard coded as", minimum_layer_size
   
           ! store to make it easy to check dim sizes in read_LUT
           LUT_dims(:,1) = [nxu,nz,ny]
  
positive.code:             loops_completed = 0
              write(*,*) "Percent Completed:"
              !$omp parallel default(shared) &
              !$omp private(i,j,k,z, u,v, layer_height, layer_height_bottom, layer_height_top, minimum_layer_size)
              ! $omp private(lt_data_m)
  
              ! initialization has to happen in each thread so each thread has its own copy
  
negative.code:         if (reverse) then
              !$omp critical (print_lock)
              write(*,*) "ERROR: reversing linear winds not set up for parallel fftw computation yet"
              write(*,*) "Use only with spatially varying wind LUTs"
              !$omp end critical (print_lock)
              stop
          endif
  