tpism_config.cdl - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch) (HTM) git clone git://src.adamsgaard.dk/pism (DIR) Log (DIR) Files (DIR) Refs (DIR) LICENSE --- tpism_config.cdl (179008B) --- 1 netcdf pism_config { 2 variables: 3 byte pism_config; 4 5 pism_config:age.enabled = "no"; 6 pism_config:age.enabled_doc = "Solve age equation (advection equation for ice age)."; 7 pism_config:age.enabled_option = "age"; 8 pism_config:age.enabled_type = "flag"; 9 10 pism_config:age.initial_value = 0.0; 11 pism_config:age.initial_value_doc = "Initial age of ice"; 12 pism_config:age.initial_value_type = "number"; 13 pism_config:age.initial_value_units = "years"; 14 15 pism_config:atmosphere.anomaly.file = ""; 16 pism_config:atmosphere.anomaly.file_doc = "Name of the file containing climate forcing fields."; 17 pism_config:atmosphere.anomaly.file_option = "atmosphere_anomaly_file"; 18 pism_config:atmosphere.anomaly.file_type = "string"; 19 20 pism_config:atmosphere.anomaly.period = 0; 21 pism_config:atmosphere.anomaly.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 22 pism_config:atmosphere.anomaly.period_option = "atmosphere_anomaly_period"; 23 pism_config:atmosphere.anomaly.period_type = "integer"; 24 pism_config:atmosphere.anomaly.period_units = "years"; 25 26 pism_config:atmosphere.anomaly.reference_year = 0; 27 pism_config:atmosphere.anomaly.reference_year_doc = "Reference year to use when ``atmosphere.anomaly.period`` is active."; 28 pism_config:atmosphere.anomaly.reference_year_option = "atmosphere_anomaly_reference_year"; 29 pism_config:atmosphere.anomaly.reference_year_type = "integer"; 30 pism_config:atmosphere.anomaly.reference_year_units = "years"; 31 32 pism_config:atmosphere.delta_P.file = ""; 33 pism_config:atmosphere.delta_P.file_doc = "Name of the file containing scalar precipitation offsets."; 34 pism_config:atmosphere.delta_P.file_option = "atmosphere_delta_P_file"; 35 pism_config:atmosphere.delta_P.file_type = "string"; 36 37 pism_config:atmosphere.delta_P.period = 0; 38 pism_config:atmosphere.delta_P.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 39 pism_config:atmosphere.delta_P.period_option = "atmosphere_delta_P_period"; 40 pism_config:atmosphere.delta_P.period_type = "integer"; 41 pism_config:atmosphere.delta_P.period_units = "years"; 42 43 pism_config:atmosphere.delta_P.reference_year = 0; 44 pism_config:atmosphere.delta_P.reference_year_doc = "Reference year to use when ``atmosphere.delta_P.period`` is active."; 45 pism_config:atmosphere.delta_P.reference_year_option = "atmosphere_delta_P_reference_year"; 46 pism_config:atmosphere.delta_P.reference_year_type = "integer"; 47 pism_config:atmosphere.delta_P.reference_year_units = "years"; 48 49 pism_config:atmosphere.delta_T.file = ""; 50 pism_config:atmosphere.delta_T.file_doc = "Name of the file containing temperature offsets."; 51 pism_config:atmosphere.delta_T.file_option = "atmosphere_delta_T_file"; 52 pism_config:atmosphere.delta_T.file_type = "string"; 53 54 pism_config:atmosphere.delta_T.period = 0; 55 pism_config:atmosphere.delta_T.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 56 pism_config:atmosphere.delta_T.period_option = "atmosphere_delta_T_period"; 57 pism_config:atmosphere.delta_T.period_type = "integer"; 58 pism_config:atmosphere.delta_T.period_units = "years"; 59 60 pism_config:atmosphere.delta_T.reference_year = 0; 61 pism_config:atmosphere.delta_T.reference_year_doc = "Reference year to use when ``atmosphere.delta_T.period`` is active."; 62 pism_config:atmosphere.delta_T.reference_year_option = "atmosphere_delta_T_reference_year"; 63 pism_config:atmosphere.delta_T.reference_year_type = "integer"; 64 pism_config:atmosphere.delta_T.reference_year_units = "years"; 65 66 pism_config:atmosphere.elevation_change.file = ""; 67 pism_config:atmosphere.elevation_change.file_doc = "Name of the file containing the reference surface elevation field (variable ``usurf``)."; 68 pism_config:atmosphere.elevation_change.file_option = "atmosphere_lapse_rate_file"; 69 pism_config:atmosphere.elevation_change.file_type = "string"; 70 71 pism_config:atmosphere.elevation_change.period = 0; 72 pism_config:atmosphere.elevation_change.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 73 pism_config:atmosphere.elevation_change.period_option = "atmosphere_lapse_rate_period"; 74 pism_config:atmosphere.elevation_change.period_type = "integer"; 75 pism_config:atmosphere.elevation_change.period_units = "years"; 76 77 pism_config:atmosphere.elevation_change.precipitation.lapse_rate = 0.0; 78 pism_config:atmosphere.elevation_change.precipitation.lapse_rate_doc = "Elevation lapse rate for the surface mass balance"; 79 pism_config:atmosphere.elevation_change.precipitation.lapse_rate_option = "precip_lapse_rate"; 80 pism_config:atmosphere.elevation_change.precipitation.lapse_rate_type = "number"; 81 pism_config:atmosphere.elevation_change.precipitation.lapse_rate_units = "(kg m-2 / year) / km"; 82 83 pism_config:atmosphere.elevation_change.precipitation.method = "shift"; 84 pism_config:atmosphere.elevation_change.precipitation.method_choices = "scale,shift"; 85 pism_config:atmosphere.elevation_change.precipitation.method_doc = "Choose the precipitation adjustment method. ``scale``: use temperature-change-dependent scaling factor. ``shift``: use the precipitation lapse rate."; 86 pism_config:atmosphere.elevation_change.precipitation.method_option = "precip_adjustment"; 87 pism_config:atmosphere.elevation_change.precipitation.method_type = "keyword"; 88 89 pism_config:atmosphere.elevation_change.reference_year = 0; 90 pism_config:atmosphere.elevation_change.reference_year_doc = "Reference year to use when ``atmosphere.elevation_change.period`` is active."; 91 pism_config:atmosphere.elevation_change.reference_year_option = "atmosphere_lapse_rate_reference_year"; 92 pism_config:atmosphere.elevation_change.reference_year_type = "integer"; 93 pism_config:atmosphere.elevation_change.reference_year_units = "years"; 94 95 pism_config:atmosphere.elevation_change.temperature_lapse_rate = 0.0; 96 pism_config:atmosphere.elevation_change.temperature_lapse_rate_doc = "Elevation lapse rate for the surface temperature"; 97 pism_config:atmosphere.elevation_change.temperature_lapse_rate_option = "temp_lapse_rate"; 98 pism_config:atmosphere.elevation_change.temperature_lapse_rate_type = "number"; 99 pism_config:atmosphere.elevation_change.temperature_lapse_rate_units = "Kelvin / km"; 100 101 pism_config:atmosphere.fausto_air_temp.c_ma = -0.7189; 102 pism_config:atmosphere.fausto_air_temp.c_ma_doc = "latitude-dependence coefficient for formula (1) in :cite:`Faustoetal2009`"; 103 pism_config:atmosphere.fausto_air_temp.c_ma_type = "number"; 104 pism_config:atmosphere.fausto_air_temp.c_ma_units = "Kelvin / degree_north"; 105 106 pism_config:atmosphere.fausto_air_temp.c_mj = -0.1585; 107 pism_config:atmosphere.fausto_air_temp.c_mj_doc = "latitude-dependence coefficient for formula (2) in :cite:`Faustoetal2009`"; 108 pism_config:atmosphere.fausto_air_temp.c_mj_type = "number"; 109 pism_config:atmosphere.fausto_air_temp.c_mj_units = "Kelvin / degree_north"; 110 111 pism_config:atmosphere.fausto_air_temp.d_ma = 314.98; 112 pism_config:atmosphere.fausto_air_temp.d_ma_doc = "41.83+273.15; base temperature for formula (1) in :cite:`Faustoetal2009`"; 113 pism_config:atmosphere.fausto_air_temp.d_ma_type = "number"; 114 pism_config:atmosphere.fausto_air_temp.d_ma_units = "Kelvin"; 115 116 pism_config:atmosphere.fausto_air_temp.d_mj = 287.85; 117 pism_config:atmosphere.fausto_air_temp.d_mj_doc = "= 14.70+273.15; base temperature for formula (2) in :cite:`Faustoetal2009`"; 118 pism_config:atmosphere.fausto_air_temp.d_mj_type = "number"; 119 pism_config:atmosphere.fausto_air_temp.d_mj_units = "Kelvin"; 120 121 pism_config:atmosphere.fausto_air_temp.gamma_ma = -0.006309; 122 pism_config:atmosphere.fausto_air_temp.gamma_ma_doc = "= -6.309 / 1km; mean slope lapse rate for formula (1) in :cite:`Faustoetal2009`"; 123 pism_config:atmosphere.fausto_air_temp.gamma_ma_type = "number"; 124 pism_config:atmosphere.fausto_air_temp.gamma_ma_units = "Kelvin / meter"; 125 126 pism_config:atmosphere.fausto_air_temp.gamma_mj = -0.005426; 127 pism_config:atmosphere.fausto_air_temp.gamma_mj_doc = "= -5.426 / 1km; mean slope lapse rate for formula (2) in :cite:`Faustoetal2009`"; 128 pism_config:atmosphere.fausto_air_temp.gamma_mj_type = "number"; 129 pism_config:atmosphere.fausto_air_temp.gamma_mj_units = "Kelvin / meter"; 130 131 pism_config:atmosphere.fausto_air_temp.kappa_ma = 0.0672; 132 pism_config:atmosphere.fausto_air_temp.kappa_ma_doc = "longitude-dependence coefficient for formula (1) in :cite:`Faustoetal2009`"; 133 pism_config:atmosphere.fausto_air_temp.kappa_ma_type = "number"; 134 pism_config:atmosphere.fausto_air_temp.kappa_ma_units = "Kelvin / degree_west"; 135 136 pism_config:atmosphere.fausto_air_temp.kappa_mj = 0.0518; 137 pism_config:atmosphere.fausto_air_temp.kappa_mj_doc = "longitude-dependence coefficient for formula (2) in :cite:`Faustoetal2009`"; 138 pism_config:atmosphere.fausto_air_temp.kappa_mj_type = "number"; 139 pism_config:atmosphere.fausto_air_temp.kappa_mj_units = "Kelvin / degree_west"; 140 141 pism_config:atmosphere.fausto_air_temp.summer_peak_day = 196; 142 pism_config:atmosphere.fausto_air_temp.summer_peak_day_doc = "day of year for July 15; used in corrected formula (4) in :cite:`Faustoetal2009`"; 143 pism_config:atmosphere.fausto_air_temp.summer_peak_day_type = "integer"; 144 pism_config:atmosphere.fausto_air_temp.summer_peak_day_units = "ordinal day number"; 145 146 pism_config:atmosphere.frac_P.file = ""; 147 pism_config:atmosphere.frac_P.file_doc = "Name of the file containing scalar precipitation scaling."; 148 pism_config:atmosphere.frac_P.file_option = "atmosphere_frac_P_file"; 149 pism_config:atmosphere.frac_P.file_type = "string"; 150 151 pism_config:atmosphere.frac_P.period = 0; 152 pism_config:atmosphere.frac_P.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 153 pism_config:atmosphere.frac_P.period_option = "atmosphere_frac_P_period"; 154 pism_config:atmosphere.frac_P.period_type = "integer"; 155 pism_config:atmosphere.frac_P.period_units = "years"; 156 157 pism_config:atmosphere.frac_P.reference_year = 0; 158 pism_config:atmosphere.frac_P.reference_year_doc = "Reference year to use when ``atmosphere.frac_P.period`` is active."; 159 pism_config:atmosphere.frac_P.reference_year_option = "atmosphere_frac_P_reference_year"; 160 pism_config:atmosphere.frac_P.reference_year_type = "integer"; 161 pism_config:atmosphere.frac_P.reference_year_units = "years"; 162 163 pism_config:atmosphere.given.file = ""; 164 pism_config:atmosphere.given.file_doc = "Name of the file containing climate forcing fields."; 165 pism_config:atmosphere.given.file_option = "atmosphere_given_file"; 166 pism_config:atmosphere.given.file_type = "string"; 167 168 pism_config:atmosphere.given.period = 0; 169 pism_config:atmosphere.given.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 170 pism_config:atmosphere.given.period_option = "atmosphere_given_period"; 171 pism_config:atmosphere.given.period_type = "integer"; 172 pism_config:atmosphere.given.period_units = "years"; 173 174 pism_config:atmosphere.given.reference_year = 0; 175 pism_config:atmosphere.given.reference_year_doc = "Reference year to use when ``atmosphere.given.period`` is active."; 176 pism_config:atmosphere.given.reference_year_option = "atmosphere_given_reference_year"; 177 pism_config:atmosphere.given.reference_year_type = "integer"; 178 pism_config:atmosphere.given.reference_year_units = "years"; 179 180 pism_config:atmosphere.models = "given"; 181 pism_config:atmosphere.models_doc = "Comma-separated list of atmosphere melt models and modifiers."; 182 pism_config:atmosphere.models_option = "atmosphere"; 183 pism_config:atmosphere.models_type = "string"; 184 185 pism_config:atmosphere.one_station.file = ""; 186 pism_config:atmosphere.one_station.file_doc = "Specifies a file containing scalar time-series 'precipitation' and 'air_temp'."; 187 pism_config:atmosphere.one_station.file_option = "atmosphere_one_station_file"; 188 pism_config:atmosphere.one_station.file_type = "string"; 189 190 pism_config:atmosphere.orographic_precipitation.background_precip_post = 0; 191 pism_config:atmosphere.orographic_precipitation.background_precip_post_doc = "Adding background precipitation after truncation"; 192 pism_config:atmosphere.orographic_precipitation.background_precip_post_option = "background_precip_post"; 193 pism_config:atmosphere.orographic_precipitation.background_precip_post_type = "number"; 194 pism_config:atmosphere.orographic_precipitation.background_precip_post_units = "mm/hr"; 195 196 pism_config:atmosphere.orographic_precipitation.background_precip_pre = 0; 197 pism_config:atmosphere.orographic_precipitation.background_precip_pre_doc = "Adding background precipitation before truncation"; 198 pism_config:atmosphere.orographic_precipitation.background_precip_pre_option = "background_precip_pre"; 199 pism_config:atmosphere.orographic_precipitation.background_precip_pre_type = "number"; 200 pism_config:atmosphere.orographic_precipitation.background_precip_pre_units = "mm/hr"; 201 202 pism_config:atmosphere.orographic_precipitation.conversion_time = 1000.0; 203 pism_config:atmosphere.orographic_precipitation.conversion_time_doc = "Cloud conversion time"; 204 pism_config:atmosphere.orographic_precipitation.conversion_time_option = "conversion_time"; 205 pism_config:atmosphere.orographic_precipitation.conversion_time_type = "number"; 206 pism_config:atmosphere.orographic_precipitation.conversion_time_units = "s"; 207 208 pism_config:atmosphere.orographic_precipitation.coriolis_latitude = 0.0; 209 pism_config:atmosphere.orographic_precipitation.coriolis_latitude_doc = "Latitude used to compute Coriolis force"; 210 pism_config:atmosphere.orographic_precipitation.coriolis_latitude_option = "coriolis_latitude"; 211 pism_config:atmosphere.orographic_precipitation.coriolis_latitude_type = "number"; 212 pism_config:atmosphere.orographic_precipitation.coriolis_latitude_units = "degrees_N"; 213 214 pism_config:atmosphere.orographic_precipitation.fallout_time = 1000.0; 215 pism_config:atmosphere.orographic_precipitation.fallout_time_doc = "Fallout time"; 216 pism_config:atmosphere.orographic_precipitation.fallout_time_option = "fallout_time"; 217 pism_config:atmosphere.orographic_precipitation.fallout_time_type = "number"; 218 pism_config:atmosphere.orographic_precipitation.fallout_time_units = "s"; 219 220 pism_config:atmosphere.orographic_precipitation.grid_size_factor = 2; 221 pism_config:atmosphere.orographic_precipitation.grid_size_factor_doc = "The spectral grid size is (Z*(grid.Mx - 1) + 1, Z*(grid.My - 1) + 1) where Z is given by this parameter."; 222 pism_config:atmosphere.orographic_precipitation.grid_size_factor_type = "integer"; 223 pism_config:atmosphere.orographic_precipitation.grid_size_factor_units = "count"; 224 225 pism_config:atmosphere.orographic_precipitation.lapse_rate = -5.8; 226 pism_config:atmosphere.orographic_precipitation.lapse_rate_doc = "Lapse rate"; 227 pism_config:atmosphere.orographic_precipitation.lapse_rate_option = "lapse_rate"; 228 pism_config:atmosphere.orographic_precipitation.lapse_rate_type = "number"; 229 pism_config:atmosphere.orographic_precipitation.lapse_rate_units = "K / km"; 230 231 pism_config:atmosphere.orographic_precipitation.moist_adiabatic_lapse_rate = -6.5; 232 pism_config:atmosphere.orographic_precipitation.moist_adiabatic_lapse_rate_doc = "Water vapor scale height"; 233 pism_config:atmosphere.orographic_precipitation.moist_adiabatic_lapse_rate_option = "moist_adiabatic_lapse_rate"; 234 pism_config:atmosphere.orographic_precipitation.moist_adiabatic_lapse_rate_type = "number"; 235 pism_config:atmosphere.orographic_precipitation.moist_adiabatic_lapse_rate_units = "K / km"; 236 237 pism_config:atmosphere.orographic_precipitation.moist_stability_frequency = 0.05; 238 pism_config:atmosphere.orographic_precipitation.moist_stability_frequency_doc = "Moist stability frequency"; 239 pism_config:atmosphere.orographic_precipitation.moist_stability_frequency_option = "moist_stability_frequency"; 240 pism_config:atmosphere.orographic_precipitation.moist_stability_frequency_type = "number"; 241 pism_config:atmosphere.orographic_precipitation.moist_stability_frequency_units = "1/s"; 242 243 pism_config:atmosphere.orographic_precipitation.reference_density = 7.4e-3; 244 pism_config:atmosphere.orographic_precipitation.reference_density_doc = "Water vapor scale height"; 245 pism_config:atmosphere.orographic_precipitation.reference_density_option = "reference_density"; 246 pism_config:atmosphere.orographic_precipitation.reference_density_type = "number"; 247 pism_config:atmosphere.orographic_precipitation.reference_density_units = "kg m-3"; 248 249 pism_config:atmosphere.orographic_precipitation.scale_factor = 1; 250 pism_config:atmosphere.orographic_precipitation.scale_factor_doc = "Precipitation scaling factor"; 251 pism_config:atmosphere.orographic_precipitation.scale_factor_option = "scale_factor"; 252 pism_config:atmosphere.orographic_precipitation.scale_factor_type = "number"; 253 pism_config:atmosphere.orographic_precipitation.scale_factor_units = "1"; 254 255 pism_config:atmosphere.orographic_precipitation.truncate = "true"; 256 pism_config:atmosphere.orographic_precipitation.truncate_doc = "Truncate precipitation at 0"; 257 pism_config:atmosphere.orographic_precipitation.truncate_option = "truncate"; 258 pism_config:atmosphere.orographic_precipitation.truncate_type = "flag"; 259 260 pism_config:atmosphere.orographic_precipitation.water_vapor_scale_height = 2500.0; 261 pism_config:atmosphere.orographic_precipitation.water_vapor_scale_height_doc = "Water vapor scale height"; 262 pism_config:atmosphere.orographic_precipitation.water_vapor_scale_height_option = "water_vapor_scale_height"; 263 pism_config:atmosphere.orographic_precipitation.water_vapor_scale_height_type = "number"; 264 pism_config:atmosphere.orographic_precipitation.water_vapor_scale_height_units = "m"; 265 266 pism_config:atmosphere.orographic_precipitation.wind_direction = 270; 267 pism_config:atmosphere.orographic_precipitation.wind_direction_doc = "The direction the wind is coming from"; 268 pism_config:atmosphere.orographic_precipitation.wind_direction_option = "wind_direction"; 269 pism_config:atmosphere.orographic_precipitation.wind_direction_type = "number"; 270 pism_config:atmosphere.orographic_precipitation.wind_direction_units = "degrees"; 271 272 pism_config:atmosphere.orographic_precipitation.wind_speed = 10; 273 pism_config:atmosphere.orographic_precipitation.wind_speed_doc = "The direction the wind is coming from"; 274 pism_config:atmosphere.orographic_precipitation.wind_speed_option = "wind_speed"; 275 pism_config:atmosphere.orographic_precipitation.wind_speed_type = "number"; 276 pism_config:atmosphere.orographic_precipitation.wind_speed_units = "m/s"; 277 278 pism_config:atmosphere.pik.file = ""; 279 pism_config:atmosphere.pik.file_doc = "Name of the file containing the reference surface elevation field (variable ``usurf``)."; 280 pism_config:atmosphere.pik.file_option = "atmosphere_pik_file"; 281 pism_config:atmosphere.pik.file_type = "string"; 282 283 pism_config:atmosphere.pik.parameterization = "martin"; 284 pism_config:atmosphere.pik.parameterization_choices = "martin,huybrechts_dewolde,martin_huybrechts_dewolde,era_interim,era_interim_sin,era_interim_lon"; 285 pism_config:atmosphere.pik.parameterization_doc = "Selects parameterizations of mean annual and mean summer near-surface air temperatures."; 286 pism_config:atmosphere.pik.parameterization_option = "atmosphere_pik"; 287 pism_config:atmosphere.pik.parameterization_type = "keyword"; 288 289 pism_config:atmosphere.precip_exponential_factor_for_temperature = 0.07041666667; 290 pism_config:atmosphere.precip_exponential_factor_for_temperature_doc = "= 0.169/2.4; in SeaRISE-Greenland formula for precipitation correction using air temperature offsets relative to present; a 7.3\% change of precipitation rate for every one degC of temperature change :cite:`Huybrechts02`"; 291 pism_config:atmosphere.precip_exponential_factor_for_temperature_type = "number"; 292 pism_config:atmosphere.precip_exponential_factor_for_temperature_units = "Kelvin-1"; 293 294 pism_config:atmosphere.precip_scaling.file = ""; 295 pism_config:atmosphere.precip_scaling.file_doc = "Name of the file containing temperature offsets to use for a precipitation correction."; 296 pism_config:atmosphere.precip_scaling.file_option = "atmosphere_precip_scaling_file"; 297 pism_config:atmosphere.precip_scaling.file_type = "string"; 298 299 pism_config:atmosphere.precip_scaling.period = 0; 300 pism_config:atmosphere.precip_scaling.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 301 pism_config:atmosphere.precip_scaling.period_option = "atmosphere_precip_scaling_period"; 302 pism_config:atmosphere.precip_scaling.period_type = "integer"; 303 pism_config:atmosphere.precip_scaling.period_units = "years"; 304 305 pism_config:atmosphere.precip_scaling.reference_year = 0; 306 pism_config:atmosphere.precip_scaling.reference_year_doc = "Reference year to use when ``atmosphere.precip_scaling.period`` is active."; 307 pism_config:atmosphere.precip_scaling.reference_year_option = "atmosphere_precip_scaling_reference_year"; 308 pism_config:atmosphere.precip_scaling.reference_year_type = "integer"; 309 pism_config:atmosphere.precip_scaling.reference_year_units = "years"; 310 311 pism_config:atmosphere.searise_greenland.file = ""; 312 pism_config:atmosphere.searise_greenland.file_doc = "Name of the file containing a precipitation field."; 313 pism_config:atmosphere.searise_greenland.file_option = "atmosphere_searise_greenland_file"; 314 pism_config:atmosphere.searise_greenland.file_type = "string"; 315 316 pism_config:atmosphere.uniform.precipitation = 1000; 317 pism_config:atmosphere.uniform.precipitation_doc = "Precipitation used by the 'uniform' atmosphere model."; 318 pism_config:atmosphere.uniform.precipitation_type = "number"; 319 pism_config:atmosphere.uniform.precipitation_units = "kg m-2 year-1"; 320 321 pism_config:atmosphere.uniform.temperature = 273.15; 322 pism_config:atmosphere.uniform.temperature_doc = "Air temperature used by the 'uniform' atmosphere model."; 323 pism_config:atmosphere.uniform.temperature_type = "number"; 324 pism_config:atmosphere.uniform.temperature_units = "Kelvin"; 325 326 pism_config:atmosphere.yearly_cycle.file = ""; 327 pism_config:atmosphere.yearly_cycle.file_doc = "Name of the file containing mean annual and mean July temperatures (`air_temp_mean_annual` and `air_temp_mean_summer`) and the `precipitation` field."; 328 pism_config:atmosphere.yearly_cycle.file_option = "atmosphere_yearly_cycle_file"; 329 pism_config:atmosphere.yearly_cycle.file_type = "string"; 330 331 pism_config:atmosphere.yearly_cycle.scaling.file = ""; 332 pism_config:atmosphere.yearly_cycle.scaling.file_doc = "Name of the file containing amplitude scaling (`amplitude_scaling`) for the near-surface air temperature."; 333 pism_config:atmosphere.yearly_cycle.scaling.file_option = "atmosphere_yearly_cycle_scaling_file"; 334 pism_config:atmosphere.yearly_cycle.scaling.file_type = "string"; 335 336 pism_config:atmosphere.yearly_cycle.scaling.period = 0; 337 pism_config:atmosphere.yearly_cycle.scaling.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 338 pism_config:atmosphere.yearly_cycle.scaling.period_type = "integer"; 339 pism_config:atmosphere.yearly_cycle.scaling.period_units = "years"; 340 341 pism_config:atmosphere.yearly_cycle.scaling.reference_year = 0; 342 pism_config:atmosphere.yearly_cycle.scaling.reference_year_doc = "Reference year to use when ``atmosphere.yearly_cycle.scaling.period`` is active."; 343 pism_config:atmosphere.yearly_cycle.scaling.reference_year_type = "integer"; 344 pism_config:atmosphere.yearly_cycle.scaling.reference_year_units = "years"; 345 346 pism_config:basal_resistance.beta_ice_free_bedrock = 1.8e9; 347 pism_config:basal_resistance.beta_ice_free_bedrock_doc = "value is for ice stream E from :cite:`HulbeMacAyeal`; thus sliding velocity, but we hope it doesn't matter much; at 100 m/year the linear sliding law gives 57040 Pa basal shear stress"; 348 pism_config:basal_resistance.beta_ice_free_bedrock_type = "number"; 349 pism_config:basal_resistance.beta_ice_free_bedrock_units = "Pascal second meter-1"; 350 351 pism_config:basal_resistance.beta_lateral_margin = 1e19; 352 pism_config:basal_resistance.beta_lateral_margin_doc = "high value of `\\beta` used to simulate drag at lateral ice margins (fjord walls, etc); the default value is chosen to disable flow in the direction along a margin"; 353 pism_config:basal_resistance.beta_lateral_margin_type = "number"; 354 pism_config:basal_resistance.beta_lateral_margin_units = "Pascal second meter-1"; 355 356 pism_config:basal_resistance.plastic.regularization = 0.01; 357 pism_config:basal_resistance.plastic.regularization_doc = "Set the value of `\\epsilon` regularization of plastic till; this is the second `\\epsilon` in formula (4.1) in :cite:`SchoofStream`"; 358 pism_config:basal_resistance.plastic.regularization_option = "plastic_reg"; 359 pism_config:basal_resistance.plastic.regularization_type = "number"; 360 pism_config:basal_resistance.plastic.regularization_units = "meter / year"; 361 362 pism_config:basal_resistance.pseudo_plastic.enabled = "no"; 363 pism_config:basal_resistance.pseudo_plastic.enabled_doc = "Use the pseudo-plastic till model (basal sliding law)."; 364 pism_config:basal_resistance.pseudo_plastic.enabled_option = "pseudo_plastic"; 365 pism_config:basal_resistance.pseudo_plastic.enabled_type = "flag"; 366 367 pism_config:basal_resistance.pseudo_plastic.q = 0.25; 368 pism_config:basal_resistance.pseudo_plastic.q_doc = "The exponent of the pseudo-plastic basal resistance model"; 369 pism_config:basal_resistance.pseudo_plastic.q_option = "pseudo_plastic_q"; 370 pism_config:basal_resistance.pseudo_plastic.q_type = "number"; 371 pism_config:basal_resistance.pseudo_plastic.q_units = "pure number"; 372 373 pism_config:basal_resistance.pseudo_plastic.sliding_scale_factor = -1.0; 374 pism_config:basal_resistance.pseudo_plastic.sliding_scale_factor_doc = "divides pseudo-plastic tauc (yield stress) by given factor; this would increase sliding by given factor in absence of membrane stresses; not used if negative or zero; not used by default"; 375 pism_config:basal_resistance.pseudo_plastic.sliding_scale_factor_option = "sliding_scale_factor_reduces_tauc"; 376 pism_config:basal_resistance.pseudo_plastic.sliding_scale_factor_type = "number"; 377 pism_config:basal_resistance.pseudo_plastic.sliding_scale_factor_units = "1"; 378 379 pism_config:basal_resistance.pseudo_plastic.u_threshold = 100.0; 380 pism_config:basal_resistance.pseudo_plastic.u_threshold_doc = "threshold velocity of the pseudo-plastic sliding law"; 381 pism_config:basal_resistance.pseudo_plastic.u_threshold_option = "pseudo_plastic_uthreshold"; 382 pism_config:basal_resistance.pseudo_plastic.u_threshold_type = "number"; 383 pism_config:basal_resistance.pseudo_plastic.u_threshold_units = "meter / year"; 384 385 pism_config:basal_yield_stress.add_transportable_water = "no"; 386 pism_config:basal_yield_stress.add_transportable_water_doc = "If 'yes' then the water amount in the transport system is added to tillwat in determining tauc (in the Mohr-Coulomb relation). Normally only the water in the till is used."; 387 pism_config:basal_yield_stress.add_transportable_water_option = "tauc_add_transportable_water"; 388 pism_config:basal_yield_stress.add_transportable_water_type = "flag"; 389 390 pism_config:basal_yield_stress.constant.value = 2e5; 391 pism_config:basal_yield_stress.constant.value_doc = "fill value for yield stress for basal till (plastic or pseudo-plastic model); note `2 \\times 10^5` Pa = 2 bar is quite strong and little sliding should occur"; 392 pism_config:basal_yield_stress.constant.value_option = "tauc"; 393 pism_config:basal_yield_stress.constant.value_type = "number"; 394 pism_config:basal_yield_stress.constant.value_units = "Pascal"; 395 396 pism_config:basal_yield_stress.ice_free_bedrock = 1e6; 397 pism_config:basal_yield_stress.ice_free_bedrock_doc = "the 'high' yield stress value used in grounded ice-free areas."; 398 pism_config:basal_yield_stress.ice_free_bedrock_option = "high_tauc"; 399 pism_config:basal_yield_stress.ice_free_bedrock_type = "number"; 400 pism_config:basal_yield_stress.ice_free_bedrock_units = "Pascal"; 401 402 pism_config:basal_yield_stress.model = "mohr_coulomb"; 403 pism_config:basal_yield_stress.model_choices = "constant,mohr_coulomb"; 404 pism_config:basal_yield_stress.model_doc = "The basal yield stress model to use when sliding is active."; 405 pism_config:basal_yield_stress.model_option = "yield_stress"; 406 pism_config:basal_yield_stress.model_type = "keyword"; 407 408 pism_config:basal_yield_stress.mohr_coulomb.delta.file = ""; 409 pism_config:basal_yield_stress.mohr_coulomb.delta.file_doc = "Name of the file containing space- and time-dependent variable `mohr_coulomb_delta` to use instead of `basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden`."; 410 pism_config:basal_yield_stress.mohr_coulomb.delta.file_option = "mohr_coulomb_delta_file"; 411 pism_config:basal_yield_stress.mohr_coulomb.delta.file_type = "string"; 412 413 pism_config:basal_yield_stress.mohr_coulomb.delta.period = 0; 414 pism_config:basal_yield_stress.mohr_coulomb.delta.period_doc = "Length of the period. Set to zero to disable."; 415 pism_config:basal_yield_stress.mohr_coulomb.delta.period_option = "mohr_coulomb_delta_period"; 416 pism_config:basal_yield_stress.mohr_coulomb.delta.period_type = "integer"; 417 pism_config:basal_yield_stress.mohr_coulomb.delta.period_units = "years"; 418 419 pism_config:basal_yield_stress.mohr_coulomb.delta.reference_year = 0; 420 pism_config:basal_yield_stress.mohr_coulomb.delta.reference_year_doc = "Reference year to use when ``basal_yield_stress.mohr_coulomb.delta.period`` is active."; 421 pism_config:basal_yield_stress.mohr_coulomb.delta.reference_year_option = "mohr_coulomb_delta_reference_year"; 422 pism_config:basal_yield_stress.mohr_coulomb.delta.reference_year_type = "integer"; 423 pism_config:basal_yield_stress.mohr_coulomb.delta.reference_year_units = "years"; 424 425 pism_config:basal_yield_stress.mohr_coulomb.tauc_to_phi.file = ""; 426 pism_config:basal_yield_stress.mohr_coulomb.tauc_to_phi.file_doc = "File containing the basal yield stress field that should be used to recover the till friction angle distribution."; 427 pism_config:basal_yield_stress.mohr_coulomb.tauc_to_phi.file_option = "tauc_to_phi"; 428 pism_config:basal_yield_stress.mohr_coulomb.tauc_to_phi.file_type = "string"; 429 430 pism_config:basal_yield_stress.mohr_coulomb.till_cohesion = 0.0; 431 pism_config:basal_yield_stress.mohr_coulomb.till_cohesion_doc = "cohesion of till; = c_0 in most references; note Schoof uses zero but Paterson pp 168--169 gives range 0--40 kPa; but Paterson notes that '... all the pairs c_0 and phi in the table would give a yield stress for Ice Stream B that exceeds the basal shear stress there...'"; 432 pism_config:basal_yield_stress.mohr_coulomb.till_cohesion_option = "till_cohesion"; 433 pism_config:basal_yield_stress.mohr_coulomb.till_cohesion_type = "number"; 434 pism_config:basal_yield_stress.mohr_coulomb.till_cohesion_units = "Pascal"; 435 436 pism_config:basal_yield_stress.mohr_coulomb.till_compressibility_coefficient = 0.12; 437 pism_config:basal_yield_stress.mohr_coulomb.till_compressibility_coefficient_doc = "coefficient of compressiblity of till; value from :cite:`Tulaczyketal2000`"; 438 pism_config:basal_yield_stress.mohr_coulomb.till_compressibility_coefficient_option = "till_compressibility_coefficient"; 439 pism_config:basal_yield_stress.mohr_coulomb.till_compressibility_coefficient_type = "number"; 440 pism_config:basal_yield_stress.mohr_coulomb.till_compressibility_coefficient_units = "pure number"; 441 442 pism_config:basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden = 0.02; 443 pism_config:basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden_doc = "`\\delta` in notes; `N_0 = \\delta P_o` where `P_o` is overburden pressure; `N_0` is reference (low) value of effective pressure (i.e. normal stress); `N_0` scales with overburden pressure unlike :cite:`Tulaczyketal2000`; default value from Greenland and Antarctic model runs"; 444 pism_config:basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden_option = "till_effective_fraction_overburden"; 445 pism_config:basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden_type = "number"; 446 pism_config:basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden_units = "pure number"; 447 448 pism_config:basal_yield_stress.mohr_coulomb.till_log_factor_transportable_water = 0.1; 449 pism_config:basal_yield_stress.mohr_coulomb.till_log_factor_transportable_water_doc = "If basal_yield_stress.add_transportable_water = yes then the water amount in the transport system is added to tillwat in determining tauc. Normally only the water in the till is used. This factor multiplies the logarithm in that formula."; 450 pism_config:basal_yield_stress.mohr_coulomb.till_log_factor_transportable_water_option = "till_log_factor_transportable_water"; 451 pism_config:basal_yield_stress.mohr_coulomb.till_log_factor_transportable_water_type = "number"; 452 pism_config:basal_yield_stress.mohr_coulomb.till_log_factor_transportable_water_units = "meters"; 453 454 pism_config:basal_yield_stress.mohr_coulomb.till_phi_default = 30.0; 455 pism_config:basal_yield_stress.mohr_coulomb.till_phi_default_doc = "fill value for till friction angle"; 456 pism_config:basal_yield_stress.mohr_coulomb.till_phi_default_option = "plastic_phi"; 457 pism_config:basal_yield_stress.mohr_coulomb.till_phi_default_type = "number"; 458 pism_config:basal_yield_stress.mohr_coulomb.till_phi_default_units = "degrees"; 459 460 pism_config:basal_yield_stress.mohr_coulomb.till_reference_effective_pressure = 1000.0; 461 pism_config:basal_yield_stress.mohr_coulomb.till_reference_effective_pressure_doc = "reference effective pressure N_0; value from :cite:`Tulaczyketal2000`"; 462 pism_config:basal_yield_stress.mohr_coulomb.till_reference_effective_pressure_type = "number"; 463 pism_config:basal_yield_stress.mohr_coulomb.till_reference_effective_pressure_units = "Pascal"; 464 465 pism_config:basal_yield_stress.mohr_coulomb.till_reference_void_ratio = 0.69; 466 pism_config:basal_yield_stress.mohr_coulomb.till_reference_void_ratio_doc = "void ratio at reference effective pressure N_0; value from :cite:`Tulaczyketal2000`"; 467 pism_config:basal_yield_stress.mohr_coulomb.till_reference_void_ratio_option = "till_reference_void_ratio"; 468 pism_config:basal_yield_stress.mohr_coulomb.till_reference_void_ratio_type = "number"; 469 pism_config:basal_yield_stress.mohr_coulomb.till_reference_void_ratio_units = "pure number"; 470 471 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.enabled = "no"; 472 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.enabled_doc = "If the option ``-topg_to_phi`` is set then this will be set to 'yes', and then MohrCoulombYieldStress will initialize the ``tillphi`` field using a piece-wise linear function of depth described by four parameters."; 473 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.enabled_type = "flag"; 474 475 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max = 15.0; 476 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max_doc = "upper value of the till friction angle; see the implementation of MohrCoulombYieldStress"; 477 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max_type = "number"; 478 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max_units = "degrees"; 479 480 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min = 5.0; 481 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min_doc = "lower value of the till friction angle; see the implementation of MohrCoulombYieldStress"; 482 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min_type = "number"; 483 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min_units = "degrees"; 484 485 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max = 1000.0; 486 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max_doc = "the elevation at which the upper value of the till friction angle is used; see the implementation of MohrCoulombYieldStress"; 487 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max_type = "number"; 488 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max_units = "meters"; 489 490 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min = -1000.0; 491 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min_doc = "the elevation at which the lower value of the till friction angle is used; see the implementation of MohrCoulombYieldStress"; 492 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min_type = "number"; 493 pism_config:basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min_units = "meters"; 494 495 pism_config:basal_yield_stress.slippery_grounding_lines = "no"; 496 pism_config:basal_yield_stress.slippery_grounding_lines_doc = "If yes, at icy grounded locations with bed elevations below sea level, within one cell of floating ice or ice-free ocean, make tauc as low as possible from the Mohr-Coulomb relation. Specifically, at such locations replace the normally-computed tauc from the Mohr-Coulomb relation, which uses the effective pressure from the modeled amount of water in the till, by the minimum value of tauc from Mohr-Coulomb, i.e. by using the effective pressure corresponding to the maximum amount of till-stored water. Does not alter the modeled or reported amount of till water, nor does this mechanism affect water conservation."; 497 pism_config:basal_yield_stress.slippery_grounding_lines_option = "tauc_slippery_grounding_lines"; 498 pism_config:basal_yield_stress.slippery_grounding_lines_type = "flag"; 499 500 pism_config:bed_deformation.bed_topography_delta_file = ""; 501 pism_config:bed_deformation.bed_topography_delta_file_doc = "The name of the file to read the topg_delta from. This field is added to the bed topography during initialization."; 502 pism_config:bed_deformation.bed_topography_delta_file_option = "topg_delta_file"; 503 pism_config:bed_deformation.bed_topography_delta_file_type = "string"; 504 505 pism_config:bed_deformation.bed_uplift_file = ""; 506 pism_config:bed_deformation.bed_uplift_file_doc = "The name of the file to read the uplift (dbdt) from. Leave empty to read it from an input file or a regridding file."; 507 pism_config:bed_deformation.bed_uplift_file_option = "uplift_file"; 508 pism_config:bed_deformation.bed_uplift_file_type = "string"; 509 510 pism_config:bed_deformation.lc.elastic_model = "yes"; 511 pism_config:bed_deformation.lc.elastic_model_doc = "Use the elastic part of the Lingle-Clark bed deformation model."; 512 pism_config:bed_deformation.lc.elastic_model_option = "bed_def_lc_elastic_model"; 513 pism_config:bed_deformation.lc.elastic_model_type = "flag"; 514 515 pism_config:bed_deformation.lc.grid_size_factor = 4; 516 pism_config:bed_deformation.lc.grid_size_factor_doc = "The spectral grid size is (Z*(grid.Mx - 1) + 1, Z*(grid.My - 1) + 1) where Z is given by this parameter. See :cite:`LingleClark`, :cite:`BLKfastearth`"; 517 pism_config:bed_deformation.lc.grid_size_factor_type = "integer"; 518 pism_config:bed_deformation.lc.grid_size_factor_units = "count"; 519 520 pism_config:bed_deformation.lc.update_interval = 10.0; 521 pism_config:bed_deformation.lc.update_interval_doc = "Interval between updates of the Lingle-Clark model"; 522 pism_config:bed_deformation.lc.update_interval_type = "number"; 523 pism_config:bed_deformation.lc.update_interval_units = "years"; 524 525 pism_config:bed_deformation.lithosphere_flexural_rigidity = 5.0e24; 526 pism_config:bed_deformation.lithosphere_flexural_rigidity_doc = "lithosphere flexural rigidity used by the bed deformation model. See :cite:`LingleClark`, :cite:`BLKfastearth`"; 527 pism_config:bed_deformation.lithosphere_flexural_rigidity_type = "number"; 528 pism_config:bed_deformation.lithosphere_flexural_rigidity_units = "Newton meter"; 529 530 pism_config:bed_deformation.mantle_density = 3300.0; 531 pism_config:bed_deformation.mantle_density_doc = "half-space (mantle) density used by the bed deformation model. See :cite:`LingleClark`, :cite:`BLKfastearth`"; 532 pism_config:bed_deformation.mantle_density_type = "number"; 533 pism_config:bed_deformation.mantle_density_units = "kg meter-3"; 534 535 pism_config:bed_deformation.mantle_viscosity = 1.0e21; 536 pism_config:bed_deformation.mantle_viscosity_doc = "half-space (mantle) viscosity used by the bed deformation model. See :cite:`LingleClark`, :cite:`BLKfastearth`"; 537 pism_config:bed_deformation.mantle_viscosity_type = "number"; 538 pism_config:bed_deformation.mantle_viscosity_units = "Pascal second"; 539 540 pism_config:bed_deformation.model = "none"; 541 pism_config:bed_deformation.model_choices = "none,iso,lc"; 542 pism_config:bed_deformation.model_doc = "Selects a bed deformation model to use. 'iso' is point-wise isostasy, 'lc' is the Lingle-Clark model (see :cite:`LingleClark`, requires FFTW3)."; 543 pism_config:bed_deformation.model_option = "bed_def"; 544 pism_config:bed_deformation.model_type = "keyword"; 545 546 pism_config:bootstrapping.defaults.bed = 1.0; 547 pism_config:bootstrapping.defaults.bed_doc = "bed elevation value to use if topg (bedrock_altitude) variable is absent in bootstrapping file"; 548 pism_config:bootstrapping.defaults.bed_type = "number"; 549 pism_config:bootstrapping.defaults.bed_units = "meters"; 550 551 pism_config:bootstrapping.defaults.bmelt = 0.0; 552 pism_config:bootstrapping.defaults.bmelt_doc = "basal melt rate value to use if variable bmelt is absent in bootstrapping file"; 553 pism_config:bootstrapping.defaults.bmelt_type = "number"; 554 pism_config:bootstrapping.defaults.bmelt_units = "meter / second"; 555 556 pism_config:bootstrapping.defaults.bwat = 0.0; 557 pism_config:bootstrapping.defaults.bwat_doc = "till water thickness value to use if variable tillwat is absent in bootstrapping file"; 558 pism_config:bootstrapping.defaults.bwat_type = "number"; 559 pism_config:bootstrapping.defaults.bwat_units = "meters"; 560 561 pism_config:bootstrapping.defaults.bwp = 0.0; 562 pism_config:bootstrapping.defaults.bwp_doc = "basal water pressure value to use if variable bwp is absent in bootstrapping file; most hydrology models do not use this value because bwp is diagnostic"; 563 pism_config:bootstrapping.defaults.bwp_type = "number"; 564 pism_config:bootstrapping.defaults.bwp_units = "Pascal"; 565 566 pism_config:bootstrapping.defaults.enwat = 0.0; 567 pism_config:bootstrapping.defaults.enwat_doc = "effective englacial water thickness value to use if variable enwat is absent in bootstrapping file"; 568 pism_config:bootstrapping.defaults.enwat_type = "number"; 569 pism_config:bootstrapping.defaults.enwat_units = "meters"; 570 571 pism_config:bootstrapping.defaults.geothermal_flux = 0.042; 572 pism_config:bootstrapping.defaults.geothermal_flux_doc = "geothermal flux value to use if bheatflx variable is absent in bootstrapping file"; 573 pism_config:bootstrapping.defaults.geothermal_flux_type = "number"; 574 pism_config:bootstrapping.defaults.geothermal_flux_units = "Watt meter-2"; 575 576 pism_config:bootstrapping.defaults.ice_thickness = 0.0; 577 pism_config:bootstrapping.defaults.ice_thickness_doc = "thickness value to use if thk (land_ice_thickness) variable is absent in bootstrapping file"; 578 pism_config:bootstrapping.defaults.ice_thickness_type = "number"; 579 pism_config:bootstrapping.defaults.ice_thickness_units = "meters"; 580 581 pism_config:bootstrapping.defaults.tillwat = 0.0; 582 pism_config:bootstrapping.defaults.tillwat_doc = "till water thickness value to use if variable tillwat is absent in bootstrapping file"; 583 pism_config:bootstrapping.defaults.tillwat_type = "number"; 584 pism_config:bootstrapping.defaults.tillwat_units = "meters"; 585 586 pism_config:bootstrapping.defaults.uplift = 0.0; 587 pism_config:bootstrapping.defaults.uplift_doc = "uplift value to use if dbdt variable is absent in bootstrapping file"; 588 pism_config:bootstrapping.defaults.uplift_type = "number"; 589 pism_config:bootstrapping.defaults.uplift_units = "meter / second"; 590 591 pism_config:bootstrapping.temperature_heuristic = "smb"; 592 pism_config:bootstrapping.temperature_heuristic_choices = "smb,quartic_guess"; 593 pism_config:bootstrapping.temperature_heuristic_doc = "The heuristic to use to initialize ice temperature during bootstrapping: 'smb' uses the surface mass balance, surface temperature, and the geothermal flux, 'quartic_guess' uses the surface temperature and the geothermal flux."; 594 pism_config:bootstrapping.temperature_heuristic_option = "boot_temperature_heuristic"; 595 pism_config:bootstrapping.temperature_heuristic_type = "keyword"; 596 597 pism_config:calving.eigen_calving.K = 0.0; 598 pism_config:calving.eigen_calving.K_doc = "Set proportionality constant to determine calving rate from strain rates. Note references :cite:`Levermannetal2012`, :cite:`Martinetal2011` use K in range `10^{9}` to `3 \\times 10^{11}` m a, that is, `3 \\times 10^{16}` to `10^{19}` m s."; 599 pism_config:calving.eigen_calving.K_option = "eigen_calving_K"; 600 pism_config:calving.eigen_calving.K_type = "number"; 601 pism_config:calving.eigen_calving.K_units = "meter second"; 602 603 pism_config:calving.float_kill.calve_near_grounding_line = "yes"; 604 pism_config:calving.float_kill.calve_near_grounding_line_doc = "Calve floating ice near the grounding line."; 605 pism_config:calving.float_kill.calve_near_grounding_line_option = "float_kill_calve_near_grounding_line"; 606 pism_config:calving.float_kill.calve_near_grounding_line_type = "flag"; 607 608 pism_config:calving.float_kill.margin_only = "no"; 609 pism_config:calving.float_kill.margin_only_doc = "Apply float_kill at ice margin cells only."; 610 pism_config:calving.float_kill.margin_only_option = "float_kill_margin_only"; 611 pism_config:calving.float_kill.margin_only_type = "flag"; 612 613 pism_config:calving.hayhurst_calving.B_tilde = 65.0; 614 pism_config:calving.hayhurst_calving.B_tilde_doc = "Effective damage rate :cite:`Mercenier2018`"; 615 pism_config:calving.hayhurst_calving.B_tilde_type = "number"; 616 pism_config:calving.hayhurst_calving.B_tilde_units = "(MPa)^r / year"; 617 618 pism_config:calving.hayhurst_calving.exponent_r = 0.43; 619 pism_config:calving.hayhurst_calving.exponent_r_doc = "Damage law exponent :cite:`Mercenier2018`"; 620 pism_config:calving.hayhurst_calving.exponent_r_type = "number"; 621 pism_config:calving.hayhurst_calving.exponent_r_units = "1"; 622 623 pism_config:calving.hayhurst_calving.sigma_threshold = 0.17; 624 pism_config:calving.hayhurst_calving.sigma_threshold_doc = "Damage threshold stress :cite:`Mercenier2018`"; 625 pism_config:calving.hayhurst_calving.sigma_threshold_type = "number"; 626 pism_config:calving.hayhurst_calving.sigma_threshold_units = "MPa"; 627 628 pism_config:calving.methods = ""; 629 pism_config:calving.methods_doc = "comma-separated list of calving methods; one or more of 'eigen_calving', 'float_kill', 'thickness_calving'"; 630 pism_config:calving.methods_option = "calving"; 631 pism_config:calving.methods_type = "string"; 632 633 pism_config:calving.thickness_calving.threshold = 50.0; 634 pism_config:calving.thickness_calving.threshold_doc = "When terminal ice thickness of floating ice shelf is less than this threshold, it will be calved off."; 635 pism_config:calving.thickness_calving.threshold_option = "thickness_calving_threshold"; 636 pism_config:calving.thickness_calving.threshold_type = "number"; 637 pism_config:calving.thickness_calving.threshold_units = "meters"; 638 639 pism_config:calving.thickness_calving.threshold_file = ""; 640 pism_config:calving.thickness_calving.threshold_file_doc = "Name of the file containing the spatially-variable thickness calving threshold."; 641 pism_config:calving.thickness_calving.threshold_file_option = "thickness_calving_threshold_file"; 642 pism_config:calving.thickness_calving.threshold_file_type = "string"; 643 644 pism_config:calving.vonmises_calving.Glen_exponent = 3.0; 645 pism_config:calving.vonmises_calving.Glen_exponent_doc = "Glen exponent in ice flow law for von Mises calving"; 646 pism_config:calving.vonmises_calving.Glen_exponent_option = "vonmises_calving_n"; 647 pism_config:calving.vonmises_calving.Glen_exponent_type = "number"; 648 pism_config:calving.vonmises_calving.Glen_exponent_units = "pure number"; 649 650 pism_config:calving.vonmises_calving.enhancement_factor = 1.0; 651 pism_config:calving.vonmises_calving.enhancement_factor_doc = "Flow enhancement factor used by the flow law in the von Mises stress computation"; 652 pism_config:calving.vonmises_calving.enhancement_factor_type = "number"; 653 pism_config:calving.vonmises_calving.enhancement_factor_units = "1"; 654 655 pism_config:calving.vonmises_calving.enhancement_factor_interglacial = 1.0; 656 pism_config:calving.vonmises_calving.enhancement_factor_interglacial_doc = "Flow enhancement factor used by the flow law in the von Mises stress computation (for ice accumulated during interglacial periods)."; 657 pism_config:calving.vonmises_calving.enhancement_factor_interglacial_type = "number"; 658 pism_config:calving.vonmises_calving.enhancement_factor_interglacial_units = "1"; 659 660 pism_config:calving.vonmises_calving.flow_law = "gpbld"; 661 pism_config:calving.vonmises_calving.flow_law_choices = "arr,arrwarm,gpbld,hooke,isothermal_glen,pb"; 662 pism_config:calving.vonmises_calving.flow_law_doc = "The custom flow law for the von Mises stress computation"; 663 pism_config:calving.vonmises_calving.flow_law_type = "keyword"; 664 665 pism_config:calving.vonmises_calving.sigma_max = 1.0e6; 666 pism_config:calving.vonmises_calving.sigma_max_doc = "Set maximum tensile stress. Note references :cite:`Morlighem2016` use 1.0e6 Pa."; 667 pism_config:calving.vonmises_calving.sigma_max_option = "vonmises_calving_calving_sigma_max"; 668 pism_config:calving.vonmises_calving.sigma_max_type = "number"; 669 pism_config:calving.vonmises_calving.sigma_max_units = "Pa"; 670 671 pism_config:calving.vonmises_calving.threshold_file = ""; 672 pism_config:calving.vonmises_calving.threshold_file_doc = "Name of the file containing the spatially-variable vonmises_calving calving threshold."; 673 pism_config:calving.vonmises_calving.threshold_file_option = "vonmises_calving_threshold_file"; 674 pism_config:calving.vonmises_calving.threshold_file_type = "string"; 675 676 pism_config:calving.vonmises_calving.use_custom_flow_law = "no"; 677 pism_config:calving.vonmises_calving.use_custom_flow_law_doc = "Use custom flow law in the von Mises stress computation"; 678 pism_config:calving.vonmises_calving.use_custom_flow_law_option = "vonmises_calving_use_custom_flow_law"; 679 pism_config:calving.vonmises_calving.use_custom_flow_law_type = "flag"; 680 681 pism_config:constants.fresh_water.density = 1000.0; 682 pism_config:constants.fresh_water.density_doc = "density of fresh water"; 683 pism_config:constants.fresh_water.density_type = "number"; 684 pism_config:constants.fresh_water.density_units = "kg meter-3"; 685 686 pism_config:constants.fresh_water.latent_heat_of_fusion = 3.34e5; 687 pism_config:constants.fresh_water.latent_heat_of_fusion_doc = "latent heat of fusion for water :cite:`AschwandenBlatter`"; 688 pism_config:constants.fresh_water.latent_heat_of_fusion_type = "number"; 689 pism_config:constants.fresh_water.latent_heat_of_fusion_units = "Joule / kg"; 690 691 pism_config:constants.fresh_water.melting_point_temperature = 273.15; 692 pism_config:constants.fresh_water.melting_point_temperature_doc = "melting point of pure water"; 693 pism_config:constants.fresh_water.melting_point_temperature_type = "number"; 694 pism_config:constants.fresh_water.melting_point_temperature_units = "Kelvin"; 695 696 pism_config:constants.fresh_water.specific_heat_capacity = 4170.0; 697 pism_config:constants.fresh_water.specific_heat_capacity_doc = "at melting point T_0 :cite:`AschwandenBlatter`"; 698 pism_config:constants.fresh_water.specific_heat_capacity_type = "number"; 699 pism_config:constants.fresh_water.specific_heat_capacity_units = "Joule / (kg Kelvin)"; 700 701 pism_config:constants.global_ocean_area = 3.625e14; 702 pism_config:constants.global_ocean_area_doc = "area of the global ocean :cite:`Cogley2011`"; 703 pism_config:constants.global_ocean_area_type = "number"; 704 pism_config:constants.global_ocean_area_units = "meter2"; 705 706 pism_config:constants.ice.beta_Clausius_Clapeyron = 7.9e-8; 707 pism_config:constants.ice.beta_Clausius_Clapeyron_doc = "Clausius-Clapeyron constant relating melting temperature and pressure: `\\beta = dT / dP` :cite:`Luethi2002`"; 708 pism_config:constants.ice.beta_Clausius_Clapeyron_type = "number"; 709 pism_config:constants.ice.beta_Clausius_Clapeyron_units = "Kelvin / Pascal"; 710 711 pism_config:constants.ice.density = 910.0; 712 pism_config:constants.ice.density_doc = "`\\rho_i`; density of ice in ice sheet"; 713 pism_config:constants.ice.density_type = "number"; 714 pism_config:constants.ice.density_units = "kg meter-3"; 715 716 pism_config:constants.ice.grain_size = 1.0; 717 pism_config:constants.ice.grain_size_doc = "Default constant ice grain size to use with the Goldsby-Kohlstedt :cite:`GoldsbyKohlstedt` flow law"; 718 pism_config:constants.ice.grain_size_option = "ice_grain_size"; 719 pism_config:constants.ice.grain_size_type = "number"; 720 pism_config:constants.ice.grain_size_units = "mm"; 721 722 pism_config:constants.ice.specific_heat_capacity = 2009.0; 723 pism_config:constants.ice.specific_heat_capacity_doc = "specific heat capacity of pure ice at melting point T_0"; 724 pism_config:constants.ice.specific_heat_capacity_type = "number"; 725 pism_config:constants.ice.specific_heat_capacity_units = "Joule / (kg Kelvin)"; 726 727 pism_config:constants.ice.thermal_conductivity = 2.10; 728 pism_config:constants.ice.thermal_conductivity_doc = "= W m-1 K-1; thermal conductivity of pure ice"; 729 pism_config:constants.ice.thermal_conductivity_type = "number"; 730 pism_config:constants.ice.thermal_conductivity_units = "Joule / (meter Kelvin second)"; 731 732 pism_config:constants.ideal_gas_constant = 8.31441; 733 pism_config:constants.ideal_gas_constant_doc = "ideal gas constant"; 734 pism_config:constants.ideal_gas_constant_type = "number"; 735 pism_config:constants.ideal_gas_constant_units = "Joule / (mol Kelvin)"; 736 737 pism_config:constants.sea_water.density = 1028.0; 738 pism_config:constants.sea_water.density_doc = "density of sea water"; 739 pism_config:constants.sea_water.density_type = "number"; 740 pism_config:constants.sea_water.density_units = "kg meter-3"; 741 742 pism_config:constants.sea_water.specific_heat_capacity = 3985.0; 743 pism_config:constants.sea_water.specific_heat_capacity_doc = "at 35 psu, value taken from `Kaye and Laby`_, section 2.7.9"; 744 pism_config:constants.sea_water.specific_heat_capacity_type = "number"; 745 pism_config:constants.sea_water.specific_heat_capacity_units = "Joule / (kg Kelvin)"; 746 747 pism_config:constants.standard_gravity = 9.81; 748 pism_config:constants.standard_gravity_doc = "acceleration due to gravity on Earth geoid"; 749 pism_config:constants.standard_gravity_type = "number"; 750 pism_config:constants.standard_gravity_units = "meter second-2"; 751 752 pism_config:energy.allow_temperature_above_melting = "no"; 753 pism_config:energy.allow_temperature_above_melting_doc = "If set to 'yes', allow temperatures above the pressure-malting point in the cold mode temperature code. Used by some verifiaction tests."; 754 pism_config:energy.allow_temperature_above_melting_type = "flag"; 755 756 pism_config:energy.basal_melt.use_grounded_cell_fraction = "true"; 757 pism_config:energy.basal_melt.use_grounded_cell_fraction_doc = "If geometry.grounded_cell_fraction is set, use the fractional floatation mask to interpolate the basal melt rate at the grounding line between grounded and floating values."; 758 pism_config:energy.basal_melt.use_grounded_cell_fraction_option = "subgl_basal_melt"; 759 pism_config:energy.basal_melt.use_grounded_cell_fraction_type = "flag"; 760 761 pism_config:energy.bedrock_thermal.conductivity = 3.0; 762 pism_config:energy.bedrock_thermal.conductivity_doc = "= W m-1 K-1; for bedrock used in thermal model :cite:`RitzEISMINT`"; 763 pism_config:energy.bedrock_thermal.conductivity_type = "number"; 764 pism_config:energy.bedrock_thermal.conductivity_units = "Joule / (meter Kelvin second)"; 765 766 pism_config:energy.bedrock_thermal.density = 3300.0; 767 pism_config:energy.bedrock_thermal.density_doc = "for bedrock used in thermal model"; 768 pism_config:energy.bedrock_thermal.density_type = "number"; 769 pism_config:energy.bedrock_thermal.density_units = "kg meter-3"; 770 771 pism_config:energy.bedrock_thermal.file = ""; 772 pism_config:energy.bedrock_thermal.file_doc = "Name of the file containing the geothermal flux field :var:`bheatflx`. Leave empty to read it from the :config:`input.file`."; 773 pism_config:energy.bedrock_thermal.file_type = "string"; 774 775 pism_config:energy.bedrock_thermal.specific_heat_capacity = 1000.0; 776 pism_config:energy.bedrock_thermal.specific_heat_capacity_doc = "for bedrock used in thermal model :cite:`RitzEISMINT`"; 777 pism_config:energy.bedrock_thermal.specific_heat_capacity_type = "number"; 778 pism_config:energy.bedrock_thermal.specific_heat_capacity_units = "Joule / (kg Kelvin)"; 779 780 pism_config:energy.ch_warming.average_channel_spacing = 20.0; 781 pism_config:energy.ch_warming.average_channel_spacing_doc = "Average spacing between elements of the cryo-hydrologic system (controls the rate of heat transfer from the CH system into the ice)."; 782 pism_config:energy.ch_warming.average_channel_spacing_type = "number"; 783 pism_config:energy.ch_warming.average_channel_spacing_units = "meters"; 784 785 pism_config:energy.ch_warming.enabled = "no"; 786 pism_config:energy.ch_warming.enabled_doc = "Enable the cryo-hydrologic warming model"; 787 pism_config:energy.ch_warming.enabled_type = "flag"; 788 789 pism_config:energy.ch_warming.residual_water_fraction = 0.005; 790 pism_config:energy.ch_warming.residual_water_fraction_doc = "residual water fraction in the cryo-hydrologic system at the end of a melt season"; 791 pism_config:energy.ch_warming.residual_water_fraction_type = "number"; 792 pism_config:energy.ch_warming.residual_water_fraction_units = "pure number"; 793 794 pism_config:energy.ch_warming.temperate_ice_thermal_conductivity_ratio = 1.0; 795 pism_config:energy.ch_warming.temperate_ice_thermal_conductivity_ratio_doc = "ratio of thermal conductivities of temperate and cold ice in the cryo-hydrologic system"; 796 pism_config:energy.ch_warming.temperate_ice_thermal_conductivity_ratio_type = "number"; 797 pism_config:energy.ch_warming.temperate_ice_thermal_conductivity_ratio_units = "pure number"; 798 799 pism_config:energy.drainage_maximum_rate = 1.58443823077064e-09; 800 pism_config:energy.drainage_maximum_rate_doc = "0.05 year-1; maximum rate at which liquid water fraction in temperate ice could possibly drain; see :cite:`AschwandenBuelerKhroulevBlatter`"; 801 pism_config:energy.drainage_maximum_rate_type = "number"; 802 pism_config:energy.drainage_maximum_rate_units = "second-1"; 803 804 pism_config:energy.drainage_target_water_fraction = 0.01; 805 pism_config:energy.drainage_target_water_fraction_doc = "liquid water fraction (omega) above which drainage occurs, but below which there is no drainage; see :cite:`AschwandenBuelerKhroulevBlatter`"; 806 pism_config:energy.drainage_target_water_fraction_type = "number"; 807 pism_config:energy.drainage_target_water_fraction_units = "1"; 808 809 pism_config:energy.enabled = "yes"; 810 pism_config:energy.enabled_doc = "Solve energy conservation equations."; 811 pism_config:energy.enabled_type = "flag"; 812 813 pism_config:energy.enthalpy.cold_bulge_max = 60270.0; 814 pism_config:energy.enthalpy.cold_bulge_max_doc = "= (2009 J kg-1 K-1) * (30 K); maximum amount by which advection can reduce the enthalpy of a column of ice below its surface enthalpy value"; 815 pism_config:energy.enthalpy.cold_bulge_max_type = "number"; 816 pism_config:energy.enthalpy.cold_bulge_max_units = "Joule / kg"; 817 818 pism_config:energy.enthalpy.temperate_ice_thermal_conductivity_ratio = 0.1; 819 pism_config:energy.enthalpy.temperate_ice_thermal_conductivity_ratio_doc = "K in cold ice is multiplied by this fraction to give K0 in :cite:`AschwandenBuelerKhroulevBlatter`"; 820 pism_config:energy.enthalpy.temperate_ice_thermal_conductivity_ratio_type = "number"; 821 pism_config:energy.enthalpy.temperate_ice_thermal_conductivity_ratio_units = "pure number"; 822 823 pism_config:energy.margin_exclude_horizontal_advection = "yes"; 824 pism_config:energy.margin_exclude_horizontal_advection_doc = "Exclude horizontal advection of energy at grid points near ice margins. See :config:`energy.margin_ice_thickness_limit`."; 825 pism_config:energy.margin_exclude_horizontal_advection_type = "flag"; 826 827 pism_config:energy.margin_exclude_strain_heating = "yes"; 828 pism_config:energy.margin_exclude_strain_heating_doc = "Exclude strain heating at grid points near ice margins. See :config:`energy.margin_ice_thickness_limit`."; 829 pism_config:energy.margin_exclude_strain_heating_type = "flag"; 830 831 pism_config:energy.margin_exclude_vertical_advection = "yes"; 832 pism_config:energy.margin_exclude_vertical_advection_doc = "Exclude vertical advection of energy at grid points near ice margins. See :config:`energy.margin_ice_thickness_limit`."; 833 pism_config:energy.margin_exclude_vertical_advection_type = "flag"; 834 835 pism_config:energy.margin_ice_thickness_limit = 100.0; 836 pism_config:energy.margin_ice_thickness_limit_doc = "use special margin treatment at grid points with a neighbor with the thickness below this limit."; 837 pism_config:energy.margin_ice_thickness_limit_type = "number"; 838 pism_config:energy.margin_ice_thickness_limit_units = "meters"; 839 840 pism_config:energy.max_low_temperature_count = 10; 841 pism_config:energy.max_low_temperature_count_doc = "Maximum number of grid points with ice temperature below energy.minimum_allowed_temperature."; 842 pism_config:energy.max_low_temperature_count_option = "max_low_temps"; 843 pism_config:energy.max_low_temperature_count_type = "integer"; 844 pism_config:energy.max_low_temperature_count_units = "count"; 845 846 pism_config:energy.minimum_allowed_temperature = 200.0; 847 pism_config:energy.minimum_allowed_temperature_doc = "Minimum allowed ice temperature"; 848 pism_config:energy.minimum_allowed_temperature_option = "low_temp"; 849 pism_config:energy.minimum_allowed_temperature_type = "number"; 850 pism_config:energy.minimum_allowed_temperature_units = "Kelvin"; 851 852 pism_config:energy.temperature_based = "no"; 853 pism_config:energy.temperature_based_doc = "Use cold ice (i.e. not polythermal) methods."; 854 pism_config:energy.temperature_based_type = "flag"; 855 856 pism_config:energy.temperature_dependent_thermal_conductivity = "no"; 857 pism_config:energy.temperature_dependent_thermal_conductivity_doc = "If yes, use varkenthSystemCtx class in the energy step. It is base on formula (4.37) in :cite:`GreveBlatter2009`. Otherwise use enthSystemCtx, which has temperature-independent thermal conductivity set by constant ice.thermal_conductivity."; 858 pism_config:energy.temperature_dependent_thermal_conductivity_option = "vark"; 859 pism_config:energy.temperature_dependent_thermal_conductivity_type = "flag"; 860 861 pism_config:enthalpy_converter.T_reference = 223.15; 862 pism_config:enthalpy_converter.T_reference_doc = "= T_0 in enthalpy formulas in :cite:`AschwandenBuelerKhroulevBlatter`"; 863 pism_config:enthalpy_converter.T_reference_type = "number"; 864 pism_config:enthalpy_converter.T_reference_units = "Kelvin"; 865 866 pism_config:enthalpy_converter.relaxed_is_temperate_tolerance = 0.001; 867 pism_config:enthalpy_converter.relaxed_is_temperate_tolerance_doc = "Tolerance within which ice is treated as temperate (cold-ice mode and diagnostics)."; 868 pism_config:enthalpy_converter.relaxed_is_temperate_tolerance_type = "number"; 869 pism_config:enthalpy_converter.relaxed_is_temperate_tolerance_units = "Kelvin"; 870 871 pism_config:flow_law.Hooke.A = 4.42165e-9; 872 pism_config:flow_law.Hooke.A_doc = "`A_{\\text{Hooke}} = (1/B_0)^n` where n=3 and B_0 = 1.928 `a^{1/3}` Pa. See :cite:`Hooke`"; 873 pism_config:flow_law.Hooke.A_type = "number"; 874 pism_config:flow_law.Hooke.A_units = "Pascal-3 second-1"; 875 876 pism_config:flow_law.Hooke.C = 0.16612; 877 pism_config:flow_law.Hooke.C_doc = "See :cite:`Hooke`"; 878 pism_config:flow_law.Hooke.C_type = "number"; 879 pism_config:flow_law.Hooke.C_units = "Kelvin^{flow_law.Hooke.k}"; 880 881 pism_config:flow_law.Hooke.Q = 7.88e4; 882 pism_config:flow_law.Hooke.Q_doc = "Activation energy, see :cite:`Hooke`"; 883 pism_config:flow_law.Hooke.Q_type = "number"; 884 pism_config:flow_law.Hooke.Q_units = "Joule / mol"; 885 886 pism_config:flow_law.Hooke.Tr = 273.39; 887 pism_config:flow_law.Hooke.Tr_doc = "See :cite:`Hooke`"; 888 pism_config:flow_law.Hooke.Tr_type = "number"; 889 pism_config:flow_law.Hooke.Tr_units = "Kelvin"; 890 891 pism_config:flow_law.Hooke.k = 1.17; 892 pism_config:flow_law.Hooke.k_doc = "See :cite:`Hooke`"; 893 pism_config:flow_law.Hooke.k_type = "number"; 894 pism_config:flow_law.Hooke.k_units = "pure number"; 895 896 pism_config:flow_law.Paterson_Budd.A_cold = 3.61e-13; 897 pism_config:flow_law.Paterson_Budd.A_cold_doc = "Paterson-Budd A_cold, see :cite:`PatersonBudd`"; 898 pism_config:flow_law.Paterson_Budd.A_cold_type = "number"; 899 pism_config:flow_law.Paterson_Budd.A_cold_units = "Pascal-3 / second"; 900 901 pism_config:flow_law.Paterson_Budd.A_warm = 1.73e3; 902 pism_config:flow_law.Paterson_Budd.A_warm_doc = "Paterson-Budd A_warm, see :cite:`PatersonBudd`"; 903 pism_config:flow_law.Paterson_Budd.A_warm_type = "number"; 904 pism_config:flow_law.Paterson_Budd.A_warm_units = "Pascal-3 / second"; 905 906 pism_config:flow_law.Paterson_Budd.Q_cold = 6.0e4; 907 pism_config:flow_law.Paterson_Budd.Q_cold_doc = "Paterson-Budd Q_cold, see :cite:`PatersonBudd`"; 908 pism_config:flow_law.Paterson_Budd.Q_cold_type = "number"; 909 pism_config:flow_law.Paterson_Budd.Q_cold_units = "Joule / mol"; 910 911 pism_config:flow_law.Paterson_Budd.Q_warm = 13.9e4; 912 pism_config:flow_law.Paterson_Budd.Q_warm_doc = "Paterson-Budd Q_warm, see :cite:`PatersonBudd`"; 913 pism_config:flow_law.Paterson_Budd.Q_warm_type = "number"; 914 pism_config:flow_law.Paterson_Budd.Q_warm_units = "Joule / mol"; 915 916 pism_config:flow_law.Paterson_Budd.T_critical = 263.15; 917 pism_config:flow_law.Paterson_Budd.T_critical_doc = "Paterson-Budd critical temperature, see :cite:`PatersonBudd`"; 918 pism_config:flow_law.Paterson_Budd.T_critical_type = "number"; 919 pism_config:flow_law.Paterson_Budd.T_critical_units = "Kelvin"; 920 921 pism_config:flow_law.Schoof_regularizing_length = 1000.0; 922 pism_config:flow_law.Schoof_regularizing_length_doc = "Regularizing length (Schoof definition)"; 923 pism_config:flow_law.Schoof_regularizing_length_type = "number"; 924 pism_config:flow_law.Schoof_regularizing_length_units = "km"; 925 926 pism_config:flow_law.Schoof_regularizing_velocity = 1.0; 927 pism_config:flow_law.Schoof_regularizing_velocity_doc = "Regularizing velocity (Schoof definition)"; 928 pism_config:flow_law.Schoof_regularizing_velocity_type = "number"; 929 pism_config:flow_law.Schoof_regularizing_velocity_units = "meter / year"; 930 931 pism_config:flow_law.gpbld.water_frac_coeff = 181.25; 932 pism_config:flow_law.gpbld.water_frac_coeff_doc = "coefficient in Glen-Paterson-Budd flow law for extra dependence of softness on liquid water fraction (omega) :cite:`GreveBlatter2009`, :cite:`LliboutryDuval1985`"; 933 pism_config:flow_law.gpbld.water_frac_coeff_type = "number"; 934 pism_config:flow_law.gpbld.water_frac_coeff_units = "pure number"; 935 936 pism_config:flow_law.gpbld.water_frac_observed_limit = 0.01; 937 pism_config:flow_law.gpbld.water_frac_observed_limit_doc = "maximum value of liquid water fraction omega for which softness values are parameterized by :cite:`LliboutryDuval1985`; used in Glen-Paterson-Budd-Lliboutry-Duval flow law; compare :cite:`AschwandenBuelerKhroulevBlatter`"; 938 pism_config:flow_law.gpbld.water_frac_observed_limit_type = "number"; 939 pism_config:flow_law.gpbld.water_frac_observed_limit_units = "1"; 940 941 pism_config:flow_law.isothermal_Glen.ice_softness = 3.1689e-24; 942 pism_config:flow_law.isothermal_Glen.ice_softness_doc = "ice softness used by IsothermalGlenIce :cite:`EISMINT96`"; 943 pism_config:flow_law.isothermal_Glen.ice_softness_type = "number"; 944 pism_config:flow_law.isothermal_Glen.ice_softness_units = "Pascal-3 second-1"; 945 946 pism_config:fracture_density.constant_fd = "no"; 947 pism_config:fracture_density.constant_fd_doc = "FIXME"; 948 pism_config:fracture_density.constant_fd_option = "constant_fd"; 949 pism_config:fracture_density.constant_fd_type = "flag"; 950 951 pism_config:fracture_density.constant_healing = "no"; 952 pism_config:fracture_density.constant_healing_doc = "Constant healing"; 953 pism_config:fracture_density.constant_healing_option = "constant_healing"; 954 pism_config:fracture_density.constant_healing_type = "flag"; 955 956 pism_config:fracture_density.enabled = "no"; 957 pism_config:fracture_density.enabled_doc = "Calculation of fracture density according to stresses and strain rate field."; 958 pism_config:fracture_density.enabled_option = "fractures"; 959 pism_config:fracture_density.enabled_type = "flag"; 960 961 pism_config:fracture_density.fd2d_scheme = "no"; 962 pism_config:fracture_density.fd2d_scheme_doc = "FIXME"; 963 pism_config:fracture_density.fd2d_scheme_option = "scheme_fd2d"; 964 pism_config:fracture_density.fd2d_scheme_type = "flag"; 965 966 pism_config:fracture_density.fracture_weighted_healing = "no"; 967 pism_config:fracture_density.fracture_weighted_healing_doc = "Fracture weighted healing"; 968 pism_config:fracture_density.fracture_weighted_healing_option = "fracture_weighted_healing"; 969 pism_config:fracture_density.fracture_weighted_healing_type = "flag"; 970 971 pism_config:fracture_density.include_grounded_ice = "no"; 972 pism_config:fracture_density.include_grounded_ice_doc = "model fracture density in grounded areas"; 973 pism_config:fracture_density.include_grounded_ice_option = "do_frac_on_grounded"; 974 pism_config:fracture_density.include_grounded_ice_type = "flag"; 975 976 pism_config:fracture_density.lefm = "no"; 977 pism_config:fracture_density.lefm_doc = "FIXME"; 978 pism_config:fracture_density.lefm_option = "lefm"; 979 pism_config:fracture_density.lefm_type = "flag"; 980 981 pism_config:fracture_density.max_shear_stress = "no"; 982 pism_config:fracture_density.max_shear_stress_doc = "Use the max. shear stress criterion."; 983 pism_config:fracture_density.max_shear_stress_option = "max_shear"; 984 pism_config:fracture_density.max_shear_stress_type = "flag"; 985 986 pism_config:fracture_density.phi0 = 0.0; 987 pism_config:fracture_density.phi0_doc = "FIXME"; 988 pism_config:fracture_density.phi0_option = "phi0"; 989 pism_config:fracture_density.phi0_type = "number"; 990 pism_config:fracture_density.phi0_units = "1"; 991 992 pism_config:fracture_density.softening_lower_limit = 1.0; 993 pism_config:fracture_density.softening_lower_limit_doc = "epsilon in equation (6) in Albrecht and Levermann, 'Fracture-induced softening for large-scale ice dynamics'"; 994 pism_config:fracture_density.softening_lower_limit_option = "fracture_softening"; 995 pism_config:fracture_density.softening_lower_limit_type = "number"; 996 pism_config:fracture_density.softening_lower_limit_units = "1"; 997 998 pism_config:frontal_melt.constant.melt_rate = 1.0; 999 pism_config:frontal_melt.constant.melt_rate_doc = "default melt rate used by the 'constant' frontal_melt model"; 1000 pism_config:frontal_melt.constant.melt_rate_option = "frontal_melt_rate"; 1001 pism_config:frontal_melt.constant.melt_rate_type = "number"; 1002 pism_config:frontal_melt.constant.melt_rate_units = "m / day"; 1003 1004 pism_config:frontal_melt.discharge_given.file = ""; 1005 pism_config:frontal_melt.discharge_given.file_doc = "Name of the file containing climate forcing fields."; 1006 pism_config:frontal_melt.discharge_given.file_option = "frontal_melt_discharge_given_file"; 1007 pism_config:frontal_melt.discharge_given.file_type = "string"; 1008 1009 pism_config:frontal_melt.discharge_given.period = 0; 1010 pism_config:frontal_melt.discharge_given.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1011 pism_config:frontal_melt.discharge_given.period_option = "frontal_melt_discharge_given_period"; 1012 pism_config:frontal_melt.discharge_given.period_type = "integer"; 1013 pism_config:frontal_melt.discharge_given.period_units = "years"; 1014 1015 pism_config:frontal_melt.discharge_given.reference_year = 0; 1016 pism_config:frontal_melt.discharge_given.reference_year_doc = "Reference year to use when ``frontal_melt.discharge_given.period`` is active."; 1017 pism_config:frontal_melt.discharge_given.reference_year_option = "frontal_melt_discharge_given_reference_year"; 1018 pism_config:frontal_melt.discharge_given.reference_year_type = "integer"; 1019 pism_config:frontal_melt.discharge_given.reference_year_units = "years"; 1020 1021 pism_config:frontal_melt.given.file = ""; 1022 pism_config:frontal_melt.given.file_doc = "Name of the file containing climate forcing fields."; 1023 pism_config:frontal_melt.given.file_option = "frontal_melt_given_file"; 1024 pism_config:frontal_melt.given.file_type = "string"; 1025 1026 pism_config:frontal_melt.given.period = 0; 1027 pism_config:frontal_melt.given.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1028 pism_config:frontal_melt.given.period_option = "frontal_melt_given_period"; 1029 pism_config:frontal_melt.given.period_type = "integer"; 1030 pism_config:frontal_melt.given.period_units = "years"; 1031 1032 pism_config:frontal_melt.given.reference_year = 0; 1033 pism_config:frontal_melt.given.reference_year_doc = "Reference year to use when ``frontal_melt.given.period`` is active."; 1034 pism_config:frontal_melt.given.reference_year_option = "frontal_melt_given_reference_year"; 1035 pism_config:frontal_melt.given.reference_year_type = "integer"; 1036 pism_config:frontal_melt.given.reference_year_units = "years"; 1037 1038 pism_config:frontal_melt.include_floating_ice = "no"; 1039 pism_config:frontal_melt.include_floating_ice_doc = "Apply frontal melt to all grid icy cells next to ocean cells"; 1040 pism_config:frontal_melt.include_floating_ice_type = "flag"; 1041 1042 pism_config:frontal_melt.models = ""; 1043 pism_config:frontal_melt.models_doc = "Comma-separated list of frontal melt models and modifiers. (Leave empty to disable.)"; 1044 pism_config:frontal_melt.models_option = "frontal_melt"; 1045 pism_config:frontal_melt.models_type = "string"; 1046 1047 pism_config:frontal_melt.routing.file = ""; 1048 pism_config:frontal_melt.routing.file_doc = "Name of the file containing climate forcing fields."; 1049 pism_config:frontal_melt.routing.file_option = "frontal_melt_routing_file"; 1050 pism_config:frontal_melt.routing.file_type = "string"; 1051 1052 pism_config:frontal_melt.routing.parameter_a = 3e-4; 1053 pism_config:frontal_melt.routing.parameter_a_doc = "parameter A in eqn. 1 in :cite:`Xu2013`"; 1054 pism_config:frontal_melt.routing.parameter_a_type = "number"; 1055 pism_config:frontal_melt.routing.parameter_a_units = "m-alpha day^(alpha-1) Celsius-beta"; 1056 1057 pism_config:frontal_melt.routing.parameter_b = 0.15; 1058 pism_config:frontal_melt.routing.parameter_b_doc = "parameter B in eqn. 1 in :cite:`Xu2013`"; 1059 pism_config:frontal_melt.routing.parameter_b_type = "number"; 1060 pism_config:frontal_melt.routing.parameter_b_units = "m day^(alpha-1) Celsius-beta"; 1061 1062 pism_config:frontal_melt.routing.period = 0; 1063 pism_config:frontal_melt.routing.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1064 pism_config:frontal_melt.routing.period_option = "frontal_melt_routing_period"; 1065 pism_config:frontal_melt.routing.period_type = "integer"; 1066 pism_config:frontal_melt.routing.period_units = "years"; 1067 1068 pism_config:frontal_melt.routing.power_alpha = 0.39; 1069 pism_config:frontal_melt.routing.power_alpha_doc = "exponent `\\alpha` in eqn. 1 in :cite:`Xu2013`"; 1070 pism_config:frontal_melt.routing.power_alpha_type = "number"; 1071 pism_config:frontal_melt.routing.power_alpha_units = "1"; 1072 1073 pism_config:frontal_melt.routing.power_beta = 1.18; 1074 pism_config:frontal_melt.routing.power_beta_doc = "exponent `\\beta` in eqn. 1 in :cite:`Xu2013`"; 1075 pism_config:frontal_melt.routing.power_beta_type = "number"; 1076 pism_config:frontal_melt.routing.power_beta_units = "1"; 1077 1078 pism_config:frontal_melt.routing.reference_year = 0; 1079 pism_config:frontal_melt.routing.reference_year_doc = "Reference year to use when ``frontal_melt.routing.period`` is active."; 1080 pism_config:frontal_melt.routing.reference_year_option = "frontal_melt_routing_reference_year"; 1081 pism_config:frontal_melt.routing.reference_year_type = "integer"; 1082 pism_config:frontal_melt.routing.reference_year_units = "years"; 1083 1084 pism_config:geometry.front_retreat.prescribed.file = ""; 1085 pism_config:geometry.front_retreat.prescribed.file_doc = "Name of the file containing the maximum ice extent mask `land_ice_area_fraction_retreat`"; 1086 pism_config:geometry.front_retreat.prescribed.file_option = "front_retreat_file"; 1087 pism_config:geometry.front_retreat.prescribed.file_type = "string"; 1088 1089 pism_config:geometry.front_retreat.prescribed.period = 0; 1090 pism_config:geometry.front_retreat.prescribed.period_doc = "Length of the period of the front retreat data. Set to zero to disable."; 1091 pism_config:geometry.front_retreat.prescribed.period_type = "integer"; 1092 pism_config:geometry.front_retreat.prescribed.period_units = "years"; 1093 1094 pism_config:geometry.front_retreat.prescribed.reference_year = 0; 1095 pism_config:geometry.front_retreat.prescribed.reference_year_doc = "Reference year to use when `geometry.front_retreat.prescribed.period` is active."; 1096 pism_config:geometry.front_retreat.prescribed.reference_year_type = "integer"; 1097 pism_config:geometry.front_retreat.prescribed.reference_year_units = "years"; 1098 1099 pism_config:geometry.front_retreat.use_cfl = "false"; 1100 pism_config:geometry.front_retreat.use_cfl_doc = "apply CFL criterion for eigen-calving rate front retreat"; 1101 pism_config:geometry.front_retreat.use_cfl_option = "front_retreat_cfl"; 1102 pism_config:geometry.front_retreat.use_cfl_type = "flag"; 1103 1104 pism_config:geometry.front_retreat.wrap_around = "false"; 1105 pism_config:geometry.front_retreat.wrap_around_doc = "If true, wrap around domain boundaries. This may be needed in some regional synthetic geometry setups."; 1106 pism_config:geometry.front_retreat.wrap_around_option = "front_retreat_wrap_around"; 1107 pism_config:geometry.front_retreat.wrap_around_type = "flag"; 1108 1109 pism_config:geometry.grounded_cell_fraction = "false"; 1110 pism_config:geometry.grounded_cell_fraction_doc = "Linear interpolation scheme ('LI' in Gladstone et al. 2010) expanded to two dimensions is used if switched on in order to evaluate the position of the grounding line on a subgrid scale."; 1111 pism_config:geometry.grounded_cell_fraction_option = "subgl"; 1112 pism_config:geometry.grounded_cell_fraction_type = "flag"; 1113 1114 pism_config:geometry.ice_free_thickness_standard = 0.01; 1115 pism_config:geometry.ice_free_thickness_standard_doc = "If ice is thinner than this standard then the mask is set to MASK_ICE_FREE_BEDROCK or MASK_ICE_FREE_OCEAN."; 1116 pism_config:geometry.ice_free_thickness_standard_type = "number"; 1117 pism_config:geometry.ice_free_thickness_standard_units = "meters"; 1118 1119 pism_config:geometry.part_grid.enabled = "no"; 1120 pism_config:geometry.part_grid.enabled_doc = "apply partially filled grid cell scheme"; 1121 pism_config:geometry.part_grid.enabled_option = "part_grid"; 1122 pism_config:geometry.part_grid.enabled_type = "flag"; 1123 1124 pism_config:geometry.part_grid.max_iterations = 10; 1125 pism_config:geometry.part_grid.max_iterations_doc = "maximum number of residual redistribution iterations"; 1126 pism_config:geometry.part_grid.max_iterations_type = "integer"; 1127 pism_config:geometry.part_grid.max_iterations_units = "count"; 1128 1129 pism_config:geometry.remove_icebergs = "no"; 1130 pism_config:geometry.remove_icebergs_doc = "identify and kill detached ice-shelf areas"; 1131 pism_config:geometry.remove_icebergs_option = "kill_icebergs"; 1132 pism_config:geometry.remove_icebergs_type = "flag"; 1133 1134 pism_config:geometry.update.enabled = "yes"; 1135 pism_config:geometry.update.enabled_doc = "Solve the mass conservation equation"; 1136 pism_config:geometry.update.enabled_option = "mass"; 1137 pism_config:geometry.update.enabled_type = "flag"; 1138 1139 pism_config:geometry.update.use_basal_melt_rate = "yes"; 1140 pism_config:geometry.update.use_basal_melt_rate_doc = "Include basal melt rate in the continuity equation"; 1141 pism_config:geometry.update.use_basal_melt_rate_option = "bmr_in_cont"; 1142 pism_config:geometry.update.use_basal_melt_rate_type = "flag"; 1143 1144 pism_config:grid.Lbz = 1000; 1145 pism_config:grid.Lbz_doc = "Thickness of the thermal bedrock layer. (Inactive if ``grid.Mbz`` < 2)"; 1146 pism_config:grid.Lbz_option = "Lbz"; 1147 pism_config:grid.Lbz_type = "number"; 1148 pism_config:grid.Lbz_units = "meters"; 1149 1150 pism_config:grid.Lx = 1500e3; 1151 pism_config:grid.Lx_doc = "Default computational box is 3000 km x 3000 km (= 2 Lx x 2 Ly) in horizontal."; 1152 pism_config:grid.Lx_type = "number"; 1153 pism_config:grid.Lx_units = "meters"; 1154 1155 pism_config:grid.Ly = 1500e3; 1156 pism_config:grid.Ly_doc = "Default computational box is 3000 km x 3000 km (= 2 Lx x 2 Ly) in horizontal."; 1157 pism_config:grid.Ly_type = "number"; 1158 pism_config:grid.Ly_units = "meters"; 1159 1160 pism_config:grid.Lz = 4000; 1161 pism_config:grid.Lz_doc = "Height of the computational domain."; 1162 pism_config:grid.Lz_option = "Lz"; 1163 pism_config:grid.Lz_type = "number"; 1164 pism_config:grid.Lz_units = "meters"; 1165 1166 pism_config:grid.Mbz = 1; 1167 pism_config:grid.Mbz_doc = "Number of thermal bedrock layers; 1 level corresponds to no bedrock."; 1168 pism_config:grid.Mbz_option = "Mbz"; 1169 pism_config:grid.Mbz_type = "integer"; 1170 pism_config:grid.Mbz_units = "count"; 1171 1172 pism_config:grid.Mx = 61; 1173 pism_config:grid.Mx_doc = "Number of grid points in the x direction."; 1174 pism_config:grid.Mx_option = "Mx"; 1175 pism_config:grid.Mx_type = "integer"; 1176 pism_config:grid.Mx_units = "count"; 1177 1178 pism_config:grid.My = 61; 1179 pism_config:grid.My_doc = "Number of grid points in the y direction."; 1180 pism_config:grid.My_option = "My"; 1181 pism_config:grid.My_type = "integer"; 1182 pism_config:grid.My_units = "count"; 1183 1184 pism_config:grid.Mz = 31; 1185 pism_config:grid.Mz_doc = "Number of vertical grid levels in the ice."; 1186 pism_config:grid.Mz_option = "Mz"; 1187 pism_config:grid.Mz_type = "integer"; 1188 pism_config:grid.Mz_units = "count"; 1189 1190 pism_config:grid.allow_extrapolation = "no"; 1191 pism_config:grid.allow_extrapolation_doc = "Allow extrapolation during regridding."; 1192 pism_config:grid.allow_extrapolation_option = "allow_extrapolation"; 1193 pism_config:grid.allow_extrapolation_type = "flag"; 1194 1195 pism_config:grid.ice_vertical_spacing = "quadratic"; 1196 pism_config:grid.ice_vertical_spacing_choices = "quadratic,equal"; 1197 pism_config:grid.ice_vertical_spacing_doc = "vertical spacing in the ice"; 1198 pism_config:grid.ice_vertical_spacing_option = "z_spacing"; 1199 pism_config:grid.ice_vertical_spacing_type = "keyword"; 1200 1201 pism_config:grid.lambda = 4.0; 1202 pism_config:grid.lambda_doc = "Vertical grid spacing parameter. Roughly equal to the factor by which the grid is coarser at an end away from the ice-bedrock interface."; 1203 pism_config:grid.lambda_type = "number"; 1204 pism_config:grid.lambda_units = "pure number"; 1205 1206 pism_config:grid.max_stencil_width = 2; 1207 pism_config:grid.max_stencil_width_doc = "Maximum width of the finite-difference stencil used in PISM."; 1208 pism_config:grid.max_stencil_width_type = "integer"; 1209 pism_config:grid.max_stencil_width_units = "count"; 1210 1211 pism_config:grid.periodicity = "xy"; 1212 pism_config:grid.periodicity_choices = "none,x,y,xy"; 1213 pism_config:grid.periodicity_doc = "horizontal grid periodicity"; 1214 pism_config:grid.periodicity_option = "periodicity"; 1215 pism_config:grid.periodicity_type = "keyword"; 1216 1217 pism_config:grid.recompute_longitude_and_latitude = "yes"; 1218 pism_config:grid.recompute_longitude_and_latitude_doc = "Re-compute longitude and latitude using grid information and provided projection parameters. Requires PROJ."; 1219 pism_config:grid.recompute_longitude_and_latitude_type = "flag"; 1220 1221 pism_config:grid.registration = "center"; 1222 pism_config:grid.registration_choices = "center,corner"; 1223 pism_config:grid.registration_doc = "horizontal grid registration"; 1224 pism_config:grid.registration_type = "keyword"; 1225 1226 pism_config:hydrology.add_water_input_to_till_storage = "yes"; 1227 pism_config:hydrology.add_water_input_to_till_storage_doc = "Add surface input to water stored in till. If no it will be added to the transportable water."; 1228 pism_config:hydrology.add_water_input_to_till_storage_type = "flag"; 1229 1230 pism_config:hydrology.cavitation_opening_coefficient = 0.5; 1231 pism_config:hydrology.cavitation_opening_coefficient_doc = "c_1 in notes; coefficient of cavitation opening term in evolution of layer thickness in hydrology::Distributed"; 1232 pism_config:hydrology.cavitation_opening_coefficient_option = "hydrology_cavitation_opening_coefficient"; 1233 pism_config:hydrology.cavitation_opening_coefficient_type = "number"; 1234 pism_config:hydrology.cavitation_opening_coefficient_units = "meter-1"; 1235 1236 pism_config:hydrology.creep_closure_coefficient = 0.04; 1237 pism_config:hydrology.creep_closure_coefficient_doc = "c_2 in notes; coefficient of creep closure term in evolution of layer thickness in hydrology::Distributed"; 1238 pism_config:hydrology.creep_closure_coefficient_option = "hydrology_creep_closure_coefficient"; 1239 pism_config:hydrology.creep_closure_coefficient_type = "number"; 1240 pism_config:hydrology.creep_closure_coefficient_units = "pure number"; 1241 1242 pism_config:hydrology.distributed.init_p_from_steady = "no"; 1243 pism_config:hydrology.distributed.init_p_from_steady_doc = "if 'yes', initialize subglacial water pressure from P(W) formula that applies in steady state"; 1244 pism_config:hydrology.distributed.init_p_from_steady_option = "hydrology_init_p_form_steady"; 1245 pism_config:hydrology.distributed.init_p_from_steady_type = "flag"; 1246 1247 pism_config:hydrology.distributed.sliding_speed_file = ""; 1248 pism_config:hydrology.distributed.sliding_speed_file_doc = "name of the file containing velbase_mag, the basal sliding speed to use with :literal:`hydrology.distributed.init_p_from_steady`"; 1249 pism_config:hydrology.distributed.sliding_speed_file_option = "hydrology_sliding_speed_file"; 1250 pism_config:hydrology.distributed.sliding_speed_file_type = "string"; 1251 1252 pism_config:hydrology.gradient_power_in_flux = 1.5; 1253 pism_config:hydrology.gradient_power_in_flux_doc = "power `\\beta` in Darcy's law `q = - k W^{\\alpha} |\\nabla \\psi|^{\\beta-2} \\nabla \\psi`, for subglacial water layer; used by hydrology::Routing and hydrology::Distributed"; 1254 pism_config:hydrology.gradient_power_in_flux_option = "hydrology_gradient_power_in_flux"; 1255 pism_config:hydrology.gradient_power_in_flux_type = "number"; 1256 pism_config:hydrology.gradient_power_in_flux_units = "pure number"; 1257 1258 pism_config:hydrology.hydraulic_conductivity = 0.001; 1259 pism_config:hydrology.hydraulic_conductivity_doc = "= k in notes; lateral conductivity, in Darcy's law, for subglacial water layer; units depend on powers alpha = hydrology.thickness_power_in_flux and beta = hydrology_potential_gradient_power_in_flux; used by hydrology::Routing and hydrology::Distributed"; 1260 pism_config:hydrology.hydraulic_conductivity_option = "hydrology_hydraulic_conductivity"; 1261 pism_config:hydrology.hydraulic_conductivity_type = "number"; 1262 pism_config:hydrology.hydraulic_conductivity_units = "`m^{2 \\beta - \\alpha} s^{2 \\beta - 3} kg^{1-\\beta}`"; 1263 1264 pism_config:hydrology.maximum_time_step = 1.0; 1265 pism_config:hydrology.maximum_time_step_doc = "maximum allowed time step length used by hydrology::Routing and hydrology::Distributed"; 1266 pism_config:hydrology.maximum_time_step_type = "number"; 1267 pism_config:hydrology.maximum_time_step_units = "years"; 1268 1269 pism_config:hydrology.model = "null"; 1270 pism_config:hydrology.model_choices = "null,routing,steady,distributed"; 1271 pism_config:hydrology.model_doc = "Basal hydrology sub-model."; 1272 pism_config:hydrology.model_option = "hydrology"; 1273 pism_config:hydrology.model_type = "keyword"; 1274 1275 pism_config:hydrology.null_diffuse_till_water = "no"; 1276 pism_config:hydrology.null_diffuse_till_water_doc = "Diffuse stored till water laterally. See equation (11) of :cite:`BBssasliding`"; 1277 pism_config:hydrology.null_diffuse_till_water_type = "flag"; 1278 1279 pism_config:hydrology.null_diffusion_distance = 2e4; 1280 pism_config:hydrology.null_diffusion_distance_doc = "diffusion distance for till water thickness; see equation (11) in :cite:`BBssasliding`; only active if hydrology.null_diffuse_till_water is set"; 1281 pism_config:hydrology.null_diffusion_distance_type = "number"; 1282 pism_config:hydrology.null_diffusion_distance_units = "meters"; 1283 1284 pism_config:hydrology.null_diffusion_time = 1000.0; 1285 pism_config:hydrology.null_diffusion_time_doc = "diffusion time for till water thickness; see equation (11) in :cite:`BBssasliding`; only active if hydrology.null_diffuse_till_water is set"; 1286 pism_config:hydrology.null_diffusion_time_type = "number"; 1287 pism_config:hydrology.null_diffusion_time_units = "years"; 1288 1289 pism_config:hydrology.null_strip_width = -1.0; 1290 pism_config:hydrology.null_strip_width_doc = "if negative then mechanism is inactive; width of strip around computational domain in which water velocity and water amount are set to zero; used by hydrology::Routing and hydrology::Distributed"; 1291 pism_config:hydrology.null_strip_width_type = "number"; 1292 pism_config:hydrology.null_strip_width_units = "meters"; 1293 1294 pism_config:hydrology.regularizing_porosity = 0.01; 1295 pism_config:hydrology.regularizing_porosity_doc = "phi_0 in notes; regularizes pressure equation by multiplying time derivative term"; 1296 pism_config:hydrology.regularizing_porosity_option = "hydrology_regularizing_porosity"; 1297 pism_config:hydrology.regularizing_porosity_type = "number"; 1298 pism_config:hydrology.regularizing_porosity_units = "pure number"; 1299 1300 pism_config:hydrology.roughness_scale = 0.1; 1301 pism_config:hydrology.roughness_scale_doc = "W_r in notes; roughness scale determining maximum amount of cavitation opening in hydrology::Distributed"; 1302 pism_config:hydrology.roughness_scale_option = "hydrology_roughness_scale"; 1303 pism_config:hydrology.roughness_scale_type = "number"; 1304 pism_config:hydrology.roughness_scale_units = "meters"; 1305 1306 pism_config:hydrology.routing.include_floating_ice = "no"; 1307 pism_config:hydrology.routing.include_floating_ice_doc = "Route subglacial water under ice shelves. This may be appropriate if a shelf is close to floatation. Note that this has no effect on ice flow."; 1308 pism_config:hydrology.routing.include_floating_ice_type = "flag"; 1309 1310 pism_config:hydrology.steady.flux_update_interval = 1.0; 1311 pism_config:hydrology.steady.flux_update_interval_doc = "interval between updates of the steady state flux"; 1312 pism_config:hydrology.steady.flux_update_interval_type = "number"; 1313 pism_config:hydrology.steady.flux_update_interval_units = "years"; 1314 1315 pism_config:hydrology.steady.input_rate_scaling = 1e7; 1316 pism_config:hydrology.steady.input_rate_scaling_doc = "input rate scaling"; 1317 pism_config:hydrology.steady.input_rate_scaling_type = "number"; 1318 pism_config:hydrology.steady.input_rate_scaling_units = "seconds"; 1319 1320 pism_config:hydrology.steady.n_iterations = 7500; 1321 pism_config:hydrology.steady.n_iterations_doc = "maxinum number of iterations to use in while estimating steady-state water flux"; 1322 pism_config:hydrology.steady.n_iterations_type = "integer"; 1323 pism_config:hydrology.steady.n_iterations_units = "count"; 1324 1325 pism_config:hydrology.steady.potential_delta = 10000.0; 1326 pism_config:hydrology.steady.potential_delta_doc = "potential adjustment used to fill sinks (smaller values require more iterations but produce fewer artifacts)"; 1327 pism_config:hydrology.steady.potential_delta_type = "number"; 1328 pism_config:hydrology.steady.potential_delta_units = "Pa"; 1329 1330 pism_config:hydrology.steady.potential_n_iterations = 1000; 1331 pism_config:hydrology.steady.potential_n_iterations_doc = "maxinum number of iterations to take while pre-processing hydraulic potential"; 1332 pism_config:hydrology.steady.potential_n_iterations_type = "integer"; 1333 pism_config:hydrology.steady.potential_n_iterations_units = "count"; 1334 1335 pism_config:hydrology.steady.volume_ratio = 0.1; 1336 pism_config:hydrology.steady.volume_ratio_doc = "water volume ratio used as the stopping criterion"; 1337 pism_config:hydrology.steady.volume_ratio_type = "number"; 1338 pism_config:hydrology.steady.volume_ratio_units = "1"; 1339 1340 pism_config:hydrology.surface_input.file = ""; 1341 pism_config:hydrology.surface_input.file_doc = "Name of the file containing ``water_input_rate``, the rate at which water from the ice surface is added to the subglacial hydrology system"; 1342 pism_config:hydrology.surface_input.file_type = "string"; 1343 1344 pism_config:hydrology.surface_input.period = 0; 1345 pism_config:hydrology.surface_input.period_doc = "Length of the period of the water input rate. Set to zero to disable."; 1346 pism_config:hydrology.surface_input.period_type = "integer"; 1347 pism_config:hydrology.surface_input.period_units = "years"; 1348 1349 pism_config:hydrology.surface_input.reference_year = 0; 1350 pism_config:hydrology.surface_input.reference_year_doc = "Reference year to use when ``hydrology.surface_input.period`` is active."; 1351 pism_config:hydrology.surface_input.reference_year_type = "integer"; 1352 pism_config:hydrology.surface_input.reference_year_units = "years"; 1353 1354 pism_config:hydrology.surface_input_from_runoff = "no"; 1355 pism_config:hydrology.surface_input_from_runoff_doc = "Use surface runoff as surface input."; 1356 pism_config:hydrology.surface_input_from_runoff_type = "flag"; 1357 1358 pism_config:hydrology.thickness_power_in_flux = 1.25; 1359 pism_config:hydrology.thickness_power_in_flux_doc = "power `\\alpha` in Darcy's law `q = - k W^{\\alpha} |\\nabla \\psi|^{\\beta-2} \\nabla \\psi`, for subglacial water layer; used by hydrology::Routing and hydrology::Distributed"; 1360 pism_config:hydrology.thickness_power_in_flux_option = "hydrology_thickness_power_in_flux"; 1361 pism_config:hydrology.thickness_power_in_flux_type = "number"; 1362 pism_config:hydrology.thickness_power_in_flux_units = "1"; 1363 1364 pism_config:hydrology.tillwat_decay_rate = 1.0; 1365 pism_config:hydrology.tillwat_decay_rate_doc = "rate at which tillwat is reduced to zero, in absence of other effects like input"; 1366 pism_config:hydrology.tillwat_decay_rate_option = "hydrology_tillwat_decay_rate"; 1367 pism_config:hydrology.tillwat_decay_rate_type = "number"; 1368 pism_config:hydrology.tillwat_decay_rate_units = "mm / year"; 1369 1370 pism_config:hydrology.tillwat_max = 2.0; 1371 pism_config:hydrology.tillwat_max_doc = "maximum effective thickness of the water stored in till"; 1372 pism_config:hydrology.tillwat_max_option = "hydrology_tillwat_max"; 1373 pism_config:hydrology.tillwat_max_type = "number"; 1374 pism_config:hydrology.tillwat_max_units = "meters"; 1375 1376 pism_config:hydrology.use_const_bmelt = "no"; 1377 pism_config:hydrology.use_const_bmelt_doc = "if 'yes', subglacial hydrology model sees basal melt rate which is constant and given by hydrology.const_bmelt"; 1378 pism_config:hydrology.use_const_bmelt_option = "hydrology_use_const_bmelt"; 1379 pism_config:hydrology.use_const_bmelt_type = "flag"; 1380 1381 pism_config:input.bootstrap = "no"; 1382 pism_config:input.bootstrap_doc = "It true, use bootstrapping heuristics when initializing PISM."; 1383 pism_config:input.bootstrap_option = "bootstrap"; 1384 pism_config:input.bootstrap_type = "flag"; 1385 1386 pism_config:input.file = ""; 1387 pism_config:input.file_doc = "Input file name"; 1388 pism_config:input.file_option = "i"; 1389 pism_config:input.file_type = "string"; 1390 1391 pism_config:input.forcing.buffer_size = 60; 1392 pism_config:input.forcing.buffer_size_doc = "number of 2D climate forcing records to keep in memory; = 5 years of monthly records"; 1393 pism_config:input.forcing.buffer_size_type = "integer"; 1394 pism_config:input.forcing.buffer_size_units = "count"; 1395 1396 pism_config:input.forcing.evaluations_per_year = 52; 1397 pism_config:input.forcing.evaluations_per_year_doc = "length of the time-series used to compute temporal averages of forcing data (such as mean annual temperature)"; 1398 pism_config:input.forcing.evaluations_per_year_type = "integer"; 1399 pism_config:input.forcing.evaluations_per_year_units = "count"; 1400 1401 pism_config:input.regrid.file = ""; 1402 pism_config:input.regrid.file_doc = "Regridding (input) file name"; 1403 pism_config:input.regrid.file_option = "regrid_file"; 1404 pism_config:input.regrid.file_type = "string"; 1405 1406 pism_config:input.regrid.vars = ""; 1407 pism_config:input.regrid.vars_doc = "Comma-separated list of variables to regrid. Leave empty to regrid all model state variables."; 1408 pism_config:input.regrid.vars_option = "regrid_vars"; 1409 pism_config:input.regrid.vars_type = "string"; 1410 1411 pism_config:inverse.design.cH1 = 0; 1412 pism_config:inverse.design.cH1_doc = "weight of derivative part of an H1 norm for inversion design variables"; 1413 pism_config:inverse.design.cH1_option = "inv_design_cH1"; 1414 pism_config:inverse.design.cH1_type = "number"; 1415 pism_config:inverse.design.cH1_units = "1"; 1416 1417 pism_config:inverse.design.cL2 = 1; 1418 pism_config:inverse.design.cL2_doc = "weight of derivative-free part of an H1 norm for inversion design variables"; 1419 pism_config:inverse.design.cL2_option = "inv_design_cL2"; 1420 pism_config:inverse.design.cL2_type = "number"; 1421 pism_config:inverse.design.cL2_units = "1"; 1422 1423 pism_config:inverse.design.func = "sobolevH1"; 1424 pism_config:inverse.design.func_choices = "sobolevH1,tv"; 1425 pism_config:inverse.design.func_doc = "functional used for inversion design variables"; 1426 pism_config:inverse.design.func_option = "inv_design_func"; 1427 pism_config:inverse.design.func_type = "keyword"; 1428 1429 pism_config:inverse.design.param = "exp"; 1430 pism_config:inverse.design.param_choices = "ident,trunc,square,exp"; 1431 pism_config:inverse.design.param_doc = "parameterization of design variables used during inversion"; 1432 pism_config:inverse.design.param_option = "inv_design_param"; 1433 pism_config:inverse.design.param_type = "keyword"; 1434 1435 pism_config:inverse.design.param_hardav_eps = 1e4; 1436 pism_config:inverse.design.param_hardav_eps_doc = "tiny vertically-averaged hardness used as a substitute for 0 in some tauc parameterizations"; 1437 pism_config:inverse.design.param_hardav_eps_type = "number"; 1438 pism_config:inverse.design.param_hardav_eps_units = "Pascal second^(1/3)"; 1439 1440 pism_config:inverse.design.param_hardav_scale = 1e8; 1441 pism_config:inverse.design.param_hardav_scale_doc = "typical size of ice hardness"; 1442 pism_config:inverse.design.param_hardav_scale_type = "number"; 1443 pism_config:inverse.design.param_hardav_scale_units = "Pascal second^(1/3)"; 1444 1445 pism_config:inverse.design.param_tauc_eps = 100; 1446 pism_config:inverse.design.param_tauc_eps_doc = "tiny yield stress used as a substitute for 0 in some tauc parameterizations"; 1447 pism_config:inverse.design.param_tauc_eps_type = "number"; 1448 pism_config:inverse.design.param_tauc_eps_units = "Pascal"; 1449 1450 pism_config:inverse.design.param_tauc_scale = 100000; 1451 pism_config:inverse.design.param_tauc_scale_doc = "typical size of yield stresses"; 1452 pism_config:inverse.design.param_tauc_scale_type = "number"; 1453 pism_config:inverse.design.param_tauc_scale_units = "Pascal"; 1454 1455 pism_config:inverse.design.param_trunc_hardav0 = 1e6; 1456 pism_config:inverse.design.param_trunc_hardav0_doc = "transition point of change to linear behaviour for design variable parameterization type 'trunc'"; 1457 pism_config:inverse.design.param_trunc_hardav0_type = "number"; 1458 pism_config:inverse.design.param_trunc_hardav0_units = "Pascal second^(1/3)"; 1459 1460 pism_config:inverse.design.param_trunc_tauc0 = 1000; 1461 pism_config:inverse.design.param_trunc_tauc0_doc = "transition point of change to linear behaviour for design variable parameterization type 'trunc'"; 1462 pism_config:inverse.design.param_trunc_tauc0_type = "number"; 1463 pism_config:inverse.design.param_trunc_tauc0_units = "Pascal"; 1464 1465 pism_config:inverse.log_ratio_scale = 10; 1466 pism_config:inverse.log_ratio_scale_doc = "Reference scale for log-ratio functionals"; 1467 pism_config:inverse.log_ratio_scale_option = "inv_log_ratio_scale"; 1468 pism_config:inverse.log_ratio_scale_type = "number"; 1469 pism_config:inverse.log_ratio_scale_units = "pure number"; 1470 1471 pism_config:inverse.max_iterations = 1000; 1472 pism_config:inverse.max_iterations_doc = "maximum iteration count"; 1473 pism_config:inverse.max_iterations_option = "inv_max_it"; 1474 pism_config:inverse.max_iterations_type = "integer"; 1475 pism_config:inverse.max_iterations_units = "count"; 1476 1477 pism_config:inverse.ssa.hardav_max = 1e10; 1478 pism_config:inverse.ssa.hardav_max_doc = "Maximum allowed value of hardav for inversions with bound constraints"; 1479 pism_config:inverse.ssa.hardav_max_type = "number"; 1480 pism_config:inverse.ssa.hardav_max_units = "Pascal second^(1/3)"; 1481 1482 pism_config:inverse.ssa.hardav_min = 0; 1483 pism_config:inverse.ssa.hardav_min_doc = "Minimum allowed value of hardav for inversions with bound constraints"; 1484 pism_config:inverse.ssa.hardav_min_type = "number"; 1485 pism_config:inverse.ssa.hardav_min_units = "Pascal second^(1/3)"; 1486 1487 pism_config:inverse.ssa.length_scale = 50e3; 1488 pism_config:inverse.ssa.length_scale_doc = "typical length scale for rescaling derivative norms"; 1489 pism_config:inverse.ssa.length_scale_type = "number"; 1490 pism_config:inverse.ssa.length_scale_units = "meters"; 1491 1492 pism_config:inverse.ssa.method = "tikhonov_lmvm"; 1493 pism_config:inverse.ssa.method_choices = "sd,nlcg,ign,tikhonov_lmvm,tikhonov_cg,tikhonov_blmvm,tikhonov_lcl,tikhonov_gn"; 1494 pism_config:inverse.ssa.method_doc = "algorithm to use for SSA inversions"; 1495 pism_config:inverse.ssa.method_option = "inv_method"; 1496 pism_config:inverse.ssa.method_type = "keyword"; 1497 1498 pism_config:inverse.ssa.tauc_max = 5e7; 1499 pism_config:inverse.ssa.tauc_max_doc = "Maximum allowed value of tauc for inversions with bound constraints"; 1500 pism_config:inverse.ssa.tauc_max_type = "number"; 1501 pism_config:inverse.ssa.tauc_max_units = "Pascal"; 1502 1503 pism_config:inverse.ssa.tauc_min = 0; 1504 pism_config:inverse.ssa.tauc_min_doc = "Minimum allowed value of tauc for inversions with bound constraints"; 1505 pism_config:inverse.ssa.tauc_min_type = "number"; 1506 pism_config:inverse.ssa.tauc_min_units = "Pascal"; 1507 1508 pism_config:inverse.ssa.tv_exponent = 1.2; 1509 pism_config:inverse.ssa.tv_exponent_doc = "Lebesgue exponent for pseudo-TV norm"; 1510 pism_config:inverse.ssa.tv_exponent_option = "inv_ssa_tv_exponent"; 1511 pism_config:inverse.ssa.tv_exponent_type = "number"; 1512 pism_config:inverse.ssa.tv_exponent_units = "pure number"; 1513 1514 pism_config:inverse.ssa.velocity_eps = 0.1; 1515 pism_config:inverse.ssa.velocity_eps_doc = "tiny size of ice velocities during inversion"; 1516 pism_config:inverse.ssa.velocity_eps_type = "number"; 1517 pism_config:inverse.ssa.velocity_eps_units = "meter / year"; 1518 1519 pism_config:inverse.ssa.velocity_scale = 100; 1520 pism_config:inverse.ssa.velocity_scale_doc = "typical size of ice velocities expected during inversion"; 1521 pism_config:inverse.ssa.velocity_scale_type = "number"; 1522 pism_config:inverse.ssa.velocity_scale_units = "meter / year"; 1523 1524 pism_config:inverse.state_func = "meansquare"; 1525 pism_config:inverse.state_func_choices = "meansquare,log_ratio,log_relative"; 1526 pism_config:inverse.state_func_doc = "functional used for inversion design variables"; 1527 pism_config:inverse.state_func_option = "inv_state_func"; 1528 pism_config:inverse.state_func_type = "keyword"; 1529 1530 pism_config:inverse.target_misfit = 100; 1531 pism_config:inverse.target_misfit_doc = "desired root misfit for SSA inversions"; 1532 pism_config:inverse.target_misfit_option = "inv_target_misfit"; 1533 pism_config:inverse.target_misfit_type = "number"; 1534 pism_config:inverse.target_misfit_units = "meter / year"; 1535 1536 pism_config:inverse.tikhonov.atol = 1e-10; 1537 pism_config:inverse.tikhonov.atol_doc = "absolute threshold for Tikhonov stopping criterion"; 1538 pism_config:inverse.tikhonov.atol_option = "tikhonov_atol"; 1539 pism_config:inverse.tikhonov.atol_type = "number"; 1540 pism_config:inverse.tikhonov.atol_units = "meter / year"; 1541 1542 pism_config:inverse.tikhonov.penalty_weight = 1; 1543 pism_config:inverse.tikhonov.penalty_weight_doc = "penalty parameter for Tikhonov inversion"; 1544 pism_config:inverse.tikhonov.penalty_weight_option = "tikhonov_penalty"; 1545 pism_config:inverse.tikhonov.penalty_weight_type = "number"; 1546 pism_config:inverse.tikhonov.penalty_weight_units = "1"; 1547 1548 pism_config:inverse.tikhonov.ptol = 0.1; 1549 pism_config:inverse.tikhonov.ptol_doc = "threshold for reaching desired misfit for adaptive Tikhonov algorithms"; 1550 pism_config:inverse.tikhonov.ptol_option = "tikhonov_ptol"; 1551 pism_config:inverse.tikhonov.ptol_type = "number"; 1552 pism_config:inverse.tikhonov.ptol_units = "pure number"; 1553 1554 pism_config:inverse.tikhonov.rtol = 5e-2; 1555 pism_config:inverse.tikhonov.rtol_doc = "relative threshold for Tikhonov stopping criterion"; 1556 pism_config:inverse.tikhonov.rtol_option = "tikhonov_rtol"; 1557 pism_config:inverse.tikhonov.rtol_type = "number"; 1558 pism_config:inverse.tikhonov.rtol_units = "1"; 1559 1560 pism_config:inverse.use_design_prior = "yes"; 1561 pism_config:inverse.use_design_prior_doc = "Use prior from inverse data file as initial guess."; 1562 pism_config:inverse.use_design_prior_option = "inv_use_design_prior"; 1563 pism_config:inverse.use_design_prior_type = "flag"; 1564 1565 pism_config:inverse.use_zeta_fixed_mask = "yes"; 1566 pism_config:inverse.use_zeta_fixed_mask_doc = "Enforce locations where the parameterized design variable should be fixed. (Automatically determined if not provided)"; 1567 pism_config:inverse.use_zeta_fixed_mask_option = "inv_use_zeta_fixed_mask"; 1568 pism_config:inverse.use_zeta_fixed_mask_type = "flag"; 1569 1570 pism_config:ocean.always_grounded = "no"; 1571 pism_config:ocean.always_grounded_doc = "Dry (ocean-less) simulation; ice is considered grounded regardless of ice thickness, bed elevation, and sea level."; 1572 pism_config:ocean.always_grounded_option = "dry"; 1573 pism_config:ocean.always_grounded_type = "flag"; 1574 1575 pism_config:ocean.anomaly.file = ""; 1576 pism_config:ocean.anomaly.file_doc = "Name of the file containing shelf basal mass flux offset fields."; 1577 pism_config:ocean.anomaly.file_option = "ocean_anomaly_file"; 1578 pism_config:ocean.anomaly.file_type = "string"; 1579 1580 pism_config:ocean.anomaly.period = 0; 1581 pism_config:ocean.anomaly.period_doc = "Length of the period of the ocean forcing data. Set to zero to disable."; 1582 pism_config:ocean.anomaly.period_option = "ocean_anomaly_period"; 1583 pism_config:ocean.anomaly.period_type = "integer"; 1584 pism_config:ocean.anomaly.period_units = "years"; 1585 1586 pism_config:ocean.anomaly.reference_year = 0; 1587 pism_config:ocean.anomaly.reference_year_doc = "Reference year to use when ``ocean.anomaly.period`` is active."; 1588 pism_config:ocean.anomaly.reference_year_option = "ocean_anomaly_reference_year"; 1589 pism_config:ocean.anomaly.reference_year_type = "integer"; 1590 pism_config:ocean.anomaly.reference_year_units = "years"; 1591 1592 pism_config:ocean.cache.update_interval = 10; 1593 pism_config:ocean.cache.update_interval_doc = "update interval of the 'cache' ocean modifier"; 1594 pism_config:ocean.cache.update_interval_option = "ocean_cache_update_interval"; 1595 pism_config:ocean.cache.update_interval_type = "integer"; 1596 pism_config:ocean.cache.update_interval_units = "years"; 1597 1598 pism_config:ocean.constant.melange_back_pressure_fraction = 0.0; 1599 pism_config:ocean.constant.melange_back_pressure_fraction_doc = "default melange back pressure fraction in the 'contant' ocean model"; 1600 pism_config:ocean.constant.melange_back_pressure_fraction_type = "number"; 1601 pism_config:ocean.constant.melange_back_pressure_fraction_units = "1"; 1602 1603 pism_config:ocean.constant.melt_rate = 0.05191419359084029; 1604 pism_config:ocean.constant.melt_rate_doc = "default melt rate used by the 'constant' ocean model (computed as `Q / (L \\rho_i)`)"; 1605 pism_config:ocean.constant.melt_rate_option = "shelf_base_melt_rate"; 1606 pism_config:ocean.constant.melt_rate_type = "number"; 1607 pism_config:ocean.constant.melt_rate_units = "m / year"; 1608 1609 pism_config:ocean.delta_T.file = ""; 1610 pism_config:ocean.delta_T.file_doc = "Name of the file containing temperature offsets."; 1611 pism_config:ocean.delta_T.file_option = "ocean_delta_T_file"; 1612 pism_config:ocean.delta_T.file_type = "string"; 1613 1614 pism_config:ocean.delta_T.period = 0; 1615 pism_config:ocean.delta_T.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1616 pism_config:ocean.delta_T.period_option = "ocean_delta_T_period"; 1617 pism_config:ocean.delta_T.period_type = "integer"; 1618 pism_config:ocean.delta_T.period_units = "years"; 1619 1620 pism_config:ocean.delta_T.reference_year = 0; 1621 pism_config:ocean.delta_T.reference_year_doc = "Reference year to use when ``ocean.delta_T.period`` is active."; 1622 pism_config:ocean.delta_T.reference_year_option = "ocean_delta_T_reference_year"; 1623 pism_config:ocean.delta_T.reference_year_type = "integer"; 1624 pism_config:ocean.delta_T.reference_year_units = "years"; 1625 1626 pism_config:ocean.delta_mass_flux.file = ""; 1627 pism_config:ocean.delta_mass_flux.file_doc = "Name of the file containing sub-shelf mass flux offsets."; 1628 pism_config:ocean.delta_mass_flux.file_option = "ocean_delta_mass_flux_file"; 1629 pism_config:ocean.delta_mass_flux.file_type = "string"; 1630 1631 pism_config:ocean.delta_mass_flux.period = 0; 1632 pism_config:ocean.delta_mass_flux.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1633 pism_config:ocean.delta_mass_flux.period_option = "ocean_delta_mass_flux_period"; 1634 pism_config:ocean.delta_mass_flux.period_type = "integer"; 1635 pism_config:ocean.delta_mass_flux.period_units = "years"; 1636 1637 pism_config:ocean.delta_mass_flux.reference_year = 0; 1638 pism_config:ocean.delta_mass_flux.reference_year_doc = "Reference year to use when ``ocean.delta_mass_flux.period`` is active."; 1639 pism_config:ocean.delta_mass_flux.reference_year_option = "ocean_delta_mass_flux_reference_year"; 1640 pism_config:ocean.delta_mass_flux.reference_year_type = "integer"; 1641 pism_config:ocean.delta_mass_flux.reference_year_units = "years"; 1642 1643 pism_config:ocean.delta_sl.file = ""; 1644 pism_config:ocean.delta_sl.file_doc = "Name of the file containing sea level offsets."; 1645 pism_config:ocean.delta_sl.file_option = "ocean_delta_sl_file"; 1646 pism_config:ocean.delta_sl.file_type = "string"; 1647 1648 pism_config:ocean.delta_sl.period = 0; 1649 pism_config:ocean.delta_sl.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1650 pism_config:ocean.delta_sl.period_option = "ocean_delta_sl_period"; 1651 pism_config:ocean.delta_sl.period_type = "integer"; 1652 pism_config:ocean.delta_sl.period_units = "years"; 1653 1654 pism_config:ocean.delta_sl.reference_year = 0; 1655 pism_config:ocean.delta_sl.reference_year_doc = "Reference year to use when ``ocean.delta_sl.period`` is active."; 1656 pism_config:ocean.delta_sl.reference_year_option = "ocean_delta_sl_reference_year"; 1657 pism_config:ocean.delta_sl.reference_year_type = "integer"; 1658 pism_config:ocean.delta_sl.reference_year_units = "years"; 1659 1660 pism_config:ocean.delta_sl_2d.file = ""; 1661 pism_config:ocean.delta_sl_2d.file_doc = "Name of the file containing climate forcing fields."; 1662 pism_config:ocean.delta_sl_2d.file_type = "string"; 1663 1664 pism_config:ocean.delta_sl_2d.period = 0; 1665 pism_config:ocean.delta_sl_2d.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1666 pism_config:ocean.delta_sl_2d.period_type = "integer"; 1667 pism_config:ocean.delta_sl_2d.period_units = "years"; 1668 1669 pism_config:ocean.delta_sl_2d.reference_year = 0; 1670 pism_config:ocean.delta_sl_2d.reference_year_doc = "Reference year to use when ``ocean.delta_sl_2d.period`` is active."; 1671 pism_config:ocean.delta_sl_2d.reference_year_type = "integer"; 1672 pism_config:ocean.delta_sl_2d.reference_year_units = "years"; 1673 1674 pism_config:ocean.frac_MBP.file = ""; 1675 pism_config:ocean.frac_MBP.file_doc = "Name of the file containing melange back-pressure scaling."; 1676 pism_config:ocean.frac_MBP.file_option = "ocean_frac_MBP_file"; 1677 pism_config:ocean.frac_MBP.file_type = "string"; 1678 1679 pism_config:ocean.frac_MBP.period = 0; 1680 pism_config:ocean.frac_MBP.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1681 pism_config:ocean.frac_MBP.period_option = "ocean_frac_MBP_period"; 1682 pism_config:ocean.frac_MBP.period_type = "integer"; 1683 pism_config:ocean.frac_MBP.period_units = "years"; 1684 1685 pism_config:ocean.frac_MBP.reference_year = 0; 1686 pism_config:ocean.frac_MBP.reference_year_doc = "Reference year to use when ``ocean.frac_MBP.period`` is active."; 1687 pism_config:ocean.frac_MBP.reference_year_option = "ocean_frac_MBP_reference_year"; 1688 pism_config:ocean.frac_MBP.reference_year_type = "integer"; 1689 pism_config:ocean.frac_MBP.reference_year_units = "years"; 1690 1691 pism_config:ocean.frac_mass_flux.file = ""; 1692 pism_config:ocean.frac_mass_flux.file_doc = "Name of the file containing sub-shelf mass flux scaling."; 1693 pism_config:ocean.frac_mass_flux.file_option = "ocean_frac_mass_flux_file"; 1694 pism_config:ocean.frac_mass_flux.file_type = "string"; 1695 1696 pism_config:ocean.frac_mass_flux.period = 0; 1697 pism_config:ocean.frac_mass_flux.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1698 pism_config:ocean.frac_mass_flux.period_option = "ocean_frac_mass_flux_period"; 1699 pism_config:ocean.frac_mass_flux.period_type = "integer"; 1700 pism_config:ocean.frac_mass_flux.period_units = "years"; 1701 1702 pism_config:ocean.frac_mass_flux.reference_year = 0; 1703 pism_config:ocean.frac_mass_flux.reference_year_doc = "Reference year to use when ``ocean.frac_mass_flux.period`` is active."; 1704 pism_config:ocean.frac_mass_flux.reference_year_option = "ocean_frac_mass_flux_reference_year"; 1705 pism_config:ocean.frac_mass_flux.reference_year_type = "integer"; 1706 pism_config:ocean.frac_mass_flux.reference_year_units = "years"; 1707 1708 pism_config:ocean.given.file = ""; 1709 pism_config:ocean.given.file_doc = "Name of the file containing climate forcing fields."; 1710 pism_config:ocean.given.file_option = "ocean_given_file"; 1711 pism_config:ocean.given.file_type = "string"; 1712 1713 pism_config:ocean.given.period = 0; 1714 pism_config:ocean.given.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1715 pism_config:ocean.given.period_option = "ocean_given_period"; 1716 pism_config:ocean.given.period_type = "integer"; 1717 pism_config:ocean.given.period_units = "years"; 1718 1719 pism_config:ocean.given.reference_year = 0; 1720 pism_config:ocean.given.reference_year_doc = "Reference year to use when ``ocean.given.period`` is active."; 1721 pism_config:ocean.given.reference_year_option = "ocean_given_reference_year"; 1722 pism_config:ocean.given.reference_year_type = "integer"; 1723 pism_config:ocean.given.reference_year_units = "years"; 1724 1725 pism_config:ocean.melange_back_pressure_fraction = 0.0; 1726 pism_config:ocean.melange_back_pressure_fraction_doc = "default melange back pressure fraction"; 1727 pism_config:ocean.melange_back_pressure_fraction_type = "number"; 1728 pism_config:ocean.melange_back_pressure_fraction_units = "1"; 1729 1730 pism_config:ocean.models = "constant"; 1731 pism_config:ocean.models_doc = "Comma-separated list of ocean models and modifiers."; 1732 pism_config:ocean.models_option = "ocean"; 1733 pism_config:ocean.models_type = "string"; 1734 1735 pism_config:ocean.pico.continental_shelf_depth = -800.0; 1736 pism_config:ocean.pico.continental_shelf_depth_doc = "Determines the edge of the continental shelf to calculate mean ocean properties"; 1737 pism_config:ocean.pico.continental_shelf_depth_option = "continental_shelf_depth"; 1738 pism_config:ocean.pico.continental_shelf_depth_type = "number"; 1739 pism_config:ocean.pico.continental_shelf_depth_units = "meters"; 1740 1741 pism_config:ocean.pico.exclude_ice_rises = "yes"; 1742 pism_config:ocean.pico.exclude_ice_rises_doc = "Exclude ice rises for box calculation in PICO"; 1743 pism_config:ocean.pico.exclude_ice_rises_option = "exclude_icerises"; 1744 pism_config:ocean.pico.exclude_ice_rises_type = "flag"; 1745 1746 pism_config:ocean.pico.file = ""; 1747 pism_config:ocean.pico.file_doc = "Name of the file containing climate forcing fields."; 1748 pism_config:ocean.pico.file_option = "ocean_pico_file"; 1749 pism_config:ocean.pico.file_type = "string"; 1750 1751 pism_config:ocean.pico.heat_exchange_coefficent = 2e-5; 1752 pism_config:ocean.pico.heat_exchange_coefficent_doc = "Turbulent heat exchange coefficient"; 1753 pism_config:ocean.pico.heat_exchange_coefficent_option = "gamma_T"; 1754 pism_config:ocean.pico.heat_exchange_coefficent_type = "number"; 1755 pism_config:ocean.pico.heat_exchange_coefficent_units = "meters second-1"; 1756 1757 pism_config:ocean.pico.maximum_ice_rise_area = 1e5; 1758 pism_config:ocean.pico.maximum_ice_rise_area_doc = "Maximum ice rise area. Patches of grounded ice that are bugger than this are treated as continental ice sheets."; 1759 pism_config:ocean.pico.maximum_ice_rise_area_type = "number"; 1760 pism_config:ocean.pico.maximum_ice_rise_area_units = "km2"; 1761 1762 pism_config:ocean.pico.number_of_boxes = 5; 1763 pism_config:ocean.pico.number_of_boxes_doc = "Default maximum number of boxes (for the largest ice shelves)"; 1764 pism_config:ocean.pico.number_of_boxes_option = "number_of_boxes"; 1765 pism_config:ocean.pico.number_of_boxes_type = "integer"; 1766 pism_config:ocean.pico.number_of_boxes_units = "pure number"; 1767 1768 pism_config:ocean.pico.overturning_coefficent = 1e6; 1769 pism_config:ocean.pico.overturning_coefficent_doc = "Overturning stregth coefficient"; 1770 pism_config:ocean.pico.overturning_coefficent_option = "overturning_coeff"; 1771 pism_config:ocean.pico.overturning_coefficent_type = "number"; 1772 pism_config:ocean.pico.overturning_coefficent_units = "meters6 seconds-1 kg-1"; 1773 1774 pism_config:ocean.pico.period = 0; 1775 pism_config:ocean.pico.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1776 pism_config:ocean.pico.period_option = "ocean_pico_period"; 1777 pism_config:ocean.pico.period_type = "integer"; 1778 pism_config:ocean.pico.period_units = "years"; 1779 1780 pism_config:ocean.pico.reference_year = 0; 1781 pism_config:ocean.pico.reference_year_doc = "Reference year to use when ``ocean.pico.period`` is active."; 1782 pism_config:ocean.pico.reference_year_option = "ocean_pico_reference_year"; 1783 pism_config:ocean.pico.reference_year_type = "integer"; 1784 pism_config:ocean.pico.reference_year_units = "years"; 1785 1786 pism_config:ocean.pik_melt_factor = 5e-3; 1787 pism_config:ocean.pik_melt_factor_doc = "dimensionless tuning parameter in the '-ocean pik' ocean heat flux parameterization; see :cite:`Martinetal2011`"; 1788 pism_config:ocean.pik_melt_factor_option = "meltfactor_pik"; 1789 pism_config:ocean.pik_melt_factor_type = "number"; 1790 pism_config:ocean.pik_melt_factor_units = "1"; 1791 1792 pism_config:ocean.sub_shelf_heat_flux_into_ice = 0.5; 1793 pism_config:ocean.sub_shelf_heat_flux_into_ice_doc = "= J meter-2 second-1; naively chosen default value for heat from ocean; see comments in pism::ocean::Constant::shelf_base_mass_flux()."; 1794 pism_config:ocean.sub_shelf_heat_flux_into_ice_type = "number"; 1795 pism_config:ocean.sub_shelf_heat_flux_into_ice_units = "W meter-2"; 1796 1797 pism_config:ocean.th.file = ""; 1798 pism_config:ocean.th.file_doc = "Name of the file containing climate forcing fields."; 1799 pism_config:ocean.th.file_option = "ocean_th_file"; 1800 pism_config:ocean.th.file_type = "string"; 1801 1802 pism_config:ocean.th.period = 0; 1803 pism_config:ocean.th.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 1804 pism_config:ocean.th.period_option = "ocean_th_period"; 1805 pism_config:ocean.th.period_type = "integer"; 1806 pism_config:ocean.th.period_units = "years"; 1807 1808 pism_config:ocean.th.reference_year = 0; 1809 pism_config:ocean.th.reference_year_doc = "Reference year to use when ``ocean.th.period`` is active."; 1810 pism_config:ocean.th.reference_year_option = "ocean_th_reference_year"; 1811 pism_config:ocean.th.reference_year_type = "integer"; 1812 pism_config:ocean.th.reference_year_units = "years"; 1813 1814 pism_config:ocean.three_equation_model_clip_salinity = "yes"; 1815 pism_config:ocean.three_equation_model_clip_salinity_doc = "Clip shelf base salinity so that it is in the range [4, 40] k/kg. See :cite:`HollandJenkins1999`."; 1816 pism_config:ocean.three_equation_model_clip_salinity_option = "clip_shelf_base_salinity"; 1817 pism_config:ocean.three_equation_model_clip_salinity_type = "flag"; 1818 1819 pism_config:output.ISMIP6 = "false"; 1820 pism_config:output.ISMIP6_doc = "Follow ISMIP6 conventions (units, variable names, \"standard names\") when writing output variables."; 1821 pism_config:output.ISMIP6_type = "flag"; 1822 1823 pism_config:output.ISMIP6_extra_variables = "lithk,orog,topg,hfgeoubed,acabf,libmassbfgr,libmassbffl,dlithkdt,velsurf,zvelsurf,velbase,zvelbase,velmean,litemptop,litempbotgr,litempbotfl,strbasemag,licalvf,lifmassbf,sftgif,sftgrf,sftflf"; 1824 pism_config:output.ISMIP6_extra_variables_doc = "Comma-separated list of fields reported by models participating in ISMIP6 simulations."; 1825 pism_config:output.ISMIP6_extra_variables_type = "string"; 1826 1827 pism_config:output.ISMIP6_ts_variables = "lim,limnsw,iareagr,iareafl,tendacabf,tendlibmassbf,tendlibmassbffl,tendlicalvf,tendlifmassbf"; 1828 pism_config:output.ISMIP6_ts_variables_doc = "Comma-separated list of scalar variables (time series) reported by models participating in ISMIP6 simulations."; 1829 pism_config:output.ISMIP6_ts_variables_type = "string"; 1830 1831 pism_config:output.backup_interval = 1.0; 1832 pism_config:output.backup_interval_doc = "wall-clock time between automatic backups"; 1833 pism_config:output.backup_interval_option = "backup_interval"; 1834 pism_config:output.backup_interval_type = "number"; 1835 pism_config:output.backup_interval_units = "hours"; 1836 1837 pism_config:output.backup_size = "small"; 1838 pism_config:output.backup_size_choices = "none,small,medium,big_2d,big"; 1839 pism_config:output.backup_size_doc = "The 'size' of a backup file. See configuration parameters output.sizes.medium, output.sizes.big_2d, output.sizes.big"; 1840 pism_config:output.backup_size_option = "backup_size"; 1841 pism_config:output.backup_size_type = "keyword"; 1842 1843 pism_config:output.extra.append = "no"; 1844 pism_config:output.extra.append_doc = "Append to an existing output file."; 1845 pism_config:output.extra.append_option = "extra_append"; 1846 pism_config:output.extra.append_type = "flag"; 1847 1848 pism_config:output.extra.file = ""; 1849 pism_config:output.extra.file_doc = "Name of the output file containing spatially-variable diagnostics."; 1850 pism_config:output.extra.file_option = "extra_file"; 1851 pism_config:output.extra.file_type = "string"; 1852 1853 pism_config:output.extra.split = "no"; 1854 pism_config:output.extra.split_doc = "Save spatially-variable diagnostics to separate files (one per time record)."; 1855 pism_config:output.extra.split_option = "extra_split"; 1856 pism_config:output.extra.split_type = "flag"; 1857 1858 pism_config:output.extra.stop_missing = "yes"; 1859 pism_config:output.extra.stop_missing_doc = "Stop if requested variable is not available instead of warning."; 1860 pism_config:output.extra.stop_missing_option = "extra_stop_missing"; 1861 pism_config:output.extra.stop_missing_type = "flag"; 1862 1863 pism_config:output.extra.times = ""; 1864 pism_config:output.extra.times_doc = "List or a range of times defining reporting intervals for spatially-variable diagnostics."; 1865 pism_config:output.extra.times_option = "extra_times"; 1866 pism_config:output.extra.times_type = "string"; 1867 1868 pism_config:output.extra.vars = ""; 1869 pism_config:output.extra.vars_doc = "Comma-separated list of spatially-variable diagnostics."; 1870 pism_config:output.extra.vars_option = "extra_vars"; 1871 pism_config:output.extra.vars_type = "string"; 1872 1873 pism_config:output.file_name = "unnamed.nc"; 1874 pism_config:output.file_name_doc = "The file to save final model results to."; 1875 pism_config:output.file_name_option = "o"; 1876 pism_config:output.file_name_type = "string"; 1877 1878 pism_config:output.fill_value = -2e9; 1879 pism_config:output.fill_value_doc = "_FillValue used when saving diagnostic quantities"; 1880 pism_config:output.fill_value_type = "number"; 1881 pism_config:output.fill_value_units = "none"; 1882 1883 pism_config:output.format = "netcdf3"; 1884 pism_config:output.format_choices = "netcdf3,netcdf4_parallel,pnetcdf,pio_pnetcdf,pio_netcdf4p,pio_netcdf4c,pio_netcdf"; 1885 pism_config:output.format_doc = "The I/O format used for spatial fields; 'netcdf3' is the default, 'netcd4_parallel' is available if PISM was built with parallel NetCDF-4, and 'pnetcdf' is available if PISM was built with PnetCDF."; 1886 pism_config:output.format_option = "o_format"; 1887 pism_config:output.format_type = "keyword"; 1888 1889 pism_config:output.ice_free_thickness_standard = 10.0; 1890 pism_config:output.ice_free_thickness_standard_doc = "If ice is thinner than this standard then a grid cell is considered ice-free for purposes of reporting glacierized area, volume, etc."; 1891 pism_config:output.ice_free_thickness_standard_type = "number"; 1892 pism_config:output.ice_free_thickness_standard_units = "meters"; 1893 1894 pism_config:output.pio.base = 0; 1895 pism_config:output.pio.base_doc = "Rank of the first I/O task"; 1896 pism_config:output.pio.base_type = "integer"; 1897 pism_config:output.pio.base_units = "count"; 1898 1899 pism_config:output.pio.n_writers = 1; 1900 pism_config:output.pio.n_writers_doc = "Number of I/O tasks to use"; 1901 pism_config:output.pio.n_writers_type = "integer"; 1902 pism_config:output.pio.n_writers_units = "count"; 1903 1904 pism_config:output.pio.stride = 1; 1905 pism_config:output.pio.stride_doc = "Offset between I/O tasks"; 1906 pism_config:output.pio.stride_type = "integer"; 1907 pism_config:output.pio.stride_units = "count"; 1908 1909 pism_config:output.runtime.area_scale_factor_log10 = 6; 1910 pism_config:output.runtime.area_scale_factor_log10_doc = "an integer; log base 10 of scale factor to use for area (in km^2) in summary line to stdout"; 1911 pism_config:output.runtime.area_scale_factor_log10_option = "summary_area_scale_factor_log10"; 1912 pism_config:output.runtime.area_scale_factor_log10_type = "integer"; 1913 pism_config:output.runtime.area_scale_factor_log10_units = "pure number"; 1914 1915 pism_config:output.runtime.time_unit_name = "year"; 1916 pism_config:output.runtime.time_unit_name_doc = "Time units used when printing model time, time step, and maximum horizontal velocity at summary to stdout. Must be valid udunits for time. (E.g. choose from year,month,day,hour,minute,second.)"; 1917 pism_config:output.runtime.time_unit_name_type = "string"; 1918 1919 pism_config:output.runtime.time_use_calendar = "yes"; 1920 pism_config:output.runtime.time_use_calendar_doc = "Whether to use the current calendar when printing model time in summary to stdout."; 1921 pism_config:output.runtime.time_use_calendar_type = "flag"; 1922 1923 pism_config:output.runtime.viewer.size = 320; 1924 pism_config:output.runtime.viewer.size_doc = "default diagnostic viewer size (number of pixels of the longer side)"; 1925 pism_config:output.runtime.viewer.size_option = "view_size"; 1926 pism_config:output.runtime.viewer.size_type = "integer"; 1927 pism_config:output.runtime.viewer.size_units = "count"; 1928 1929 pism_config:output.runtime.viewer.variables = ""; 1930 pism_config:output.runtime.viewer.variables_doc = "comma-separated list of map-plane diagnostic quantities to view at runtime"; 1931 pism_config:output.runtime.viewer.variables_option = "view"; 1932 pism_config:output.runtime.viewer.variables_type = "string"; 1933 1934 pism_config:output.runtime.volume_scale_factor_log10 = 6; 1935 pism_config:output.runtime.volume_scale_factor_log10_doc = "an integer; log base 10 of scale factor to use for volume (in km^3) in summary line to stdout"; 1936 pism_config:output.runtime.volume_scale_factor_log10_option = "summary_vol_scale_factor_log10"; 1937 pism_config:output.runtime.volume_scale_factor_log10_type = "integer"; 1938 pism_config:output.runtime.volume_scale_factor_log10_units = "pure number"; 1939 1940 pism_config:output.size = "medium"; 1941 pism_config:output.size_choices = "none,small,medium,big_2d,big"; 1942 pism_config:output.size_doc = "The 'size' of an output file. See configuration parameters output.sizes.medium, output.sizes.big_2d, output.sizes.big"; 1943 pism_config:output.size_option = "o_size"; 1944 pism_config:output.size_type = "keyword"; 1945 1946 pism_config:output.sizes.big = "cts,liqfrac,temp,temp_pa,uvel,vvel,wvel,wvel_rel"; 1947 pism_config:output.sizes.big_doc = "Comma-separated list of variables to write to the output (in addition to model_state variables and variables listed in output.sizes.medium and output.sizes.big_2d) if 'big' output size is selected. Does not include fields written by sub-models."; 1948 pism_config:output.sizes.big_type = "string"; 1949 1950 pism_config:output.sizes.big_2d = "age,bfrict,bheatflx,bmelt,bwp,bwprel,dbdt,effbwp,enthalpybase,enthalpysurf,flux_divergence,hardav,hydroinput,lat,litho_temp,lon,nuH,rank,tempbase,tempicethk,tempicethk_basal,temppabase,tempsurf,thk,thksmooth,tillphi,topg,velbar,velbase,wallmelt,wvelbase"; 1951 pism_config:output.sizes.big_2d_doc = "Comma-separated list of variables to write to the output (in addition to model_state variables and variables listed in output.sizes.medium) if 'big_2d' output size is selected. Does not include fields written by boundary models."; 1952 pism_config:output.sizes.big_2d_type = "string"; 1953 1954 pism_config:output.sizes.medium = "bwat,bwatvel,climatic_mass_balance,diffusivity,enthalpy,flux,flux_mag,ice_surface_temp,liqfrac,mask,schoofs_theta,strain_rates,taub_mag,tauc,taud_mag,temp_pa,tillwat,topgsmooth,usurf,velbar_mag,velbase_mag,velsurf,velsurf_mag,wvelsurf"; 1955 pism_config:output.sizes.medium_doc = "Comma-separated list of variables to write to the output (in addition to model_state variables) if 'medium' output size (the default) is selected. Does not include fields written by sub-models."; 1956 pism_config:output.sizes.medium_type = "string"; 1957 1958 pism_config:output.snapshot.file = ""; 1959 pism_config:output.snapshot.file_doc = "Snapshot (output) file name (or prefix, if saving to individual files)."; 1960 pism_config:output.snapshot.file_option = "save_file"; 1961 pism_config:output.snapshot.file_type = "string"; 1962 1963 pism_config:output.snapshot.size = "small"; 1964 pism_config:output.snapshot.size_choices = "none,small,medium,big_2d,big"; 1965 pism_config:output.snapshot.size_doc = "The 'size' of a snapshot file. See configuration parameters output.sizes.medium, output.sizes.big_2d, output.sizes.big"; 1966 pism_config:output.snapshot.size_option = "save_size"; 1967 pism_config:output.snapshot.size_type = "keyword"; 1968 1969 pism_config:output.snapshot.split = "no"; 1970 pism_config:output.snapshot.split_doc = "Save model state snapshots to separate files (one per time record)."; 1971 pism_config:output.snapshot.split_option = "save_split"; 1972 pism_config:output.snapshot.split_type = "flag"; 1973 1974 pism_config:output.snapshot.times = ""; 1975 pism_config:output.snapshot.times_doc = "List or a range of times to save model state snapshots at."; 1976 pism_config:output.snapshot.times_option = "save_times"; 1977 pism_config:output.snapshot.times_type = "string"; 1978 1979 pism_config:output.timeseries.append = "false"; 1980 pism_config:output.timeseries.append_doc = "If true, append to the scalar time series output file."; 1981 pism_config:output.timeseries.append_option = "ts_append"; 1982 pism_config:output.timeseries.append_type = "flag"; 1983 1984 pism_config:output.timeseries.buffer_size = 10000; 1985 pism_config:output.timeseries.buffer_size_doc = "Number of scalar diagnostic time-series records to hold in memory before writing to disk. (PISM writes this many time-series records to reduce I/O costs.) Send the USR2 signal to flush time-series."; 1986 pism_config:output.timeseries.buffer_size_type = "integer"; 1987 pism_config:output.timeseries.buffer_size_units = "count"; 1988 1989 pism_config:output.timeseries.filename = ""; 1990 pism_config:output.timeseries.filename_doc = "Name of the file to save scalar time series to. Leave empty to disable reporting scalar time-series."; 1991 pism_config:output.timeseries.filename_option = "ts_file"; 1992 pism_config:output.timeseries.filename_type = "string"; 1993 1994 pism_config:output.timeseries.times = ""; 1995 pism_config:output.timeseries.times_doc = "List or range of times defining reporting time intervals."; 1996 pism_config:output.timeseries.times_option = "ts_times"; 1997 pism_config:output.timeseries.times_type = "string"; 1998 1999 pism_config:output.timeseries.variables = ""; 2000 pism_config:output.timeseries.variables_doc = "Requested scalar (time-series) diagnostics. Leave empty to save all available diagnostics."; 2001 pism_config:output.timeseries.variables_option = "ts_vars"; 2002 pism_config:output.timeseries.variables_type = "string"; 2003 2004 pism_config:output.use_MKS = "false"; 2005 pism_config:output.use_MKS_doc = "Use MKS units in output files."; 2006 pism_config:output.use_MKS_type = "flag"; 2007 2008 pism_config:regional.no_model_strip = 5.0; 2009 pism_config:regional.no_model_strip_doc = "Default width of the 'no model strip' in regional setups."; 2010 pism_config:regional.no_model_strip_option = "no_model_strip"; 2011 pism_config:regional.no_model_strip_type = "number"; 2012 pism_config:regional.no_model_strip_units = "km"; 2013 2014 pism_config:regional.no_model_yield_stress = 1000.0; 2015 pism_config:regional.no_model_yield_stress_doc = "High yield stress used in the 'no_model_mask' area in the regional mode."; 2016 pism_config:regional.no_model_yield_stress_type = "number"; 2017 pism_config:regional.no_model_yield_stress_units = "kPa"; 2018 2019 pism_config:regional.zero_gradient = "false"; 2020 pism_config:regional.zero_gradient_doc = "Use zero ice thickness and ice surface gradient in the no_model_mask area."; 2021 pism_config:regional.zero_gradient_option = "zero_grad_where_no_model"; 2022 pism_config:regional.zero_gradient_type = "flag"; 2023 2024 pism_config:run_info.institution = ""; 2025 pism_config:run_info.institution_doc = "Institution name. This string is written to output files as the 'institution' global attribute."; 2026 pism_config:run_info.institution_option = "institution"; 2027 pism_config:run_info.institution_type = "string"; 2028 2029 pism_config:run_info.title = ""; 2030 pism_config:run_info.title_doc = "Free-form string containing a concise description of the current run. This string is written to output files as the 'title' global attribute."; 2031 pism_config:run_info.title_option = "title"; 2032 pism_config:run_info.title_type = "string"; 2033 2034 pism_config:sea_level.models = "constant"; 2035 pism_config:sea_level.models_doc = "Comma-separated list of sea level models and modifiers."; 2036 pism_config:sea_level.models_option = "sea_level"; 2037 pism_config:sea_level.models_type = "string"; 2038 2039 pism_config:stress_balance.calving_front_stress_bc = "no"; 2040 pism_config:stress_balance.calving_front_stress_bc_doc = "Apply CFBC condition as in :cite:`Albrechtetal2011`, :cite:`Winkelmannetal2011`. May only apply to some stress balances; e.g. SSAFD as of May 2011. If not set then a strength-extension is used, as in :cite:`BBssasliding`."; 2041 pism_config:stress_balance.calving_front_stress_bc_option = "cfbc"; 2042 pism_config:stress_balance.calving_front_stress_bc_type = "flag"; 2043 2044 pism_config:stress_balance.ice_free_thickness_standard = 10.0; 2045 pism_config:stress_balance.ice_free_thickness_standard_doc = "If ice is thinner than this standard then a cell is considered ice-free for purposes of computing ice velocity distribution."; 2046 pism_config:stress_balance.ice_free_thickness_standard_type = "number"; 2047 pism_config:stress_balance.ice_free_thickness_standard_units = "meters"; 2048 2049 pism_config:stress_balance.model = "sia"; 2050 pism_config:stress_balance.model_choices = "none,prescribed_sliding,weertman_sliding,sia,ssa,prescribed_sliding+sia,weertman_sliding+sia,ssa+sia"; 2051 pism_config:stress_balance.model_doc = "Stress balance model"; 2052 pism_config:stress_balance.model_option = "stress_balance"; 2053 pism_config:stress_balance.model_type = "keyword"; 2054 2055 pism_config:stress_balance.prescribed_sliding.file = ""; 2056 pism_config:stress_balance.prescribed_sliding.file_doc = "The name of the file containing prescribed sliding velocity (variable names: `ubar`, `vbar`)."; 2057 pism_config:stress_balance.prescribed_sliding.file_type = "string"; 2058 2059 pism_config:stress_balance.sia.Glen_exponent = 3.0; 2060 pism_config:stress_balance.sia.Glen_exponent_doc = "Glen exponent in ice flow law for SIA"; 2061 pism_config:stress_balance.sia.Glen_exponent_option = "sia_n"; 2062 pism_config:stress_balance.sia.Glen_exponent_type = "number"; 2063 pism_config:stress_balance.sia.Glen_exponent_units = "pure number"; 2064 2065 pism_config:stress_balance.sia.bed_smoother.range = 5.0e3; 2066 pism_config:stress_balance.sia.bed_smoother.range_doc = "half-width of smoothing domain for stressbalance::BedSmoother, in implementing :cite:`Schoofbasaltopg2003` bed roughness parameterization for SIA; set value to zero to turn off mechanism"; 2067 pism_config:stress_balance.sia.bed_smoother.range_option = "bed_smoother_range"; 2068 pism_config:stress_balance.sia.bed_smoother.range_type = "number"; 2069 pism_config:stress_balance.sia.bed_smoother.range_units = "meters"; 2070 2071 pism_config:stress_balance.sia.bed_smoother.theta_min = 0.0; 2072 pism_config:stress_balance.sia.bed_smoother.theta_min_doc = "minimum value of `\\theta` in :cite:`Schoofbasaltopg2003` bed roughness parameterization for SIA"; 2073 pism_config:stress_balance.sia.bed_smoother.theta_min_type = "number"; 2074 pism_config:stress_balance.sia.bed_smoother.theta_min_units = "1"; 2075 2076 pism_config:stress_balance.sia.e_age_coupling = "no"; 2077 pism_config:stress_balance.sia.e_age_coupling_doc = "Couple the SIA enhancement factor to age as in :cite:`Greve`."; 2078 pism_config:stress_balance.sia.e_age_coupling_option = "e_age_coupling"; 2079 pism_config:stress_balance.sia.e_age_coupling_type = "flag"; 2080 2081 pism_config:stress_balance.sia.enhancement_factor = 1.0; 2082 pism_config:stress_balance.sia.enhancement_factor_doc = "Flow enhancement factor for SIA"; 2083 pism_config:stress_balance.sia.enhancement_factor_option = "sia_e"; 2084 pism_config:stress_balance.sia.enhancement_factor_type = "number"; 2085 pism_config:stress_balance.sia.enhancement_factor_units = "1"; 2086 2087 pism_config:stress_balance.sia.enhancement_factor_interglacial = 1.0; 2088 pism_config:stress_balance.sia.enhancement_factor_interglacial_doc = "Flow enhancement factor for SIA; used for ice accumulated during interglacial periods."; 2089 pism_config:stress_balance.sia.enhancement_factor_interglacial_option = "sia_e_interglacial"; 2090 pism_config:stress_balance.sia.enhancement_factor_interglacial_type = "number"; 2091 pism_config:stress_balance.sia.enhancement_factor_interglacial_units = "1"; 2092 2093 pism_config:stress_balance.sia.flow_law = "gpbld"; 2094 pism_config:stress_balance.sia.flow_law_choices = "arr,arrwarm,gk,gpbld,hooke,isothermal_glen,pb"; 2095 pism_config:stress_balance.sia.flow_law_doc = "The SIA flow law."; 2096 pism_config:stress_balance.sia.flow_law_option = "sia_flow_law"; 2097 pism_config:stress_balance.sia.flow_law_type = "keyword"; 2098 2099 pism_config:stress_balance.sia.grain_size_age_coupling = "no"; 2100 pism_config:stress_balance.sia.grain_size_age_coupling_doc = "Use age of the ice to compute grain size to use with the Goldsby-Kohlstedt :cite:`GoldsbyKohlstedt` flow law"; 2101 pism_config:stress_balance.sia.grain_size_age_coupling_option = "grain_size_age_coupling"; 2102 pism_config:stress_balance.sia.grain_size_age_coupling_type = "flag"; 2103 2104 pism_config:stress_balance.sia.limit_diffusivity = "no"; 2105 pism_config:stress_balance.sia.limit_diffusivity_doc = "Limit SIA diffusivity by `stress_balance.sia.max_diffusivity`."; 2106 pism_config:stress_balance.sia.limit_diffusivity_option = "limit_sia_diffusivity"; 2107 pism_config:stress_balance.sia.limit_diffusivity_type = "flag"; 2108 2109 pism_config:stress_balance.sia.max_diffusivity = 100.0; 2110 pism_config:stress_balance.sia.max_diffusivity_doc = "Maximum allowed diffusivity of the SIA flow. PISM stops with an error message if the SIA diffusivity exceeds this limit."; 2111 pism_config:stress_balance.sia.max_diffusivity_type = "number"; 2112 pism_config:stress_balance.sia.max_diffusivity_units = "m2 s-1"; 2113 2114 pism_config:stress_balance.sia.surface_gradient_method = "haseloff"; 2115 pism_config:stress_balance.sia.surface_gradient_method_choices = "eta,haseloff,mahaffy"; 2116 pism_config:stress_balance.sia.surface_gradient_method_doc = "method used for surface gradient calculation at staggered grid points"; 2117 pism_config:stress_balance.sia.surface_gradient_method_option = "gradient"; 2118 pism_config:stress_balance.sia.surface_gradient_method_type = "keyword"; 2119 2120 pism_config:stress_balance.ssa.Glen_exponent = 3.0; 2121 pism_config:stress_balance.ssa.Glen_exponent_doc = "Glen exponent in ice flow law for SSA"; 2122 pism_config:stress_balance.ssa.Glen_exponent_option = "ssa_n"; 2123 pism_config:stress_balance.ssa.Glen_exponent_type = "number"; 2124 pism_config:stress_balance.ssa.Glen_exponent_units = "pure number"; 2125 2126 pism_config:stress_balance.ssa.compute_surface_gradient_inward = "no"; 2127 pism_config:stress_balance.ssa.compute_surface_gradient_inward_doc = "If yes then use inward first-order differencing in computing surface gradient in the SSA objects."; 2128 pism_config:stress_balance.ssa.compute_surface_gradient_inward_type = "flag"; 2129 2130 pism_config:stress_balance.ssa.dirichlet_bc = "no"; 2131 pism_config:stress_balance.ssa.dirichlet_bc_doc = "apply SSA velocity Dirichlet boundary condition"; 2132 pism_config:stress_balance.ssa.dirichlet_bc_option = "ssa_dirichlet_bc"; 2133 pism_config:stress_balance.ssa.dirichlet_bc_type = "flag"; 2134 2135 pism_config:stress_balance.ssa.enhancement_factor = 1.0; 2136 pism_config:stress_balance.ssa.enhancement_factor_doc = "Flow enhancement factor for SSA"; 2137 pism_config:stress_balance.ssa.enhancement_factor_option = "ssa_e"; 2138 pism_config:stress_balance.ssa.enhancement_factor_type = "number"; 2139 pism_config:stress_balance.ssa.enhancement_factor_units = "1"; 2140 2141 pism_config:stress_balance.ssa.enhancement_factor_interglacial = 1.0; 2142 pism_config:stress_balance.ssa.enhancement_factor_interglacial_doc = "Flow enhancement factor for SSA; used for ice accumulated during interglacial periods."; 2143 pism_config:stress_balance.ssa.enhancement_factor_interglacial_option = "ssa_e_interglacial"; 2144 pism_config:stress_balance.ssa.enhancement_factor_interglacial_type = "number"; 2145 pism_config:stress_balance.ssa.enhancement_factor_interglacial_units = "1"; 2146 2147 pism_config:stress_balance.ssa.epsilon = 1.0e13; 2148 pism_config:stress_balance.ssa.epsilon_doc = "Initial amount of regularization in computation of product of effective viscosity and thickness (`\\nu H`). This default value for `\\nu H` comes e.g. from a hardness for the Ross ice shelf (`\\bar B`) = 1.9e8 Pa `s^{1/3}` :cite:`MacAyealetal` and a typical strain rate of 0.001 1/year for the Ross ice shelf, giving `\\nu = (\\bar B) / (2 \\cdot 0.001^{2/3})` = 9.49e+14 Pa s ~ 30 MPa year, the value in :cite:`Ritzetal2001`, but with a tiny thickness `H` of about 1 cm."; 2149 pism_config:stress_balance.ssa.epsilon_option = "ssa_eps"; 2150 pism_config:stress_balance.ssa.epsilon_type = "number"; 2151 pism_config:stress_balance.ssa.epsilon_units = "Pascal second meter"; 2152 2153 pism_config:stress_balance.ssa.fd.brutal_sliding = "false"; 2154 pism_config:stress_balance.ssa.fd.brutal_sliding_doc = "Enhance sliding speed brutally."; 2155 pism_config:stress_balance.ssa.fd.brutal_sliding_option = "brutal_sliding"; 2156 pism_config:stress_balance.ssa.fd.brutal_sliding_type = "flag"; 2157 2158 pism_config:stress_balance.ssa.fd.brutal_sliding_scale = 1.0; 2159 pism_config:stress_balance.ssa.fd.brutal_sliding_scale_doc = "Brutal SSA Sliding Scale"; 2160 pism_config:stress_balance.ssa.fd.brutal_sliding_scale_option = "brutal_sliding_scale"; 2161 pism_config:stress_balance.ssa.fd.brutal_sliding_scale_type = "number"; 2162 pism_config:stress_balance.ssa.fd.brutal_sliding_scale_units = "1"; 2163 2164 pism_config:stress_balance.ssa.fd.lateral_drag.enabled = "false"; 2165 pism_config:stress_balance.ssa.fd.lateral_drag.enabled_doc = "set viscosity at ice shelf margin next to ice free bedrock as friction parameterization"; 2166 pism_config:stress_balance.ssa.fd.lateral_drag.enabled_type = "flag"; 2167 2168 pism_config:stress_balance.ssa.fd.lateral_drag.viscosity = 5.0e15; 2169 pism_config:stress_balance.ssa.fd.lateral_drag.viscosity_doc = "Staggered Viscosity used as side friction parameterization."; 2170 pism_config:stress_balance.ssa.fd.lateral_drag.viscosity_option = "nu_bedrock"; 2171 pism_config:stress_balance.ssa.fd.lateral_drag.viscosity_type = "number"; 2172 pism_config:stress_balance.ssa.fd.lateral_drag.viscosity_units = "Pascal second"; 2173 2174 pism_config:stress_balance.ssa.fd.max_iterations = 300; 2175 pism_config:stress_balance.ssa.fd.max_iterations_doc = "Maximum number of Picard iterations for the ice viscosity computation, in the SSAFD object"; 2176 pism_config:stress_balance.ssa.fd.max_iterations_option = "ssafd_picard_maxi"; 2177 pism_config:stress_balance.ssa.fd.max_iterations_type = "integer"; 2178 pism_config:stress_balance.ssa.fd.max_iterations_units = "count"; 2179 2180 pism_config:stress_balance.ssa.fd.max_speed = 300000; 2181 pism_config:stress_balance.ssa.fd.max_speed_doc = "Upper bound for the ice speed computed by the SSAFD solver."; 2182 pism_config:stress_balance.ssa.fd.max_speed_option = "ssafd_max_speed"; 2183 pism_config:stress_balance.ssa.fd.max_speed_type = "number"; 2184 pism_config:stress_balance.ssa.fd.max_speed_units = "km s-1"; 2185 2186 pism_config:stress_balance.ssa.fd.nuH_iter_failure_underrelaxation = 0.8; 2187 pism_config:stress_balance.ssa.fd.nuH_iter_failure_underrelaxation_doc = "In event of 'Effective viscosity not converged' failure, use outer iteration rule nuH <- nuH + f (nuH - nuH_old), where f is this parameter."; 2188 pism_config:stress_balance.ssa.fd.nuH_iter_failure_underrelaxation_option = "ssafd_nuH_iter_failure_underrelaxation"; 2189 pism_config:stress_balance.ssa.fd.nuH_iter_failure_underrelaxation_type = "number"; 2190 pism_config:stress_balance.ssa.fd.nuH_iter_failure_underrelaxation_units = "pure number"; 2191 2192 pism_config:stress_balance.ssa.fd.relative_convergence = 1.0e-4; 2193 pism_config:stress_balance.ssa.fd.relative_convergence_doc = "Relative change tolerance for the effective viscosity in the SSAFD object"; 2194 pism_config:stress_balance.ssa.fd.relative_convergence_option = "ssafd_picard_rtol"; 2195 pism_config:stress_balance.ssa.fd.relative_convergence_type = "number"; 2196 pism_config:stress_balance.ssa.fd.relative_convergence_units = "1"; 2197 2198 pism_config:stress_balance.ssa.fd.replace_zero_diagonal_entries = "yes"; 2199 pism_config:stress_balance.ssa.fd.replace_zero_diagonal_entries_doc = "Replace zero diagonal entries in the SSAFD matrix with basal_resistance.beta_ice_free_bedrock to avoid solver failures."; 2200 pism_config:stress_balance.ssa.fd.replace_zero_diagonal_entries_type = "flag"; 2201 2202 pism_config:stress_balance.ssa.flow_law = "gpbld"; 2203 pism_config:stress_balance.ssa.flow_law_choices = "arr,arrwarm,gpbld,hooke,isothermal_glen,pb"; 2204 pism_config:stress_balance.ssa.flow_law_doc = "The SSA flow law."; 2205 pism_config:stress_balance.ssa.flow_law_option = "ssa_flow_law"; 2206 pism_config:stress_balance.ssa.flow_law_type = "keyword"; 2207 2208 pism_config:stress_balance.ssa.method = "fd"; 2209 pism_config:stress_balance.ssa.method_choices = "fd,fem"; 2210 pism_config:stress_balance.ssa.method_doc = "Algorithm for computing the SSA solution."; 2211 pism_config:stress_balance.ssa.method_option = "ssa_method"; 2212 pism_config:stress_balance.ssa.method_type = "keyword"; 2213 2214 pism_config:stress_balance.ssa.read_initial_guess = "yes"; 2215 pism_config:stress_balance.ssa.read_initial_guess_doc = "Read the initial guess from the input file when re-starting."; 2216 pism_config:stress_balance.ssa.read_initial_guess_option = "ssa_read_initial_guess"; 2217 pism_config:stress_balance.ssa.read_initial_guess_type = "flag"; 2218 2219 pism_config:stress_balance.ssa.strength_extension.constant_nu = 9.48680701906572e+14; 2220 pism_config:stress_balance.ssa.strength_extension.constant_nu_doc = "The SSA is made elliptic by use of a constant value for the product of viscosity (nu) and thickness (H). This value for nu comes from hardness (bar B)=1.9e8 `Pa s^{1/3}` :cite:`MacAyealetal` and a typical strain rate of 0.001 year-1: `\\nu = (\\bar B) / (2 \\cdot 0.001^{2/3})`. Compare the value of 9.45e14 Pa s = 30 MPa year in :cite:`Ritzetal2001`."; 2221 pism_config:stress_balance.ssa.strength_extension.constant_nu_type = "number"; 2222 pism_config:stress_balance.ssa.strength_extension.constant_nu_units = "Pascal second"; 2223 2224 pism_config:stress_balance.ssa.strength_extension.min_thickness = 50.0; 2225 pism_config:stress_balance.ssa.strength_extension.min_thickness_doc = "The SSA is made elliptic by use of a constant value for the product of viscosity (nu) and thickness (H). At ice thicknesses below this value the product nu*H switches from the normal vertical integral to a constant value. The geometry itself is not affected by this value."; 2226 pism_config:stress_balance.ssa.strength_extension.min_thickness_type = "number"; 2227 pism_config:stress_balance.ssa.strength_extension.min_thickness_units = "meters"; 2228 2229 pism_config:stress_balance.vertical_velocity_approximation = "centered"; 2230 pism_config:stress_balance.vertical_velocity_approximation_choices = "centered,upstream"; 2231 pism_config:stress_balance.vertical_velocity_approximation_doc = "Vertical velocity FD approximation. \"Upstream\" uses first-order finite difference to compute u_x and v_y. Uses basal velocity to make decisions."; 2232 pism_config:stress_balance.vertical_velocity_approximation_option = "vertical_velocity_approximation"; 2233 pism_config:stress_balance.vertical_velocity_approximation_type = "keyword"; 2234 2235 pism_config:stress_balance.weertman_sliding.A = 1.8e-16; 2236 pism_config:stress_balance.weertman_sliding.A_doc = "Sliding parameter in the Weertman-style sliding parameterization :cite:`Tomkin2007`"; 2237 pism_config:stress_balance.weertman_sliding.A_type = "number"; 2238 pism_config:stress_balance.weertman_sliding.A_units = "Pa-3 year-1 m-2"; 2239 2240 pism_config:stress_balance.weertman_sliding.k = 0.2; 2241 pism_config:stress_balance.weertman_sliding.k_doc = "The ratio of the basal water pressure and the ice overburden pressure in the Weertman-style sliding parameterization."; 2242 pism_config:stress_balance.weertman_sliding.k_type = "number"; 2243 pism_config:stress_balance.weertman_sliding.k_units = "1"; 2244 2245 pism_config:surface.anomaly.file = ""; 2246 pism_config:surface.anomaly.file_doc = "Name of the file containing climate forcing fields."; 2247 pism_config:surface.anomaly.file_option = "surface_anomaly_file"; 2248 pism_config:surface.anomaly.file_type = "string"; 2249 2250 pism_config:surface.anomaly.period = 0; 2251 pism_config:surface.anomaly.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2252 pism_config:surface.anomaly.period_option = "surface_anomaly_period"; 2253 pism_config:surface.anomaly.period_type = "integer"; 2254 pism_config:surface.anomaly.period_units = "years"; 2255 2256 pism_config:surface.anomaly.reference_year = 0; 2257 pism_config:surface.anomaly.reference_year_doc = "Reference year to use when ``surface.anomaly.period`` is active."; 2258 pism_config:surface.anomaly.reference_year_option = "surface_anomaly_reference_year"; 2259 pism_config:surface.anomaly.reference_year_type = "integer"; 2260 pism_config:surface.anomaly.reference_year_units = "years"; 2261 2262 pism_config:surface.cache.update_interval = 10; 2263 pism_config:surface.cache.update_interval_doc = "Update interval (in years) for the `-surface cache` modifier."; 2264 pism_config:surface.cache.update_interval_type = "integer"; 2265 pism_config:surface.cache.update_interval_units = "year"; 2266 2267 pism_config:surface.delta_T.file = ""; 2268 pism_config:surface.delta_T.file_doc = "Name of the file containing temperature offsets."; 2269 pism_config:surface.delta_T.file_option = "surface_delta_T_file"; 2270 pism_config:surface.delta_T.file_type = "string"; 2271 2272 pism_config:surface.delta_T.period = 0; 2273 pism_config:surface.delta_T.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2274 pism_config:surface.delta_T.period_option = "surface_delta_T_period"; 2275 pism_config:surface.delta_T.period_type = "integer"; 2276 pism_config:surface.delta_T.period_units = "years"; 2277 2278 pism_config:surface.delta_T.reference_year = 0; 2279 pism_config:surface.delta_T.reference_year_doc = "Reference year to use when ``surface.delta_T.period`` is active."; 2280 pism_config:surface.delta_T.reference_year_option = "surface_delta_T_reference_year"; 2281 pism_config:surface.delta_T.reference_year_type = "integer"; 2282 pism_config:surface.delta_T.reference_year_units = "years"; 2283 2284 pism_config:surface.elevation_change.file = ""; 2285 pism_config:surface.elevation_change.file_doc = "Name of the file containing the reference surface elevation field (variable ``usurf``)."; 2286 pism_config:surface.elevation_change.file_option = "surface_elevation_change_file"; 2287 pism_config:surface.elevation_change.file_type = "string"; 2288 2289 pism_config:surface.elevation_change.period = 0; 2290 pism_config:surface.elevation_change.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2291 pism_config:surface.elevation_change.period_option = "surface_elevation_change_period"; 2292 pism_config:surface.elevation_change.period_type = "integer"; 2293 pism_config:surface.elevation_change.period_units = "years"; 2294 2295 pism_config:surface.elevation_change.reference_year = 0; 2296 pism_config:surface.elevation_change.reference_year_doc = "Reference year to use when ``surface.elevation_change.period`` is active."; 2297 pism_config:surface.elevation_change.reference_year_option = "surface_elevation_change_reference_year"; 2298 pism_config:surface.elevation_change.reference_year_type = "integer"; 2299 pism_config:surface.elevation_change.reference_year_units = "years"; 2300 2301 pism_config:surface.elevation_change.smb.exp_factor = 0; 2302 pism_config:surface.elevation_change.smb.exp_factor_doc = "Exponential for the surface mass balance."; 2303 pism_config:surface.elevation_change.smb.exp_factor_option = "smb_exp_factor"; 2304 pism_config:surface.elevation_change.smb.exp_factor_type = "number"; 2305 pism_config:surface.elevation_change.smb.exp_factor_units = "Kelvin-1"; 2306 2307 pism_config:surface.elevation_change.smb.lapse_rate = 0; 2308 pism_config:surface.elevation_change.smb.lapse_rate_doc = "Lapse rate for the surface mass balance."; 2309 pism_config:surface.elevation_change.smb.lapse_rate_option = "smb_lapse_rate"; 2310 pism_config:surface.elevation_change.smb.lapse_rate_type = "number"; 2311 pism_config:surface.elevation_change.smb.lapse_rate_units = "(m / year) / km"; 2312 2313 pism_config:surface.elevation_change.smb.method = "shift"; 2314 pism_config:surface.elevation_change.smb.method_choices = "scale,shift"; 2315 pism_config:surface.elevation_change.smb.method_doc = "Choose the SMB adjustment method. ``scale``: use temperature-change-dependent scaling factor. ``shift``: use the SMB lapse rate."; 2316 pism_config:surface.elevation_change.smb.method_option = "smb_adjustment"; 2317 pism_config:surface.elevation_change.smb.method_type = "keyword"; 2318 2319 pism_config:surface.elevation_change.temperature_lapse_rate = 0; 2320 pism_config:surface.elevation_change.temperature_lapse_rate_doc = "Lapse rate for the temperature at the top of the ice."; 2321 pism_config:surface.elevation_change.temperature_lapse_rate_option = "temp_lapse_rate"; 2322 pism_config:surface.elevation_change.temperature_lapse_rate_type = "number"; 2323 pism_config:surface.elevation_change.temperature_lapse_rate_units = "K / km"; 2324 2325 pism_config:surface.force_to_thickness.alpha = 0.01; 2326 pism_config:surface.force_to_thickness.alpha_doc = "exponential coefficient in force-to-thickness mechanism"; 2327 pism_config:surface.force_to_thickness.alpha_option = "force_to_thickness_alpha"; 2328 pism_config:surface.force_to_thickness.alpha_type = "number"; 2329 pism_config:surface.force_to_thickness.alpha_units = "year-1"; 2330 2331 pism_config:surface.force_to_thickness.ice_free_alpha_factor = 1.0; 2332 pism_config:surface.force_to_thickness.ice_free_alpha_factor_doc = "surface.force_to_thickness.alpha is multiplied by this factor in areas that are ice-free according to the target ice thickness and surface.force_to_thickness.ice_free_thickness_threshold"; 2333 pism_config:surface.force_to_thickness.ice_free_alpha_factor_option = "force_to_thickness_ice_free_alpha_factor"; 2334 pism_config:surface.force_to_thickness.ice_free_alpha_factor_type = "number"; 2335 pism_config:surface.force_to_thickness.ice_free_alpha_factor_units = "1"; 2336 2337 pism_config:surface.force_to_thickness.ice_free_thickness_threshold = 1.0; 2338 pism_config:surface.force_to_thickness.ice_free_thickness_threshold_doc = "threshold of ice thickness in the force-to-thickness target field. Used to determine whether to use surface.force_to_thickness.ice_free_alpha_factor."; 2339 pism_config:surface.force_to_thickness.ice_free_thickness_threshold_option = "force_to_thickness_ice_free_thickness_threshold"; 2340 pism_config:surface.force_to_thickness.ice_free_thickness_threshold_type = "number"; 2341 pism_config:surface.force_to_thickness.ice_free_thickness_threshold_units = "meters"; 2342 2343 pism_config:surface.force_to_thickness.start_time = -4.54e9; 2344 pism_config:surface.force_to_thickness.start_time_doc = "Starting time for the \"force to thickness\" modifier; the default is \"start from the creation of the Earth.\""; 2345 pism_config:surface.force_to_thickness.start_time_type = "number"; 2346 pism_config:surface.force_to_thickness.start_time_units = "years"; 2347 2348 pism_config:surface.force_to_thickness_file = ""; 2349 pism_config:surface.force_to_thickness_file_doc = "The name of the file to read the target ice thickness from."; 2350 pism_config:surface.force_to_thickness_file_option = "force_to_thickness_file"; 2351 pism_config:surface.force_to_thickness_file_type = "string"; 2352 2353 pism_config:surface.given.file = ""; 2354 pism_config:surface.given.file_doc = "Name of the file containing climate forcing fields."; 2355 pism_config:surface.given.file_option = "surface_given_file"; 2356 pism_config:surface.given.file_type = "string"; 2357 2358 pism_config:surface.given.period = 0; 2359 pism_config:surface.given.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2360 pism_config:surface.given.period_option = "surface_given_period"; 2361 pism_config:surface.given.period_type = "integer"; 2362 pism_config:surface.given.period_units = "years"; 2363 2364 pism_config:surface.given.reference_year = 0; 2365 pism_config:surface.given.reference_year_doc = "Reference year to use when ``surface.given.period`` is active."; 2366 pism_config:surface.given.reference_year_option = "surface_given_reference_year"; 2367 pism_config:surface.given.reference_year_type = "integer"; 2368 pism_config:surface.given.reference_year_units = "years"; 2369 2370 pism_config:surface.given.smb_max = 91000; 2371 pism_config:surface.given.smb_max_doc = "Maximum climatic mass balance value (used to check input data). Corresponds to 100 m/year ice equivalent."; 2372 pism_config:surface.given.smb_max_type = "number"; 2373 pism_config:surface.given.smb_max_units = "kg m-2 year-1"; 2374 2375 pism_config:surface.ismip6.file = ""; 2376 pism_config:surface.ismip6.file_doc = "Name of the file containing climate forcing anomaly fields."; 2377 pism_config:surface.ismip6.file_option = "surface_ismip6_file"; 2378 pism_config:surface.ismip6.file_type = "string"; 2379 2380 pism_config:surface.ismip6.period = 0; 2381 pism_config:surface.ismip6.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2382 pism_config:surface.ismip6.period_option = "surface_ismip6_period"; 2383 pism_config:surface.ismip6.period_type = "integer"; 2384 pism_config:surface.ismip6.period_units = "years"; 2385 2386 pism_config:surface.ismip6.reference_file = ""; 2387 pism_config:surface.ismip6.reference_file_doc = "Name of the file containing reference climate forcing fields."; 2388 pism_config:surface.ismip6.reference_file_option = "surface_ismip6_reference_file"; 2389 pism_config:surface.ismip6.reference_file_type = "string"; 2390 2391 pism_config:surface.ismip6.reference_year = 0; 2392 pism_config:surface.ismip6.reference_year_doc = "Reference year to use when ``surface.ismip6.period`` is active."; 2393 pism_config:surface.ismip6.reference_year_type = "integer"; 2394 pism_config:surface.ismip6.reference_year_units = "years"; 2395 2396 pism_config:surface.models = "given"; 2397 pism_config:surface.models_doc = "Comma-separated list of surface models and modifiers."; 2398 pism_config:surface.models_option = "surface"; 2399 pism_config:surface.models_type = "string"; 2400 2401 pism_config:surface.pdd.air_temp_all_precip_as_rain = 275.15; 2402 pism_config:surface.pdd.air_temp_all_precip_as_rain_doc = "threshold temperature above which all precipitation is rain; must exceed surface.pdd.air_temp_all_precip_as_snow to avoid division by zero, because difference is in a denominator"; 2403 pism_config:surface.pdd.air_temp_all_precip_as_rain_type = "number"; 2404 pism_config:surface.pdd.air_temp_all_precip_as_rain_units = "Kelvin"; 2405 2406 pism_config:surface.pdd.air_temp_all_precip_as_snow = 273.15; 2407 pism_config:surface.pdd.air_temp_all_precip_as_snow_doc = "threshold temperature below which all precipitation is snow"; 2408 pism_config:surface.pdd.air_temp_all_precip_as_snow_type = "number"; 2409 pism_config:surface.pdd.air_temp_all_precip_as_snow_units = "Kelvin"; 2410 2411 pism_config:surface.pdd.balance_year_start_day = 274; 2412 pism_config:surface.pdd.balance_year_start_day_doc = "day of year for October 1st, beginning of the balance year in northern latitudes."; 2413 pism_config:surface.pdd.balance_year_start_day_type = "integer"; 2414 pism_config:surface.pdd.balance_year_start_day_units = "ordinal day number"; 2415 2416 pism_config:surface.pdd.factor_ice = 0.00879120879120879; 2417 pism_config:surface.pdd.factor_ice_doc = "EISMINT-Greenland value :cite:`RitzEISMINT`; = (8 mm liquid-water-equivalent) / (pos degree day)"; 2418 pism_config:surface.pdd.factor_ice_type = "number"; 2419 pism_config:surface.pdd.factor_ice_units = "meter / (Kelvin day)"; 2420 2421 pism_config:surface.pdd.factor_snow = 0.0032967032967033; 2422 pism_config:surface.pdd.factor_snow_doc = "EISMINT-Greenland value :cite:`RitzEISMINT`; = (3 mm liquid-water-equivalent) / (pos degree day)"; 2423 pism_config:surface.pdd.factor_snow_type = "number"; 2424 pism_config:surface.pdd.factor_snow_units = "meter / (Kelvin day)"; 2425 2426 pism_config:surface.pdd.fausto.T_c = 272.15; 2427 pism_config:surface.pdd.fausto.T_c_doc = "= -1 + 273.15; for formula (6) in :cite:`Faustoetal2009`"; 2428 pism_config:surface.pdd.fausto.T_c_type = "number"; 2429 pism_config:surface.pdd.fausto.T_c_units = "Kelvin"; 2430 2431 pism_config:surface.pdd.fausto.T_w = 283.15; 2432 pism_config:surface.pdd.fausto.T_w_doc = "= 10 + 273.15; for formula (6) in :cite:`Faustoetal2009`"; 2433 pism_config:surface.pdd.fausto.T_w_type = "number"; 2434 pism_config:surface.pdd.fausto.T_w_units = "Kelvin"; 2435 2436 pism_config:surface.pdd.fausto.beta_ice_c = 0.015; 2437 pism_config:surface.pdd.fausto.beta_ice_c_doc = "water-equivalent thickness; for formula (6) in :cite:`Faustoetal2009`"; 2438 pism_config:surface.pdd.fausto.beta_ice_c_type = "number"; 2439 pism_config:surface.pdd.fausto.beta_ice_c_units = "meter / (Kelvin day)"; 2440 2441 pism_config:surface.pdd.fausto.beta_ice_w = 0.007; 2442 pism_config:surface.pdd.fausto.beta_ice_w_doc = "water-equivalent thickness; for formula (6) in :cite:`Faustoetal2009`"; 2443 pism_config:surface.pdd.fausto.beta_ice_w_type = "number"; 2444 pism_config:surface.pdd.fausto.beta_ice_w_units = "meter / (Kelvin day)"; 2445 2446 pism_config:surface.pdd.fausto.beta_snow_c = 0.003; 2447 pism_config:surface.pdd.fausto.beta_snow_c_doc = "water-equivalent thickness; for formula (6) in :cite:`Faustoetal2009`"; 2448 pism_config:surface.pdd.fausto.beta_snow_c_type = "number"; 2449 pism_config:surface.pdd.fausto.beta_snow_c_units = "meter / (Kelvin day)"; 2450 2451 pism_config:surface.pdd.fausto.beta_snow_w = 0.003; 2452 pism_config:surface.pdd.fausto.beta_snow_w_doc = "water-equivalent thickness; for formula (6) in :cite:`Faustoetal2009`"; 2453 pism_config:surface.pdd.fausto.beta_snow_w_type = "number"; 2454 pism_config:surface.pdd.fausto.beta_snow_w_units = "meter / (Kelvin day)"; 2455 2456 pism_config:surface.pdd.fausto.enabled = "false"; 2457 pism_config:surface.pdd.fausto.enabled_doc = "Set PDD parameters using formulas (6) and (7) in :cite:`Faustoetal2009`"; 2458 pism_config:surface.pdd.fausto.enabled_option = "pdd_fausto"; 2459 pism_config:surface.pdd.fausto.enabled_type = "flag"; 2460 2461 pism_config:surface.pdd.fausto.latitude_beta_w = 72.0; 2462 pism_config:surface.pdd.fausto.latitude_beta_w_doc = "latitude below which to use warm case, in formula (6) in :cite:`Faustoetal2009`"; 2463 pism_config:surface.pdd.fausto.latitude_beta_w_type = "number"; 2464 pism_config:surface.pdd.fausto.latitude_beta_w_units = "degree_north"; 2465 2466 pism_config:surface.pdd.firn_compaction_to_accumulation_ratio = 0.75; 2467 pism_config:surface.pdd.firn_compaction_to_accumulation_ratio_doc = "How much firn as a fraction of accumulation is turned into ice"; 2468 pism_config:surface.pdd.firn_compaction_to_accumulation_ratio_type = "number"; 2469 pism_config:surface.pdd.firn_compaction_to_accumulation_ratio_units = "1"; 2470 2471 pism_config:surface.pdd.firn_depth_file = ""; 2472 pism_config:surface.pdd.firn_depth_file_doc = "The name of the file to read the firn_depth from."; 2473 pism_config:surface.pdd.firn_depth_file_option = "pdd_firn_depth_file"; 2474 pism_config:surface.pdd.firn_depth_file_type = "string"; 2475 2476 pism_config:surface.pdd.interpret_precip_as_snow = "no"; 2477 pism_config:surface.pdd.interpret_precip_as_snow_doc = "Interpret precipitation as snow fall."; 2478 pism_config:surface.pdd.interpret_precip_as_snow_type = "flag"; 2479 2480 pism_config:surface.pdd.max_evals_per_year = 52; 2481 pism_config:surface.pdd.max_evals_per_year_doc = "maximum number of times the PDD scheme will ask for air temperature and precipitation to build location-dependent time series for computing (expected) number of positive degree days and snow accumulation; the default means the PDD uses weekly samples of the annual cycle; see also surface.pdd.std_dev"; 2482 pism_config:surface.pdd.max_evals_per_year_type = "integer"; 2483 pism_config:surface.pdd.max_evals_per_year_units = "count"; 2484 2485 pism_config:surface.pdd.method = "expectation_integral"; 2486 pism_config:surface.pdd.method_choices = "expectation_integral,repeatable_random_process,random_process"; 2487 pism_config:surface.pdd.method_doc = "PDD implementation method"; 2488 pism_config:surface.pdd.method_option = "pdd_method"; 2489 pism_config:surface.pdd.method_type = "keyword"; 2490 2491 pism_config:surface.pdd.positive_threshold_temp = 273.15; 2492 pism_config:surface.pdd.positive_threshold_temp_doc = "temperature used to determine meaning of 'positive' degree day"; 2493 pism_config:surface.pdd.positive_threshold_temp_type = "number"; 2494 pism_config:surface.pdd.positive_threshold_temp_units = "Kelvin"; 2495 2496 pism_config:surface.pdd.refreeze = 0.6; 2497 pism_config:surface.pdd.refreeze_doc = "EISMINT-Greenland value :cite:`RitzEISMINT`"; 2498 pism_config:surface.pdd.refreeze_type = "number"; 2499 pism_config:surface.pdd.refreeze_units = "1"; 2500 2501 pism_config:surface.pdd.refreeze_ice_melt = "yes"; 2502 pism_config:surface.pdd.refreeze_ice_melt_doc = "If set to 'yes', refreeze surface.pdd.refreeze fraction of melted ice, otherwise all of the melted ice runs off."; 2503 pism_config:surface.pdd.refreeze_ice_melt_type = "flag"; 2504 2505 pism_config:surface.pdd.std_dev = 5.0; 2506 pism_config:surface.pdd.std_dev_doc = "std dev of daily temp variation; = EISMINT-Greenland value :cite:`RitzEISMINT`"; 2507 pism_config:surface.pdd.std_dev_type = "number"; 2508 pism_config:surface.pdd.std_dev_units = "Kelvin"; 2509 2510 pism_config:surface.pdd.std_dev.file = ""; 2511 pism_config:surface.pdd.std_dev.file_doc = "The name of the file to read `air_temp_sd` (standard deviation of air temperature) from."; 2512 pism_config:surface.pdd.std_dev.file_option = "pdd_sd_file"; 2513 pism_config:surface.pdd.std_dev.file_type = "string"; 2514 2515 pism_config:surface.pdd.std_dev.period = 0; 2516 pism_config:surface.pdd.std_dev.period_doc = "Length of the period of the climate forcing data. Set to zero to disable."; 2517 pism_config:surface.pdd.std_dev.period_option = "pdd_sd_period"; 2518 pism_config:surface.pdd.std_dev.period_type = "integer"; 2519 pism_config:surface.pdd.std_dev.period_units = "years"; 2520 2521 pism_config:surface.pdd.std_dev.reference_year = 0; 2522 pism_config:surface.pdd.std_dev.reference_year_doc = "Reference year to use when ``surface.pdd.std_dev.period`` is active."; 2523 pism_config:surface.pdd.std_dev.reference_year_option = "pdd_sd_reference_year"; 2524 pism_config:surface.pdd.std_dev.reference_year_type = "integer"; 2525 pism_config:surface.pdd.std_dev.reference_year_units = "years"; 2526 2527 pism_config:surface.pdd.std_dev_lapse_lat_base = 72.0; 2528 pism_config:surface.pdd.std_dev_lapse_lat_base_doc = "std_dev is a function of latitude, with value surface.pdd.std_dev at this latitude; this value only active if surface.pdd.std_dev_lapse_lat_rate is nonzero"; 2529 pism_config:surface.pdd.std_dev_lapse_lat_base_type = "number"; 2530 pism_config:surface.pdd.std_dev_lapse_lat_base_units = "degree_north"; 2531 2532 pism_config:surface.pdd.std_dev_lapse_lat_rate = 0.0; 2533 pism_config:surface.pdd.std_dev_lapse_lat_rate_doc = "std_dev is a function of latitude, with rate of change with respect to latitude given by this constant"; 2534 pism_config:surface.pdd.std_dev_lapse_lat_rate_type = "number"; 2535 pism_config:surface.pdd.std_dev_lapse_lat_rate_units = "Kelvin / degree_north"; 2536 2537 pism_config:surface.pdd.std_dev_param_a = -0.15; 2538 pism_config:surface.pdd.std_dev_param_a_doc = "Parameter a in Sigma = a*T + b, with T in degrees C. Used only if surface.pdd.std_dev_use_param is set to yes."; 2539 pism_config:surface.pdd.std_dev_param_a_type = "number"; 2540 pism_config:surface.pdd.std_dev_param_a_units = "pure number"; 2541 2542 pism_config:surface.pdd.std_dev_param_b = 0.66; 2543 pism_config:surface.pdd.std_dev_param_b_doc = "Parameter b in Sigma = a*T + b, with T in degrees C. Used only if surface.pdd.std_dev_use_param is set to yes."; 2544 pism_config:surface.pdd.std_dev_param_b_type = "number"; 2545 pism_config:surface.pdd.std_dev_param_b_units = "Kelvin"; 2546 2547 pism_config:surface.pdd.std_dev_use_param = "no"; 2548 pism_config:surface.pdd.std_dev_use_param_doc = "Parameterize standard deviation as a linear function of air temperature over ice-covered grid cells. The region of application is controlled by geometry.ice_free_thickness_standard."; 2549 pism_config:surface.pdd.std_dev_use_param_type = "flag"; 2550 2551 pism_config:surface.pressure = 0.0; 2552 pism_config:surface.pressure_doc = "atmospheric pressure; = pressure at ice surface"; 2553 pism_config:surface.pressure_type = "number"; 2554 pism_config:surface.pressure_units = "Pascal"; 2555 2556 pism_config:time.calendar = "365_day"; 2557 pism_config:time.calendar_choices = "standard,gregorian,proleptic_gregorian,noleap,365_day,360_day,julian,none"; 2558 pism_config:time.calendar_doc = "The calendar to use."; 2559 pism_config:time.calendar_option = "calendar"; 2560 pism_config:time.calendar_type = "keyword"; 2561 2562 pism_config:time.dimension_name = "time"; 2563 pism_config:time.dimension_name_doc = "The name of the time dimension in PISM output files."; 2564 pism_config:time.dimension_name_type = "string"; 2565 2566 pism_config:time.eemian_end = -114500.0; 2567 pism_config:time.eemian_end_doc = "End of the Eemian interglacial period. See :cite:`Greve97Greenland`."; 2568 pism_config:time.eemian_end_type = "number"; 2569 pism_config:time.eemian_end_units = "years"; 2570 2571 pism_config:time.eemian_start = -132000.0; 2572 pism_config:time.eemian_start_doc = "Start of the Eemian interglacial period. See :cite:`Greve97Greenland`."; 2573 pism_config:time.eemian_start_type = "number"; 2574 pism_config:time.eemian_start_units = "years"; 2575 2576 pism_config:time.holocene_start = -11000.0; 2577 pism_config:time.holocene_start_doc = "Start of the Holocene interglacial period. See :cite:`Greve97Greenland`."; 2578 pism_config:time.holocene_start_type = "number"; 2579 pism_config:time.holocene_start_units = "years"; 2580 2581 pism_config:time.reference_date = "1-1-1"; 2582 pism_config:time.reference_date_doc = "year-month-day; reference date used for calendar computations and in PISM output files"; 2583 pism_config:time.reference_date_type = "string"; 2584 2585 pism_config:time.run_length = 1000; 2586 pism_config:time.run_length_doc = "Default run length"; 2587 pism_config:time.run_length_type = "number"; 2588 pism_config:time.run_length_units = "years"; 2589 2590 pism_config:time.start_year = 0; 2591 pism_config:time.start_year_doc = "Start year."; 2592 pism_config:time.start_year_type = "number"; 2593 pism_config:time.start_year_units = "years"; 2594 2595 pism_config:time_stepping.adaptive_ratio = 0.12; 2596 pism_config:time_stepping.adaptive_ratio_doc = "Adaptive time stepping ratio for the explicit scheme for the mass balance equation; :cite:`BBL`, inequality (25)"; 2597 pism_config:time_stepping.adaptive_ratio_option = "adapt_ratio"; 2598 pism_config:time_stepping.adaptive_ratio_type = "number"; 2599 pism_config:time_stepping.adaptive_ratio_units = "1"; 2600 2601 pism_config:time_stepping.count_steps = "no"; 2602 pism_config:time_stepping.count_steps_doc = "If yes, IceModel::run() will count the number of time steps it took. Sometimes useful for performance evaluation. Counts all steps, regardless of whether processes (mass continuity, energy, velocity, ...) occurred within the step."; 2603 pism_config:time_stepping.count_steps_option = "count_steps"; 2604 pism_config:time_stepping.count_steps_type = "flag"; 2605 2606 pism_config:time_stepping.hit_extra_times = "yes"; 2607 pism_config:time_stepping.hit_extra_times_doc = "Modify the time-stepping mechanism to hit times requested using -extra_times."; 2608 pism_config:time_stepping.hit_extra_times_option = "extra_force_output_times"; 2609 pism_config:time_stepping.hit_extra_times_type = "flag"; 2610 2611 pism_config:time_stepping.hit_multiples = 0.0; 2612 pism_config:time_stepping.hit_multiples_doc = "Hit every X years, where X is specified using this parameter. Use 0 to disable"; 2613 pism_config:time_stepping.hit_multiples_option = "timestep_hit_multiples"; 2614 pism_config:time_stepping.hit_multiples_type = "number"; 2615 pism_config:time_stepping.hit_multiples_units = "years"; 2616 2617 pism_config:time_stepping.hit_save_times = "no"; 2618 pism_config:time_stepping.hit_save_times_doc = "Modify the time-stepping mechanism to hit times requested using -save_times."; 2619 pism_config:time_stepping.hit_save_times_option = "save_force_output_times"; 2620 pism_config:time_stepping.hit_save_times_type = "flag"; 2621 2622 pism_config:time_stepping.hit_ts_times = "no"; 2623 pism_config:time_stepping.hit_ts_times_doc = "Modify the time-stepping mechanism to hit times requested using -ts_times."; 2624 pism_config:time_stepping.hit_ts_times_type = "flag"; 2625 2626 pism_config:time_stepping.maximum_time_step = 60.0; 2627 pism_config:time_stepping.maximum_time_step_doc = "Maximum allowed time step length"; 2628 pism_config:time_stepping.maximum_time_step_option = "max_dt"; 2629 pism_config:time_stepping.maximum_time_step_type = "number"; 2630 pism_config:time_stepping.maximum_time_step_units = "years"; 2631 2632 pism_config:time_stepping.skip.enabled = "no"; 2633 pism_config:time_stepping.skip.enabled_doc = "Use the temperature, age, and SSA stress balance computation skipping mechanism."; 2634 pism_config:time_stepping.skip.enabled_option = "skip"; 2635 pism_config:time_stepping.skip.enabled_type = "flag"; 2636 2637 pism_config:time_stepping.skip.max = 10; 2638 pism_config:time_stepping.skip.max_doc = "Number of mass-balance steps, including SIA diffusivity updates, to perform before a the temperature, age, and SSA stress balance computations are done"; 2639 pism_config:time_stepping.skip.max_option = "skip_max"; 2640 pism_config:time_stepping.skip.max_type = "integer"; 2641 pism_config:time_stepping.skip.max_units = "count"; 2642 2643 pism_config:long_name = "PISM configuration flags and parameters."; 2644 pism_config:long_name_doc = "The 'long_name' attribute is required by CF conventions. It is not used by PISM itself."; 2645 }