def standardTillPhiVec(grid,name='tillphi'):
tillphi = IceModelVec2S()
tillphi.create(grid,name,True,WIDE_STENCIL)
# // ghosted to allow the "redundant" computation of tauc
# // PROPOSED standard_name = land_ice_basal_material_friction_angle
tillphi.set_attrs("climate_steady", "friction angle for till under grounded ice sheet",
"degrees", "")
tillphi.time_independent = True
return tillphi
def standardBasalMeltRateVec(grid,name='bmelt'):
bmr = IceModelVec2S()
bmr.create(grid, "bmelt", True, WIDE_STENCIL);
bmr.set_attrs("model_state",
"ice basal melt rate in ice thickness per time",
"m s-1", "land_ice_basal_melt_rate")
bmr.set_glaciological_units("m year-1")
bmr.write_in_glaciological_units = True;
bmr.set_attr("comment", "positive basal melt rate corresponds to ice loss");
return bmr
def standardCBar(grid,name="cbar"):
cbar = IceModelVec2S()
cbar.create(grid, name, False);
cbar.set_attrs("diagnostic",
"magnitude of vertically-integrated horizontal velocity of ice",
"m s-1", "")
cbar.set_glaciological_units("m year-1")
cbar.set_attr("valid_min", 0.0);
cbar.write_in_glaciological_units = True
return cbar
def standardBasalWaterVec(grid,name='bwat'):
bwat = IceModelVec2S()
bwat.create(grid,name,True,WIDE_STENCIL)
bwat.set_attrs("model_state", "effective thickness of subglacial melt water",
"m", "")
#// NB! Effective thickness of subglacial melt water *does* vary from 0 to hmelt_max meters only.
bwat.set_attr("valid_min", 0.0)
valid_max = PISM.global_config().get("hmelt_max")
bwat.set_attr("valid_max", valid_max )
return bwat
def standardYieldStressVec(grid,name='tauc'):
# yield stress for basal till (plastic or pseudo-plastic model)
tauc = IceModelVec2S()
tauc.create(grid, name, True, WIDE_STENCIL);
tauc.set_attrs("diagnostic", "yield stress for basal till (plastic or pseudo-plastic model)", "Pa", "");
return tauc;
def standardBedrockElevationVec(grid,name='topg'):
bed = IceModelVec2S()
bed.create(grid, name, True, WIDE_STENCIL);
bed.set_attrs("model_state", "bedrock surface elevation", "m", "bedrock_altitude");
return bed
def standardIceThicknessVec(grid,name='thk'):
thickness = IceModelVec2S();
thickness.create(grid, name, True, WIDE_STENCIL);
thickness.set_attrs("model_state", "land ice thickness", "m", "land_ice_thickness");
thickness.set_attr("valid_min", 0.0);
return thickness
def standardIceSurfaceVec(grid,name='usurf'):
surface = IceModelVec2S();
surface.create(grid, name, True, WIDE_STENCIL)
surface.set_attrs("diagnostic", "ice upper surface elevation", "m", "surface_altitude");
return surface;