Changeset 21626 for issm/trunk-jpl/src/m/classes/calvingdev.m

Timestamp:

03/23/17 15:47:16 (8 years ago)

Author:

youngmc3

Message:

CHG: added calving stress thresholds and calving max parameters

File:

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

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

@@ -6 +6 @@
 classdef calvingdev
         properties (SetAccess=public) 
-                coeff         = NaN;
+                stress_threshold_groundedice = 0.;
+                stress_threshold_floatingice = 0.;
                 meltingrate   = NaN;
         end
@@ -29 +30 @@
                 end % }}}
                 function self = extrude(self,md) % {{{
-                        self.coeff=project3d(md,'vector',self.coeff,'type','node');
                         self.meltingrate=project3d(md,'vector',self.meltingrate,'type','node');
                 end % }}}
@@ -35 +35 @@
 
                         %Proportionality coefficient in Pi model
-                        self.coeff=2e13;
+                        %self.coeff=2e13;
+
+                        %Default sigma max
+                        self.stress_threshold_groundedice = 1e6;
+                        self.stress_threshold_floatingice = 150e3;
                 end % }}}
                 function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -41 +45 @@
                         if (~strcmp(solution,'TransientSolution') | md.transient.ismovingfront==0), return; end
 
-                        md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
+                        md = checkfield(md,'fieldname','calving.stress_threshold_groundedice','>',0,'numel',1,'nan',1,'Inf',1);
+                        md = checkfield(md,'fieldname','calving.stress_threshold_floatingice','>',0,'numel',1,'nan',1,'Inf',1);
                         md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'timeseries',1,'>=',0);
                 end % }}}
                 function disp(self) % {{{
                         disp(sprintf('   Calving Pi parameters:'));
-                        fielddisplay(self,'coeff','proportionality coefficient in Pi model');
+                        fielddisplay(self,'stress_threshold_groundedice','sigma_max applied to grounded ice only [Pa]');
+                        fielddisplay(self,'stress_threshold_floatingice','sigma_max applied to floating ice only [Pa]');
                         fielddisplay(self,'meltingrate','melting rate at given location [m/a]');
 
@@ -53 +59 @@
                         yts=md.constants.yts;
                         WriteData(fid,prefix,'name','md.calving.law','data',2,'format','Integer');
-                        WriteData(fid,prefix,'object',self,'fieldname','coeff','format','DoubleMat','mattype',1);
-                        WriteData(fid,prefix,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1,'scale',1./yts);
+                        WriteData(fid,prefix,'object',self,'fieldname','stress_threshold_groundedice','format','DoubleMat','mattype',1);
+                        WriteData(fid,prefix,'object',self,'fieldname','stress_threshold_floatingice','format','DoubleMat','mattype',1);
+                        WriteData(fid,prefix,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts);
                 end % }}}
         end