SICOPOLIS V3.3
|
Declarations of global variables for SICOPOLIS. More...
Go to the source code of this file.
Modules | |
module | sico_variables_m |
Declarations of global variables for SICOPOLIS. | |
Variables | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::maske |
maske(j,i): Ice-land-ocean mask. 0: grounded ice, 1: ice-free land, 2: ocean, 3: floating ice More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::maske_old |
maske_old(j,i): Old value of maske (at the previous time step) More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::maske_neu |
maske_neu(j,i): New value of maske computed during an integration step More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::mask_run |
mask_run(j,i): mask indicating melt type. 2: visible (ocean, for later developments), 1: visible (grounded ice), -1: hidden on land, -2: hidden in ocean More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::n_cts |
n_cts(j,i): Mask for thermal conditions. -1: cold ice base, 0: temperate ice base with cold ice above, 1: temperate ice base with temperate ice layer above (only for POLY) More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::n_cts_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::kc_cts |
kc_cts(j,i): Position kc of the CTS (for COLD, ENTC, ENTM) More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::kc_cts_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_inner_point |
flag_inner_point(j,i): Inner-point flag. .true.: inner point, .false.: margin point More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_grounding_line_1 |
flag_grounding_line_1(j,i): Grounding line flag. .true.: grounding line point (grounded ice point with at least one floating ice neighbour), .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_grounding_line_2 |
flag_grounding_line_2(j,i): Grounding line flag. .true.: grounding line point (floating ice point with at least one grounded ice neighbour), .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_calving_front_1 |
flag_calving_front_1(j,i): Calving front flag. .true.: calving front point (floating ice point with at least one ocean neighbour), .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_calving_front_2 |
flag_calving_front_2(j,i): Calving front flag. .true.: calving front point (ocean point with at least one floating ice neighbour), .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_shelfy_stream_x |
flag_shelfy_stream_x(j,i): Shelfy stream flag in x-direction, at (i+1/2,j). .true.: shelfy stream point .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_shelfy_stream_y |
flag_shelfy_stream_y(j,i): Shelfy stream flag in y-direction, at (i,j+1/2). .true.: shelfy stream point .false.: otherwise More... | |
logical, dimension(0:jmax, 0:imax) | sico_variables_m::flag_shelfy_stream |
flag_shelfy_stream(j,i): Shelfy stream flag on the main grid. .true.: grounded ice, and at least one neighbour on the staggered grid is a shelfy stream point .false.: otherwise More... | |
real(dp), dimension(0:imax) | sico_variables_m::xi |
xi(i): Coordinate xi (= x) of grid point i More... | |
real(dp), dimension(0:jmax) | sico_variables_m::eta |
eta(j): Coordinate eta (= y) of grid point j More... | |
real(dp), dimension(0:kcmax) | sico_variables_m::zeta_c |
zeta_c(kc): Sigma coordinate zeta_c of grid point kc More... | |
real(dp), dimension(0:ktmax) | sico_variables_m::zeta_t |
zeta_t(kt): Sigma coordinate zeta_t of grid point kt More... | |
real(dp), dimension(0:krmax) | sico_variables_m::zeta_r |
zeta_r(kr): Sigma coordinate zeta_r of grid point kr More... | |
real(dp) | sico_variables_m::aa |
aa: Exponential stretch parameter of the non-equidistant vertical grid in the upper (kc) ice domain More... | |
logical | sico_variables_m::flag_aa_nonzero |
flag_aa_nonzero: Flag for the exponential stretch parameter aa. .true.: aa greater than zero (non-equidistant grid) .false.: aa equal to zero (equidistant grid) More... | |
real(dp) | sico_variables_m::ea |
ea: Abbreviation for exp(aa) More... | |
real(dp), dimension(0:kcmax) | sico_variables_m::eaz_c |
eaz_c(kc): Abbreviation for exp(aa*zeta(kc)) More... | |
real(dp), dimension(0:kcmax) | sico_variables_m::eaz_c_quotient |
eaz_c_quotient(kc): Abbreviation for (eaz_c(kc)-1.0)/(ea-1.0) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::lambda |
lambda(j,i): Geographic longitude of grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::phi |
phi(j,i): Geographic latitude of grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::area |
area(j,i): Area of grid cell associated with grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g11_g |
sq_g11_g(j,i): Square root of the coefficient g11 of the metric tensor on grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g22_g |
sq_g22_g(j,i): Square root of the coefficient g22 of the metric tensor on grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::insq_g11_g |
insq_g11_g(j,i): Inverse square root of g11 on grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::insq_g22_g |
insq_g22_g(j,i): Inverse square root of g22 on grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g11_sgx |
sq_g11_sgx(j,i): Square root of g11, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g11_sgy |
sq_g11_sgy(j,i): Square root of g11, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g22_sgx |
sq_g22_sgx(j,i): Square root of g22, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::sq_g22_sgy |
sq_g22_sgy(j,i): Square root of g22, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::insq_g11_sgx |
insq_g11_sgx(j,i): Inverse square root of g11, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::insq_g22_sgy |
insq_g22_sgy(j,i): Inverse square root of g22, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zs |
zs(j,i): Coordinate z of the surface topography More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zm |
zm(j,i): Coordinate z of the bottom of the upper (kc) ice domain = top of the lower (kt) ice domain (position of the CTS for POLY, equal to zb for ISOT, COLD, ENTC, ENTM) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zb |
zb(j,i): Coordinate z of the ice base More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zl |
zl(j,i): Coordinate z of the lithosphere surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zl0 |
zl0(j,i): zl for isostatically relaxed ice-free conditions More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::wss |
wss(j,i): Isostatic steady-state displacement of the lithosphere More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::flex_rig_lith |
flex_rig_lith(j,i): Flexural rigidity of the lithosphere More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::time_lag_asth |
time_lag_asth(j,i): Time lag of the relaxing asthenosphere More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_c |
H_c(j,i): Thickness of ice in the upper (kc) domain (thickness of the cold-ice layer for POLY, entire ice thickness for ISOT, COLD, ENTC, ENTM) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_t |
H_t(j,i): Thickness of ice in the lower (kt) domain (thickness of the temperate layer for POLY, redundant and thus set to zero for ISOT, COLD, ENTC, ENTM) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzs_dxi |
dzs_dxi(j,i): Derivative of zs by xi (at i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzm_dxi |
dzm_dxi(j,i): Derivative of zm by xi (at i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzb_dxi |
dzb_dxi(j,i): Derivative of zb by xi (at i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_c_dxi |
dH_c_dxi(j,i): Derivative of H_c by xi (at i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_t_dxi |
dH_t_dxi(j,i): Derivative of H_t by xi (at i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzs_deta |
dzs_deta(j,i): Derivative of zs by eta (at i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzm_deta |
dzm_deta(j,i): Derivative of zm by eta (at i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzb_deta |
dzb_deta(j,i): Derivative of zb by eta (at i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_c_deta |
dH_c_deta(j,i): Derivative of H_c by eta (at i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_t_deta |
dH_t_deta(j,i): Derivative of H_t by eta (at i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzs_dxi_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzm_dxi_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzb_dxi_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_c_dxi_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_t_dxi_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzs_deta_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzm_deta_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzb_deta_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_c_deta_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_t_deta_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzs_dtau |
dzs_dtau(j,i): Derivative of zs by tau (time) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzm_dtau |
dzm_dtau(j,i): Derivative of zm by tau (time) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzb_dtau |
dzb_dtau(j,i): Derivative of zb by tau (time) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dzl_dtau |
dzl_dtau(j,i): Derivative of zl by tau (time) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_c_dtau |
dH_c_dtau(j,i): Derivative of H_c by tau (time) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::dh_t_dtau |
dH_t_dtau(j,i): Derivative of H_t by tau (time) More... | |
integer(i4b), dimension(0:jmax, 0:imax) | sico_variables_m::p_weert |
p_weert(j,i): Weertman exponent for the basal shear stress More... | |
integer(i4b), dimension(0:jmax, 0:imax) | sico_variables_m::q_weert |
q_weert(j,i): Weertman exponent for the basal pressure More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::p_weert_inv |
p_weert_inv(j,i): Inverse of p_weert More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::c_slide |
c_slide(j,i): Basal sliding coefficient More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::d_help_b |
d_help_b(j,i): Auxiliary quantity for the computation of vx_b and vy_b More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::c_drag |
c_drag(j,i): Auxiliary quantity for the computation of the basal drag More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::p_b_w |
p_b_w(j,i): Basal water pressure More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vx_b |
vx_b(j,i): Velocity in x-direction at the ice base, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vy_b |
vy_b(j,i): Velocity in y-direction at the ice base, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vx_m |
vx_m(j,i): Mean (depth-averaged) velocity in x-direction, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vy_m |
vy_m(j,i): Mean (depth-averaged) velocity in y-direction, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::ratio_sl_x |
ratio_sl_x(j,i): Ratio of basal to surface velocity (slip ratio) in x-direction, at (i+1/2,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::ratio_sl_y |
ratio_sl_y(j,i): Ratio of basal to surface velocity (slip ratio) in y-direction, at (i,j+1/2) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vx_b_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vy_b_g |
(.)_g(j,i): Staggered-grid quantity (.) interpolated to grid point (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vz_b |
vz_b(j,i): Velocity in z-direction at the ice base, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vz_m |
vz_m(j,i): Velocity in z-direction at the position z=zm (interface between the upper (kc) and the lower (kt) domain), at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vx_s_g |
vx_s_g(j,i): Velocity in x-direction at the ice surface, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vy_s_g |
vy_s_g(j,i): Velocity in x-direction at the ice surface, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vz_s |
vz_s(j,i): Velocity in z-direction at the ice surface, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::flui_ave_sia |
flui_ave_sia(j,i): Depth-averaged fluidity of the SIA More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_diff |
h_diff(j,i): Diffusivity of the SIA evolution equation of the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::qx |
qx(j,i): Volume flux in x-direction (depth-integrated vx, at (i+1/2,j)) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::qy |
qy(j,i): Volume flux in y-direction (depth-integrated vy, at (i,j+1/2)) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::q_gl_g |
q_gl_g(j,i): Volume flux across the grounding line, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::q_geo |
q_geo(j,i): Geothermal heat flux More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_b |
temp_b(j,i): Basal temperature More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temph_b |
temph_b(j,i): Basal temperature relative to the pressure melting point More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::q_bm |
Q_bm(j,i): Basal melting rate. More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::q_tld |
Q_tld(j,i): Water drainage rate from the temperate layer. More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::q_b_tot |
Q_b_tot(j,i): Sum of Q_bm and Q_tld. More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_w |
H_w(j,i): Thickness of the water column under the ice base. More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::accum |
accum(j,i): Accumulation rate at the ice surface (includes liquid precipitation = rainfall!) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::evap |
evap(j,i): Evaporation rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::runoff |
runoff(j,i): Runoff rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::as_perp |
as_perp(j,i): Accumulation-ablation function at the ice surface (SMB) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::as_perp_apl |
as_perp_apl(j,i): Applied accumulation-ablation function (SMB) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::mb_source_apl |
mb_source_apl(j,i): Applied mass balance source (SMB, BMB, calving) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::runoff_top |
runoff_top(j,i): Runoff rate (accounting at ice surface) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::runoff_bot |
runoff_bot(j,i): Runoff rate (accounting at ice bottom) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::disc_top |
disc_top(j,i): Ice discharge rate (accounting at ice surface) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::calv_grounded |
calv_grounded(j,i): Calving rate of grounded ice More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_s |
temp_s(j,i): Ice surface temperature More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::am_perp |
am_perp(j,i): Ice volume flux across the z=zm interface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::am_perp_st |
am_perp_st(j,i): Steady-state part of am_perp (without contribution of dzm_dtau) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zs_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zm_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zb_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zl_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_c_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_t_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zs_ref |
zs_ref(j,i): Reference elevation for precip_present, temp_ma_present and temp_mj_present More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::accum_present |
accum_present(j,i): Present-day accumulation rate at the ice surface (for EISMINT, ISMIP HEINO and the north and south polar caps of Mars) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::precip_ma_present |
precip_ma_present(j,i): Present-day mean annual precipitation rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::precip_ma_lgm_anom |
precip_ma_lgm_anom(j,i): LGM anomaly (ratio LGM/present) of the mean annual precipitation rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_ma_present |
temp_ma_present(j,i): Present-day mean annual surface temperature More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_mj_present |
temp_mj_present(j,i): Present-day mean summer (northern hemisphere: July, southern hemisphere: January) surface temperature More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_ma_lgm_anom |
temp_ma_lgm_anom(j,i): LGM anomaly (difference LGM - present) of the mean annual surface temperature More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::temp_mj_lgm_anom |
temp_mj_lgm_anom(j,i): LGM anomaly (difference LGM - present) of the mean summer (northern hemisphere: July, southern hemisphere: January) surface temperature More... | |
real(dp), dimension(-jmax:jmax,-imax:imax) | sico_variables_m::dist_dxdy |
dist_dxdy(jr,ir): Distance between grid points with delta_i=ir, delta_j=jr More... | |
real(dp), dimension(0:jmax, 0:imax, 12) | sico_variables_m::precip_present |
precip_present(j,i,n): Present-day mean monthly precipitation rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax, 12) | sico_variables_m::precip_lgm_anom |
precip_lgm_anom(j,i,n): LGM anomaly (ratio LGM/present) of the mean monthly precipitation rate at the ice surface More... | |
real(dp), dimension(0:jmax, 0:imax, 12) | sico_variables_m::gamma_precip_lgm_anom |
gamma_precip_lgm_anom(j,i,n): negative natural logarithm of precip_lgm_anom(j,i,n) More... | |
real(dp), dimension(0:jmax, 0:imax, 12) | sico_variables_m::temp_mm_present |
temp_mm_present(j,i,n): Present-day mean monthly surface temperature More... | |
real(dp), dimension(0:jmax, 0:imax, 12) | sico_variables_m::temp_mm_lgm_anom |
temp_mm_lgm_anom(j,i,n): LGM anomaly (difference LGM - present) of the mean monthly surface temperature More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::d_help_c |
d_help_c(kc,j,i): Auxiliary quantity for the computation of vx, vy und zs More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::vx_c |
vx_c(kc,j,i): Velocity in x-direction in the upper (kc) ice domain (at (i+1/2,j,kc)) More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::vy_c |
vy_c(kc,j,i): Velocity in y-direction in the upper (kc) ice domain (at (i,j+1/2,kc)) More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::vz_c |
vz_c(kc,j,i): Velocity in z-direction in the upper (kc) ice domain (at (i,j,kc+1/2)) More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::temp_c |
temp_c(kc,j,i): Temperature in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::temp_c_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::temp_c_m |
temp_c_m(kc,j,i): Melting temperature in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::age_c |
age_c(kc,j,i): Age in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::age_c_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::txz_c |
txz_c(kc,j,i): Shear stress txz in the upper (kc) ice domain (at (i+1/2,j,kc)) More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::tyz_c |
tyz_c(kc,j,i): Shear stress tyz in the upper (kc) ice domain (at (i,j+1/2,kc)) More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::sigma_c |
sigma_c(kc,j,i): Effective stress in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::enh_c |
enh_c(kc,j,i): Flow enhancement factor in the upper (kc) ice domain More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::de_ssa |
de_ssa(j,i): Effective strain rate of the SSA, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vis_int_g |
vis_int_g(j,i): Depth-integrated viscosity of the SSA, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vx_g |
vx_g(j,i): Velocity in x-direction of the SSA, at (i,j) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::vy_g |
vy_g(j,i): Velocity in y-direction of the SSA, at (i,j) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::d_help_t |
d_help_t(kt,j,i): Auxiliary quantity for the computation of vx, vy und zs More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::vx_t |
vx_t(kt,j,i): Velocity in x-direction in the lower (kt) ice domain (at (i+1/2,j,kt)) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::vy_t |
vy_t(kt,j,i): Velocity in y-direction in the lower (kt) ice domain (at (i,j+1/2,kt)) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::vz_t |
vz_t(kt,j,i): Velocity in z-direction in the lower (kt) ice domain (at (i,j,kt+1/2)) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::omega_t |
omega_t(kt,j,i): Water content in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::omega_t_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::temp_t_m |
temp_t_m(kt,j,i): Melting temperature in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::age_t |
age_t(kt,j,i): Age in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::age_t_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::txz_t |
txz_t(kt,j,i): Shear stress txz in the lower (kt) ice domain (at (i+1/2,j,kt)) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::tyz_t |
tyz_t(kt,j,i): Shear stress tyz in the lower (kt) ice domain (at (i,j+1/2,kt)) More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::sigma_t |
sigma_t(kt,j,i): Effective stress in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::enh_t |
enh_t(kt,j,i): Flow enhancement factor in the lower (kt) ice domain More... | |
real(dp), dimension(0:krmax, 0:jmax, 0:imax) | sico_variables_m::temp_r |
temp_r(kr,j,i): Temperature in the bedrock More... | |
real(dp), dimension(0:krmax, 0:jmax, 0:imax) | sico_variables_m::temp_r_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::enth_c |
enth_c(kc,j,i): Enthalpy in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::enth_c_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::omega_c |
omega_c(kc,j,i): Water content in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::omega_c_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::enth_t |
enth_t(kt,j,i): Enthalpy in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::enth_t_neu |
(.)_neu: New value of quantity (.) computed during an integration step More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::dxx_c |
dxx_c(kc,j,i): Strain rate dxx in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::dyy_c |
dyy_c(kc,j,i): Strain rate dyy in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::dxy_c |
dxy_c(kc,j,i): Strain rate dxy in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::dxz_c |
dxz_c(kc,j,i): Strain rate dxz in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::dyz_c |
dyz_c(kc,j,i): Strain rate dyz in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::de_c |
de_c(kc,j,i): Full effective strain rate in the upper (kc) ice domain More... | |
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) | sico_variables_m::lambda_shear_c |
lambda_shear_c(kc,j,i): Shear fraction in the upper (kc) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::dxx_t |
dxx_t(kt,j,i): Strain rate dxx in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::dyy_t |
dyy_t(kt,j,i): Strain rate dyy in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::dxy_t |
dxy_t(kt,j,i): Strain rate dxy in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::dxz_t |
dxz_t(kt,j,i): Strain rate dxz in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::dyz_t |
dyz_t(kt,j,i): Strain rate dyz in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::de_t |
de_t(kt,j,i): Full effective strain rate in the lower (kt) ice domain More... | |
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) | sico_variables_m::lambda_shear_t |
lambda_shear_t(kt,j,i): Shear fraction in the lower (kt) ice domain More... | |
real(dp) | sico_variables_m::rho |
RHO: Density of ice. More... | |
real(dp) | sico_variables_m::rho_w |
RHO_W: Density of pure water. More... | |
real(dp) | sico_variables_m::rho_sw |
RHO_SW: Density of sea water. More... | |
real(dp) | sico_variables_m::l |
L: Latent heat of ice. More... | |
real(dp) | sico_variables_m::g |
G: Acceleration due to gravity. More... | |
real(dp) | sico_variables_m::nue |
NUE: Water diffusivity in ice. More... | |
real(dp) | sico_variables_m::beta |
BETA: Clausius-Clapeyron gradient of ice. More... | |
real(dp) | sico_variables_m::delta_tm_sw |
DELTA_TM_SW: Melting point depression of sea water due to its average salinity. More... | |
real(dp) | sico_variables_m::omega_max |
OMEGA_MAX: Threshold value for the water content of temperate ice. More... | |
real(dp) | sico_variables_m::h_r |
H_R: Thickness of the modelled lithosphere layer. More... | |
real(dp) | sico_variables_m::rho_c_r |
RHO_C_R: Density times specific heat of the lithosphere. More... | |
real(dp) | sico_variables_m::kappa_r |
KAPPA_R: Heat conductivity of the lithosphere. More... | |
real(dp) | sico_variables_m::rho_a |
RHO_A: Density of the asthenosphere. More... | |
real(dp) | sico_variables_m::r_t |
R_T: Coefficient of the water-content dependence in the rate factor for temperate ice. More... | |
real(dp) | sico_variables_m::r |
R: Radius of the planet. More... | |
real(dp) | sico_variables_m::a |
A: Semi-major axis of the planet. More... | |
real(dp) | sico_variables_m::b |
B: Semi-minor axis of the planet. More... | |
real(dp) | sico_variables_m::lambda0 |
LAMBDA0: Reference longitude (central meridian) of the stereographic projection. More... | |
real(dp) | sico_variables_m::phi0 |
PHI0: Standard parallel of the stereographic projection. More... | |
real(dp) | sico_variables_m::s_stat_0 |
S_STAT_0: Standard deviation of the air termperature for the degree-day model. More... | |
real(dp) | sico_variables_m::beta1_0 |
BETA1_0: Degree-day factor for snow. More... | |
real(dp) | sico_variables_m::beta1_lt_0 |
BETA1_LT_0: Degree-day factor for snow at low summer temperatures. More... | |
real(dp) | sico_variables_m::beta1_ht_0 |
BETA1_HT_0: Degree-day factor for snow at high summer temperatures. More... | |
real(dp) | sico_variables_m::beta2_0 |
BETA2_0: Degree-day factor for ice. More... | |
real(dp) | sico_variables_m::beta2_lt_0 |
BETA2_LT_0: Degree-day factor for ice at low summer temperatures. More... | |
real(dp) | sico_variables_m::beta2_ht_0 |
BETA2_HT_0: Degree-day factor for ice at high summer temperatures. More... | |
real(dp) | sico_variables_m::phi_sep_0 |
PHI_SEP_0: Separation latitude for the computation of the degree-day factors beta1 and beta2: Equatorward of phi_sep, only the high-temperature values are used, whereas poleward of phi_sep, beta1 and beta2 are temperature-dependent. More... | |
real(dp) | sico_variables_m::pmax_0 |
PMAX_0: Saturation factor for the formation of superimposed ice. More... | |
real(dp) | sico_variables_m::mu_0 |
MU_0: Firn-warming correction. More... | |
real(dp), dimension(-190:10) | sico_variables_m::rf |
RF(n): Tabulated values for the rate factor of cold ice. More... | |
real(dp), dimension(-190:10) | sico_variables_m::kappa |
KAPPA(n): Tabulated values for the heat conductivity of ice. More... | |
real(dp), dimension(-190:10) | sico_variables_m::c |
C(n): Tabulated values for the specific heat of ice. More... | |
real(dp), parameter | sico_variables_m::pi = 3.141592653589793_dp |
pi: Constant pi More... | |
real(dp), parameter | sico_variables_m::pi_inv = 1.0_dp/pi |
pi_inv: Inverse of pi More... | |
real(dp), parameter | sico_variables_m::pi_180 = pi/180.0_dp |
pi_180: pi divided by 180 (-> deg to rad) More... | |
real(dp), parameter | sico_variables_m::pi_180_inv = 180.0_dp/pi |
pi_180_inv: 180 divided by pi (-> rad to deg) More... | |
real(dp), parameter | sico_variables_m::euler = 2.718281828459045_dp |
euler: Euler number More... | |
real(dp), parameter | sico_variables_m::eps = 1.0e-05_dp |
eps: Small number More... | |
real(dp), parameter | sico_variables_m::epsi = 1.0e-12_dp |
epsi: Very small number More... | |
real(sp), parameter | sico_variables_m::eps_sp = epsilon(1.0_sp) |
eps_sp: Small number to single-precision accuracy More... | |
real(dp), parameter | sico_variables_m::eps_sp_dp = real(eps_sp, dp) |
eps_sp_dp: Small number to single-precision accuracy in double precision More... | |
real(dp), parameter | sico_variables_m::eps_dp = epsilon(1.0_dp) |
eps_dp: Small number to double-precision accuracy More... | |
real(dp) | sico_variables_m::year_zero |
year_zero: SICOPOLIS year zero in astronomical year numbering [ = signed year CE (AD) ] More... | |
character(len=64) | sico_variables_m::ch_domain_long |
ch_domain_long: Long name of the computational domain More... | |
character(len=16) | sico_variables_m::ch_domain_short |
ch_domain_short: Short name of the computational domain More... | |
integer(i2b) | sico_variables_m::forcing_flag |
forcing_flag: Flag for the forcing type. 1: forcing by a spatially constant surface temperature anomaly (delta_ts), 2: forcing by a glacial index (glac_index), 3: forcing by time-dependent surface temperature and precipitation data. More... | |
integer(i4b) | sico_variables_m::n_core |
n_core: Number of positions to be considered in the time-series file for deep boreholes More... | |
real(dp), dimension(:), allocatable | sico_variables_m::lambda_core |
lambda_core(n): Geographical longitude of the prescribed borehole positions More... | |
real(dp), dimension(:), allocatable | sico_variables_m::phi_core |
phi_core(n): Geographical latitude of the prescribed borehole positions More... | |
real(dp), dimension(:), allocatable | sico_variables_m::x_core |
x_core(n): Coordinate xi (= x) of the prescribed borehole positions More... | |
real(dp), dimension(:), allocatable | sico_variables_m::y_core |
y_core(n): Coordinate eta (= y) of the prescribed borehole positions More... | |
character(len=16), dimension(:), allocatable | sico_variables_m::ch_core |
ch_core(n): Names of the prescribed borehole positions More... | |
integer(i4b) | sico_variables_m::grip_time_min |
grip_time_min: Minimum time of the data values for the surface temperature anomaly More... | |
integer(i4b) | sico_variables_m::grip_time_stp |
grip_time_stp: Time step of the data values for the surface temperature anomaly More... | |
integer(i4b) | sico_variables_m::grip_time_max |
grip_time_max: Maximum time of the data values for the surface temperature anomaly More... | |
integer(i4b) | sico_variables_m::ndata_grip |
ndata_grip: Number of data values for the surface temperature anomaly More... | |
real(dp), dimension(:), allocatable | sico_variables_m::griptemp |
griptemp(n): Data values for the surface temperature anomaly More... | |
integer(i4b) | sico_variables_m::gi_time_min |
gi_time_min: Minimum time of the data values for the glacial index More... | |
integer(i4b) | sico_variables_m::gi_time_stp |
gi_time_stp: Time step of the data values for the glacial index More... | |
integer(i4b) | sico_variables_m::gi_time_max |
gi_time_max: Maximum time of the data values for the glacial index More... | |
integer(i4b) | sico_variables_m::ndata_gi |
ndata_gi: Number of data values for the glacial index More... | |
real(dp), dimension(:), allocatable | sico_variables_m::glacial_index |
glacial_index(n): Data values for the glacial index More... | |
integer(i4b) | sico_variables_m::specmap_time_min |
specmap_time_min: Minimum time of the data values for the sea level More... | |
integer(i4b) | sico_variables_m::specmap_time_stp |
specmap_time_stp: Time step of the data values for the sea level More... | |
integer(i4b) | sico_variables_m::specmap_time_max |
specmap_time_max: Maximum time of the data values for the sea level More... | |
integer(i4b) | sico_variables_m::ndata_specmap |
ndata_specmap: Number of data values for the sea level More... | |
real(dp), dimension(:), allocatable | sico_variables_m::specmap_zsl |
specmap_zsl(n): Data values for the sea level More... | |
real(dp) | sico_variables_m::time_target_topo_init |
time_target_topo_init: Initial time for target-topography adjustment More... | |
real(dp) | sico_variables_m::time_target_topo_final |
time_target_topo_final: Final time for target-topography adjustment More... | |
real(dp) | sico_variables_m::target_topo_tau |
target_topo_tau: Relaxation time for target-topography adjustment More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::maske_target |
maske_target(j,i): Target topography (ice-land-ocean mask) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zs_target |
zs_target(j,i): Target topography (ice surface) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zb_target |
zb_target(j,i): Target topography (ice base) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::zl_target |
zl_target(j,i): Target topography (lithosphere surface) More... | |
real(dp), dimension(0:jmax, 0:imax) | sico_variables_m::h_target |
H_target(j,i): Target topography (ice thickness) More... | |
integer(i2b), dimension(0:jmax, 0:imax) | sico_variables_m::maske_maxextent |
maske_maxextent(j,i): Maximum ice extent mask. 0: not allowed to glaciate, 1: allowed to glaciate. More... | |
integer(i4b) | sico_variables_m::ncid_ser |
ncid_ser: ID of the NetCDF time-series output file More... | |
integer(i4b) | sico_variables_m::ncid_core |
ncid_core: ID of the NetCDF time-series output file for the deep ice cores More... | |
integer(i4b) | sico_variables_m::ncid_temp_precip |
ncid_temp_precip: ID of the NetCDF file containing the surface-temperature and precipitation data as functions of time More... | |
integer(i4b) | sico_variables_m::ndata_temp_precip |
ndata_temp_precip: Number of surface-temperature and precipitation data More... | |
real(dp) | sico_variables_m::temp_precip_time_min |
temp_precip_time_min: Minimum time of the surface-temperature and precipitation data More... | |
real(dp) | sico_variables_m::temp_precip_time_stp |
temp_precip_time_stp: Time step of the surface-temperature and precipitation data More... | |
real(dp) | sico_variables_m::temp_precip_time_max |
temp_precip_time_max: Maximum time of the surface-temperature and precipitation data More... | |
real(dp), dimension(:), allocatable | sico_variables_m::temp_precip_time |
temp_precip_time(n): Times of the surface-temperature and precipitation data More... | |
real(dp), dimension(-10000:10000) | sico_variables_m::kei |
kei(n): Tabulated values of the kei function (Kelvin function of zero order) More... | |
integer(i4b) | sico_variables_m::n_data_kei |
n_data_kei: Number of tabulated values of the kei function More... | |
real(dp) | sico_variables_m::kei_r_max |
kei_r_max: Maximum value of the argument r of the tabulated kei function More... | |
real(dp) | sico_variables_m::kei_r_incr |
kei_r_incr: Increment of the argument r of the tabulated kei function More... | |
integer(i4b), parameter | sico_variables_m::rcl1 = 3*8*(IMAX+1) |
rcl1: Maximum length of record for input files (with factor 3 safety margin) More... | |
integer(i4b), parameter | sico_variables_m::rcl2 = 3 *(IMAX+1) |
rcl2: Maximum length of record for input mask files (with factor 3 safety margin) More... | |
real(dp), parameter | sico_variables_m::no_value_pos_1 = 1.11e+11_dp |
no_value_pos_1: Positive no-value parameter More... | |
real(dp), parameter | sico_variables_m::no_value_pos_2 = 9.999e+03_dp |
no_value_pos_2: Positive no-value parameter More... | |
real(dp), parameter | sico_variables_m::no_value_neg_1 = -1.11e+11_dp |
no_value_neg_1: Negative no-value parameter More... | |
real(dp), parameter | sico_variables_m::no_value_neg_2 = -9.999e+03_dp |
no_value_neg_2: Negative no-value parameter More... | |
Declarations of global variables for SICOPOLIS.
Copyright 2009-2017 Ralf Greve
This file is part of SICOPOLIS.
SICOPOLIS 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.
SICOPOLIS 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 SICOPOLIS. If not, see http://www.gnu.org/licenses/.
Definition in file sico_variables_m.F90.