anchor.diff:      integer :: n_spd_values=10
   
       complex,parameter :: j= (0,1)
  +    
  +    real, parameter :: SMALL_VALUE = 1e-15
   
   contains
   
  
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: 
      logical :: variable_N
      logical :: smooth_nsq
      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
  