A_bed A_surface aspect_schoof B_bed B_schoof blatter_coarsening_hook blatter_restriction_hook C_bed C_surface compute_2D_principal_strain_rates compute_2D_stresses copy_staggered_vec create D2 diff_centered diff_uphill disable_ew enable_ew H0_schoof H_exact is_marginal ISMIPHOMTest L_schoof m_schoof ssa_test_cfbc_grid ssa_test_const_grid ssa_test_i_grid ssa_test_j_grid ssa_test_linear_grid ssa_test_plug_grid SSAFactory SSAFDSNESConvergenceTest u_exact u_exterior weight zero_above_ice

is_marginal() static bool pism::stressbalance::is_marginal ( int  i, int  j, const array::CellType1 &  cell_type, bool  ssa_dirichlet_bc  ) static Checks if a cell is near or at the ice front. You need to create array::AccessScope object and add cell_type to it. Note that a cell is a CFBC location of one of four direct neighbors is ice-free. If one of the diagonal neighbors is ice-free we don't use the CFBC, but we do need to compute weights used in the SSA discretization (see assemble_matrix()) to avoid differencing across interfaces between icy and ice-free cells. This method ensures that checks in assemble_rhs() and assemble_matrix() are consistent. Definition at line 270 of file SSAFDBase.cc. References pism::array::CellType1::box_int(), pism::mask::ice_free(), pism::mask::ice_free_ocean(), and pism::mask::icy(). Referenced by pism::stressbalance::SSAFDBase::assemble_rhs(), and pism::stressbalance::SSAFDBase::fd_operator().