anchor.diff:          !$omp end do
   
           ! finally deallocate all temporary arrays that were created... chould be a datastructure and a subroutine...
  -!       if (.not.save_intermediate) then
           deallocate(k,l,kl,sig,denom,m,ineta,msq,mimag)
   
           !$omp end parallel
  
positive.code: 
      end function calc_direction
  
      ! calculate the strength of the u wind field given a direction [0-2*pi] and magnitude
      pure function calc_speed(u, v) result(speed)
          implicit none
          real, intent(in) :: u,v
  
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
  