59 integer(i4b) :: i, j, kc
66 #if (defined(NMARS) || defined(SMARS)) /* Polar caps of Mars */ 72 #else /* all other domains */ 101 integer(i4b) :: i, j, kc
135 integer(i4b) :: i, j, kc
136 real(dp) :: kappa_const_val
137 real(dp) :: temp_ice_base
141 #if (defined(NMARS) || defined(SMARS)) /* Polar caps of Mars */ 143 #else /* all other domains */ 150 if (
maske(j,i)<=2_i2b)
then 155 + (
q_geo(j,i)/kappa_const_val) &
167 + (temp_ice_base-
temp_s(j,i)) &
179 + (temp_ice_base-
temp_s(j,i)) &
207 integer(i4b) :: i, j, kc
208 real(dp) :: kappa_const_val, c_const_val
209 real(dp) :: as_val, H_val, qgeo_val, K, z_above_base
210 real(dp) :: erf_val_1, erf_val_2
211 real(dp) :: temp_ice_base, temp_scale_factor
215 #if (defined(NMARS) || defined(SMARS)) /* Polar caps of Mars */ 217 c_const_val =
c_val(-100.0_dp)
218 #else /* all other domains */ 220 c_const_val =
c_val(-10.0_dp)
226 if (
maske(j,i)<=2_i2b)
then 232 h_val = max(
h_c(j,i) ,
eps)
235 k = sqrt( (kappa_const_val/(
rho*c_const_val)) * (h_val/as_val) )
237 erf_val_1 = erf(
h_c(j,i)/(sqrt(2.0_dp)*k))
241 erf_val_2 = erf(z_above_base/(sqrt(2.0_dp)*k))
243 + (qgeo_val/kappa_const_val) &
244 * sqrt(0.5_dp*
pi)*k*(erf_val_1-erf_val_2)
249 temp_scale_factor = (temp_ice_base-
temp_s(j,i)) &
251 else if (
maske(j,i) == 3_i2b)
then 253 temp_scale_factor = (temp_ice_base-
temp_s(j,i)) &
256 temp_scale_factor = 1.0_dp
283 integer(i4b) :: i, j, kc
308 subroutine init_temp_r()
312 integer(i4b) :: i, j, kr
328 end subroutine init_temp_r
333 subroutine init_water()
340 end subroutine init_water
345 subroutine init_age()
349 #if (defined(ASF)) /* Austfonna */ 351 age_c = 3500.0_dp*year_sec
352 age_t = 3500.0_dp*year_sec
354 #elif (defined(NMARS) || defined(SMARS)) /* Polar caps of Mars */ 356 age_c = 1.0e+06_dp*year_sec
357 age_t = 1.0e+06_dp*year_sec
359 #else /* all other domains */ 361 age_c = 15000.0_dp*year_sec
362 age_t = 15000.0_dp*year_sec
366 end subroutine init_age
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) temp_c
temp_c(kc,j,i): Temperature in the upper (kc) ice domain
real(dp) kappa_r
KAPPA_R: Heat conductivity of the lithosphere.
real(dp) delta_tm_sw
DELTA_TM_SW: Melting point depression of sea water due to its average salinity.
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) age_c
age_c(kc,j,i): Age in the upper (kc) ice domain
real(dp), parameter epsi
epsi: Very small number
real(dp), dimension(0:jmax, 0:imax) q_geo
q_geo(j,i): Geothermal heat flux
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) age_t
age_t(kt,j,i): Age in the lower (kt) ice domain
real(dp), dimension(0:krmax, 0:jmax, 0:imax) temp_r
temp_r(kr,j,i): Temperature in the bedrock
real(dp) h_r
H_R: Thickness of the modelled lithosphere layer.
Material properties of ice: Rate factor, heat conductivity, specific heat (heat capacity), creep function, viscosity.
real(dp), parameter eps
eps: Small number
subroutine, public init_temp_water_age_2()
Initial temperature, water content and age (case ANF_DAT==2: ice-free conditions with relaxed bedrock...
real(dp) function, public kappa_val(temp_val)
Heat conductivity of ice: Linear interpolation of tabulated values in KAPPA(.).
integer(i2b), dimension(0:jmax, 0:imax) maske
maske(j,i): Ice-land-ocean mask. 0: grounded ice, 1: ice-free land, 2: ocean, 3: floating ice ...
real(dp), dimension(0:jmax, 0:imax) as_perp
as_perp(j,i): Accumulation-ablation function at the ice surface (SMB)
subroutine, public init_temp_water_age_1_2()
Initial temperature, water content and age (case ANF_DAT==1, TEMP_INIT==2: present-day initial topogr...
real(dp), dimension(0:jmax, 0:imax) temp_s
temp_s(j,i): Ice surface temperature
Declarations of global variables for SICOPOLIS (for the ANT domain).
Initial temperature, water content and age.
Declarations of kind types for SICOPOLIS.
real(dp), dimension(0:jmax, 0:imax) h_c
H_c(j,i): Thickness of ice in the upper (kc) domain (thickness of the cold-ice layer for POLY...
real(dp), dimension(0:krmax) zeta_r
zeta_r(kr): Sigma coordinate zeta_r of grid point kr
real(dp), dimension(0:kcmax, 0:jmax, 0:imax) omega_c
omega_c(kc,j,i): Water content in the upper (kc) ice domain
real(dp) function, public c_val(temp_val)
Specific heat of ice: Linear interpolation of tabulated values in C(.).
real(dp), dimension(0:kcmax) eaz_c_quotient
eaz_c_quotient(kc): Abbreviation for (eaz_c(kc)-1.0)/(ea-1.0)
subroutine, public init_temp_water_age_1_3()
Initial temperature, water content and age (case ANF_DAT==1, TEMP_INIT==3: present-day initial topogr...
real(dp) beta
BETA: Clausius-Clapeyron gradient of ice.
real(dp), parameter pi
pi: Constant pi
real(dp), dimension(0:ktmax, 0:jmax, 0:imax) omega_t
omega_t(kt,j,i): Water content in the lower (kt) ice domain
real(dp) rho
RHO: Density of ice.
subroutine, public init_temp_water_age_1_4()
Initial temperature, water content and age (case ANF_DAT==1, TEMP_INIT==4: present-day initial topogr...
subroutine, public init_temp_water_age_1_1()
Initial temperature, water content and age (case ANF_DAT==1, TEMP_INIT==1: present-day initial topogr...
Declarations of global variables for SICOPOLIS.