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 }