quadric_mod.f

Go to the documentation of this file.

      MODULE quadric

!     Maximum of the number of quadrics that can be read
      INTEGER, PARAMETER          :: dim_quadric = 500
!     Nnumber of quadrics
      INTEGER                     :: n_quadric
!     Current Quadric
      INTEGER :: quadric_id
!     form of quadric : 'normal' or one of the pre-defined quadrics
      CHARACTER (LEN=12), DIMENSION(DIM_QUADRIC) :: quadric_form
!     Scale factor for quadrics
      DOUBLE PRECISION :: quadric_scale
!     Characteristic values of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: lambda_x,lambda_y,lambda_z
!     d - coefficient of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: dquadric
!     Translation components of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: t_x,t_y,t_z
!     Rotation angles (Deg) of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: theta_x,theta_y,theta_z
!     Radius for either Spere or Cylinder (pre-defined quadrics)
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: radius
!     Radii for Torus
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: torus_r1, torus_r2
!     Radii for U-coil
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: ucoil_r1, ucoil_r2
      !     Y-location of bends for U-coil
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: ucoil_y1, ucoil_y2
!     Radii for Bend
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: bend_r1, bend_r2
!     Angles for Bend
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: bend_theta1, bend_theta2
!     Y-locations of cylinder-cone-cylinder
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: c2c_y1, c2c_y2
!     Radii of cylinder-cone-cylinder
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: c2c_r1, c2c_r2
!     Half-angle for cone (pre-defined quadrics)
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: half_angle
!     Reactor 1 lower, upper cylinder radii
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: reactor1_r1,reactor1_r2
!     Reactor 1 lower, upper conical transition between cylinders
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: reactor1_y1,reactor1_y2
!     Reactor 1 lower, upper rounding locations
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: reactor1_yr1,reactor1_yr2
!     Reactor 1 lower, upper rounding radii
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: reactor1_rr1,reactor1_rr2
!     Reactor 1 lower, upper rounding angles
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: reactor1_theta1,reactor1_theta2
!     Normal vector components for plane (pre-defined quadrics)
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: n_x,n_y,n_z
!     A-matrices of the quadrics
      DOUBLE PRECISION, DIMENSION(3,3,DIM_QUADRIC) :: a_quadric
!     Translation-matrices of the quadrics
      DOUBLE PRECISION, DIMENSION(1,3,DIM_QUADRIC) :: t_quadric
!     Clipping range  of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: clip_xmin,clip_xmax,clip_ymin,clip_ymax,clip_zmin,clip_zmax
!     Piecewise range  of the quadrics
      DOUBLE PRECISION, DIMENSION(DIM_QUADRIC) :: piece_xmin,piece_xmax,piece_ymin,piece_ymax,piece_zmin,piece_zmax
!     Clip flag
      LOGICAL, DIMENSION(DIM_QUADRIC) :: fluid_in_clipped_region
!     Boundary condition ID
      INTEGER :: bc_id_q(dim_quadric)
!     Maximum number of groups
      INTEGER, PARAMETER :: dim_group = 50
!     Number of groups
      INTEGER :: n_group
!     Number of quadric in each group
      INTEGER,DIMENSION(DIM_GROUP) :: group_size
!     Quadric ID list in each group
      INTEGER,DIMENSION(DIM_GROUP,DIM_QUADRIC) :: group_q
!     Quadric relation in each group
      CHARACTER(LEN=9),DIMENSION(DIM_GROUP) :: group_relation
!     Relation between groups
      CHARACTER(LEN=9),DIMENSION(DIM_GROUP) :: relation_with_previous
!     Tolerance for intersection between quadrics and planes
      DOUBLE PRECISION :: tol_f
!     Maximum number of iterations while finding intersection between geometry and grid
      INTEGER :: itermax_int

    CONTAINS

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Function name: CROSS_PRODUCT                                        C
!  Purpose: Performs the cross product between two vectors             C
!           C = A x B                                                  C
!                                                                      C
!  Author: Jeff Dietiker                              Date: 21-Feb-08  C
!  Reviewer:                                          Date:            C
!                                                                      C
!  Revision Number #                                  Date: ##-###-##  C
!  Author: #                                                           C
!  Purpose: #                                                          C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      FUNCTION cross_product(A,B)

        IMPLICIT NONE
        DOUBLE PRECISION, DIMENSION(3) :: CROSS_PRODUCT
        DOUBLE PRECISION, INTENT(IN), DIMENSION(3) :: A,B

        cross_product(1) = a(2) * b(3) - a(3) * b(2)
        cross_product(2) = a(3) * b(1) - a(1) * b(3)
        cross_product(3) = a(1) * b(2) - a(2) * b(1)

      END FUNCTION cross_product

    END MODULE quadric