Define some much-used constants for global use (all defined as PARAMETER). More...

Functions/Subroutines
subroutine	tabu_fw
	Estimate friction coefficients in oscillatory boundary layers using tabulation on Kelvin functions. More...

subroutine	kzeone (X, Y, RE0, IM0, RE1, IM1)
	June 1999 adaptation to CRESTb, all tests on range of (x,y) have been bypassed, we implicitly expect X to be positive or \|x,y\| non zero. More...

subroutine	kerkei (X, KER, KEI)
	Computes the values of the zeroth order Kelvin function Ker and Kei. More...

Variables
logical, parameter	tstout = .FALSE.
	TSTOUT Flag for generation of test files. More...

real, parameter	grav = 9.806
	GRAV Acc. More...

real, parameter	dwat = 1000.
	DWAT Density of water (kg/m3). More...

real, parameter	dair = 1.225
	DAIR Density of air (kg/m3). More...

real, parameter	nu_air = 1.4E-5
	NU_AIR Kinematic viscosity of air (m2/s). More...

real, parameter	nu_water = 1.31E-6
	NU_WATER Kinematic viscosity of water (m2/s). More...

real, parameter	sed_sg = 2.65
	SED_SG Specific gravity of sediments (N.D.). More...

real, parameter	kappa = 0.40
	KAPPA von Karman's constant (N.D.). More...

real, parameter	pi = 3.141592653589793
	PI Value of Pi. More...

real, parameter	tpi = 2.0 * PI
	TPI 2*Pi. More...

real, parameter	hpi = 0.5 * PI
	HPI 1/2*Pi. More...

real, parameter	tpiinv = 1. / TPI
	TPIINV Inverse of 2*Pi. More...

real, parameter	hpiinv = 1. / HPI
	HPIINV Inverse of 1/2*Pi. More...

real, parameter	rade = 180. / PI
	RADE Conversion factor from radians to degrees. More...

real, parameter	dera = PI / 180.
	DERA Conversion factor from degrees to radians. More...

real, parameter	radius = 4.E7 * TPIINV
	RADIUS Radius of the earth (m). More...

real, parameter	g2pi3i = 1. / ( GRAV*2 TPI**3 )
	G2PI3I Inverse of gravity^2 * (2*Pi)^3. More...

real, parameter	g1pi1i = 1. / ( GRAV * TPI )
	G1PI1I Inverse of gravity * 2 * Pi. More...

real	undef = -999.9
	UNDEF Value for undefined variable in output. More...

character(*), parameter	file_endian = ENDIANNESS
	FILE_ENDIAN Filled by preprocessor with 'big_endian', 'little_endian', or 'native'. More...

integer, parameter	sizefwtable =300
	SIZEFWTABLE. More...

real, dimension(0:sizefwtable)	fwtable
	FWTABLE. More...

real	delab
	DELAB. More...

real, parameter	abmin = -1.
	ABMIN. More...

integer, parameter	srce_direct = 0
	srce_direct More...

integer, parameter	srce_imp_post = 1
	srce_imp_post More...

integer, parameter	srce_imp_pre = 2
	srce_imp_pre More...

integer, parameter	debug_node = 1014
	DEBUG_NODE Node number used for debugging. More...

integer, parameter	debug_element = 50
	DEBUG_ELEMENT Element number used for debug. More...

logical	lpdlib = .FALSE.
	LPDLIB Logical for using the PDLIB library. More...

logical	lsetup = .FALSE.
	LSETUP Logical LSETUP is not used. More...

logical	is_esmf_component = .FALSE.
	IS_ESMF_COMPONENT Flag for model invoked via ESMF. More...

Detailed Description

Define some much-used constants for global use (all defined as PARAMETER).

Author: H. L. Tolman

Date: 05-Jun-2018

Function/Subroutine Documentation

◆ kerkei()

subroutine constants::kerkei	(	real	X,
		real	KER,
		real	KEI
	)

Computes the values of the zeroth order Kelvin function Ker and Kei.

These functions are used to determine the friction factor fw as a function of the bottom roughness length assuming a linear profile of eddy viscosity (See Grant and Madsen, 1979).

Parameters

X
KER
KEI

Author: N/A

Date: N/A

Definition at line 424 of file constants.F90.

     !**********************************************************************
     ! Computes the values of the zeroth order Kelvin function Ker and Kei
     ! These functions are used to determine the friction factor fw as a
     ! function of the bottom roughness length assuming a linear profile
     ! of eddy viscosity (See Grant and Madsen, 1979)
     !**********************************************************************
     IMPLICIT NONE
  
     DOUBLE PRECISION ZR,ZI,CYR,CYI,CYR1,CYI1
     REAL X,KER,KEI
  
     zr=x*.50d0*sqrt(2.0d0)
     zi=zr
     CALL kzeone(zr, zi, cyr, cyi,cyr1,cyi1)
     ker=cyr/exp(zr)
     kei=cyi/exp(zr)

References kzeone().

Referenced by tabu_fw().

◆ kzeone()

subroutine constants::kzeone	(	double precision	X,
		double precision	Y,
		double precision	RE0,
		double precision	IM0,
		double precision	RE1,
		double precision	IM1
	)

June 1999 adaptation to CRESTb, all tests on range of (x,y) have been bypassed, we implicitly expect X to be positive or |x,y| non zero.

