Provides a derived type which specifies the boundary conditions for the ice shelf model used by Dallaston et al. (2015). These are Dirichlet conditions at the lower bound of the first condition as well as, for thickness, the upper bound.
Constructs a boundary condition object for an ice shelf based on the conditions used in Dallaston et al. (2015).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=r8), | intent(in) | :: | thickness | The ice thickness at the inflowing ice shelf boundary |
||
| real(kind=r8), | intent(in) | :: | velocity | The longitudinal ice velocity at the inflowing ice shelf boundary |
||
| real(kind=r8), | intent(in) | :: | chi | The dimensionless ratio $\chi \equiv \frac{\rho_igh_0x_x}{2\eta_0u_0}$ |
A type with procedures for getting the boundary conditions of the ice shelf model used by Dallaston et al. (2015). These are Dirichlet conditions at the lower bound of the first condition as well as, for thickness, the upper bound.
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=r8), | private | :: | thickness | = | 1.0_r8 | The thickness of the glacier at the inflowing boundary |
|
| real(kind=r8), | private | :: | velocity | = | 1.0_r8 | The velocity of the glacier at the inflowing boundary |
|
| real(kind=r8), | private | :: | chi | = | 1.0_r8 | The dimensionless ratio $\chi \equiv \frac{\rho_igh_0x_x}{2\eta_0u_0}$ |
| private pure function constructor(thickness, velocity, chi) | Constructs a boundary condition object for an ice shelf based on the conditions used in Dallaston et al. (2015). |
| procedure, public :: thickness_upper_bound => bound_array | Returns a 1D array which should be passed as the
|
| procedure, public :: thickness_upper_type => bound_type | Returns an array indicating what type of boundary conditions apply for thickness at the upper boundary of each dimension. The types are specified using the parameters in [boundary_types_mod]. |
| procedure, public :: thickness_lower_bound => dallaston2015_lower_bound | Returns a 1D array which should be passed as the
|
| procedure, public :: velocity_lower_bound => dallaston2015_lower_bound | Returns a 1D array which should be passed as the
|
| procedure, public :: velocity_upper_bound => dallaston2015_upper_bound | Returns a 1D array which should be passed as the
|
| procedure, public :: thickness_lower_type => dallaston2015_lower_type | Returns an array indicating what type of boundary conditions apply for thickness at the lower boundary of each dimension. The types are specified using the parameters in [boundary_types_mod]. |
| procedure, public :: velocity_lower_type => dallaston2015_lower_type | Returns an array indicating what type of boundary conditions apply for velocity at the lower boundary of each dimension. The types are specified using the parameters in [boundary_types_mod]. |
| procedure, public :: velocity_upper_type => dallaston2015_upper_type | Returns an array indicating what type of boundary conditions apply for velocity at the upper boundary of each dimension. The types are specified using the parameters in [boundary_types_mod]. |
| procedure, public :: boundary_residuals => dallaston2015_residuals | Returns an array consisting of the difference between the required boundary values and those which actually exist. This can then be appended to a glacier's state vector. The order in which these are listed is as follows: lower thickness boundary, upper thickness boundary, lower velocity boundary, and upper velocity boundary. |
Constructs a boundary condition object for an ice shelf based on the conditions used in Dallaston et al. (2015).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=r8), | intent(in) | :: | thickness | The ice thickness at the inflowing ice shelf boundary |
||
| real(kind=r8), | intent(in) | :: | velocity | The longitudinal ice velocity at the inflowing ice shelf boundary |
||
| real(kind=r8), | intent(in) | :: | chi | The dimensionless ratio $\chi \equiv \frac{\rho_igh_0x_x}{2\eta_0u_0}$ |
Indicates that one layer of cells at the lower boundary in the first dimension should be omitted. This is appropriate for Dirichlet boundary conditions.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(dallaston2015_glacier_boundary), | intent(in) | :: | this |
Indicates that one layer of cells at the upper boundary in the first dimension should be omitted for thickness. This is appropriate for Dirichlet boundary conditions.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(dallaston2015_glacier_boundary), | intent(in) | :: | this |
Specifies that the lower boundary in the first dimension has Dirichlet boundary conditions.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(dallaston2015_glacier_boundary), | intent(in) | :: | this |
Specifies that the upper boundary in the first dimension has Neumann boundary conditions.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(dallaston2015_glacier_boundary), | intent(in) | :: | this |
Returns the difference between the glacier conditions of the plume and the Dirichlet conditions prescribed in the model of Dallaston et al. (2015)
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(dallaston2015_glacier_boundary), | intent(in) | :: | this | |||
| class(scalar_field), | intent(in) | :: | thickness | A field containing the thickness of the glacier |
||
| class(vector_field), | intent(in) | :: | velocity | A field containing the flow velocity of the glacier |
||
| class(scalar_field), | intent(in) | :: | viscosity | A field containing the viscosity of the ice in the glacier. |
||
| real(kind=r8), | intent(in) | :: | t | The time at which the boundary conditions are to be calculated. |
An array containing the difference between the required boundary values and those which are actually present. They are stored in the order: lower thickness boundary, upper thickness boundary, lower velocity boundary, and upper velocity boundary.