anchor.diff:                                                                  +domain%qrain(i,j,k)+domain%qsnow(i,j,k))
                           domain%nsquared(i,j,k) = max(min(domain%nsquared(i,j,k) * nsq_calibration(i,k), max_stability), min_stability)
                       else
  -                        ! Low-res boundary condition variables will be in a different array format.  It should be 
  +                        ! Low-res boundary condition variables will be in a different array format.  It should be
                           ! easy enough to call calc_stability after e.g. transposing z and y dimension, but some
  -                        ! e.g. pii will not be set in the forcing data, so this may need a little thought. 
  +                        ! e.g. pii will not be set in the forcing data, so this may need a little thought.
                           domain%nsquared(i,j,k) = 3e-6
                       endif
                   end do
  
positive.code:                             npos = step
                          endif
                      end do
                      
                      ! Calculate the weights and the "next" u/v position
                      ! "next" usually = pos+1 but for edge cases next = 1 or n
                      dweight = calc_weight(dir_values, dpos, nextd, curdir)
                      sweight = calc_weight(spd_values, spos, nexts, curspd)
                      nweight = calc_weight(nsq_values, npos, nextn, curnsq)
                      
                      ! perform linear interpolation between LUT values
                      if (k<=ny) then
                          wind_first =      nweight  * (dweight * u_LUT(spos, dpos,npos, i,j,k) + (1-dweight) * u_LUT(spos, nextd,npos, i,j,k))   &
  
negative.code:         else
              rev=.False.
          endif
  
          if (present(useDensity)) then
              useD=useDensity
          else
              useD=.False.
          endif
  
          ! this is a little trickier, because it does have to be domain dependant... could at least be stored in the domain though...
          if (rev) then
              linear_contribution = options%lt_options%rm_linear_contribution
  