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 — Method

Params(; <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
sec_per_year: 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)
accumulation_rate: uniform accumulation_rate (m/yr)
glen_a_activation_energy: activation energy in glen b calculation
glen_a_ref: array of glen a reference values used in glen b calculation
glen_temperature_ref: reference temperature using in glen b calculation
glen_n: exponent in glen b calculation
glen_reg_strain_rate: strain rate regularization value
weertman_c: basal sliding field of coefficients
weertman_m: sliding law exponent
weertman_reg_speed: regularization speed, used to prevent bed speed going to zero
sea_level_wrt_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

source