Index: /issm/trunk-jpl/src/m/classes/groundingline.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/groundingline.js	(revision 28117)
+++ /issm/trunk-jpl/src/m/classes/groundingline.js	(revision 28118)
@@ -11,4 +11,5 @@
 		this.friction_interpolation	= 'SubelementFriction1';
 		this.melt_interpolation		= 'NoMeltOnPartiallyFloating';
+		this.intrusion_distance 	= 0;
 
 	}// }}}
@@ -17,5 +18,6 @@
 		fielddisplay(this,'migration',"type of grounding line migration: 'SubelementMigration','SoftMigration','AggressiveMigration','Contact' or 'None'");
 		fielddisplay(this,'friction_interpolation',"type of friction interpolation forpartially floating elements: 'NoFrictionOnPartiallyFloating','SubelementFriction1' or 'SubelementFriction2'");
-		fielddisplay(this,'melt_interpolation',"type of melt interpolation forpartially floating elements: 'NoMeltOnPartiallyFloating','FullMeltOnPartiallyFloating','SubelementMelt1' or 'SubelementMelt2'");
+		fielddisplay(this,'melt_interpolation',"type of melt interpolation forpartially floating elements: 'NoMeltOnPartiallyFloating','FullMeltOnPartiallyFloating','SubelementMelt1','SubelementMelt2' or 'IntrusionMelt'");
+		fielddisplay(this,'intrusion_distance','distance of seawater intrusion from grounding line [m]');
 
 	}// }}}
@@ -27,5 +29,6 @@
 			checkfield(md,'fieldname','groundingline.migration','values',['None', 'AggressiveMigration', 'SoftMigration', 'Contact', 'GroundingOnly']);
 			checkfield(md,'fieldname','groundingline.friction_interpolation','values',['NoFrictionOnPartiallyFloating', 'SubelementFriction1', 'SubelementFriction2']);
-			checkfield(md,'fieldname','groundingline.melt_interpolation','values',['NoMeltOnPartiallyFloating', 'SubelementMelt1', 'SubelementMelt2', 'FullMeltOnPartiallyFloating']);
+			checkfield(md,'fieldname','groundingline.melt_interpolation','values',['NoMeltOnPartiallyFloating', 'SubelementMelt1', 'SubelementMelt2', 'IntrusionMelt', 'FullMeltOnPartiallyFloating']);
+			checkfield(md,'fieldname','groundingline.intrusion_distance','>=',0)
 
 			if (this.migration !='None'){
@@ -49,4 +52,6 @@
 			WriteData(fid,prefix,'data',this.friction_interpolation,'name','md.groundingline.friction_interpolation','format','String');
 			WriteData(fid,prefix,'data',this.melt_interpolation,'name','md.groundingline.melt_interpolation','format','String');
+			WriteData(fid,prefix,'object',self,'fieldname','intrusion_distance','format','Double');
+
 		}//}}}
 		this.fix=function() { //{{{
@@ -57,4 +62,5 @@
 	this.friction_interpolation = '';
 	this.melt_interpolation     = '';
+	this.intrusion_distance		= 0;
 	this.setdefaultparameters();
 	//}}}
Index: /issm/trunk-jpl/src/m/classes/groundingline.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/groundingline.m	(revision 28117)
+++ /issm/trunk-jpl/src/m/classes/groundingline.m	(revision 28118)
@@ -89,4 +89,5 @@
 			writejsstring(fid,[modelname '.groundingline.friction_interpolation'],self.friction_interpolation);
 			writejsstring(fid,[modelname '.groundingline.melt_interpolation'],self.melt_interpolation);
+			writejsdouble(fif,[modelname '.groundingline.intrusion_distance'],self.intrusion_distance);
 			writejscellstring(fid,[modelname '.groundingline.requested_outputs'],self.requested_outputs);
 
Index: /issm/trunk-jpl/src/m/classes/groundingline.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/groundingline.py	(revision 28117)
+++ /issm/trunk-jpl/src/m/classes/groundingline.py	(revision 28118)
@@ -18,4 +18,5 @@
         self.friction_interpolation = ''
         self.melt_interpolation = ''
+        self.intrusion_distance = 0
         self.requested_outptuts = []
 
@@ -29,5 +30,5 @@
         s += '{}\n'.format(fielddisplay(self, 'migration', 'type of grounding line migration: \'SoftMigration\', \'SubelementMigration\', \'AggressiveMigration\', \'Contact\', \'None\''))
         s += '{}\n'.format(fielddisplay(self, 'migration', 'type of friction interpolation on partially floating elements: ''SubelementFriction1'', ''SubelementFriction2'', ''NoFrictionOnPartiallyFloating'''))
-        s += '{}\n'.format(fielddisplay(self, 'migration', 'type of melt interpolation on partially floating elements: \'SubelementMelt1\', \'SubelementMelt2\', \'NoMeltOnPartiallyFloating\', \'FullMeltOnPartiallyFloating\''))
+        s += '{}\n'.format(fielddisplay(self, 'migration', 'type of melt interpolation on partially floating elements: \'SubelementMelt1\', \'SubelementMelt2\', \'IntrusionMelt\', \'NoMeltOnPartiallyFloating\', \'FullMeltOnPartiallyFloating\''))
         s += '{}\n'.format(fielddisplay(self, 'requested_outputs', 'additional outputs requested'))
         return s
@@ -44,4 +45,5 @@
         self.friction_interpolation = 'SubelementFriction1'
         self.melt_interpolation = 'NoMeltOnPartiallyFloating'
+        self.intrusion_distance = 0
         # Default output
         self.requested_outputs = ['default']
@@ -53,5 +55,6 @@
         md = checkfield(md, 'fieldname', 'groundingline.migration', 'values', ['None', 'SubelementMigration', 'AggressiveMigration', 'SoftMigration', 'Contact', 'GroundingOnly'])
         md = checkfield(md, 'fieldname', 'groundingline.friction_interpolation', 'values', ['SubelementFriction1', 'SubelementFriction2', 'NoFrictionOnPartiallyFloating'])
-        md = checkfield(md, 'fieldname', 'groundingline.melt_interpolation', 'values', ['SubelementMelt1', 'SubelementMelt2', 'NoMeltOnPartiallyFloating', 'FullMeltOnPartiallyFloating'])
+        md = checkfield(md, 'fieldname', 'groundingline.melt_interpolation', 'values', ['SubelementMelt1', 'SubelementMelt2', 'IntrusionMelt', 'NoMeltOnPartiallyFloating', 'FullMeltOnPartiallyFloating'])
+        md = checkfield(md, 'fieldname', 'groundingline.intrusion_distance','>=',0)
         md = checkfield(md, 'fieldname', 'groundingline.requested_outputs', 'stringrow', 1)
 
@@ -72,4 +75,5 @@
         WriteData(fid, prefix, 'data', self.friction_interpolation, 'name', 'md.groundingline.friction_interpolation', 'format', 'String')
         WriteData(fid, prefix, 'data', self.melt_interpolation, 'name', 'md.groundingline.melt_interpolation', 'format', 'String')
+        WriteData(fid, prefix, 'object', self, 'fieldname', 'intrusion_distance', 'format', 'Double')
         
         # Process requested outputs