anchor.diff:  
       logical :: variable_N
       logical :: smooth_nsq
  -    real :: N_squared
  +    real    :: N_squared
       !! unfortunately these have to be allocated every call because we could be calling on both the high-res
       !! domain and the low res domain (to "remove" the linear winds)
       real,allocatable,dimension(:,:)::k,l,kl,sig
  
positive.code:         real :: BV_freq
  
          BV_freq = gravity * (log(th_top)-log(th_bot)) / (z_top - z_bot)
  
      end function calc_dry_stability
  
      !>----------------------------------------------------------
  
negative.code:     integer :: n_spd_values=10
  
      complex,parameter :: j= (0,1)
      
      real, parameter :: SMALL_VALUE = 1e-15
  
  contains
  
  