The variables X and Y are the real and imaginary parts of the argument of the first two modified bessel functions of the second kind,k0 and k1. Re0,im0,re1 and im1 give the real and imaginary parts of exp(x)*k0 and exp(x)*k1, respectively. Although the real notation used in this subroutine may seem inelegant when compared with the complex notation that fortran allows, this version runs about 30 percent faster than one written using complex variables.

Copyright: This subroutine is copyright by ACM, see http://www.acm.org/pubs/copyright_policy/softwareCRnotice.html. ACM declines any responsibility of any kind.

Parameters

X	Real part of argument to modified Bessel functions.
Y	Imaginary part of argument to modified Bessel functions.
RE0	Real part of exp(x)*k0.
IM0	Imaginary part of exp(x)*k0.
RE1	Real part of exp(x)*k1.
IM1	Imaginary part of exp(x)*k1.

Author: N/A

Date: N/A

Definition at line 244 of file constants.F90.

     !  June 1999 adaptation to CRESTb, all tests on range of (x,y) have been
     !  bypassed, we implicitly expect X to be positive or |x,y| non zero
     !
     ! This subroutine is copyright by ACM
     ! see http://www.acm.org/pubs/copyright_policy/softwareCRnotice.html
     ! ACM declines any responsibility of any kind
     !
     ! THE VARIABLES X AND Y ARE THE REAL AND IMAGINARY PARTS OF
     ! THE ARGUMENT OF THE FIRST TWO MODIFIED BESSEL FUNCTIONS
     ! OF THE SECOND KIND,K0 AND K1.  RE0,IM0,RE1 AND IM1 GIVE
     ! THE REAL AND IMAGINARY PARTS OF EXP(X)*K0 AND EXP(X)*K1,
     ! RESPECTIVELY.  ALTHOUGH THE REAL NOTATION USED IN THIS
     ! SUBROUTINE MAY SEEM INELEGANT WHEN COMPARED WITH THE
     ! COMPLEX NOTATION THAT FORTRAN ALLOWS, THIS VERSION RUNS
     ! ABOUT 30 PERCENT FASTER THAN ONE WRITTEN USING COMPLEX
     ! VARIABLES.
     ! ACM Libraries
     ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     IMPLICIT NONE
     DOUBLE PRECISION X, Y, X2, Y2, RE0, IM0, RE1, IM1, &
          R1, R2, T1, T2, P1, P2, RTERM, ITERM, L
     DOUBLE PRECISION , PARAMETER, DIMENSION(8) :: EXSQ = &
          (/ 0.5641003087264d0,0.4120286874989d0,0.1584889157959d0, &
          0.3078003387255d-1,0.2778068842913d-2,0.1000044412325d-3, &
          0.1059115547711d-5,0.1522475804254d-8 /)
     DOUBLE PRECISION , PARAMETER, DIMENSION(8) :: TSQ = &
          (/ 0.0d0,3.19303633920635d-1,1.29075862295915d0, &
          2.95837445869665d0,5.40903159724444d0,8.80407957805676d0, &
          1.34685357432515d1,2.02499163658709d1 /)
     INTEGER N,M,K
     ! THE ARRAYS TSQ AND EXSQ CONTAIN THE SQUARE OF THE
     ! ABSCISSAS AND THE WEIGHT FACTORS USED IN THE GAUSS-
     ! HERMITE QUADRATURE.
     r2 = x*x + y*y
     IF (r2.GE.1.96d2) GO TO 50
     IF (r2.GE.1.849d1) GO TO 30
     ! THIS SECTION CALCULATES THE FUNCTIONS USING THE SERIES
     ! EXPANSIONS
     x2 = x/2.0d0
     y2 = y/2.0d0
     p1 = x2*x2
     p2 = y2*y2
     t1 = -(dlog(p1+p2)/2.0d0+0.5772156649015329d0)
     ! THE CONSTANT IN THE PRECEDING STATEMENT IS EULER*S
     ! CONSTANT
     t2 = -datan2(y,x)
     x2 = p1 - p2
     y2 = x*y2
     rterm = 1.0d0
     iterm = 0.0d0
     re0 = t1
     im0 = t2
     t1 = t1 + 0.5d0
     re1 = t1
     im1 = t2
     p2 = dsqrt(r2)
     l = 2.106d0*p2 + 4.4d0
     IF (p2.LT.8.0d-1) l = 2.129d0*p2 + 4.0d0
     DO n=1,int(l)
       p1 = n
       p2 = n*n
       r1 = rterm
       rterm = (r1*x2-iterm*y2)/p2
       iterm = (r1*y2+iterm*x2)/p2
       t1 = t1 + 0.5d0/p1
       re0 = re0 + t1*rterm - t2*iterm
       im0 = im0 + t1*iterm + t2*rterm
       p1 = p1 + 1.0d0
       t1 = t1 + 0.5d0/p1
       re1 = re1 + (t1*rterm-t2*iterm)/p1
       im1 = im1 + (t1*iterm+t2*rterm)/p1
     END DO
     r1 = x/r2 - 0.5d0*(x*re1-y*im1)
     r2 = -y/r2 - 0.5d0*(x*im1+y*re1)
     p1 = dexp(x)
     re0 = p1*re0
     im0 = p1*im0
     re1 = p1*r1
     im1 = p1*r2
     RETURN
     ! THIS SECTION CALCULATES THE FUNCTIONS USING THE INTEGRAL
     ! REPRESENTATION, EQN 3, EVALUATED WITH 15 POINT GAUSS-
     ! HERMITE QUADRATURE
 30  x2 = 2.0d0*x
     y2 = 2.0d0*y
     r1 = y2*y2
     p1 = dsqrt(x2*x2+r1)
     p2 = dsqrt(p1+x2)
     t1 = exsq(1)/(2.0d0*p1)
     re0 = t1*p2
     im0 = t1/p2
     re1 = 0.0d0
     im1 = 0.0d0
     DO n=2,8
       t2 = x2 + tsq(n)
       p1 = dsqrt(t2*t2+r1)
       p2 = dsqrt(p1+t2)
       t1 = exsq(n)/p1
       re0 = re0 + t1*p2
       im0 = im0 + t1/p2
       t1 = exsq(n)*tsq(n)
       re1 = re1 + t1*p2
       im1 = im1 + t1/p2
     END DO
     t2 = -y2*im0
     re1 = re1/r2
     r2 = y2*im1/r2
     rterm = 1.41421356237309d0*dcos(y)
     iterm = -1.41421356237309d0*dsin(y)
     ! THE CONSTANT IN THE PREVIOUS STATEMENTS IS,OF COURSE,
     ! SQRT(2.0).
     im0 = re0*iterm + t2*rterm
     re0 = re0*rterm - t2*iterm
     t1 = re1*rterm - r2*iterm
     t2 = re1*iterm + r2*rterm
     re1 = t1*x + t2*y
     im1 = -t1*y + t2*x
     RETURN
     ! THIS SECTION CALCULATES THE FUNCTIONS USING THE
     ! ASYMPTOTIC EXPANSIONS
 50  rterm = 1.0d0
     iterm = 0.0d0
     re0 = 1.0d0
     im0 = 0.0d0
     re1 = 1.0d0
     im1 = 0.0d0
     p1 = 8.0d0*r2
     p2 = dsqrt(r2)
     l = 3.91d0+8.12d1/p2
     r1 = 1.0d0
     r2 = 1.0d0
     m = -8
     k = 3
     DO n=1,int(l)
       m = m + 8
       k = k - m
       r1 = float(k-4)*r1
       r2 = float(k)*r2
       t1 = float(n)*p1
       t2 = rterm
       rterm = (t2*x+iterm*y)/t1
       iterm = (-t2*y+iterm*x)/t1
       re0 = re0 + r1*rterm
       im0 = im0 + r1*iterm
       re1 = re1 + r2*rterm
       im1 = im1 + r2*iterm
     END DO
     t1 = dsqrt(p2+x)
     t2 = -y/t1
     p1 = 8.86226925452758d-1/p2
     ! THIS CONSTANT IS SQRT(PI)/2.0, WITH PI=3.14159...
     rterm = p1*dcos(y)
     iterm = -p1*dsin(y)
     r1 = re0*rterm - im0*iterm
     r2 = re0*iterm + im0*rterm
     re0 = t1*r1 - t2*r2
     im0 = t1*r2 + t2*r1
     r1 = re1*rterm - im1*iterm
     r2 = re1*iterm + im1*rterm
     re1 = t1*r1 - t2*r2
     im1 = t1*r2 + t2*r1
     RETURN

