anchor.diff:      use io_routines,               only : io_read, io_write
       use string,                    only : str
       use linear_theory_lut_disk_io, only : read_LUT, write_LUT
  -    use mod_atm_utilities
  +    use mod_atm_utilities,         only : calc_stability, calc_u, calc_v, calc_speed, calc_direction, blocking_fraction
       use array_utilities,           only : smooth_array, calc_weight, linear_space
   
       implicit none
  
positive.code:                 ! $omp end parallel
              endif
          else
                  domain%nsquared = options%lt_options%N_squared
          endif
  
          ! smooth array has it's own parallelization, so this probably can't go in a critical section
  
negative.code:         nz  = size(domain%u,2)
          nxu = size(domain%u,1)
          nyv = size(domain%v,3)
  
          ! Here we just set externalN based on two options.
          ! if NfromForcing is True and nvar is defined then N is to be
          ! read from the forcing data set instead of being calculated by ICAR,
  