Changeset 24666 for issm/trunk-jpl/src/m/classes/frictionjosh.m

Timestamp:

03/16/20 02:01:18 (5 years ago)

Author:

rueckamp

Message:

CHG: add a limit for effective pressure

File:

• issm/trunk-jpl/src/m/classes/frictionjosh.m (modified) (5 diffs)

issm/trunk-jpl/src/m/classes/frictionjosh.m

@@ -9 +9 @@
                 pressure_adjusted_temperature = NaN;
                 gamma      = 0.;
+                effective_pressure_limit = 0;
         end
         methods
@@ -30 +31 @@
                         self.gamma = 1.;
 
+                        %Default 0
+                        self.effective_pressure_limit = 0;
+
                 end % }}}
                 function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -39 +43 @@
                         md = checkfield(md,'fieldname','friction.pressure_adjusted_temperature','NaN',1,'Inf',1);
                         md = checkfield(md,'fieldname','friction.gamma','numel',1,'NaN',1,'Inf',1,'>',0.);
+                        md = checkfield(md,'fieldname','friction.effective_pressure_limit','numel',[1],'>=',0);
 
                         %Check that temperature is provided
@@ -48 +53 @@
                         fielddisplay(self,'pressure_adjusted_temperature','friction pressure_adjusted_temperature (T - Tpmp) [K]');
                         fielddisplay(self,'gamma','(T - Tpmp)/gamma [K]');
+                        fielddisplay(self,'effective_pressure_limit','Neff do not allow to fall below a certain limit: effective_pressure_limit*rho_ice*g*thickness (default 0)');
                 end % }}}
                 function marshall(self,prefix,md,fid) % {{{
@@ -55 +61 @@
                         WriteData(fid,prefix,'class','friction','object',self,'fieldname','pressure_adjusted_temperature','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
                         WriteData(fid,prefix,'class','friction','object',self,'fieldname','gamma','format','Double');
+                        WriteData(fid,prefix,'object',self,'class','friction','fieldname','effective_pressure_limit','format','Double');
                 end % }}}
         end