Referenced by kerkei().

◆ tabu_fw()

subroutine constants::tabu_fw

Estimate friction coefficients in oscillatory boundary layers using tabulation on Kelvin functions.

Author: F. Ardhuin

Date: 28-Feb-2013

Definition at line 120 of file constants.F90.

     !/
     !/                  +-----------------------------------+
     !/                  | WAVEWATCH III           NOAA/NCEP |
     !/                  |            F. Ardhuin             |
     !/                  |                        FORTRAN 90 |
     !/                  | Last update :         28-Feb-2013 |
     !/                  +-----------------------------------+
     !/
     !/    19-Oct-2007 : Origination.                        ( version 3.13 )
     !/    28-Feb-2013 : Caps the friction factor to 0.5     ( version 4.08 )
     !/
     !  1. Purpose :
     !     TO estimate friction coefficients in oscillatory boundary layers
     !     METHOD.
     !      tabulation on Kelvin functions
     !
     !  2. Method :
     !
     !  3. Parameters :
     !
     !     Parameter list
     !     ----------------------------------------------------------------
     !     ----------------------------------------------------------------
     !
     !  4. Subroutines used :
     !
     !      Name      Type  Module   Description
     !     ----------------------------------------------------------------
     !      STRACE    Subr. W3SERVMD Subroutine tracing.
     !     ----------------------------------------------------------------
     !
     !  5. Called by :
     !
     !      Name      Type  Module   Description
     !     ----------------------------------------------------------------
     !      WW3_GRID  Prog. WW3_GRID Model grid initialization
     !     ----------------------------------------------------------------
     !
     !  6. Error messages :
     !
     !       None.
     !
     !  7. Remarks :
     !
     !  8. Structure :
     !
     !     See source code.
     !
     !  9. Switches :
     !
     !     !/S  Enable subroutine tracing.
     !
     ! 10. Source code :
     !
     !/ ------------------------------------------------------------------- /
     IMPLICIT NONE
     INTEGER, PARAMETER      :: NITER=100
     REAL   , PARAMETER      :: XM=0.50
     REAL   , PARAMETER      :: EPS1=0.00001
     ! ----------------------------------------------------------------------
     INTEGER I,ITER
     REAL KER, KEI
     REAL ABR,ABRLOG,L10,FACT,FSUBW,FSUBWMEMO,dzeta0,dzeta0memo
     !
     delab   = (abmax-abmin)/real(sizefwtable)
     l10=alog(10.)
     DO i=0,sizefwtable
       !
       !  index I in this table corresponds to a normalized roughness z0/ABR = 10^ABMIN+REAL(I)*DELAB
       !
       abrlog=abmin+real(i)*delab
       abr=exp(abrlog*l10)
       fact=1/abr/(21.2*kappa)
       fsubw=0.05
       dzeta0=0.
       DO iter=1,niter
         fsubwmemo=fsubw
         dzeta0memo=dzeta0
         dzeta0=fact*fsubw**(-.5)
         CALL kerkei(2.*sqrt(dzeta0),ker,kei)
         fsubw=.08/(ker**2+kei**2)
         fsubw=.5*(fsubwmemo+fsubw)
         dzeta0=.5*(dzeta0memo+dzeta0)
       END DO
       !
       ! Maximum value of 0.5 for fe is based on field
       ! and lab experiment by Lowe et al. JGR 2005, 2007
       !
       fwtable(i)  = min(fsubw,0.5)
       !           WRITE(994,*) 'Friction factor:',I,ABR,FWTABLE(I)
     END DO
     RETURN

