def __repr__(self): # {{{
string = "Weertman sliding law parameters: Sigma_b = C^(-1/m) * |u_b|^(1/m-1) * u_b"
string = "%s\n%s" % (
string, fielddisplay(self, "C", "friction coefficient [SI]"))
string = "%s\n%s" % (string, fielddisplay(self, "m", "m exponent"))
return string
def __repr__(self): # {{{
string=" debug parameters:"
string="%s\n%s"%(string,fielddisplay(self,"valgrind","use Valgrind to debug (0 or 1)"))
string="%s\n%s"%(string,fielddisplay(self,"gprof","use gnu-profiler to find out where the time is spent"))
string="%s\n%s"%(string,fielddisplay(self,'profiling','enables profiling (memory, flops, time)'))
return string
def __repr__(self): # {{{
string = " MISMIP+ basal melt parameterization\n"
string = "%s\n%s" % (
string,
fielddisplay(self, "groundedice_melting_rate",
"basal melting rate (positive if melting) [m/yr]"))
string = "%s\n%s" % (
string,
fielddisplay(
self, "meltrate_factor",
"Melt-rate rate factor [1/yr] (sign is opposite to MISMIP+ benchmark to remain consistent with ISSM convention of positive values for melting)"
))
string = "%s\n%s" % (
string,
fielddisplay(
self, "threshold_thickness",
"Threshold thickness for saturation of basal melting [m]"))
string = "%s\n%s" % (string,
fielddisplay(
self, "upperdepth_melt",
"Depth above which melt rate is zero [m]"))
string = "%s\n%s" % (string,
fielddisplay(self, "geothermalflux",
"Geothermal heat flux [W/m^2]"))
return string
def __repr__(self): # {{{
s = " automatic differentiation parameters:\n"
s += "%s\n" % fielddisplay(
self, 'isautodiff',
"indicates if the automatic differentiation is activated")
s += "%s\n" % fielddisplay(self, 'dependents',
"list of dependent variables")
s += "%s\n" % fielddisplay(self, 'independents',
"list of independent variables")
s += "%s\n" % fielddisplay(
self, 'driver', "ADOLC driver ('fos_forward' or 'fov_forward')")
s += "%s\n" % fielddisplay(
self, 'obufsize',
"Number of operations per buffer (==OBUFSIZE in usrparms.h)")
s += "%s\n" % fielddisplay(
self, 'lbufsize',
"Number of locations per buffer (==LBUFSIZE in usrparms.h)")
s += "%s\n" % fielddisplay(
self, 'cbufsize',
"Number of values per buffer (==CBUFSIZE in usrparms.h)")
s += "%s\n" % fielddisplay(
self, 'tbufsize',
"Number of taylors per buffer (<=TBUFSIZE in usrparms.h)")
s += "%s\n" % fielddisplay(
self, 'gcTriggerRatio',
"free location block sorting/consolidation triggered if the ratio between allocated and used locations exceeds gcTriggerRatio"
)
s += "%s\n" % fielddisplay(
self, 'gcTriggerMaxSize',
"free location block sorting/consolidation triggered if the allocated locations exceed gcTriggerMaxSize)"
)
return s
def __repr__(self): # {{{
string = " surface forcings parameters:"
string = "%s\n%s" % (
string,
fielddisplay(
self, 'isd18opd',
'is delta18o parametrisation from present day temperature and precipitation activated (0 or 1, default is 0)'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'desfac',
'desertification elevation factor (between 0 and 1, default is 0.5) [m]'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 's0p',
'should be set to elevation from precip source (between 0 and a few 1000s m, default is 0) [m]'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 's0t',
'should be set to elevation from temperature source (between 0 and a few 1000s m, default is 0) [m]'
))
string = "%s\n%s" % (
string,
fielddisplay(self, 'rlaps', 'present day lapse rate [degree/km]'))
if self.isd18opd:
string = "%s\n%s" % (
string,
fielddisplay(
self, 'temperatures_presentday',
'monthly present day surface temperatures [K], required if delta18o/mungsm is activated'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'precipitations_presentday',
'monthly surface precipitation [m/yr water eq], required if delta18o or mungsm is activated'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'delta18o',
'delta18o [per mil], required if pdd is activated and delta18o activated'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'dpermil',
'degree per mil, required if d18opd is activated'))
string = "%s\n%s" % (string,
fielddisplay(self, 'requested_outputs',
'additional outputs requested'))
return string
def __repr__(self): # {{{
string=' private parameters: do not change'
string="%s\n%s"%(string,fielddisplay(self,'isconsistent','is model self consistent'))
string="%s\n%s"%(string,fielddisplay(self,'runtimename','name of the run launched'))
string="%s\n%s"%(string,fielddisplay(self,'bamg','structure with mesh properties constructed if bamg is used to mesh the domain'))
string="%s\n%s"%(string,fielddisplay(self,'solution','type of solution launched'))
return string
def __repr__(self): # {{{
string=' Calving parameters:'
string="%s\n%s"%(string,fielddisplay(self,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding'))
string="%s\n%s"%(string,fielddisplay(self,'spclevelset','levelset constraints (NaN means no constraint)'))
string="%s\n%s"%(string,fielddisplay(self,'calvingrate','calving rate at given location [m/a]'))
string="%s\n%s"%(string,fielddisplay(self,'meltingrate','melting rate at given location [m/a]'))
return string
def __repr__(self): # {{{
string = ' radaroverlay parameters:'
string = "%s\n%s" % (
string, fielddisplay(self, 'pwr', 'radar power image (matrix)'))
string = "%s\n%s" % (
string, fielddisplay(self, 'x', 'corresponding x coordinates [m]'))
string = "%s\n%s" % (
string, fielddisplay(self, 'y', 'corresponding y coordinates [m]'))
return string
def __repr__(self): # {{{
string = "Basal shear stress parameters: Sigma_b = coefficient^2 * Neff ^r * |u_b|^(s-1) * u_b,\n(effective stress Neff=rho_ice*g*thickness+rho_water*g*bed, r=q/p and s=1/p)"
string = "%s\n%s" % (string,
fielddisplay(self, "coefficient",
"friction coefficient [SI]"))
string = "%s\n%s" % (string, fielddisplay(self, "p", "p exponent"))
string = "%s\n%s" % (string, fielddisplay(self, "q", "q exponent"))
return string
def __repr__(self): # {{{
string = ' miscellaneous parameters:'
string = "%s\n%s" % (
string, fielddisplay(self, 'notes', 'notes in a cell of strings'))
string = "%s\n%s" % (string, fielddisplay(self, 'name', 'model name'))
string = "%s\n%s" % (string,
fielddisplay(self, 'dummy',
'empty field to store some data'))
return string
def __repr__(self): # {{{
string = " surface forcings parameters:"
string = "%s\n%s" % (string,
fielddisplay(
self, 'mass_balance',
'surface mass balance [m/yr ice eq]'))
string = "%s\n%s" % (string,
fielddisplay(self, 'requested_outputs',
'additional outputs requested'))
return string
def __repr__(self): # {{{
string = ' rifts parameters:'
string = "%s\n%s" % (
string,
fielddisplay(
self, 'riftstruct',
'structure containing all rift information (vertices coordinates, segments, type of melange, ...)'
))
string = "%s\n%s" % (string, fielddisplay(self, 'riftproperties', ''))
return string
def __repr__(self): # {{{
s =" independent variable:\n"
s+="%s\n" % fielddisplay(self,'name',"variable name (must match corresponding Enum)")
s+="%s\n" % fielddisplay(self,'type',"type of variable ('vertex' or 'scalar')")
if not numpy.isnan(self.fos_forward_index):
s+="%s\n" % fielddisplay(self,'fos_forward_index',"index for fos_foward driver of ADOLC")
if numpy.any(numpy.logical_not(numpy.isnan(self.fov_forward_indices))):
s+="%s\n" % fielddisplay(self,'fov_forward_indices',"indices for fov_foward driver of ADOLC")
return s
def __repr__(self): # {{{
string = ' Calving Levermann parameters:'
string = "%s\n%s" % (
string,
fielddisplay(self, 'coeff',
'proportionality coefficient in Levermann model'))
string = "%s\n%s" % (string,
fielddisplay(
self, 'meltingrate',
'melting rate at given location [m/a]'))
return string
def __repr__(self): # {{{
string = ' steadystate solution parameters:'
string = "%s\n%s" % (string,
fielddisplay(self, 'reltol',
'relative tolerance criterion'))
string = "%s\n%s" % (string,
fielddisplay(self, 'maxiter',
'maximum number of iterations'))
string = "%s\n%s" % (string,
fielddisplay(self, 'requested_outputs',
'additional requested outputs'))
return string
def __repr__(self): # {{{
string = ' Calving parameters:'
string = "%s\n%s" % (string,
fielddisplay(
self, 'calvingrate',
'calving rate at given location [m/a]'))
string = "%s\n%s" % (string,
fielddisplay(
self, 'meltingrate',
'melting rate at given location [m/a]'))
return string
def __repr__(self): # {{{
string = " geometry parameters:"
string = "%s\n%s" % (string,
fielddisplay(self, 'surface',
'ice upper surface elevation [m]'))
string = "%s\n%s" % (
string, fielddisplay(self, 'thickness', 'ice thickness [m]'))
string = "%s\n%s" % (
string, fielddisplay(self, 'base', 'ice base elevation [m]'))
string = "%s\n%s" % (string,
fielddisplay(self, 'bed', 'bed elevation [m]'))
return string
def __repr__(self): # {{{
string = "Basal shear stress parameters: Sigma_b = min(coefficient^2 * Neff ^r * |u_b|^(s-1) * u_b,\n coefficientcoulomb^2 * rho_i * g * (h-h_f)), (effective stress Neff=rho_ice*g*thickness+rho_water*g*bed, r=q/p and s=1/p)."
string = "%s\n%s" % (
string,
fielddisplay(self, "coefficient",
"power law (Weertman) friction coefficient [SI]"))
string = "%s\n%s" % (string,
fielddisplay(self, "coefficientcoulomb",
"Coulomb friction coefficient [SI]"))
string = "%s\n%s" % (string, fielddisplay(self, "p", "p exponent"))
string = "%s\n%s" % (string, fielddisplay(self, "q", "q exponent"))
return string
def __repr__(self): # {{{
s = " Model results:\n"
if 'step' in self.__dict__:
s += "%s\n" % fielddisplay(self, 'step', "step number")
if 'time' in self.__dict__:
s += "%s\n" % fielddisplay(self, 'time', "time value")
if 'SolutionType' in self.__dict__:
s += "%s\n" % fielddisplay(self, 'SolutionType', "solution type")
for name in self.__dict__.iterkeys():
if name not in [
'step', 'time', 'SolutionType', 'errlog', 'outlog'
]:
if isinstance(getattr(self, name), list):
s += "%s\n" % fielddisplay(self, name,
"model results list")
elif isinstance(getattr(self, name), results):
s += "%s\n" % fielddisplay(self, name,
"model results case")
else:
s += "%s\n" % fielddisplay(self, name, "")
if 'errlog' in self.__dict__:
s += "%s\n" % fielddisplay(self, 'errlog', "error log file")
if 'outlog' in self.__dict__:
s += "%s\n" % fielddisplay(self, 'outlog', "output log file")
return s
def __repr__(self): # {{{
string = ' hydrologyshreve solution parameters:'
string = "%s\n%s" % (
string,
fielddisplay(
self, 'spcwatercolumn',
'water thickness constraints (NaN means no constraint) [m]'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'stabilization',
'artificial diffusivity (default is 1). can be more than 1 to increase diffusivity.'
))
return string
def __repr__(self): # {{{
string = " basal forcings parameters:"
string = "%s\n%s" % (
string,
fielddisplay(self, "groundedice_melting_rate",
"basal melting rate (positive if melting) [m/yr]"))
string = "%s\n%s" % (
string,
fielddisplay(self, "floatingice_melting_rate",
"basal melting rate (positive if melting) [m/yr]"))
string = "%s\n%s" % (string,
fielddisplay(self, "geothermalflux",
"geothermal heat flux [W/m^2]"))
return string
def __repr__(self): # {{{
string = " constants parameters:"
string = "%s\n%s" % (
string,
fielddisplay(self, "g", "gravitational acceleration [m/s^2]"))
string = "%s\n%s" % (
string,
fielddisplay(self, "yts", "number of seconds in a year [s/yr]"))
string = "%s\n%s" % (
string,
fielddisplay(
self, "referencetemperature",
"reference temperature used in the enthalpy model [K]"))
return string
def __repr__(self): # {{{
string = " masks:"
string = "%s\n%s" % (
string,
fielddisplay(
self, "groundedice_levelset",
"is ice grounded ? grounded ice if > 0, grounding line position if = 0, floating ice if < 0"
))
string = "%s\n%s" % (
string,
fielddisplay(
self, "ice_levelset",
"presence of ice if < 0, icefront position if = 0, no ice if > 0"
))
return string
def __repr__(self): # {{{
print('enter repr')
s = "List of toolkits options per analysis:\n\n"
for analysis in list(vars(self).keys()):
s += "%s\n" % fielddisplay(self, analysis, '')
return s
def __repr__(self): # {{{
s = " dependent variable:\n"
s += "%s\n" % fielddisplay(
self, 'name', "variable name (must match corresponding String)")
s += "%s\n" % fielddisplay(self, 'type',
"type of variable ('vertex' or 'scalar')")
if not np.isnan(self.fos_reverse_index):
s += "%s\n" % fielddisplay(
self, 'fos_reverse_index',
"index for fos_reverse driver of ADOLC")
if self.exp:
s += "%s\n" % fielddisplay(
self, 'exp', "file needed to compute dependent variable")
s += "%s\n" % fielddisplay(self, 'segments', "mass flux segments")
return s
def __repr__(self): # {{{
string = " surface forcings parameters (SMB=accumulation-runoff-evaporation) :"
string = "%s\n%s" % (string,
fielddisplay(self, 'accumulation',
'accumulated snow [m/yr ice eq]'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'runoff',
'amount of ice melt lost from the ice column [m/yr ice eq]'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'evaporation',
'mount of ice lost to evaporative processes [m/yr ice eq]'))
string = "%s\n%s" % (string,
fielddisplay(self, 'requested_outputs',
'additional outputs requested'))
return string
def __repr__(self): # {{{
string = " Outputdefinitions:"
string = "%s\n%s" % (
string,
fielddisplay(
self, "definitions",
"list of potential outputs that can be requested, but which need additional data to be defined"
))
return string
def __repr__(self): # {{{
string = " Massfluxatgate:"
string = "%s\n%s" % (
string,
fielddisplay(self, 'name',
'identifier for this massfluxatgate response'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'definitionstring',
'string that identifies this output definition uniquely, from Outputdefinition[1-10]'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'profilename',
'name of file (shapefile or argus file) defining a profile (or gate)'
))
return string
def __repr__(self): # {{{
string = ' Thermal solution parameters:'
string = "%s\n%s" % (
string,
fielddisplay(
self, 'spctemperature',
'temperature constraints (NaN means no constraint) [K]'))
string = "%s\n%s" % (string,
fielddisplay(
self, 'stabilization',
'0: no, 1: artificial_diffusivity, 2: SUPG'))
string = "%s\n%s" % (string,
fielddisplay(
self, 'maxiter',
'maximum number of non linear iterations'))
string = "%s\n%s" % (string,
fielddisplay(self, 'reltol',
'relative tolerance criterion'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'penalty_lock',
'stabilize unstable thermal constraints that keep zigzagging after n iteration (default is 0, no stabilization)'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'penalty_threshold',
'threshold to declare convergence of thermal solution (default is 0)'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'isenthalpy',
'use an enthalpy formulation to include temperate ice (default is 0)'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'isdynamicbasalspc',
'enable dynamic setting of basal forcing. required for enthalpy formulation (default is 0)'
))
string = "%s\n%s" % (string,
fielddisplay(self, 'requested_outputs',
'additional outputs requested'))
return string
def __repr__(self): # {{{
string = ' balance thickness solution parameters:'
string = "%s\n%s" % (
string,
fielddisplay(
self, 'spcthickness',
'thickness constraints (NaN means no constraint) [m]'))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'thickening_rate',
'ice thickening rate used in the mass conservation (dh/dt) [m/yr]'
))
string = "%s\n%s" % (
string,
fielddisplay(
self, 'stabilization',
"0: None, 1: SU, 2: SSA's artificial diffusivity, 3:DG"))
return string
