CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches
module_nst_parameters.f90
1
4
12
13 use machine, only : kind_phys
14 !
15 ! air constants and coefficients from the atmospehric model
16 use physcons, only: &
17 eps => con_eps & !< con_rd/con_rv (nd)
18 ,cp_a => con_cp &
19 ,epsm1 => con_epsm1 &
20 ,hvap => con_hvap &
21 ,sigma_r => con_sbc &
22 ,grav => con_g &
23 ,omega => con_omega &
24 ,rvrdm1 => con_fvirt &
25 ,rd => con_rd &
26 ,rocp => con_rocp &
27 ,pi => con_pi
28
29 implicit none
30
31 private
32
33 public :: sigma_r
34 public :: zero, one, half
35 public :: niter_conv, niter_z_w, niter_sfs
36 public :: z_w_max, z_w_min, z_w_ini, z_c_max, z_c_ini, eps_z_w, eps_conv, eps_sfs
37 public :: ri_c, ri_g, omg_m, omg_sh, tc_w, visw, cp_w, t0k, ustar_a_min, delz, exp_const
38 public :: rad2deg, const_rot, tw_max, sst_max, solar_time_6am, tau_min, wd_max
39
40 real(kind_phys), parameter :: zero = 0.0_kind_phys, one = 1.0_kind_phys, half = 0.5_kind_phys
41 !
42 ! note: take timestep from here later
43 integer :: &
44 niter_conv = 5, &
45 niter_z_w = 5, &
46 niter_sfs = 5
47 !
48 ! general constants
49 real (kind=kind_phys), parameter :: &
50 sec_in_day = 86400. &
51 ,sec_in_hour = 3600. &
52 ,solar_time_6am = 21600.0 &
53 ,const_rot = 0.000073 &
54 ,ri_c = 0.65 &
55 ,ri_g = 0.25 &
56 ,eps_z_w = 0.01 &
57 ,eps_conv = 0.01 &
58 ,eps_sfs = 0.01 &
59 ,z_w_max = 20.0 &
60 ,z_w_min = 0.2 &
61 ,z_w_ini = 0.2 &
62 ,z_c_max = 0.01 &
63 ,z_c_ini = 0.001 &
64 ,ustar_a_min = 0.031 &
65 ,tau_min = 0.005 &
66 ,exp_const = 9.5 &
67 ,delz = 0.1 &
68 ,von = 0.4 &
69 ,t0k = 273.16 &
70 ,gray = 0.97 &
71 ,sst_max = 308.16 &
72 ,tw_max = 5.0 &
73 ,wd_max = 2.0 &
74 ,omg_m = 1.0 &
75 ,omg_rot = 1.0 &
76 ,omg_sh = 1.0 &
77 ,visw = 1.e-6 &
78 ,novalue = 0 &
79 ,smallnumber = 1.e-6 &
80 ,timestep_oc = sec_in_day/8. &
81 ,radian = 2.*pi/180. &
82 ,rad2deg = 180./pi &
83 ,cp_w = 4000. &
84 ,rho0_w = 1022.0 &
85 ,vis_w = 1.e-6 &
86 ,tc_w = 0.6 &
87 ,capa_w = 3950.0 &
88 ,thref = 1.0e-3
89
90!!$!============================================
91!!$
92!!$ ,lvapor=2.453e6 & ! latent heat of vaporization note: make it function of t ????? note the same as hvap
93!!$ ,alpha=1 ! thermal expansion coefficient
94!!$ ,beta ! saline contraction coefficient
95!!$ ,cp=1 !=1 specific heat of sea water
96!!$ ,g=1 ! acceleration due to gravity
97!!$ ,kw=1 ! thermal conductivity of water
98!!$ ,nu=1 !kinematic wiscosity
99!!$ ,rho_w=1 !water density
100!!$ ,rho_a=1 !air density
101!!$ ,l_vapr=2.453e6
102!!$ ,novalue=--1.0e+10
103!!$
104!!$c factors
105!!$ beta=1.2 !given as 1.25 in fairall et al.(1996)
106!!$ von=0.4 ! von karman's "constant"
107!!$c fdg=1.00 ! fairall's lkb rr to von karman adjustment
108!!$ fdg=1.00 !based on results from flux workshop august 1995
109!!$ tok=273.16 ! celsius to kelvin
110!!$ twopi=3.14159*2.
111!!$
112!!$c air constants and coefficients
113!!$ rgas=287.1 !j/kg/k gas const. dry air
114!!$ xlv=(2.501-0.00237*ts)*1e+6 !j/kg latent heat of vaporization at ts
115!!$ cpa=1004.67 !j/kg/k specific heat of dry air (businger 1982)
116!!$ cpv=cpa*(1+0.84*q) !moist air - currently not used (businger 1982)
117!!$ rhoa=p*100./(rgas*(t+tok)*(1.+.61*q)) !kg/m3 moist air density ( " )
118!!$ visa=1.326e-5*(1+6.542e-3*t+8.301e-6*t*t-4.84e-9*t*t*t) !m2/s
119!!$ !kinematic viscosity of dry air - andreas (1989) crrel rep. 89-11
120!!$c
121!!$c cool skin constants
122!!$ al=2.1e-5*(ts+3.2)**0.79 !water thermal expansion coefft.
123!!$ be=0.026 !salinity expansion coefft.
124!!$ cpw=4000. !j/kg/k specific heat water
125!!$ rhow=1022. !kg/m3 density water
126!!$ visw=1.e-6 !m2/s kinematic viscosity water
127!!$ tcw=0.6 !w/m/k thermal conductivity water
128!!$ bigc=16.*grav*cpw*(rhow*visw)**3/(tcw*tcw*rhoa*rhoa)
129!!$ wetc=0.622*xlv*qs/(rgas*(ts+tok)**2) !correction for dq;slope of sat. vap.
130!!$
131!!$!
132!!$! functions
133!!$
134!!$
135!!$ real, parameter :: timestep=86400. !integration time step, second
136!!$
137!!$ real, parameter :: grav =9.81 !gravity, kg/m/s^2
138!!$ real, parameter :: capa =3950.0 !heat capacity of sea water
139!!$ real, parameter :: rhoref = 1024.438 !sea water reference density, kg/m^3
140!!$ real , parameter :: hslab=50.0 !slab ocean depth
141!!$ real , parameter :: bad=-1.0e+10
142!!$ real , parameter :: tmin=2.68e+02
143!!$ real , parameter :: tmax=3.11e+02
144!!$
145!!$ real, parameter :: grav =9.81 !gravity, kg/m/s^2
146!!$ real, parameter :: capa =3950.0 !heat capacity of sea water
147!!$ real, parameter :: rhoref = 1024.438 !sea water reference density, kg/m^3
148!!$ real, parameter :: tmin=2.68e+02 !normal minimal temp
149!!$ real, parameter :: tmax=3.11e+02 !normal max temp
150!!$ real, parameter :: smin=1.0 !normal minimal salt
151!!$ real, parameter :: smax=50. !normal maximum salt
152!!$ real, parameter :: visct=1.e-5 !viscocity for temperature diffusion
153!!$ real, parameter :: viscs=1.e-5 !viscocity for salt diffusion
154!!$
155!!$
156end module module_nst_parameters