mom_ice_shelf_dynamics::ice_shelf_dyn_cs integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_u_shelf id_u_shelf = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_v_shelf id_v_shelf = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_t_shelf id_t_shelf = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_ground_frac id_ground_frac = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_col_thick id_col_thick = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_od_av id_od_av = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_u_mask id_u_mask = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_v_mask id_v_mask = -1 Diagnostic handles. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::id_t_mask id_t_mask = -1 Diagnostic handles. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::u_shelf u_shelf => NULL() the zonal velocity of the ice shelf/sheet on q-points (B grid) [L T-1 ~> m s-1] real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::v_shelf v_shelf => NULL() the meridional velocity of the ice shelf/sheet on q-points (B grid) [L T-1 ~> m s-1] real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::u_face_mask u_face_mask => NULL() mask for velocity boundary conditions on the C-grid u-face - this is because the FEM cares about FACES THAT GET INTEGRATED OVER, not vertices. Will represent boundary conditions on computational boundary (or permanent boundary between fast-moving and near-stagnant ice FOR NOW: 1=interior bdry, 0=no-flow boundary, 2=stress bdry condition, 3=inhomogeneous Dirichlet boundary, 4=flux boundary: at these faces a flux will be specified which will override velocities; a homogeneous velocity condition will be specified (this seems to give the solver less difficulty) real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::v_face_mask v_face_mask => NULL() A mask for velocity boundary conditions on the C-grid v-face, with valued defined similarly to u_face_mask. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::u_face_mask_bdry u_face_mask_bdry => NULL() A duplicate copy of u_face_mask? real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::v_face_mask_bdry v_face_mask_bdry => NULL() A duplicate copy of v_face_mask? real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::u_flux_bdry_val u_flux_bdry_val => NULL() The ice volume flux per unit face length into the cell through open boundary u-faces (where u_face_mask=4) [Z L T-1 ~> m2 s-1]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::v_flux_bdry_val v_flux_bdry_val => NULL() The ice volume flux per unit face length into the cell through open boundary v-faces (where v_face_mask=4) [Z L T-1 ~> m2 s-1]?? real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::umask umask => NULL() u-mask on the actual degrees of freedom (B grid) 1=normal node, 3=inhomogeneous boundary node, 0 - no flow node (will also get ice-free nodes) real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::vmask vmask => NULL() v-mask on the actual degrees of freedom (B grid) 1=normal node, 3=inhomogeneous boundary node, 0 - no flow node (will also get ice-free nodes) real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::calve_mask calve_mask => NULL() a mask to prevent the ice shelf front from advancing past its initial position (but it may retreat) real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::t_shelf t_shelf => NULL() Vertically integrated temperature in the ice shelf/stream, on corner-points (B grid) [degC]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::tmask tmask => NULL() A mask on tracer points that is 1 where there is ice. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::ice_visc ice_visc => NULL() Glen's law ice viscosity, often in [R L4 Z T-1 ~> kg m2 s-1]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::thickness_bdry_val thickness_bdry_val => NULL() The ice thickness at an inflowing boundary [Z ~> m]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::u_bdry_val u_bdry_val => NULL() The zonal ice velocity at inflowing boundaries [L yr-1 ~> m yr-1]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::v_bdry_val v_bdry_val => NULL() The meridional ice velocity at inflowing boundaries [L yr-1 ~> m yr-1]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::h_bdry_val h_bdry_val => NULL() The ice thickness at inflowing boundaries [m]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::t_bdry_val t_bdry_val => NULL() The ice temperature at inflowing boundaries [degC]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::basal_traction basal_traction => NULL() The area integrated nonlinear part of "linearized" basal stress [R Z L2 T-1 ~> kg s-1]. The exact form depends on basal law exponent and/or whether flow is "hybridized" a la Goldberg 2011. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::od_rt od_rt => NULL() A running total for calculating OD_av. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::ground_frac_rt ground_frac_rt => NULL() A running total for calculating ground_frac. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::od_av od_av => NULL() The time average open ocean depth [Z ~> m]. real, dimension(:,:), pointer real, dimension(:,:), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::ground_frac ground_frac => NULL() Fraction of the time a cell is "exposed", i.e. the column thickness is below a threshold and interacting with the rock [nondim]. When this is 1, the ice-shelf is grounded. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::od_rt_counter od_rt_counter = 0 A counter of the number of contributions to OD_rt. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::velocity_update_time_step velocity_update_time_step The time interval over which to update the ice shelf velocity using the nonlinear elliptic equation, or 0 to update every timestep [T ~> s]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::elapsed_velocity_time elapsed_velocity_time The elapsed time since the ice velocities were last updated [T ~> s]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::g_earth g_earth The gravitational acceleration [L2 Z-1 T-2 ~> m s-2]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::density_ice density_ice A typical density of ice [R ~> kg m-3]. logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::gl_regularize gl_regularize Specifies whether to regularize the floatation condition at the grounding line as in Goldberg Holland Schoof 2009. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::n_sub_regularize n_sub_regularize partition of cell over which to integrate for interpolated grounding line the (rectangular) is divided into nxn equally-sized rectangles, over which basal contribution is integrated (iterative quadrature) logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::gl_couple gl_couple whether to let the floatation condition be determined by ocean column thickness means update_OD_ffrac will be called (note: GL_regularize and GL_couple should be exclusive) real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::cfl_factor cfl_factor A factor used to limit subcycled advective timestep in uncoupled runs i.e. dt <= CFL_factor * min(dx / u) real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::a_glen_isothermal a_glen_isothermal Ice viscosity parameter in Glen's Law, [Pa-3 s-1]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::n_glen n_glen Nonlinearity exponent in Glen's Law. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::eps_glen_min eps_glen_min Min. strain rate to avoid infinite Glen's law viscosity, [year-1]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::c_basal_friction c_basal_friction Coefficient in sliding law tau_b = C u^(n_basal_fric), in units= Pa (m yr-1)-(n_basal_fric) real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::n_basal_fric n_basal_fric Exponent in sliding law tau_b = C u^(m_slide) real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::density_ocean_avg density_ocean_avg A typical ocean density [R ~> kg m-3]. This does not affect ocean circulation or thermodynamics. It is used to estimate the gravitational driving force at the shelf front (until we think of a better way to do it, but any difference will be negligible). real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::thresh_float_col_depth thresh_float_col_depth The water column depth over which the shelf if considered to be floating. logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::moving_shelf_front moving_shelf_front Specify whether to advance shelf front (and calve). logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::calve_to_mask calve_to_mask If true, calve off the ice shelf when it passes the edge of a mask. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::min_thickness_simple_calve min_thickness_simple_calve min. ice shelf thickness criteria for calving [Z ~> m]. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::cg_tolerance cg_tolerance The tolerance in the CG solver, relative to initial residual, that determines when to stop the conjugate gradient iterations. real real mom_ice_shelf_dynamics::ice_shelf_dyn_cs::nonlinear_tolerance nonlinear_tolerance The fractional nonlinear tolerance, relative to the initial error, that sets when to stop the iterative velocity solver. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::cg_max_iterations cg_max_iterations The maximum number of iterations that can be used in the CG solver. integer integer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::nonlin_solve_err_mode nonlin_solve_err_mode 1: exit vel solve based on nonlin residual 2: exit based on "fixed point" metric (|u - u_last| / |u| < tol) where | | is infty-norm logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::debug debug If true, write verbose checksums for debugging purposes and use reproducible sums. logical logical mom_ice_shelf_dynamics::ice_shelf_dyn_cs::module_is_initialized module_is_initialized = .false. True if this module has been initialized. type(diag_ctrl), pointer type(diag_ctrl), pointer mom_ice_shelf_dynamics::ice_shelf_dyn_cs::diag diag => NULL() A structure that is used to control diagnostic output. The control structure for the ice shelf dynamics. real logical diag_ctrl mom_ice_shelf_dynamics::ice_shelf_dyn_cs moving_shelf_front calve_to_mask gl_couple gl_regularize debug module_is_initialized diag id_ground_frac id_od_av cg_max_iterations id_t_shelf id_col_thick id_u_shelf id_v_shelf nonlin_solve_err_mode od_rt_counter id_v_mask ... v_face_mask ground_frac u_face_mask_bdry a_glen_isothermal thresh_float_col_depth n_basal_fric basal_traction u_bdry_val u_face_mask od_rt ... integer mom_ice_shelf_dynamics::ice_shelf_dyn_csa_glen_isothermal mom_ice_shelf_dynamics::ice_shelf_dyn_csbasal_traction mom_ice_shelf_dynamics::ice_shelf_dyn_csc_basal_friction mom_ice_shelf_dynamics::ice_shelf_dyn_cscalve_mask mom_ice_shelf_dynamics::ice_shelf_dyn_cscalve_to_mask mom_ice_shelf_dynamics::ice_shelf_dyn_cscfl_factor mom_ice_shelf_dynamics::ice_shelf_dyn_cscg_max_iterations mom_ice_shelf_dynamics::ice_shelf_dyn_cscg_tolerance mom_ice_shelf_dynamics::ice_shelf_dyn_csdebug mom_ice_shelf_dynamics::ice_shelf_dyn_csdensity_ice mom_ice_shelf_dynamics::ice_shelf_dyn_csdensity_ocean_avg mom_ice_shelf_dynamics::ice_shelf_dyn_csdiag mom_ice_shelf_dynamics::ice_shelf_dyn_cselapsed_velocity_time mom_ice_shelf_dynamics::ice_shelf_dyn_cseps_glen_min mom_ice_shelf_dynamics::ice_shelf_dyn_csg_earth mom_ice_shelf_dynamics::ice_shelf_dyn_csgl_couple mom_ice_shelf_dynamics::ice_shelf_dyn_csgl_regularize mom_ice_shelf_dynamics::ice_shelf_dyn_csground_frac mom_ice_shelf_dynamics::ice_shelf_dyn_csground_frac_rt mom_ice_shelf_dynamics::ice_shelf_dyn_csh_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csice_visc mom_ice_shelf_dynamics::ice_shelf_dyn_csid_col_thick mom_ice_shelf_dynamics::ice_shelf_dyn_csid_ground_frac mom_ice_shelf_dynamics::ice_shelf_dyn_csid_od_av mom_ice_shelf_dynamics::ice_shelf_dyn_csid_t_mask mom_ice_shelf_dynamics::ice_shelf_dyn_csid_t_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csid_u_mask mom_ice_shelf_dynamics::ice_shelf_dyn_csid_u_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csid_v_mask mom_ice_shelf_dynamics::ice_shelf_dyn_csid_v_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csmin_thickness_simple_calve mom_ice_shelf_dynamics::ice_shelf_dyn_csmodule_is_initialized mom_ice_shelf_dynamics::ice_shelf_dyn_csmoving_shelf_front mom_ice_shelf_dynamics::ice_shelf_dyn_csn_basal_fric mom_ice_shelf_dynamics::ice_shelf_dyn_csn_glen mom_ice_shelf_dynamics::ice_shelf_dyn_csn_sub_regularize mom_ice_shelf_dynamics::ice_shelf_dyn_csnonlin_solve_err_mode mom_ice_shelf_dynamics::ice_shelf_dyn_csnonlinear_tolerance mom_ice_shelf_dynamics::ice_shelf_dyn_csod_av mom_ice_shelf_dynamics::ice_shelf_dyn_csod_rt mom_ice_shelf_dynamics::ice_shelf_dyn_csod_rt_counter mom_ice_shelf_dynamics::ice_shelf_dyn_cst_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_cst_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csthickness_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csthresh_float_col_depth mom_ice_shelf_dynamics::ice_shelf_dyn_cstmask mom_ice_shelf_dynamics::ice_shelf_dyn_csu_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csu_face_mask mom_ice_shelf_dynamics::ice_shelf_dyn_csu_face_mask_bdry mom_ice_shelf_dynamics::ice_shelf_dyn_csu_flux_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csu_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csumask mom_ice_shelf_dynamics::ice_shelf_dyn_csv_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csv_face_mask mom_ice_shelf_dynamics::ice_shelf_dyn_csv_face_mask_bdry mom_ice_shelf_dynamics::ice_shelf_dyn_csv_flux_bdry_val mom_ice_shelf_dynamics::ice_shelf_dyn_csv_shelf mom_ice_shelf_dynamics::ice_shelf_dyn_csvelocity_update_time_step mom_ice_shelf_dynamics::ice_shelf_dyn_csvmask