References abmin, delab, fwtable, kappa, kerkei(), and sizefwtable.

Referenced by w3gridmd::w3grid().

Variable Documentation

◆ abmin

real, parameter constants::abmin = -1.

ABMIN.

Definition at line 94 of file constants.F90.

94 REAL, PARAMETER :: ABMIN = -1.

Referenced by tabu_fw(), w3sbt4md::w3sbt4(), w3sic2md::w3sic2(), and w3sic3md::w3sic3().

◆ dair

real, parameter constants::dair = 1.225

DAIR Density of air (kg/m3).

Definition at line 63 of file constants.F90.

63 REAL, PARAMETER :: DAIR = 1.225

Referenced by w3wdatmd::w3dimw(), w3exnc(), w3expo(), w3src6md::w3sds6(), w3wavemd::w3wave(), and wmupdtmd::wmupd2().

◆ debug_element

integer, parameter constants::debug_element = 50

DEBUG_ELEMENT Element number used for debug.

Definition at line 100 of file constants.F90.

100 INTEGER, PARAMETER :: DEBUG_ELEMENT = 50

◆ debug_node

integer, parameter constants::debug_node = 1014

DEBUG_NODE Node number used for debugging.

Definition at line 99 of file constants.F90.

99 INTEGER, PARAMETER :: DEBUG_NODE = 1014

Referenced by w3sdb1md::w3sdb1(), w3src4md::w3sds4(), and w3wavemd::w3wave().

◆ delab

real constants::delab

DELAB.

Definition at line 93 of file constants.F90.

93 REAL :: DELAB

Referenced by tabu_fw(), w3sbt4md::w3sbt4(), w3sic2md::w3sic2(), and w3sic3md::w3sic3().

◆ dera

real, parameter constants::dera = PI / 180.

DERA Conversion factor from degrees to radians.

Definition at line 77 of file constants.F90.

77 REAL, PARAMETER :: DERA = pi / 180.

Referenced by w3smcomd::calc_interp(), w3oacpmd::cpl_oasis_grid(), w3servmd::dist_sphere(), w3snl3md::insnl3(), w3partmd::ptmean(), w3psmcmd::smcdcxy(), w3psmcmd::smcdhxy(), w3psmcmd::smcgradn(), w3ogcmmd::snd_fields_to_ocean(), w3triamd::ug_gradients(), w3updtmd::w3dzxy(), w3exnc(), w3expo(), w3iopomd::w3iope(), w3iogomd::w3outg(), w3partmd::w3part(), w3prtide(), w3psmcmd::w3psmc(), w3gdatmd::w3setref(), w3ref1md::w3sref(), w3strt(), w3servmd::w3thrtn(), w3updtmd::w3ucur(), w3updtmd::w3uini(), w3updtmd::w3ulev(), w3updtmd::w3utau(), w3updtmd::w3utrn(), w3updtmd::w3uwnd(), w3wavemd::w3wave(), w3pro1md::w3xyp1(), w3pro2md::w3xyp2(), w3pro3md::w3xyp3(), w3servmd::w3xyrtn(), wmgridmd::wmgeql(), and wmgridmd::wmghgh().

◆ dwat

real, parameter constants::dwat = 1000.

DWAT Density of water (kg/m3).

Definition at line 62 of file constants.F90.

62 REAL, PARAMETER :: DWAT = 1000.

