anchor.diff:  
       logical :: variable_N
       logical :: smooth_nsq
  +    logical :: using_blocked_flow
       real    :: N_squared
   
       !! unfortunately these have to be allocated every call because we could be calling on both the high-res
  
positive.code:         use_spatial_linear_fields = options%lt_options%spatial_linear_fields    ! use a spatially varying linear wind perturbation
          use_linear_mask           = options%lt_options%linear_mask              ! use a spatial mask for the linear wind field
          use_nsq_calibration       = options%lt_options%nsq_calibration          ! use a spatial mask to calibrate the nsquared (brunt vaisala frequency) field
  
          ! Look up table generation parameters, range for each parameter, and number of steps to cover that range
          dirmax = options%lt_options%dirmax
  
negative.code:             do j=1, nz
                  do i=1, nxu
  
                      if (using_blocked_flow) then
                          blocked = blocking_fraction(domain%froude(min(i,nx),min(k,ny)))
                      else
                          blocked = 0
                      endif
                      if (blocked<1) then
  
                          uk = min(k,ny)
  