This module contains constants and parameters used in GFS near surface sea temperature scheme.

Variables
real(kind_phys), parameter, public	zero = 0.0_kind_phys

real(kind_phys), parameter, public	one = 1.0_kind_phys

real(kind_phys), parameter, public	half = 0.5_kind_phys

integer, public	niter_conv = 5

integer, public	niter_z_w = 5

integer, public	niter_sfs = 5

real(kind=kind_phys), parameter	sec_in_day = 86400.

real(kind=kind_phys), parameter	sec_in_hour = 3600.

real(kind=kind_phys), parameter, public	solar_time_6am = 21600.0

real(kind=kind_phys), parameter, public	const_rot = 0.000073
	constant to calculate corioli force

real(kind=kind_phys), parameter, public	ri_c = 0.65

real(kind=kind_phys), parameter, public	ri_g = 0.25

real(kind=kind_phys), parameter, public	eps_z_w = 0.01
	criteria to finish iterations for z_w

real(kind=kind_phys), parameter, public	eps_conv = 0.01
	criteria to finish iterations for d_conv

real(kind=kind_phys), parameter, public	eps_sfs = 0.01
	criteria to finish iterations for d_sfs

real(kind=kind_phys), parameter, public	z_w_max = 20.0
	max warm layer thickness

real(kind=kind_phys), parameter, public	z_w_min = 0.2
	min warm layer thickness

real(kind=kind_phys), parameter, public	z_w_ini = 0.2
	initial warm layer thickness in dtl_onset

real(kind=kind_phys), parameter, public	z_c_max = 0.01
	maximum of sub-layer thickness (m)

real(kind=kind_phys), parameter, public	z_c_ini = 0.001
	initial value of z_c

real(kind=kind_phys), parameter, public	ustar_a_min = 0.031
	minimum of friction wind speed (m/s): 0.031 ~ 1m/s at 10 m hight

real(kind=kind_phys), parameter, public	tau_min = 0.005
	minimum of wind stress for dtm

real(kind=kind_phys), parameter, public	exp_const = 9.5
	coefficient in exponet profile

real(kind=kind_phys), parameter, public	delz = 0.1
	vertical increment for integral calculation (m)

real(kind=kind_phys), parameter	von = 0.4
	von karman's "constant"

real(kind=kind_phys), parameter, public	t0k = 273.16
	celsius to kelvin

real(kind=kind_phys), parameter	gray = 0.97

real(kind=kind_phys), parameter, public	sst_max = 308.16

real(kind=kind_phys), parameter, public	tw_max = 5.0

real(kind=kind_phys), parameter, public	wd_max = 2.0

real(kind=kind_phys), parameter, public	omg_m = 1.0
	trace factor to apply salinity effect

real(kind=kind_phys), parameter	omg_rot = 1.0
	trace factor to apply rotation effect

real(kind=kind_phys), parameter, public	omg_sh = 1.0
	trace factor to apply sensible heat due to rainfall effect

real(kind=kind_phys), parameter, public	visw = 1.e-6
	m2/s kinematic viscosity water

real(kind=kind_phys), parameter	novalue = 0

real(kind=kind_phys), parameter	smallnumber = 1.e-6

real(kind=kind_phys), parameter	timestep_oc = sec_in_day/8.
	time step in the ocean model (3 hours)

real(kind=kind_phys), parameter	radian = 2.*pi/180.

real(kind=kind_phys), parameter, public	rad2deg = 180./pi

real(kind=kind_phys), parameter, public	cp_w = 4000.
	specific heat water (j/kg/k )

real(kind=kind_phys), parameter	rho0_w = 1022.0
	density water (kg/m3 ) (or 1024.438)

real(kind=kind_phys), parameter	vis_w = 1.e-6
	kinematic viscosity water (m2/s )

real(kind=kind_phys), parameter, public	tc_w = 0.6
	thermal conductivity water (w/m/k )

real(kind=kind_phys), parameter	capa_w = 3950.0
	heat capacity of sea water !

real(kind=kind_phys), parameter	thref = 1.0e-3
	reference value of specific volume (m**3/kg)

Variables