Referenced by w3sic3md::drfun_dble_cheng(), w3sic3md::drfun_quad_cheng(), wmesmfmd::fieldindex(), w3dispmd::liu_reverse_dispersion(), w3src6md::tau_wave_atmos(), w3wavset::trig_compute_lh_stress(), w3exnc(), w3expo(), w3fld1md::w3fld1(), w3fld2md::w3fld2(), w3sic3md::w3ic3wncg_v1(), w3iogomd::w3outg(), w3sbt8md::w3sbt8(), w3sbt9md::w3sbt9(), w3src4md::w3sds4(), w3sic3md::w3sic3(), w3src6md::w3sin6(), w3sis2md::w3sis2(), w3srcemd::w3srce(), and w3swldmd::w3swl4().

◆ file_endian

character(*), parameter constants::file_endian = ENDIANNESS

FILE_ENDIAN Filled by preprocessor with 'big_endian', 'little_endian', or 'native'.

Definition at line 86 of file constants.F90.

86 CHARACTER(*), PARAMETER :: FILE_ENDIAN = endianness

Referenced by w3snl4md::insnl4(), w3gkemd::prepkernelio(), w3smcomd::read_smcint(), w3src4md::tabu_tauhf2(), w3bounc(), w3bound(), w3exgo(), w3fldsmd::w3fldo(), w3iobcmd::w3iobc(), w3iogomd::w3iogo(), w3iogrmd::w3iogr(), w3iopomd::w3iopo(), w3iorsmd::w3iors(), w3iosfmd::w3iosf(), w3iotrmd::w3iotr(), w3outp(), w3prep(), w3prnc(), w3trck(), w3trnc(), w3strkmd::wavetracking_nws_v2(), and m_fileio::z_fileio().

◆ fwtable

real, dimension(0:sizefwtable) constants::fwtable

FWTABLE.

Definition at line 92 of file constants.F90.

92 REAL :: FWTABLE(0:SIZEFWTABLE)

Referenced by tabu_fw(), w3sbt4md::w3sbt4(), w3sic2md::w3sic2(), w3sic3md::w3sic3(), and w3src4md::w3sin4().

◆ g1pi1i

real, parameter constants::g1pi1i = 1. / ( GRAV * TPI )

G1PI1I Inverse of gravity * 2 * Pi.

Definition at line 82 of file constants.F90.

82 REAL, PARAMETER :: G1PI1I = 1. / ( grav * tpi )

Referenced by w3src2md::w3sds2().

◆ g2pi3i

real, parameter constants::g2pi3i = 1. / ( GRAV**2 * TPI**3 )

G2PI3I Inverse of gravity^2 * (2*Pi)^3.

Definition at line 81 of file constants.F90.

81 REAL, PARAMETER :: G2PI3I = 1. / ( grav**2 * tpi**3 )

Referenced by w3src2md::w3sds2().

◆ grav

real, parameter constants::grav = 9.806

GRAV Acc.

of gravity (m/s2).

Definition at line 61 of file constants.F90.

61 REAL, PARAMETER :: GRAV = 9.806

Referenced by pdlib_w3profsmd::action_limiter_local(), w3iogomd::calc_u3stokes(), w3src3md::calc_ustar(), w3src4md::calc_ustar(), pdlib_w3profsmd::calcarray_jacobi_source_1(), pdlib_w3profsmd::calcarray_jacobi_source_2(), w3snl1md::couple(), w3gig1md::df1f2theta(), w3dispmd::distab(), w3sic3md::drfun_dble_cheng(), w3sic3md::drfun_quad_cheng(), wmesmfmd::fieldindex(), w3sic5md::fsdisp(), w3sic3md::ic3table_cheng(), w3src4md::insin4(), w3snl2md::insnl2(), w3snl3md::insnl3(), w3snl1md::insnlgqm(), w3snlsmd::insnls(), w3dispmd::liu_reverse_dispersion(), pdlib_w3profsmd::pdlib_init(), pdlib_w3profsmd::pdlib_jacobi_gauss_seidel_block(), w3iogomd::secondhh(), w3fld1md::sig2wn(), w3iogomd::skewness(), w3src3md::tabu_stress(), w3src4md::tabu_stress(), w3src3md::tabu_tauhf(), w3src4md::tabu_tauhf(), w3src4md::tabu_tauhf2(), w3src6md::tau_wave_atmos(), w3wavset::trig_compute_lh_stress(), w3canomd::w3add2ndorder(), w3exnc(), w3expo(), w3fld1md::w3fld1(), w3fld2md::w3fld2(), w3flx2md::w3flx2(), w3flx3md::w3flx3(), w3flx5md::w3flx5(), w3sic3md::w3ic3wncg_v1(), w3iogomd::w3outg(), w3psmcmd::w3psmc(), w3sis2md::w3rpwnice(), w3sbs1md::w3sbs1(), w3sbt4md::w3sbt4(), w3sbt8md::w3sbt8(), w3sbt9md::w3sbt9(), w3src2md::w3sds2(), w3src3md::w3sds3(), w3src4md::w3sds4(), w3src6md::w3sds6(), w3sic2md::w3sic2(), w3sic3md::w3sic3(), w3src2md::w3sin2(), w3src3md::w3sin3(), w3src4md::w3sin4(), w3src6md::w3sin6(), w3sis2md::w3sis2(), w3sln1md::w3sln1(), w3snl3md::w3snl3(), w3snl5md::w3snl5(), w3snlsmd::w3snls(), w3src2md::w3spr2(), w3src4md::w3spr4(), w3srcemd::w3srce(), w3ref1md::w3sref(), w3str1md::w3str1(), w3swldmd::w3swl4(), w3swldmd::w3swl6(), w3updtmd::w3uini(), w3updtmd::w3uwnd(), w3pro1md::w3xyp1(), w3pro2md::w3xyp2(), w3pro3md::w3xyp3(), w3dispmd::wavnu1(), w3dispmd::wavnu2(), w3dispmd::wavnu3(), wmgridmd::wmgeql(), and wmgridmd::wmghgh().

