-
Notifications
You must be signed in to change notification settings - Fork 1
/
NumberKinds.f90
109 lines (95 loc) · 5.26 KB
/
NumberKinds.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
module NumberKindsModule
!------------------------------------------------------------------------------
! Lagrangian Particle / Panel Method - Spherical Model
!------------------------------------------------------------------------------
!
!> @author
!> Peter Bosler, Sandia National Laboratories Center for Computing Research
!
!> @defgroup NumberKinds Number kinds module
!> Defines global numerical constants for use by all LPPM modules and executables.
!>
!
!
! DESCRIPTION:
!> @file
!> Defines global numerical constants for use by all LPPM modules and executables.
!
!------------------------------------------------------------------------------
implicit none
public
!> @defgroup TypeConstants Type Constants
!> @brief sets compiler-dependent kind variables for primitive data types
integer, parameter :: KREAL = kind(0.d0) !<< compiler generated kind for double precision reals @ingroup TypeConstants
integer, parameter :: KINT = kind(1) !<< compiler generated kind for default integer @ingroup TypeConstants
integer, parameter :: KLOG = kind(.TRUE.) !<< compiler generated kind logical @ingroup TypeConstants
real(KREAL), parameter :: ZERO_TOL = 1.0d-10 !<< zero tolerance for real numbers @ingroup TypeConstants
!> @defgroup PhysicalConstants Physical Constants
!> @brief Physical constants referenced from @cite CAM3.0
real(KREAL), parameter :: PI = 3.1415926535897932384626433832795027975_KREAL !<< Pi @ingroup PhysicalConstants
real(KREAL), parameter :: RAD_2_DEG = 180.0_kreal/PI !<< convert radians to degrees @ingroup PhysicalConstants
real(KREAL), parameter :: DEG_2_RAD = PI / 180.0_kreal !<< convert degress to radians @ingroup PhysicalConstants
real(KREAL), parameter :: GRAV = 9.80616 !<< acceleration due to gravity {m s^(-2)} @ingroup PhysicalConstants
real(KREAL), parameter :: ONE_DAY = 86140.0_kreal !<< Earth's sidereal day {s} @ingroup PhysicalConstants
real(KREAL), parameter :: EARTH_RADIUS = 6371220_kreal !<< mean radius of the Earth { m } @ingroup PhysicalConstants
! real(KREAL), save :: OMEGA = 2.0_kreal*PI / ONE_DAY !<< rotation rate of sphere {s^(-1)} @ingroup PhysicalConstants
real(KREAL), parameter :: EARTH_SURFACE_AREA = 4.0_kreal * PI * EARTH_RADIUS * EARTH_RADIUS
!> @defgroup IOConstants I/O Constants
!> @brief constants for file input and output
integer(KINT), parameter :: STD_ERR = 0
integer(KINT), parameter :: STD_IN = 5
integer(KINT), parameter :: STD_OUT = 6
integer(KINT), parameter :: READ_UNIT = 11
integer(KINT), parameter :: WRITE_UNIT_1 = 12
integer(KINT), parameter :: WRITE_UNIT_2 = 13
integer(KINT), parameter :: WRITE_UNIT_3 = 14
integer(KINT), parameter :: MAX_STRING_LENGTH = 256
!> @defgroup MeshConstants Mesh and application constants
!> @brief constants for mesh and problem definition
integer(KINT), parameter :: QUAD_PANEL = 4 !<< panelKind parameter for quadrilateral panels @ingroup MeshConstants
integer(KINT), parameter :: TRI_PANEL = 3 !<< panelKind parameter for triangular panels @ingroup MeshConstants
integer(KINT), parameter :: PLANAR_GEOM = 81
integer(KINT), parameter :: SPHERE_GEOM = 82
integer(KINT), parameter :: EUCLIDEAN_3D = 83
integer(KINT), parameter :: MAX_VERTEX_DEGREE = 10 !<< panelKind parameter for triangular panels @ingroup MeshConstants
integer(KINT), parameter :: ADVECTION_SOLVER = 90 !<< problemKind parameter for spherical advection @ingroup MeshConstants
integer(KINT), parameter :: BVE_SOLVER = 91 !<< problemKind parameter for quadrilateral panels @ingroup MeshConstants
integer(KINT), parameter :: SWE_SOLVER = 92 !<< problemKind parameter for quadrilateral panels @ingroup MeshConstants
integer(KINT), parameter :: PLANE_SOLVER = 99 !<< problemKind parameter for quadrilateral panels @ingroup MeshConstants
integer(KINT), parameter :: MAX_POLYGON_SIDES = 20 !<< defines maximum allowable polygon size in a Voronoi mesh @ingroup MeshConstants
integer(KINT), parameter :: FREE_BOUNDARIES = 101 !<< boundary condition parameter for planar meshes @ingroup MeshConstants
integer(KINT), parameter :: PERIODIC_BOUNDARIES = 102 !<< boundary condition parameter for planar meshes @ingroup MeshConstants
integer(KINT), parameter :: TRI_HEX_SEED = 201
integer(KINT), parameter :: QUAD_RECT_SEED = 202
integer(KINT), parameter :: POLAR_DISC_SEED = 203
integer(KINT), parameter :: QUAD_RECT_PERIODIC_SEED = 204
integer(KINT), parameter :: ICOS_TRI_SPHERE_SEED = 205
integer(KINT), parameter :: CUBED_SPHERE_SEED = 206
real(KREAL), save :: SphereRadius = 1.0_kreal
!> @defgroup MPIVariables MPI Variables
!> @brief variables to define MPI computing environment
integer(KINT), save :: numProcs = 1 !<< number of processess @ingroup MPIVariables
integer(KINT), save :: procRank = 0 !<< process rank @ingroup MPIVariables
interface
real(8) function scalarFnOf2DSpace( x, y)
real(8), intent(in) :: x, y
end function
end interface
interface
real(8) function scalarFnOf3DSpace( x, y, z)
real(8), intent(in) :: x, y, z
end function
end interface
interface
function vectorFnOf2DSpace( x, y )
real(8), dimension(2) :: vectorFnOf2DSpace
real(8), intent(in) :: x, y
end function
end interface
interface
function vectorFnOf3DSpace( x, y, z)
real(8), dimension(3) :: vectorFnOf3DSpace
real(8), intent(in) :: x, y, z
end function
end interface
end module NumberKindsModule