Index: ../trunk-jpl/src/m/classes/seaice.m
===================================================================
--- ../trunk-jpl/src/m/classes/seaice.m	(revision 18493)
+++ ../trunk-jpl/src/m/classes/seaice.m	(revision 18494)
@@ -20,6 +20,7 @@
 			end
 		end % }}}
 		function md = checkconsistency(obj,md,solution,analyses) % {{{
+			if solution~=SeaiceSolutionEnum(), return; end
 			md = checkfield(md,'fieldname','seaice.thickness','size',[md.mesh.numberofelements 1],'NaN',1,'>=',0);
 			md = checkfield(md,'fieldname','seaice.concentration','size',[md.mesh.numberofelements 1],'NaN',1,'>=',0);
 			md = checkfield(md,'fieldname','seaice.spcvx','size',[md.mesh.numberofvertices 1]);
Index: ../trunk-jpl/src/m/classes/matseaice.m
===================================================================
--- ../trunk-jpl/src/m/classes/matseaice.m	(revision 18493)
+++ ../trunk-jpl/src/m/classes/matseaice.m	(revision 18494)
@@ -5,6 +5,7 @@
 
 classdef matseaice
 	properties (SetAccess=public) 
+		rho_ice                = 0.;
 		poisson                = 0.;
 		young_modulus          = 0.;
 		damage                 = NaN;
@@ -35,6 +36,9 @@
 		end % }}}
 		function obj = setdefaultparameters(obj) % {{{
 
+			%Ice density
+			obj.rho_ice = 917.; %kg/m^3
+
 			%poisson modulus
 			obj.poisson=0.3;
 
@@ -58,6 +62,7 @@
 
 		end % }}}
 		function md = checkconsistency(obj,md,solution,analyses) % {{{
+			md = checkfield(md,'fieldname','materials.rho_ice','NaN',1,'>',0,'numel',1);
 			md = checkfield(md,'fieldname','materials.poisson','NaN',1,'>',0,'<',.5,'numel',1);
 			md = checkfield(md,'fieldname','materials.young_modulus','NaN',1,'>',0,'numel',1);
 			md = checkfield(md,'fieldname','materials.damage','NaN',1,'>=',0,'<',1,'size',[md.mesh.numberofelements 1]);
@@ -70,6 +75,7 @@
 		function disp(obj) % {{{
 			disp(sprintf('   Sea Ice Material:'));
 
+			fielddisplay(obj,'rho_ice','ice density [kg/m^3]');
 			fielddisplay(obj,'poisson','poisson ratio for undamaged ice [no unit]');
 			fielddisplay(obj,'young_modulus','Young modulus for undamaged ice [Pa]');
 			fielddisplay(obj,'damage','damage, between 0 (no damage) and 1 (fully damaged)');
@@ -81,6 +87,7 @@
 		end % }}}
 		function marshall(obj,md,fid) % {{{
 			WriteData(fid,'enum',MaterialsEnum(),'data',MatseaiceEnum(),'format','Integer');
+			WriteData(fid,'object',obj,'class','materials','fieldname','rho_ice','format','Double');
 			WriteData(fid,'object',obj,'class','materials','fieldname','poisson','format','Double');
 			WriteData(fid,'object',obj,'class','materials','fieldname','young_modulus','format','Double');
 			WriteData(fid,'object',obj,'class','materials','fieldname','damage','format','DoubleMat','mattype',2);