◆ hpi

real, parameter constants::hpi = 0.5 * PI

HPI 1/2*Pi.

Definition at line 73 of file constants.F90.

73 REAL, PARAMETER :: HPI = 0.5 * pi

◆ hpiinv

real, parameter constants::hpiinv = 1. / HPI

HPIINV Inverse of 1/2*Pi.

Definition at line 75 of file constants.F90.

75 REAL, PARAMETER :: HPIINV = 1. / hpi

◆ is_esmf_component

logical constants::is_esmf_component = .FALSE.

IS_ESMF_COMPONENT Flag for model invoked via ESMF.

Definition at line 109 of file constants.F90.

109 LOGICAL :: IS_ESMF_COMPONENT = .false.

Referenced by wmesmfmd::setservices(), wminitmd::wminit(), and wminitmd::wminitnml().

◆ kappa

real, parameter constants::kappa = 0.40

KAPPA von Karman's constant (N.D.).

Definition at line 69 of file constants.F90.

69 REAL, PARAMETER :: KAPPA = 0.40

Referenced by w3src4md::calc_ustar(), tabu_fw(), w3src4md::tabu_stress(), w3src4md::tabu_tauhf(), w3src4md::tabu_tauhf2(), w3fld1md::w3fld1(), w3fld2md::w3fld2(), w3flx5md::w3flx5(), and w3src4md::w3sin4().

◆ lpdlib

logical constants::lpdlib = .FALSE.

LPDLIB Logical for using the PDLIB library.

Definition at line 101 of file constants.F90.

101 LOGICAL :: LPDLIB = .false.

Referenced by pdlib_w3profsmd::block_solver_init(), wmesmfmd::createexpmesh(), wmesmfmd::fieldindex(), w3parall::get_jsea_ibelong(), w3parall::init_get_isea(), w3parall::init_get_jsea_isproc(), wmmdatmd::init_get_jsea_isproc_glob(), w3parall::prop_freq_shift(), w3parall::prop_freq_shift_m2(), w3parall::prop_refraction_pr3(), w3triamd::readmsh(), w3triamd::readmsh_iobp(), pdlib_w3profsmd::set_iobdp_pdlib(), pdlib_w3profsmd::set_iobpa_pdlib(), w3parall::set_up_nseal_nsealm(), w3profsmd::setdepth(), wmesmfmd::setexport(), wmesmfmd::setupimpbmsk(), w3triamd::ug_gradients(), w3adatmd::w3dima(), w3wdatmd::w3dimw(), w3initmd::w3init(), w3iopomd::w3iope(), w3iorsmd::w3iors(), w3initmd::w3mpii(), w3initmd::w3mpio(), w3wavemd::w3scat(), w3shel(), w3ref1md::w3sref(), and w3wavemd::w3wave().

◆ lsetup

logical constants::lsetup = .FALSE.

LSETUP Logical LSETUP is not used.

Definition at line 102 of file constants.F90.

102 LOGICAL :: LSETUP = .false.

◆ nu_air

real, parameter constants::nu_air = 1.4E-5

NU_AIR Kinematic viscosity of air (m2/s).

Definition at line 64 of file constants.F90.

64 REAL, PARAMETER :: NU_AIR = 1.4e-5

Referenced by w3src4md::calc_ustar(), w3flx5md::w3flx5(), w3src4md::w3sin4(), and w3src4md::w3spr4().

◆ nu_water

real, parameter constants::nu_water = 1.31E-6

NU_WATER Kinematic viscosity of water (m2/s).

Definition at line 67 of file constants.F90.

67 REAL, PARAMETER :: NU_WATER = 1.31e-6

Referenced by w3sbt8md::w3sbt8(), and w3sbt9md::w3sbt9().

◆ pi

real, parameter constants::pi = 3.141592653589793

PI Value of Pi.

Definition at line 71 of file constants.F90.

71 REAL, PARAMETER :: PI = 3.141592653589793

Referenced by w3fld1md::appendtail(), w3strkmd::mean_anglev2(), w3strkmd::mean_anglev3(), w3tidemd::setvuf(), w3tidemd::setvuf_fast(), w3psmcmd::smcdhxy(), w3ogcmmd::snd_fields_to_ocean(), w3uostmd::uost_srctrmcompute(), w3gig1md::w3addig(), w3bounc(), w3bound(), w3bullmd::w3bull(), w3servmd::w3eqtoll(), w3fld1md::w3fld1(), w3fld2md::w3fld2(), w3servmd::w3lltoeq(), w3iogomd::w3outg(), w3outp(), w3psmcmd::w3psmc(), w3sbs1md::w3sbs1(), w3sbt8md::w3sbt8(), w3sbt9md::w3sbt9(), w3sdb1md::w3sdb1(), w3src4md::w3sds4(), w3gdatmd::w3setref(), w3sis2md::w3sis2(), w3str1md::w3str1(), w3updtmd::w3ucur(), w3updtmd::w3utau(), w3updtmd::w3uwnd(), w3dispmd::wavnu3(), and ww3_systrk().

