|
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)
|
|