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:         lt_data%u_accumulator = 0
          lt_data%v_accumulator = 0
  
          current_z = step_size / 2
          do i=1,n_steps
              call linear_perturbation_at_height(U,V,Nsq,current_z, fourier_terrain, lt_data)
              lt_data%u_accumulator = lt_data%u_accumulator + lt_data%u_perturb
  
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
  