SAMSIM
Variables
mo_parameters Module Reference

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(wpmax_flux_plate = 50.0
 Maximal heating rate of a heating plate. More...
 
real(wpk_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(wpk_styropor = 0.8_wp
 Niels, 2017 add: heat conduction of styropor (empirical value to fit measurement data) More...
 

Detailed Description

Module determines physical constants to be used by the SAMSIM Seaice model.

Many values are taken from Notz 2005, Table 5.2.

Author
Philipp Griewank

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/.

Revision History
Started by Philipp Griewank 2010-07-08
add parameters: heat conductivity of styropor under special sea ice lab conditions and ratio of penetrating melt water by Niels Fuchs, MPIMET (2017-03-01)

Variable Documentation

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::kappa_l = k_l/rho_l/c_l

heat diffusivity of water

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].