valib_mod.F90

 
 
!*******************************************************************************
! Edit to switch to single-precision version (double is the default):
! double precision version
#define VA_REAL_F c_double
! single precision version
!#define VA_REAL_F c_float
!*******************************************************************************
 
module valib_mod
   use iso_c_binding
   implicit none
   interface
      subroutine valib_contrarotation(A, S_RAVG) &
         bind(c, name='valib_contrarotation')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: S_RAVG(9)
      end subroutine valib_contrarotation
 
      subroutine valib_corotation(A, omega_RAVG) &
         bind(c, name='valib_corotation')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: omega_RAVG(3)
      end subroutine valib_corotation
 
      subroutine valib_delta(A, delta) &
         bind(c, name='valib_delta')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: delta
      end subroutine valib_delta
 
      subroutine valib_lambda2(A, lambda_2) &
         bind(c, name='valib_lambda2')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: lambda_2
      end subroutine valib_lambda2
 
      subroutine valib_qcriterion(A, qcriterion) &
         bind(c, name='valib_qcriterion')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: qcriterion
      end subroutine valib_qcriterion
 
      subroutine valib_qdcriterion(A, qdcriterion) &
         bind(c, name='valib_qdcriterion')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: qdcriterion
      end subroutine valib_qdcriterion
 
      subroutine valib_qmcriterion(A, qmcriterion) &
         bind(c, name='valib_qmcriterion')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: qmcriterion
      end subroutine valib_qmcriterion
 
      subroutine valib_qwcriterion(A, qwcriterion) &
         bind(c, name='valib_qwcriterion')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: qwcriterion
      end subroutine valib_qwcriterion
 
      subroutine valib_rortex(A, rortex) &
         bind(c, name='valib_rortex')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: rortex(3)
      end subroutine valib_rortex
 
      subroutine valib_swirling_strength(A, lambda_ci, lambda_cr) &
         bind(c, name='valib_swirling_strength')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         real(VA_REAL_F), intent(out) :: lambda_ci
         real(VA_REAL_F), intent(out) :: lambda_cr
      end subroutine valib_swirling_strength
 
      subroutine valib_triple_decomposition(A, num_intervals, &
                                            residual_vorticity, &
                                            residual_strain, &
                                            shear) &
         bind(c, name='valib_triple_decomposition')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)    :: A(9)
         integer(c_int), intent(in) :: num_intervals
         real(VA_REAL_F), intent(out) :: residual_vorticity
         real(VA_REAL_F), intent(out) :: residual_strain
         real(VA_REAL_F), intent(out) :: shear
      end subroutine valib_triple_decomposition
 
      subroutine valib_triple_decomposition_4norms(A, num_intervals, &
                                                   residual_vorticity, &
                                                   residual_strain, &
                                                   shear_vorticity, &
                                                   shear_strain) &
         bind(c, name='valib_triple_decomposition_4norms')
         use iso_c_binding
         implicit none
         real(VA_REAL_F), intent(in)  :: A(9)
         integer(c_int), intent(in) :: num_intervals
         real(VA_REAL_F), intent(out) :: residual_vorticity
         real(VA_REAL_F), intent(out) :: residual_strain
         real(VA_REAL_F), intent(out) :: shear_vorticity
         real(VA_REAL_F), intent(out) :: shear_strain
      end subroutine valib_triple_decomposition_4norms
   end interface
 
end module valib_mod