Index: /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.m	(revision 26443)
+++ /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.m	(revision 26444)
@@ -6,5 +6,6 @@
 classdef calvingcrevassedepth
 	properties (SetAccess=public) 
-		crevasse_opening_stress=1.;
+		crevasse_opening_stress=1;
+		crevasse_threshold     =1.;
 		water_height = 0.;
 	end
@@ -32,6 +33,7 @@
 		function self = setdefaultparameters(self) % {{{
 			
-			self.water_height = 0.;
-			crevasse_opening_stress=1.;
+			crevasse_threshold      = 1.;
+			crevasse_opening_stress = 1;
+         self.water_height       = 0.;
 		end % }}}
 		function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -40,4 +42,5 @@
 
 			md = checkfield(md,'fieldname','calving.crevasse_opening_stress','numel',[1],'values',[0,1]);
+         md = checkfield(md,'fieldname','calving.crevasse_threshold','numel',[1],'values',[0,1]);
 			md = checkfield(md,'fieldname','calving.water_height','NaN',1,'Inf',1,'timeseries',1,'>=',0);
 		end % }}}
@@ -45,4 +48,5 @@
 			disp(sprintf('   Calving Pi parameters:'));
 			fielddisplay(self,'crevasse_opening_stress','0: stress only in the ice-flow direction, 1: max principal');
+			fielddisplay(self,'crevasse_threshold','ratio of full thickness to calve (e.g. 0.75 is for 75% of the total ice thickness)');
 			fielddisplay(self,'water_height','water height in the crevasse [m]');
 
@@ -52,4 +56,5 @@
 			WriteData(fid,prefix,'name','md.calving.law','data',6,'format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','crevasse_opening_stress','format','Integer');
+         WriteData(fid,prefix,'object',self,'fieldname','crevasse_threshold','format','Double');
 			WriteData(fid,prefix,'object',self,'fieldname','water_height','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
 		end % }}}
Index: /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.py	(revision 26443)
+++ /issm/trunk-jpl/src/m/classes/calvingcrevassedepth.py	(revision 26444)
@@ -14,24 +14,22 @@
     def __init__(self):  # {{{
         self.crevasse_opening_stress = 1
-        self.water_height = 0
+        self.crevasse_threshold      = 1.
+        self.water_height            = 0.
 
         #self.setdefaultparameters() # Uncomment if/when setdefaultparameters is used
     #}}}
-
     def __repr__(self):  # {{{
         s = '   Calving Pi parameters:'
         s += '{}\n'.format(fielddisplay(self, 'crevasse_opening_stress', '0: stress only in the ice-flow direction, 1: max principal'))
+        s += '{}\n'.format(fielddisplay(self, 'crevasse_threshold', 'ratio of full thickness to calve (e.g. 0.75 is for 75% of the total ice thickness)'))
         s += '{}\n'.format(fielddisplay(self, 'water_height', 'water height in the crevasse [m]'))
         return s
     #}}}
-
     def setdefaultparameters(self):  # {{{
         return self
     #}}}
-
     def extrude(self, md):  # {{{
         return self
     #}}}
-
     def checkconsistency(self, md, solution, analyses):  # {{{
         #Early return
@@ -40,13 +38,14 @@
 
         md = checkfield(md, 'fieldname', 'calving.crevasse_opening_stress', 'numel', [1], 'values', [0,1])
+        md = checkfield(md, 'fieldname', 'calving.crevasse_threshold', 'numel', [1], '>', 0., '<=', 1.)
         md = checkfield(md, 'fieldname', 'calving.water_height', 'NaN', 1, 'Inf', 1, 'timeseries', 1, '>=', 0) 
 
         return md
     # }}}
-
     def marshall(self, prefix, md, fid):  # {{{
         yts = md.constants.yts
         WriteData(fid, prefix, 'name', 'md.calving.law', 'data', 6, 'format', 'Integer')
         WriteData(fid, prefix, 'object', self, 'fieldname', 'crevasse_opening_stress', 'format', 'Integer')
+        WriteData(fid, prefix, 'object', self, 'fieldname', 'crevasse_threshold', 'format', 'Double')
         WriteData(fid, prefix, 'object', self, 'fieldname', 'water_height', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', md.constants.yts)
     # }}}