19 #ifndef PART_GRID_THRESHOLD_THICKNESS_H
20 #define PART_GRID_THRESHOLD_THICKNESS_H
22 #include "pism/util/stencils.hh"
27 stencils::Star<double> ice_thickness,
28 stencils::Star<double> surface_elevation,
29 double bed_elevation);
double part_grid_threshold_thickness(stencils::Star< int > cell_type, stencils::Star< double > ice_thickness, stencils::Star< double > surface_elevation, double bed_elevation)
Compute threshold thickness used when deciding if a partially-filled cell should be considered 'full'...