anchor.diff:  
       end function calc_stability
   
  -    ! calculate a single brunt vaisala frequency effectively averaged across the entire domain
  +    !>----------------------------------------------------------
  +    !! Calculate a single brunt vaisala frequency effectively averaged across the entire domain
  +    !!
  +    !!----------------------------------------------------------
       pure function calc_domain_stability(domain) result(BV_freq)
           implicit none
           class(linearizable_type),intent(in)::domain
  
positive.code: 
                          debug=.False. ! after the first time through set debug to false
  
                          u_LUT(k,i,j,:,:,:)=(domain%u-calc_u(dir_values(i),spd_values(k)))
                          v_LUT(k,i,j,:,:,:)=(domain%v-calc_v(dir_values(i),spd_values(k)))
                      end do
                  end do
              end do
              write(*,*) char(10),"--------  Linear wind look up table generation complete ---------"
              domain%u=savedU
              domain%v=savedV
          end if
  
          deallocate(savedU,savedV)
  
negative.code: 
          !$omp do
          do z=1,nz
              U = U_layers(z)
              V = V_layers(z)
              if ((abs(U)+abs(V))>0.5) then
                  sig  = U*k+V*l
                  where(sig==0.0) sig=1e-15
  