Index: /issm/trunk-jpl/src/c/EnumDefinitions/EnumDefinitions.h
===================================================================
--- /issm/trunk-jpl/src/c/EnumDefinitions/EnumDefinitions.h	(revision 11191)
+++ /issm/trunk-jpl/src/c/EnumDefinitions/EnumDefinitions.h	(revision 11192)
@@ -126,4 +126,5 @@
 	MiscellaneousNameEnum, //FIXME: only used by qmu, should not be marshalled (already in queueing script)
 	PrognosticHydrostaticAdjustmentEnum,
+	PrognosticMinThicknessEnum,
 	PrognosticPenaltyFactorEnum,
 	PrognosticSpcthicknessEnum,
Index: /issm/trunk-jpl/src/c/modules/EnumToStringx/EnumToStringx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/EnumToStringx/EnumToStringx.cpp	(revision 11191)
+++ /issm/trunk-jpl/src/c/modules/EnumToStringx/EnumToStringx.cpp	(revision 11192)
@@ -130,4 +130,5 @@
 		case MiscellaneousNameEnum : return "MiscellaneousName";
 		case PrognosticHydrostaticAdjustmentEnum : return "PrognosticHydrostaticAdjustment";
+		case PrognosticMinThicknessEnum : return "PrognosticMinThickness";
 		case PrognosticPenaltyFactorEnum : return "PrognosticPenaltyFactor";
 		case PrognosticSpcthicknessEnum : return "PrognosticSpcthickness";
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 11191)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 11192)
@@ -51,4 +51,5 @@
 	parameters->AddObject(iomodel->CopyConstantObject(PrognosticStabilizationEnum));
 	parameters->AddObject(iomodel->CopyConstantObject(DiagnosticPenaltyFactorEnum));
+	parameters->AddObject(iomodel->CopyConstantObject(PrognosticMinThicknessEnum));
 	parameters->AddObject(iomodel->CopyConstantObject(PrognosticPenaltyFactorEnum));
 	parameters->AddObject(iomodel->CopyConstantObject(ThermalPenaltyFactorEnum));
Index: /issm/trunk-jpl/src/c/modules/StringToEnumx/StringToEnumx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/StringToEnumx/StringToEnumx.cpp	(revision 11191)
+++ /issm/trunk-jpl/src/c/modules/StringToEnumx/StringToEnumx.cpp	(revision 11192)
@@ -128,4 +128,5 @@
 	else if (strcmp(name,"MiscellaneousName")==0) return MiscellaneousNameEnum;
 	else if (strcmp(name,"PrognosticHydrostaticAdjustment")==0) return PrognosticHydrostaticAdjustmentEnum;
+	else if (strcmp(name,"PrognosticMinThickness")==0) return PrognosticMinThicknessEnum;
 	else if (strcmp(name,"PrognosticPenaltyFactor")==0) return PrognosticPenaltyFactorEnum;
 	else if (strcmp(name,"PrognosticSpcthickness")==0) return PrognosticSpcthicknessEnum;
Index: /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11191)
+++ /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11192)
@@ -1785,5 +1785,5 @@
 	int    i,hydroadjustment;
 	int*   doflist = NULL;
-	double rho_ice,rho_water;
+	double rho_ice,rho_water,minthickness;
 	double newthickness[numdof];
 	double newbed[numdof];
@@ -1801,9 +1801,10 @@
 
 	/*Use the dof list to index into the solution vector and extrude it */
+	this->parameters->FindParam(&minthickness,PrognosticMinThicknessEnum);
 	for(i=0;i<numdof2d;i++){
 		newthickness[i]=solution[doflist[i]];
 		if(isnan(newthickness[i])) _error_("NaN found in solution vector");
 		/*Constrain thickness to be at least 1m*/
-		if(newthickness[i]<1) newthickness[i]=1;
+		if(newthickness[i]<minthickness) newthickness[i]=minthickness;
 		newthickness[i+numdof2d]=newthickness[i];
 	}
