module_nst_parameters Module Reference

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)