anchor.diff:                          v_perturbation(i,j,k)=v_perturbation(i,j,k) * (1-linear_update_fraction) &
                                       + linear_update_fraction * (sweight*wind_first + (1-sweight)*wind_second)
   
  -                        domain%v(i,j,k) = domain%v(i,j,k) + v_perturbation(i,j,k)
  +                        domain%v(i,j,k) = domain%v(i,j,k) + v_perturbation(i,j,k) * linear_contribution
                       endif
                   end do
               end do
  
positive.code:                         u_perturbation(i,j,k)=u_perturbation(i,j,k) * (1-linear_update_fraction) &
                                      + linear_update_fraction * (sweight*wind_first + (1-sweight)*wind_second)
  
                          domain%u(i,j,k) = domain%u(i,j,k) + u_perturbation(i,j,k)
                      endif
                      if (i<=nx) then
                          wind_first =      nweight  * (dweight * v_LUT(spos,dpos,npos,i,j,k)  + (1-dweight) * v_LUT(spos,nextd,npos,i,j,k))    &
  
negative.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, linear_contribution, 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
  