◆ rade

real, parameter constants::rade = 180. / PI

RADE Conversion factor from radians to degrees.

Definition at line 76 of file constants.F90.

76 REAL, PARAMETER :: RADE = 180. / pi

Referenced by w3servmd::dist_sphere(), w3src4md::insin4(), w3snl2md::insnl2(), w3snl3md::insnl3(), w3snlsmd::insnls(), w3partmd::ptmean(), s2grid(), w3servmd::uv_to_mag_dir(), uv_to_mag_dir(), w3iosfmd::w3cprt(), w3cspcmd::w3cspc(), w3gdatmd::w3dims(), w3exgb(), w3exgo(), w3exnc(), w3expo(), w3iopomd::w3iope(), w3iosfmd::w3iosf(), w3iogomd::w3outg(), w3smcomd::w3s2xy_smc(), w3smcomd::w3s2xy_smcnn(), w3smcomd::w3s2xy_smcnn_int(), w3smcomd::w3s2xy_smcrg_dir(), w3src4md::w3sds4(), w3gdatmd::w3setref(), w3src1md::w3sin1(), w3src2md::w3sin2(), w3src3md::w3sin3(), w3sln1md::w3sln1(), and w3srcemd::w3srce().

◆ radius

real, parameter constants::radius = 4.E7 * TPIINV

RADIUS Radius of the earth (m).

Definition at line 79 of file constants.F90.

79 REAL, PARAMETER :: RADIUS = 4.e7 * tpiinv

Referenced by w3smcomd::calc_interp(), w3oacpmd::cpl_oasis_grid(), w3psmcmd::smcgradn(), w3triamd::ug_gradients(), w3updtmd::w3dzxy(), w3exnc(), w3expo(), w3psmcmd::w3psmc(), w3ref1md::w3sref(), w3updtmd::w3uini(), w3updtmd::w3utrn(), w3wavemd::w3wave(), w3pro1md::w3xyp1(), w3pro2md::w3xyp2(), w3pro3md::w3xyp3(), wmgridmd::wmgeql(), and wmgridmd::wmghgh().

◆ sed_sg

real, parameter constants::sed_sg = 2.65

SED_SG Specific gravity of sediments (N.D.).

Definition at line 68 of file constants.F90.

68 REAL, PARAMETER :: SED_SG = 2.65

Referenced by w3sbt4md::w3sbt4().

◆ sizefwtable

integer, parameter constants::sizefwtable =300

SIZEFWTABLE.

Definition at line 91 of file constants.F90.

91 INTEGER, PARAMETER :: SIZEFWTABLE=300

Referenced by tabu_fw(), w3sbt4md::w3sbt4(), w3sic2md::w3sic2(), w3sic3md::w3sic3(), and w3src4md::w3sin4().

◆ srce_direct

integer, parameter constants::srce_direct = 0

srce_direct

Definition at line 96 of file constants.F90.

96 INTEGER, PARAMETER :: srce_direct = 0

Referenced by w3srcemd::w3srce(), and w3wavemd::w3wave().

◆ srce_imp_post

integer, parameter constants::srce_imp_post = 1

srce_imp_post

Definition at line 97 of file constants.F90.

97 INTEGER, PARAMETER :: srce_imp_post = 1

Referenced by w3srcemd::w3srce(), and w3wavemd::w3wave().

◆ srce_imp_pre

integer, parameter constants::srce_imp_pre = 2

srce_imp_pre

Definition at line 98 of file constants.F90.

98 INTEGER, PARAMETER :: srce_imp_pre = 2

Referenced by w3srcemd::w3srce(), and w3wavemd::w3wave().

◆ tpi

real, parameter constants::tpi = 2.0 * PI

TPI 2*Pi.

Definition at line 72 of file constants.F90.

72 REAL, PARAMETER :: TPI = 2.0 * pi

