Parameters
A Params
object is a WAVI.jl object that stores parameters related to physical parameters involved in the problem. The following table summarises these parameters:
Keyword Argument | Description | Default Value |
---|---|---|
g | Gravitional acceleration | 9.81 m/s^2 |
density_ice | Ice density | 918.0 kg/m^3 |
density_ocean | Ocean density | 1028.0 kg/m^3 |
gas_const | Gas constant used in Arrhenius relation in determining Glen flow parameter | 8.314 |
sec_per_year | Seconds per year | 3.16e7 s |
default_thickness | Default ice thickness choice, value set everywhere if ice thickness is not set | 100 m |
default_viscosity | Default ice viscosity choice, value set everywhere if ice thickness is not set | 1e7 Pa a^{1/3} |
default_temperature | Default ice temperature choice, value set everywhere if ice thickness is not set | 265.7 K |
default_damage | Default dimensionless ice damage choice, value set everywhere if ice thickness is not set | 0 |
accumulation_rate | Accumulation rate. Can be a scalar (same value everywhere), or an array of the same size as the numerical grid. | 0.3 m/a |
glen_a_activation_energy | Activation energy used in Glen flow law | 5.8631e4 J |
glen_a_ref | Reference value used in Glen flow law | 1.5453e-17 |
glen_temperature_ref | Reference temperature used in Glen flow law | 263.15 K |
glen_n | Exponent in Glen flow law | 3 |
glen_reg_strain_rate | Strain rate regularization value, sets a lower bound on the strain rate | 1e-5 |
weertman_c | Sliding coefficient in Weertman sliding law. Can be a scalar (same value everywhere), or an array of the same size as the numerical grid. | 1e4 Pa a^{1/3}/m^{1/3} |
weertman_m | Exponent in Weertman sliding law | 3 |
weertman_reg_speed | Minimum basal sliding speed | 1e-5 m/a |
sea_level_wrt_geoid | Reference sea level | 0 m |
minimum_thickness | Minimum ice thickness (on grid points in ice domain) | 50 m |
WAVI.Params
— MethodParams(; <kwargs>)
Construct a WAVI.jl parameters object for holding physical parameters.
Keyword arguments
- 'dt': model timestep (NB: simulation timestep set in timestepping params, this value is updated when model embedded to the value specified in timestepping_params when passed to simulation)
- 'g': gravitational acceleration (m^2 / s)
- 'density_ice': ice density (kg / m^3)
- 'density_ocean' ocean water density (kg / m^3)
- 'gas_const': gas constant in glen b calculation
- 'secperyear': seconds per year (s)
- 'default_thickness': thickness value reverted to if no initial thickness passed (m)
- 'default_viscosity': viscosity value reverted to if no initial thickness passed (Pa s)
- 'default_temperature': temperature value reverted to if no initial thickness passed (K)
- 'default_damage': damage value reverted to if no initial thickness passed (dimensionless)
- 'accumulationrate': uniform accumulationrate (m/yr)
- 'glenaactivation_energy': activation energy in glen b calculation
- 'glenaref': array of glen a reference values used in glen b calculation
- 'glentemperatureref': reference temperature using in glen b calculation
- 'glen_n': exponent in glen b calculation
- 'glenregstrain_rate': strain rate regularization value
- 'weertman_c': basal sliding field of coefficients
- 'weertman_m': sliding law exponent
- 'weertmanregspeed': regularization speed, used to prevent bed speed going to zero
- 'sealevelwrt_geoid': reference sea level
- 'minimum_thickness': minimum ice thickness on model domain
- 'evolveShelves': flag for turning on and off the evolution of the shelves in the forward run_simulation
- 'smallHAF': small value of HAF used within update_thickness when not evolving shelves