Index: /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp	(revision 11191)
+++ /issm/trunk-jpl/src/c/objects/Elements/Tria.cpp	(revision 11192)
@@ -1603,5 +1603,5 @@
 	int       i,hydroadjustment;
 	int*      doflist=NULL;
-	double    rho_ice,rho_water;
+	double    rho_ice,rho_water,minthickness;
 	double    newthickness[numdof];
 	double    newbed[numdof];
@@ -1615,9 +1615,10 @@
 
 	/*Use the dof list to index into the solution vector: */
+	this->parameters->FindParam(&minthickness,PrognosticMinThicknessEnum);
 	for(i=0;i<numdof;i++){
 		newthickness[i]=solution[doflist[i]];
 		if(isnan(newthickness[i])) _error_("NaN found in solution vector");
 		/*Constrain thickness to be at least 1m*/
-		if(newthickness[i]<1) newthickness[i]=1;
+		if(newthickness[i]<minthickness) newthickness[i]=minthickness;
 	}
 
Index: /issm/trunk-jpl/src/m/classes/prognostic.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/prognostic.m	(revision 11191)
+++ /issm/trunk-jpl/src/m/classes/prognostic.m	(revision 11192)
@@ -7,4 +7,5 @@
 	properties (SetAccess=public) 
 		 spcthickness           = NaN;
+		 min_thickness          = 0;
 		 hydrostatic_adjustment = 0;
 		 stabilization          = 0;
@@ -36,4 +37,7 @@
 			obj.penalty_factor=3;
 
+			%Minimum ice thickness that can be used
+			obj.min_thickness=1;
+
 			%Hydrostatic adjustment
 			obj.hydrostatic_adjustment='Absolute';
@@ -47,4 +51,5 @@
 			checkfield(md,'prognostic.hydrostatic_adjustment','values',{'Absolute' 'Incremental'});
 			checkfield(md,'prognostic.stabilization','values',[0 1 3]);
+			checkfield(md,'prognostic.min_thickness','>',0);
 
 		end % }}}
@@ -52,4 +57,5 @@
 			disp(sprintf('   Prognostic solution parameters:'));
 			fielddisplay(obj,'spcthickness','thickness constraints (NaN means no constraint)');
+			fielddisplay(obj,'min_thickness','minimum ice thickness allowed');
 			fielddisplay(obj,'hydrostatic_adjustment','adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' ');
 			fielddisplay(obj,'stabilization','0->no, 1->artificial_diffusivity, 3->discontinuous Galerkin');
@@ -62,4 +68,5 @@
 		function marshall(obj,fid) % {{{
 			WriteData(fid,'object',obj,'fieldname','spcthickness','format','DoubleMat','mattype',1);
+			WriteData(fid,'object',obj,'fieldname','min_thickness','format','Double');
 			WriteData(fid,'data',StringToEnum(obj.hydrostatic_adjustment),'format','Integer','enum',PrognosticHydrostaticAdjustmentEnum);
 			WriteData(fid,'object',obj,'fieldname','stabilization','format','Integer');
Index: /issm/trunk-jpl/src/m/enum/PrognosticMinThicknessEnum.m
===================================================================
--- /issm/trunk-jpl/src/m/enum/PrognosticMinThicknessEnum.m	(revision 11192)
+++ /issm/trunk-jpl/src/m/enum/PrognosticMinThicknessEnum.m	(revision 11192)
@@ -0,0 +1,11 @@
+function macro=PrognosticMinThicknessEnum()
+%PROGNOSTICMINTHICKNESSENUM - Enum of PrognosticMinThickness
+%
+%   WARNING: DO NOT MODIFY THIS FILE
+%            this file has been automatically generated by src/c/EnumDefinitions/Synchronize.sh
+%            Please read src/c/EnumDefinitions/README for more information
+%
+%   Usage:
+%      macro=PrognosticMinThicknessEnum()
+
+macro=StringToEnum('PrognosticMinThickness');