Referenced by pdlib_w3profsmd::action_limiter_local(), w3fld1md::appendtail(), pdlib_w3profsmd::block_solver_init(), w3iogomd::calc_u3stokes(), pdlib_w3profsmd::calcarray_jacobi_source_1(), pdlib_w3profsmd::calcarray_jacobi_source_2(), expand(), wmesmfmd::fieldindex(), w3sic5md::fsdisp(), w3sis2md::insis2(), w3snl5md::insnl5(), w3dispmd::liu_forward_dispersion(), w3dispmd::liu_reverse_dispersion(), pdlib_w3profsmd::pdlib_init(), pdlib_w3profsmd::pdlib_jacobi_gauss_seidel_block(), w3partmd::ptmean(), w3iogomd::secondhh(), w3iogomd::skewness(), w3src3md::tabu_tauhf(), w3src6md::tau_wave_atmos(), uvtocart(), w3gig1md::w3addig(), w3bounc(), w3bound(), w3bullmd::w3bull(), w3iosfmd::w3cprt(), w3cspcmd::w3cspc(), w3exnc(), w3expo(), w3fld1md::w3fld1(), w3fld2md::w3fld2(), w3flx2md::w3flx2(), w3flx3md::w3flx3(), w3sic3md::w3ic3wncg_v1(), w3sic5md::w3ic5wncg(), w3initmd::w3init(), w3iopomd::w3iope(), w3pro1md::w3ktp1(), w3pro2md::w3ktp2(), w3pro3md::w3ktp3(), w3iogomd::w3outg(), w3outp(), w3partmd::w3part(), w3sbt4md::w3sbt4(), w3src3md::w3sds3(), w3src4md::w3sds4(), w3src6md::w3sds6(), w3gdatmd::w3setref(), w3sic1md::w3sic1(), w3sic2md::w3sic2(), w3sic3md::w3sic3(), w3sic4md::w3sic4(), w3sic5md::w3sic5(), w3src2md::w3sin2(), w3src3md::w3sin3(), w3src4md::w3sin4(), w3src6md::w3sin6(), w3sis2md::w3sis2(), w3snl1md::w3snl1(), w3snl2md::w3snl2(), w3snl3md::w3snl3(), w3snl4md::w3snl4(), w3snl5md::w3snl5(), w3snl1md::w3snlgqm(), w3snlsmd::w3snls(), w3srcemd::w3srce(), w3ref1md::w3sref(), w3str1md::w3str1(), w3strt(), w3servmd::w3thrtn(), w3updtmd::w3ucur(), w3updtmd::w3ulev(), w3updtmd::w3utau(), w3updtmd::w3uwnd(), w3servmd::w3xyrtn(), and w3sic3md::wn_precalc_cheng().

◆ tpiinv

real, parameter constants::tpiinv = 1. / TPI

TPIINV Inverse of 2*Pi.

Definition at line 74 of file constants.F90.

74 REAL, PARAMETER :: TPIINV = 1. / tpi

Referenced by pdlib_w3profsmd::block_solver_init(), w3iogomd::calc_u3stokes(), pdlib_w3profsmd::calcarray_jacobi_source_1(), pdlib_w3profsmd::calcarray_jacobi_source_2(), wmesmfmd::calcstokes3d(), pdlib_w3profsmd::compute_mean_param(), w3src3md::insin3(), w3src4md::insin4(), w3sis2md::insis2(), w3snl1md::insnl1(), w3snl2md::insnl2(), w3snl3md::insnl3(), w3snlsmd::insnls(), pdlib_w3profsmd::pdlib_jacobi_gauss_seidel_block(), w3partmd::ptmean(), w3iogomd::skewness(), w3bounc(), w3bound(), w3crnc(), w3exnc(), w3expo(), w3iogomd::w3outg(), w3outp(), w3sdb1md::w3sdb1(), w3src2md::w3sds2(), w3sis2md::w3sis2(), w3snl1md::w3snl1(), w3snl3md::w3snl3(), w3snlsmd::w3snls(), w3src0md::w3spr0(), w3src1md::w3spr1(), w3src2md::w3spr2(), w3src3md::w3spr3(), w3src4md::w3spr4(), w3src6md::w3spr6(), w3srcemd::w3srce(), w3ref1md::w3sref(), w3strt(), w3updtmd::w3uini(), and w3wavemd::w3wave().

◆ tstout

logical, parameter constants::tstout = .FALSE.

TSTOUT Flag for generation of test files.

Definition at line 56 of file constants.F90.

56 LOGICAL, PARAMETER :: TSTOUT = .false.

Referenced by w3initmd::w3init(), wminitmd::wminit(), and wminitmd::wminitnml().

◆ undef

real constants::undef = -999.9

UNDEF Value for undefined variable in output.

Definition at line 84 of file constants.F90.

84 REAL :: UNDEF = -999.9

Referenced by gxoutf(), w3partmd::ptmean(), w3nmlounfmd::read_field_nml(), s2grid(), w3servmd::uv_to_mag_dir(), uv_to_mag_dir(), w3exgb(), w3exgi(), w3exgo(), w3exnc(), w3grib(), w3grid_interp(), w3initmd::w3init(), w3iogomd::w3iogo(), w3ounf(), w3ounp(), w3outf(), w3iogomd::w3outg(), w3outp(), w3smcomd::w3s2xy_smc(), w3smcomd::w3s2xy_smcnn(), w3smcomd::w3s2xy_smcnn_int(), w3smcomd::w3s2xy_smcrg(), w3smcomd::w3s2xy_smcrg_dir(), w3servmd::w3thrtn(), w3wavemd::w3wave(), and w3servmd::w3xyrtn().

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ kerkei()

◆ kzeone()

◆ tabu_fw()

Variable Documentation

◆ abmin

◆ dair

◆ debug_element

◆ debug_node

◆ delab

◆ dera

◆ dwat

◆ file_endian

◆ fwtable

◆ g1pi1i

◆ g2pi3i

◆ grav

◆ hpi

◆ hpiinv

◆ is_esmf_component

◆ kappa

◆ lpdlib

◆ lsetup

◆ nu_air

◆ nu_water

◆ pi

◆ rade

◆ radius

◆ sed_sg

◆ sizefwtable

◆ srce_direct

◆ srce_imp_post

◆ srce_imp_pre

◆ tpi

◆ tpiinv

◆ tstout

◆ undef