SAMSIM
Variables
mo_data Module Reference

Sets data and contains all flag descriptions. More...

Variables

real(wp), dimension(:), allocatable h
 Enthalpy [J]. More...
 
real(wp), dimension(:), allocatable h_abs
 specific Enthalpy [J/kg] More...
 
real(wp), dimension(:), allocatable q
 Heat in layer [J]. More...
 
real(wp), dimension(:), allocatable fl_q
 Heat flux between layers [J/s]. More...
 
real(wp), dimension(:), allocatable t
 Temperature [C]. More...
 
real(wp), dimension(:), allocatable s_bu
 Bulk Salinity [g/kg]. More...
 
real(wp), dimension(:), allocatable fl_s
 Salinity flux [(g/s]. More...
 
real(wp), dimension(:), allocatable s_abs
 Absolute Salinity [g]. More...
 
real(wp), dimension(:), allocatable s_br
 Brine salinity [g/kg]. More...
 
real(wp), dimension(:), allocatable thick
 Layer thickness [m]. More...
 
real(wp), dimension(:), allocatable m
 Mass [kg]. More...
 
real(wp), dimension(:), allocatable fl_m
 Mass fluxes between layers [kg]. More...
 
real(wp), dimension(:), allocatable v_s
 Volume [m^3] of solid. More...
 
real(wp), dimension(:), allocatable v_l
 Volume [m^3] of liquid. More...
 
real(wp), dimension(:), allocatable v_g
 Volume [m^3] of gas. More...
 
real(wp), dimension(:), allocatable v_ex
 Volume of brine due expelled due to freezing [m^3] of solid, gas & liquid. More...
 
real(wp), dimension(:), allocatable phi
 Solid mass fraction. More...
 
real(wp), dimension(:), allocatable psi_s
 Solid volume fraction. More...
 
real(wp), dimension(:), allocatable psi_l
 Liquid volume fraction. More...
 
real(wp), dimension(:), allocatable psi_g
 Gas volume fraction. More...
 
real(wp), dimension(:), allocatable ray
 Rayleigh number of each layer. More...
 
real(wp), dimension(:), allocatable perm
 
real(wp), dimension(:), allocatable flush_v
 
real(wp), dimension(:), allocatable flush_h
 
real(wp), dimension(:), allocatable flush_v_old
 
real(wp), dimension(:), allocatable flush_h_old
 Permeability [?]. More...
 
real(wp) dt
 Timestep [s]. More...
 
real(wp) thick_0
 Initial layer thickness [m]. More...
 
real(wp) time
 Time [s]. More...
 
real(wp) freeboard
 Height of ice surface above (or below) waterlevel [m]. More...
 
real(wp) t_freeze
 Freezing temperature [C]. More...
 
integer nlayer
 Number of layers. More...
 
integer n_bottom
 Number of bottom layers. More...
 
integer n_middle
 Number of middle layers. More...
 
integer n_top
 Number of top layers. More...
 
integer n_active
 Number of Layers active in the present. More...
 
integer i
 Index, normally used for time. More...
 
integer k
 Index, normally used for layer. More...
 
integer styropor_flag
 
real(wp) time_out
 Time between outputs [s]. More...
 
real(wp) time_total
 Time of simulation [s]. More...
 
integer i_time
 Number of timesteps. More...
 
integer i_time_out
 Number of timesteps between each output. More...
 
integer n_time_out
 Counts number of timesteps between output. More...
 
character *12000 format_t
 
character *12000 format_psi
 
character *12000 format_thick
 
character *12000 format_snow
 
character *12000 format_integer
 
character *12000 format_t2m_top
 
character *12000 format_bgc
 
character *12000 format_melt
 Format strings for output. Niels(2017) add: melt output. More...
 
character *12000 format_perm
 Niels(2017) add: permeability output. More...
 
real(wp) t_bottom
 Temperature of water beneath the ice [C]. More...
 
real(wp) t_top
 Temperature at the surface [C]. More...
 
real(wp) s_bu_bottom
 Salinity beneath the ice [g/kg]. More...
 
real(wp) t2m
 Two meter Temperature [C]. More...
 
real(wp) fl_q_bottom
 Bottom heat flux [J*s]. More...
 
real(wp) psi_s_snow
 Solid volume fraction of snow layer. More...
 
real(wp) psi_l_snow
 Liquid volume fraction of snow layer. More...
 
real(wp) psi_g_snow
 Gas volume fraction of snow layer. More...
 
real(wp) phi_s
 Solid mass fraction of snow layer. More...
 
real(wp) s_abs_snow
 Absolute salinity of snow layer [g]. More...
 
real(wp) h_abs_snow
 Absolute enthalpy of snow layer [J]. More...
 
real(wp) m_snow
 Mass of snow layer [kg]. More...
 
real(wp) t_snow
 Temperature of snow layer [C]. More...
 
real(wp) thick_snow
 
real(wp) test
 Thickness of snow layer [m]. More...
 
real(wp) liquid_precip
 Liquid precip, [meter of water/s]. More...
 
real(wp) solid_precip
 Solid precip, [meter of water /s]. More...
 
real(wp) fl_q_snow
 flow of heat into the snow layer More...
 
real(wp) energy_stored
 Total amount of energy stored, control is freezing point temperature of S_bu_bottom [J]. More...
 
real(wp) total_resist
 Thermal resistance of the whole column []. More...
 
real(wp) surface_water
 Percentage of water fraction in the top 5cm [%]. More...
 
real(wp) freshwater
 Meters of freshwater stored in column [m]. More...
 
real(wp) thickness
 Meters of ice [m]. More...
 
real(wp) bulk_salin
 Salt/Mass [ppt]. More...
 
real(wp) thick_min
 Parameter for snow, determines when snow is in thermal equilibrium with the ice and when it is totally neglected. More...
 
real(wp), save t_test
 First guess for getT subroutine. More...
 
real(wp) albedo
 Amount of short wave radiation which is reflected at the top surface. More...
 
real(wp) fl_sw
 Incoming shortwave radiation [W/m**2]. More...
 
real(wp) fl_lw
 Incoming longwave radiation [W/m**2]. More...
 
real(wp) fl_sen
 Sensitive heat flux [W/m**2]. More...
 
real(wp) fl_lat
 Latent heat flux [W/m**2]. More...
 
real(wp) fl_rest
 Bundled longwave,sensitive and latent heat flux [W/m**2]. More...
 
real(wp), dimension(:), allocatable fl_rad
 Energy flux of absorbed sw radiation of each layer [J/s]. More...
 
real(wp) grav_drain
 brine flux of gravity drainage between two outputs [kg/s] More...
 
real(wp) grav_salt
 salt flux moved by gravity drainage between two outputs [kg*ppt/s] More...
 
real(wp) grav_temp
 average temperature of gravity drainage brine between two outputs [T] More...
 
real(wp) melt_thick
 thickness of fully liquid part of top layer [m] More...
 
real(wp) melt_thick_snow
 
real(wp) melt_thick_snow_old
 Niels(2017) add: thickness of excess fully liquid part from snow_melt_processes [m]. More...
 
real(wp), dimension(3) melt_thick_output
 Niels, 2017 add: output field of surface liquid meltwater sizes. More...
 
real(wp) alpha_flux_instable
 Proportionality constant which determines energy flux by the temperature difference T_top>T2m [W/C]. More...
 
real(wp) alpha_flux_stable
 Proportionality constant which determines energy flux by the temperature difference T_top<T2m [W/C]. More...
 
integer atmoflux_flag
 1: Use mean climatology of Notz, 2: Use imported reanalysis data, 3: use fixed values defined in mo_init More...
 
integer grav_flag
 1: no gravity drainage, 2: Gravity drainage, 3: Simple Drainage More...
 
integer prescribe_flag
 1: nothing happens, 2: prescribed Salinity profile is prescribed at each timestep (does not disable brine dynamics, just overwrites the salinity!) More...
 
integer grav_heat_flag
 1: nothing happens, 2: compensates heatfluxes in grav_flag = 2 More...
 
integer flush_heat_flag
 1: nothing happens, 2: compensates heatfluxes in flush_flag = 5 More...
 
integer turb_flag
 1: No bottom turbulence, 2: Bottom mixing More...
 
integer salt_flag
 1: Sea salt, 2: NaCL More...
 
integer boundflux_flag
 1: top and bottom cooling plate, 2:top Notz fluxes, bottom cooling plate 3: top flux=a*(T-T_s) More...
 
integer flush_flag
 1: no flushing, 4:meltwater is removed artificially, 5:vert and horiz flushing, 6: simplified More...
 
integer flood_flag
 1: no flooding, 2:normal flooding, 3:simple flooding More...
 
integer bottom_flag
 1: nothing changes, 2: deactivates all bottom layer dynamics, useful for some debugging and idealized tests More...
 
integer debug_flag
 1: no raw layer output, 2: each layer is output at every timestep (warning, file size can be very large) More...
 
integer precip_flag
 0: solid and liquid precipitation, 1:phase determined by T2m More...
 
integer harmonic_flag
 1: minimal permeability is used to calculate Rayleigh number, 2:harmonic mean is used for Rayleigh number More...
 
integer tank_flag
 1: nothing, 2: S_bu_bottom and bgc_bottom are calculated as if the experiment is conducted in a tank More...
 
integer albedo_flag
 1: simple albedo, 2: normal albedo, see func_albedo for details More...
 
integer lab_snow_flag
 Niels, 2017 add: 0: lab setup without snow covers, 1: lab setup include snow influence on heat fluxes. More...
 
integer freeboard_snow_flag
 Niels, 2017 add: 0: respect the mass of snow in the freeboard calculation, 1: don't. More...
 
integer snow_flush_flag
 Niels, 2017 add: 0: all meltwater from snow forms slush, 1: meltwater partly leads to flushing, ratio defined by "k_snow_flush". More...
 
integer snow_precip_flag
 Niels, 2017 add: 0: all precipitation is set to zero, 1: physical behaviour. More...
 
integer length_input
 Sets the input length for atmoflux_flag==2, common value of 13169. More...
 
real(wp), dimension(:), allocatable tinput
 Niels, 2017 add: used to read in top temperature for field experiment tests, dimension needs to be set in the code. More...
 
real(wp), dimension(:), allocatable precipinput
 Niels, 2017 add: used to read in precipation for field experiment tests, dimension needs to be set in the code. More...
 
real(wp), dimension(:), allocatable ocean_t_input
 Niels, 2017 add: used to read in ocean temperature for field experiment tests, dimension needs to be set in the code. More...
 
real(wp), dimension(:), allocatable ocean_flux_input
 Niels, 2017 add: used to read in oceanic heat flux for field experiment tests, dimension needs to be set in the code. More...
 
real(wp), dimension(:), allocatable styropor_input
 Niels, 2017 add: if styropor is used in the lab on top of the ice to simulate snow heat fluxes. More...
 
real(wp), dimension(:), allocatable ttop_input
 Niels, 2017 add: used for testcase 111, comparison with greenland harp data, uppermost harp temperature is seen as Ttop. More...
 
real(wp), dimension(:), allocatable fl_sw_input
 Used to read in sw fluxes from ERA for atmoflux_flag==2. More...
 
real(wp), dimension(:), allocatable fl_lw_input
 Used to read in lw fluxes from ERA for atmoflux_flag==2. More...
 
real(wp), dimension(:), allocatable t2m_input
 Used to read in 2Tm from ERA for atmoflux_flag==2. More...
 
real(wp), dimension(:), allocatable precip_input
 Used to read in precipitation from ERA for atmoflux_flag==2. More...
 
real(wp), dimension(:), allocatable time_input
 Used to read in time from ERA for atmoflux_flag==2. More...
 
integer time_counter
 Keeps track of input data. More...
 
integer bgc_flag
 1: no bgc, 2:bgc More...
 
integer n_bgc
 Number of chemicals. More...
 
real(wp), dimension(:,:), allocatable fl_brine_bgc
 Brine fluxes in a matrix, [kg/s], first index is the layer of origin, and the second index is the layer of arrival. More...
 
real(wp), dimension(:,:), allocatable bgc_abs
 Absolute amount of chemicals [kmol] for each tracer. More...
 
real(wp), dimension(:,:), allocatable bgc_bu
 Bulk amounts of chemicals [kmol/kg]. More...
 
real(wp), dimension(:,:), allocatable bgc_br
 Brine concentrations of chems [kmol/kg]. More...
 
real(wp), dimension(:), allocatable bgc_bottom
 Bulk concentrations of chems below the ice [kmol/kg]. More...
 
real(wp), dimension(:), allocatable bgc_total
 Total of chems, for lab experiments with a fixed total amount. More...
 
real(wp) m_total
 Total initial water mass, for lab experiments with a fixed total amount. More...
 
real(wp) s_total
 Total initial salt mass, for lab experiments with a fixed total amount. More...
 
real(wp) tank_depth
 water depth in meters, used to calculate concentrations below ice for tank experiments More...
 
character *3 flush_question ='No!'
 Niels, 2017 add: used to indicate in stdout wether flushing occurs at this moment or not. More...
 
real(wp) melt_err =0._wp
 Niels, 2017 add: used to check how much meltwater vanishes in flushing routine. More...
 
integer length_input_lab
 Niels, 2017 add: used to allocate lab testcase input arrays in mo_init, set value in testcases. More...
 

Detailed Description

Sets data and contains all flag descriptions.

All data needed by mo_grotz are set in this module. Most arrays are allocated after the needed dimension is specified for each testcase in mo_init.f90.

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
Initialized by Philipp Griewank, IMPRS (2010-07-14)
Add several variables by Niels Fuchs, MPIMET (2017-03-01)

Variable Documentation

real(wp) mo_data::albedo

Amount of short wave radiation which is reflected at the top surface.

integer mo_data::albedo_flag

1: simple albedo, 2: normal albedo, see func_albedo for details

real(wp) mo_data::alpha_flux_instable

Proportionality constant which determines energy flux by the temperature difference T_top>T2m [W/C].

real(wp) mo_data::alpha_flux_stable

Proportionality constant which determines energy flux by the temperature difference T_top<T2m [W/C].

integer mo_data::atmoflux_flag

1: Use mean climatology of Notz, 2: Use imported reanalysis data, 3: use fixed values defined in mo_init

real(wp), dimension(:,:), allocatable mo_data::bgc_abs

Absolute amount of chemicals [kmol] for each tracer.

real(wp), dimension(:), allocatable mo_data::bgc_bottom

Bulk concentrations of chems below the ice [kmol/kg].

real(wp), dimension(:,:), allocatable mo_data::bgc_br

Brine concentrations of chems [kmol/kg].

real(wp), dimension(:,:), allocatable mo_data::bgc_bu

Bulk amounts of chemicals [kmol/kg].

integer mo_data::bgc_flag

1: no bgc, 2:bgc

real(wp), dimension(:), allocatable mo_data::bgc_total

Total of chems, for lab experiments with a fixed total amount.

integer mo_data::bottom_flag

1: nothing changes, 2: deactivates all bottom layer dynamics, useful for some debugging and idealized tests

integer mo_data::boundflux_flag

1: top and bottom cooling plate, 2:top Notz fluxes, bottom cooling plate 3: top flux=a*(T-T_s)

real(wp) mo_data::bulk_salin

Salt/Mass [ppt].

integer mo_data::debug_flag

1: no raw layer output, 2: each layer is output at every timestep (warning, file size can be very large)

real(wp) mo_data::dt

Timestep [s].

real(wp) mo_data::energy_stored

Total amount of energy stored, control is freezing point temperature of S_bu_bottom [J].

real(wp), dimension(:,:), allocatable mo_data::fl_brine_bgc

Brine fluxes in a matrix, [kg/s], first index is the layer of origin, and the second index is the layer of arrival.

real(wp) mo_data::fl_lat

Latent heat flux [W/m**2].

real(wp) mo_data::fl_lw

Incoming longwave radiation [W/m**2].

real(wp), dimension(:), allocatable mo_data::fl_lw_input

Used to read in lw fluxes from ERA for atmoflux_flag==2.

real(wp), dimension(:), allocatable mo_data::fl_m

Mass fluxes between layers [kg].

real(wp), dimension(:), allocatable mo_data::fl_q

Heat flux between layers [J/s].

real(wp) mo_data::fl_q_bottom

Bottom heat flux [J*s].

real(wp) mo_data::fl_q_snow

flow of heat into the snow layer

real(wp), dimension(:), allocatable mo_data::fl_rad

Energy flux of absorbed sw radiation of each layer [J/s].

real(wp) mo_data::fl_rest

Bundled longwave,sensitive and latent heat flux [W/m**2].

real(wp), dimension(:), allocatable mo_data::fl_s

Salinity flux [(g/s].

real(wp) mo_data::fl_sen

Sensitive heat flux [W/m**2].

real(wp) mo_data::fl_sw

Incoming shortwave radiation [W/m**2].

real(wp), dimension(:), allocatable mo_data::fl_sw_input

Used to read in sw fluxes from ERA for atmoflux_flag==2.

integer mo_data::flood_flag

1: no flooding, 2:normal flooding, 3:simple flooding

integer mo_data::flush_flag

1: no flushing, 4:meltwater is removed artificially, 5:vert and horiz flushing, 6: simplified

real(wp), dimension(:), allocatable mo_data::flush_h
real(wp), dimension(:), allocatable mo_data::flush_h_old

Permeability [?].

integer mo_data::flush_heat_flag

1: nothing happens, 2: compensates heatfluxes in flush_flag = 5

character*3 mo_data::flush_question ='No!'

Niels, 2017 add: used to indicate in stdout wether flushing occurs at this moment or not.

real(wp), dimension(:), allocatable mo_data::flush_v
real(wp), dimension(:), allocatable mo_data::flush_v_old
character*12000 mo_data::format_bgc
character*12000 mo_data::format_integer
character*12000 mo_data::format_melt

Format strings for output. Niels(2017) add: melt output.

character*12000 mo_data::format_perm

Niels(2017) add: permeability output.

character*12000 mo_data::format_psi
character*12000 mo_data::format_snow
character*12000 mo_data::format_t
character*12000 mo_data::format_t2m_top
character*12000 mo_data::format_thick
real(wp) mo_data::freeboard

Height of ice surface above (or below) waterlevel [m].

integer mo_data::freeboard_snow_flag

Niels, 2017 add: 0: respect the mass of snow in the freeboard calculation, 1: don't.

real(wp) mo_data::freshwater

Meters of freshwater stored in column [m].

real(wp) mo_data::grav_drain

brine flux of gravity drainage between two outputs [kg/s]

integer mo_data::grav_flag

1: no gravity drainage, 2: Gravity drainage, 3: Simple Drainage

integer mo_data::grav_heat_flag

1: nothing happens, 2: compensates heatfluxes in grav_flag = 2

real(wp) mo_data::grav_salt

salt flux moved by gravity drainage between two outputs [kg*ppt/s]

real(wp) mo_data::grav_temp

average temperature of gravity drainage brine between two outputs [T]

real(wp), dimension(:), allocatable mo_data::h

Enthalpy [J].

real(wp), dimension(:), allocatable mo_data::h_abs

specific Enthalpy [J/kg]

real(wp) mo_data::h_abs_snow

Absolute enthalpy of snow layer [J].

integer mo_data::harmonic_flag

1: minimal permeability is used to calculate Rayleigh number, 2:harmonic mean is used for Rayleigh number

integer mo_data::i

Index, normally used for time.

integer mo_data::i_time

Number of timesteps.

integer mo_data::i_time_out

Number of timesteps between each output.

integer mo_data::k

Index, normally used for layer.

integer mo_data::lab_snow_flag

Niels, 2017 add: 0: lab setup without snow covers, 1: lab setup include snow influence on heat fluxes.

integer mo_data::length_input

Sets the input length for atmoflux_flag==2, common value of 13169.

integer mo_data::length_input_lab

Niels, 2017 add: used to allocate lab testcase input arrays in mo_init, set value in testcases.

real(wp) mo_data::liquid_precip

Liquid precip, [meter of water/s].

real(wp), dimension(:), allocatable mo_data::m

Mass [kg].

real(wp) mo_data::m_snow

Mass of snow layer [kg].

real(wp) mo_data::m_total

Total initial water mass, for lab experiments with a fixed total amount.

real(wp) mo_data::melt_err =0._wp

Niels, 2017 add: used to check how much meltwater vanishes in flushing routine.

real(wp) mo_data::melt_thick

thickness of fully liquid part of top layer [m]

real(wp), dimension(3) mo_data::melt_thick_output

Niels, 2017 add: output field of surface liquid meltwater sizes.

real(wp) mo_data::melt_thick_snow
real(wp) mo_data::melt_thick_snow_old

Niels(2017) add: thickness of excess fully liquid part from snow_melt_processes [m].

integer mo_data::n_active

Number of Layers active in the present.

integer mo_data::n_bgc

Number of chemicals.

integer mo_data::n_bottom

Number of bottom layers.

integer mo_data::n_middle

Number of middle layers.

integer mo_data::n_time_out

Counts number of timesteps between output.

integer mo_data::n_top

Number of top layers.

integer mo_data::nlayer

Number of layers.

real(wp), dimension(:), allocatable mo_data::ocean_flux_input

Niels, 2017 add: used to read in oceanic heat flux for field experiment tests, dimension needs to be set in the code.

real(wp), dimension(:), allocatable mo_data::ocean_t_input

Niels, 2017 add: used to read in ocean temperature for field experiment tests, dimension needs to be set in the code.

real(wp), dimension(:), allocatable mo_data::perm
real(wp), dimension(:), allocatable mo_data::phi

Solid mass fraction.

real(wp) mo_data::phi_s

Solid mass fraction of snow layer.

integer mo_data::precip_flag

0: solid and liquid precipitation, 1:phase determined by T2m

real(wp), dimension(:), allocatable mo_data::precip_input

Used to read in precipitation from ERA for atmoflux_flag==2.

real(wp), dimension(:), allocatable mo_data::precipinput

Niels, 2017 add: used to read in precipation for field experiment tests, dimension needs to be set in the code.

integer mo_data::prescribe_flag

1: nothing happens, 2: prescribed Salinity profile is prescribed at each timestep (does not disable brine dynamics, just overwrites the salinity!)

real(wp), dimension(:), allocatable mo_data::psi_g

Gas volume fraction.

real(wp) mo_data::psi_g_snow

Gas volume fraction of snow layer.

real(wp), dimension(:), allocatable mo_data::psi_l

Liquid volume fraction.

real(wp) mo_data::psi_l_snow

Liquid volume fraction of snow layer.

real(wp), dimension(:), allocatable mo_data::psi_s

Solid volume fraction.

real(wp) mo_data::psi_s_snow

Solid volume fraction of snow layer.

real(wp), dimension(:), allocatable mo_data::q

Heat in layer [J].

real(wp), dimension(:), allocatable mo_data::ray

Rayleigh number of each layer.

real(wp), dimension(:), allocatable mo_data::s_abs

Absolute Salinity [g].

real(wp) mo_data::s_abs_snow

Absolute salinity of snow layer [g].

real(wp), dimension(:), allocatable mo_data::s_br

Brine salinity [g/kg].

real(wp), dimension(:), allocatable mo_data::s_bu

Bulk Salinity [g/kg].

real(wp) mo_data::s_bu_bottom

Salinity beneath the ice [g/kg].

real(wp) mo_data::s_total

Total initial salt mass, for lab experiments with a fixed total amount.

integer mo_data::salt_flag

1: Sea salt, 2: NaCL

integer mo_data::snow_flush_flag

Niels, 2017 add: 0: all meltwater from snow forms slush, 1: meltwater partly leads to flushing, ratio defined by "k_snow_flush".

integer mo_data::snow_precip_flag

Niels, 2017 add: 0: all precipitation is set to zero, 1: physical behaviour.

real(wp) mo_data::solid_precip

Solid precip, [meter of water /s].

integer mo_data::styropor_flag
real(wp), dimension(:), allocatable mo_data::styropor_input

Niels, 2017 add: if styropor is used in the lab on top of the ice to simulate snow heat fluxes.

real(wp) mo_data::surface_water

Percentage of water fraction in the top 5cm [%].

real(wp), dimension(:), allocatable mo_data::t

Temperature [C].

real(wp) mo_data::t2m

Two meter Temperature [C].

real(wp), dimension(:), allocatable mo_data::t2m_input

Used to read in 2Tm from ERA for atmoflux_flag==2.

real(wp) mo_data::t_bottom

Temperature of water beneath the ice [C].

real(wp) mo_data::t_freeze

Freezing temperature [C].

real(wp) mo_data::t_snow

Temperature of snow layer [C].

real(wp), save mo_data::t_test

First guess for getT subroutine.

real(wp) mo_data::t_top

Temperature at the surface [C].

real(wp) mo_data::tank_depth

water depth in meters, used to calculate concentrations below ice for tank experiments

integer mo_data::tank_flag

1: nothing, 2: S_bu_bottom and bgc_bottom are calculated as if the experiment is conducted in a tank

real(wp) mo_data::test

Thickness of snow layer [m].

real(wp), dimension(:), allocatable mo_data::thick

Layer thickness [m].

real(wp) mo_data::thick_0

Initial layer thickness [m].

real(wp) mo_data::thick_min

Parameter for snow, determines when snow is in thermal equilibrium with the ice and when it is totally neglected.

real(wp) mo_data::thick_snow
real(wp) mo_data::thickness

Meters of ice [m].

real(wp) mo_data::time

Time [s].

integer mo_data::time_counter

Keeps track of input data.

real(wp), dimension(:), allocatable mo_data::time_input

Used to read in time from ERA for atmoflux_flag==2.

real(wp) mo_data::time_out

Time between outputs [s].

real(wp) mo_data::time_total

Time of simulation [s].

real(wp), dimension(:), allocatable mo_data::tinput

Niels, 2017 add: used to read in top temperature for field experiment tests, dimension needs to be set in the code.

real(wp) mo_data::total_resist

Thermal resistance of the whole column [].

real(wp), dimension(:), allocatable mo_data::ttop_input

Niels, 2017 add: used for testcase 111, comparison with greenland harp data, uppermost harp temperature is seen as Ttop.

integer mo_data::turb_flag

1: No bottom turbulence, 2: Bottom mixing

real(wp), dimension(:), allocatable mo_data::v_ex

Volume of brine due expelled due to freezing [m^3] of solid, gas & liquid.

real(wp), dimension(:), allocatable mo_data::v_g

Volume [m^3] of gas.

real(wp), dimension(:), allocatable mo_data::v_l

Volume [m^3] of liquid.

real(wp), dimension(:), allocatable mo_data::v_s

Volume [m^3] of solid.