SAMSIM
|
Module determines physical constants to be used by the SAMSIM Seaice model. More...
Variables | |
integer, parameter | wp = SELECTED_REAL_KIND(12, 307) |
set working precision _wp More... | |
real, parameter | pi = 3.1415_wp |
real, parameter | grav = 9.8061_wp |
gravitational constant [m/s^2] More... | |
real(wp), parameter | k_s = 2.2_wp |
solid heat conductivity [J / m s K] 2.2 More... | |
real(wp), parameter | k_l = 0.523_wp |
liquid heat conductivity [J / m s K] 0.523 More... | |
real(wp), parameter | c_s = 2020.0_wp |
solid heat capacity [J/ kg K] More... | |
real(wp), parameter | c_s_beta = 7.6973_wp |
linear solid heat capacity approximation [J/ kg K^2] c_s = c_s+c_s_beta*T More... | |
real(wp), parameter | c_l = 3400._wp |
liquid heat capacity [J/ kg K] More... | |
real(wp), parameter | rho_s = 920._wp |
density of solid [kg / m^3] More... | |
real(wp), parameter | rho_l = 1028.0_wp |
density of liquid [kg / m^3] More... | |
real(wp), parameter | latent_heat = 333500._wp |
latent heat release [J/kg] More... | |
real(wp), parameter | zerok = 273.15_wp |
Zero degrees Celsius in Kelvin [K]. More... | |
real(wp), parameter | bbeta = 0.8_wp*1e-3 |
concentration expansion coefficient [kg / (m^3 ppt)] More... | |
real(wp), parameter | mu = 2.55_wp*1e-3 |
dynamic viscosity [kg /m s] More... | |
real(wp), parameter | kappa_l = k_l/rho_l/c_l |
heat diffusivity of water More... | |
real(wp), parameter | sigma = 5.6704_wp*1e-8 |
Stefan Boltzmann constant [W/(m^2*K^4)]. More... | |
real(wp), parameter | psi_s_min = 0.05_wp |
The amount of ice that the lowest layer can have before it counts as an ice layer. More... | |
real(wp), parameter | neg_free = -0.05_wp |
The distance the freeboard can be below 0 before water starts flooding through cracks. More... | |
real(wp), parameter | x_grav = 0.000584_wp |
real(wp), parameter | ray_crit = 4.89_wp |
real(wp), parameter | para_flush_horiz = 1.0_wp |
determines relationship of horizontal flow distance in during flushing (guess 1) More... | |
real(wp), parameter | para_flush_gamma = 0.9_wp |
Strength of desalination per timestep (guess) More... | |
real(wp), parameter | psi_s_top_min = 0.40_wp |
if psi_s is below this value meltwater forms (guess) 0.4 More... | |
real(wp), parameter | ratio_flood = 1.50_wp |
Ratio of flooded to dissolve snow, plays an important role in subroutine flood. More... | |
real(wp), parameter | ref_salinity = 34._wp |
Reference salinity [g/kg] used to calculate freshwater column. More... | |
real(wp), parameter | rho_snow = 330._wp |
density of new snow [kg/m**3], !< Niels, 2017 add: can be adjusted to lab values if they are measured More... | |
real(wp), parameter | gas_snow_ice = 0.10_wp |
volume of gas percentage in new snow ice due to flooding, no longer used More... | |
real(wp), parameter | gas_snow_ice2 = 0.20_wp |
volume of gas percentage in new snow ice due to snow melting (Eicken 95) More... | |
real(wp), parameter | emissivity_ice = 0.95_wp |
Emissivity of water and ice. More... | |
real(wp), parameter | emissivity_snow = 1.00_wp |
Emissivity of Snow. More... | |
real(wp), parameter | penetr = 0.30_wp |
Amount of penetrating sw radiation. More... | |
real(wp), parameter | extinc = 2.00_wp |
Extinction coefficient of ice. More... | |
real(wp), parameter | turb_a = 0.1_wp*0.05_wp*rho_l/86400._wp |
Standard turbulence [kg/s] WARNING no source, just set so that 5cm of water are overturned each day. More... | |
real(wp), parameter | turb_b = 0.05_wp |
Exponential turbulence slope [m**3/kg] WARNING no source, simple guess. More... | |
real(wp) | max_flux_plate = 50.0 |
Maximal heating rate of a heating plate. More... | |
real(wp) | k_snow_flush = 0.75_wp |
Niels, 2017 add: Percentage of excess liquid water content in the snow that is used for flushing instead of forming slush. More... | |
real(wp) | k_styropor = 0.8_wp |
Niels, 2017 add: heat conduction of styropor (empirical value to fit measurement data) More... | |
Module determines physical constants to be used by the SAMSIM Seaice model.
Many values are taken from Notz 2005, Table 5.2.
COPYRIGHT
This file is part of SAMSIM.
SAMSIM is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
SAMSIM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with SAMSIM. If not, see http://www.gnu.org/licenses/.
real(wp), parameter mo_parameters::bbeta = 0.8_wp*1e-3 |
concentration expansion coefficient [kg / (m^3 ppt)]
real(wp), parameter mo_parameters::c_l = 3400._wp |
liquid heat capacity [J/ kg K]
real(wp), parameter mo_parameters::c_s = 2020.0_wp |
solid heat capacity [J/ kg K]
real(wp), parameter mo_parameters::c_s_beta = 7.6973_wp |
linear solid heat capacity approximation [J/ kg K^2] c_s = c_s+c_s_beta*T
real(wp), parameter mo_parameters::emissivity_ice = 0.95_wp |
Emissivity of water and ice.
real(wp), parameter mo_parameters::emissivity_snow = 1.00_wp |
Emissivity of Snow.
real(wp), parameter mo_parameters::extinc = 2.00_wp |
Extinction coefficient of ice.
real(wp), parameter mo_parameters::gas_snow_ice = 0.10_wp |
volume of gas percentage in new snow ice due to flooding, no longer used
real(wp), parameter mo_parameters::gas_snow_ice2 = 0.20_wp |
volume of gas percentage in new snow ice due to snow melting (Eicken 95)
real, parameter mo_parameters::grav = 9.8061_wp |
gravitational constant [m/s^2]
real(wp), parameter mo_parameters::k_l = 0.523_wp |
liquid heat conductivity [J / m s K] 0.523
real(wp), parameter mo_parameters::k_s = 2.2_wp |
solid heat conductivity [J / m s K] 2.2
real(wp) mo_parameters::k_snow_flush = 0.75_wp |
Niels, 2017 add: Percentage of excess liquid water content in the snow that is used for flushing instead of forming slush.
real(wp) mo_parameters::k_styropor = 0.8_wp |
Niels, 2017 add: heat conduction of styropor (empirical value to fit measurement data)
real(wp), parameter mo_parameters::latent_heat = 333500._wp |
latent heat release [J/kg]
real(wp) mo_parameters::max_flux_plate = 50.0 |
Maximal heating rate of a heating plate.
real(wp), parameter mo_parameters::mu = 2.55_wp*1e-3 |
dynamic viscosity [kg /m s]
real(wp), parameter mo_parameters::neg_free = -0.05_wp |
The distance the freeboard can be below 0 before water starts flooding through cracks.
real(wp), parameter mo_parameters::para_flush_gamma = 0.9_wp |
Strength of desalination per timestep (guess)
real(wp), parameter mo_parameters::para_flush_horiz = 1.0_wp |
determines relationship of horizontal flow distance in during flushing (guess 1)
real(wp), parameter mo_parameters::penetr = 0.30_wp |
Amount of penetrating sw radiation.
real, parameter mo_parameters::pi = 3.1415_wp |
real(wp), parameter mo_parameters::psi_s_min = 0.05_wp |
The amount of ice that the lowest layer can have before it counts as an ice layer.
real(wp), parameter mo_parameters::psi_s_top_min = 0.40_wp |
if psi_s is below this value meltwater forms (guess) 0.4
real(wp), parameter mo_parameters::ratio_flood = 1.50_wp |
Ratio of flooded to dissolve snow, plays an important role in subroutine flood.
real(wp), parameter mo_parameters::ray_crit = 4.89_wp |
real(wp), parameter mo_parameters::ref_salinity = 34._wp |
Reference salinity [g/kg] used to calculate freshwater column.
real(wp), parameter mo_parameters::rho_l = 1028.0_wp |
density of liquid [kg / m^3]
real(wp), parameter mo_parameters::rho_s = 920._wp |
density of solid [kg / m^3]
real(wp), parameter mo_parameters::rho_snow = 330._wp |
density of new snow [kg/m**3], !< Niels, 2017 add: can be adjusted to lab values if they are measured
real(wp), parameter mo_parameters::sigma = 5.6704_wp*1e-8 |
Stefan Boltzmann constant [W/(m^2*K^4)].
real(wp), parameter mo_parameters::turb_a = 0.1_wp*0.05_wp*rho_l/86400._wp |
Standard turbulence [kg/s] WARNING no source, just set so that 5cm of water are overturned each day.
real(wp), parameter mo_parameters::turb_b = 0.05_wp |
Exponential turbulence slope [m**3/kg] WARNING no source, simple guess.
integer, parameter mo_parameters::wp = SELECTED_REAL_KIND(12, 307) |
set working precision _wp
real(wp), parameter mo_parameters::x_grav = 0.000584_wp |
real(wp), parameter mo_parameters::zerok = 273.15_wp |
Zero degrees Celsius in Kelvin [K].