anchor.diff:          original_buffer = options%lt_options%buffer
           variable_N = options%lt_options%variable_N
           smooth_nsq = options%lt_options%smooth_nsq
  -        linear_contribution = options%lt_options%linear_contribution
   
           stability_window_size = options%lt_options%stability_window_size        ! window to average nsq over
           max_stability = options%lt_options%max_stability                        ! limits on the calculated Brunt Vaisala Frequency
  
positive.code:         !$omp private(spos, dpos, npos, nexts,nextd, nextn,n, smoothz), &
          !$omp private(wind_first, wind_second, curspd, curdir, curnsq, sweight,dweight, nweight), &
          !$omp shared(domain, spd_values, dir_values, nsq_values, u_LUT, v_LUT), &
          !$omp shared(u_perturbation, v_perturbation, linear_update_fraction, nsq_calibration), &
          !$omp shared(min_stability, max_stability, n_dir_values, n_spd_values, n_nsq_values, smooth_nsq)
          !$omp do
          do k=1,ny
  
negative.code:         complex(C_DOUBLE_COMPLEX),allocatable,dimension(:,:)::complex_terrain
          type(C_PTR) :: plan
          integer::nx,ny,nz, save_buffer
          real :: saved_linear_contribution ! temporary variable to save the linear contribution state so we can restore
  
          ! store module level variables so we don't have to pass options through everytime
          ! lots of these things probably need to be moved to the linearizable class so they
  