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) * linear_contribution
  +                        domain%v(i,j,k) = domain%v(i,j,k) + v_perturbation(i,j,k)
                       endif
                   end do
               end do
  
positive.code:                 
                  linear_contribution = saved_linear_contribution
                  linear_mask = linear_mask * linear_contribution
              ! else
              !     write(*,*) "Skipping spatial wind field for presumed domain repeat"
              endif
  
negative.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))    &
  