Index: /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/matice.js
===================================================================
--- /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/matice.js	(revision 20817)
+++ /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/matice.js	(revision 20817)
@@ -0,0 +1,170 @@
+//MATICE class definition
+//
+//   Usage:
+//      matice=matice();
+
+function matice(){
+	//methods
+		this.setdefaultparameters = function(){ // {{{
+
+			//ice density (kg/m^3)
+			this.rho_ice=917.;
+
+			//ocean water density (kg/m^3)
+			this.rho_water=1023.;
+
+			//fresh water density (kg/m^3)
+			this.rho_freshwater=1000.;
+
+			//water viscosity (N.s/m^2)
+			this.mu_water=0.001787;  
+
+			//ice heat capacity cp (J/kg/K)
+			this.heatcapacity=2093.;
+
+			//ice latent heat of fusion L (J/kg)
+			this.latentheat=3.34*10^5;
+
+			//ice thermal conductivity (W/m/K)
+			this.thermalconductivity=2.4;
+			
+			//wet ice thermal conductivity (W/m/K)
+			this.temperateiceconductivity=.24;
+
+			//the melting point of ice at 1 atmosphere of pressure in K
+			this.meltingpoint=273.15;
+
+			//rate of change of melting point with pressure (K/Pa)
+			this.beta=9.8*10^-8;
+
+			//mixed layer (ice-water interface) heat capacity (J/kg/K)
+			this.mixed_layer_capacity=3974.;
+
+			//thermal exchange velocity (ice-water interface) (m/s)
+			this.thermal_exchange_velocity=1.00*10^-4;
+
+			//Rheology law: what is the temperature dependence of B with T
+			//available: none, paterson and arrhenius
+			this.rheology_law='Paterson';
+
+			// GIA:
+			this.lithosphere_shear_modulus  = 6.7*Math.pow(10,10);  // (Pa)
+			this.lithosphere_density        = 3.32;       // (g/cm^-3)
+			this.mantle_shear_modulus       = 1.45*Math.pow(10,11); // (Pa)
+			this.mantle_density             = 3.34;       // (g/cm^-3)
+			
+			//SLR
+			this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
+
+
+		} //}}}
+		this.disp = function() {// {{{
+			console.log(sprintf('   Materials:'));
+
+			fielddisplay(this,'rho_ice','ice density [kg/m^3]');
+			fielddisplay(this,'rho_water','ocean water density [kg/m^3]');
+			fielddisplay(this,'rho_freshwater','fresh water density [kg/m^3]');
+			fielddisplay(this,'mu_water','water viscosity [N s/m^2]');
+			fielddisplay(this,'heatcapacity','heat capacity [J/kg/K]');
+			fielddisplay(this,'thermalconductivity','ice thermal conductivity [W/m/K]');
+			fielddisplay(this,'temperateiceconductivity','temperate ice thermal conductivity [W/m/K]');
+			fielddisplay(this,'meltingpoint','melting point of ice at 1atm in K');
+			fielddisplay(this,'latentheat','latent heat of fusion [J/kg]');
+			fielddisplay(this,'beta','rate of change of melting point with pressure [K/Pa]');
+			fielddisplay(this,'mixed_layer_capacity','mixed layer capacity [W/kg/K]');
+			fielddisplay(this,'thermal_exchange_velocity','thermal exchange velocity [m/s]');
+			fielddisplay(this,'rheology_B','flow law parameter [Pa/s^(1/n)]');
+			fielddisplay(this,'rheology_n',"Glen's flow law exponent");
+			fielddisplay(this,'rheology_law',"law for the temperature dependance of the rheology: 'None', 'Cuffey', 'CuffeyTemperate', 'Paterson', 'Arrhenius' or 'LliboutryDuval'");
+			fielddisplay(this,'lithosphere_shear_modulus','Lithosphere shear modulus [Pa]');
+			fielddisplay(this,'lithosphere_density','Lithosphere density [g/cm^-3]');
+			fielddisplay(this,'mantle_shear_modulus','Mantle shear modulus [Pa]');
+			fielddisplay(this,'mantle_density','Mantle density [g/cm^-3]');
+			fielddisplay(this,'earth_density','Mantle density [kg/m^-3]');
+
+		} // }}}
+        this.extrude = function(md) {
+            this.rheology_B=project3d(md,'vector',this.rheology_B,'type','node');
+            this.rheology_n=project3d(md,'vector',this.rheology_n,'type','element');
+            return this;
+        }
+		this.classname = function() {// {{{
+			return "matice";
+		} // }}}
+		this.checkconsistency = function(md,solution,analyses) { // {{{
+			checkfield(md,'fieldname','materials.rho_ice','>',0);
+			checkfield(md,'fieldname','materials.rho_water','>',0);
+			checkfield(md,'fieldname','materials.rho_freshwater','>',0);
+			checkfield(md,'fieldname','materials.mu_water','>',0);
+			checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1,'Inf',1);
+			checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements,1]);
+			checkfield(md,'fieldname','materials.rheology_law','values',['None','Cuffey','CuffeyTemperate','Paterson','Arrhenius','LliboutryDuval']);
+
+			if(ArrayAnyEqual(ArrayIsMember(GiaAnalysisEnum(),analyses),1)){
+				checkfield(md,'fieldname','materials.lithosphere_shear_modulus','>',0,'numel',1);
+				checkfield(md,'fieldname','materials.lithosphere_density','>',0,'numel',1);
+				checkfield(md,'fieldname','materials.mantle_shear_modulus','>',0,'numel',1);
+				checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',1);
+			}
+			if (ArrayAnyEqual(ArrayIsMember(SealevelriseAnalysisEnum(),analyses),1)){
+				checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
+			}
+
+
+		} // }}}
+		this.marshall=function(md,prefix,fid) { //{{{
+			WriteData(fid,prefix,'name','md.materials.type','data',MaticeEnum(),'format','Integer');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','rho_ice','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','rho_water','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','rho_freshwater','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','mu_water','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','heatcapacity','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','latentheat','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','thermalconductivity','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','temperateiceconductivity','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','meltingpoint','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','beta','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','mixed_layer_capacity','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','thermal_exchange_velocity','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1);
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','rheology_n','format','DoubleMat','mattype',2);
+			WriteData(fid,prefix,'data',StringToEnum(this.rheology_law),'name','md.materials.rheology_law','format','Integer');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','lithosphere_shear_modulus','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','lithosphere_density','format','Double','scale',10^3);
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','mantle_shear_modulus','format','Double');
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','mantle_density','format','Double','scale',10^3);
+			WriteData(fid,prefix,'object',this,'class','materials','fieldname','earth_density','format','Double');
+
+		}//}}}
+		this.fix=function() { //{{{
+		}//}}}
+	//properties 
+	// {{{
+		this.rho_ice                    = 0.;
+		this.rho_water                  = 0.;
+		this.rho_freshwater             = 0.;
+		this.mu_water                   = 0.;
+		this.heatcapacity               = 0.;
+		this.latentheat                 = 0.;
+		this.thermalconductivity        = 0.;
+		this.temperateiceconductivity   = 0.;
+		this.meltingpoint               = 0.;
+		this.beta                       = 0.;
+		this.mixed_layer_capacity       = 0.;
+		this.thermal_exchange_velocity  = 0.;
+		this.rheology_B   = NaN;
+		this.rheology_n   = NaN;
+		this.rheology_law = '';
+
+		//gia: 
+		this.lithosphere_shear_modulus  = 0.;
+		this.lithosphere_density        = 0.;
+		this.mantle_shear_modulus       = 0.;
+		this.mantle_density             = 0.;
+
+		//SLR
+		this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
+
+		this.setdefaultparameters();
+		//}}}
+}
Index: /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/model.js
===================================================================
--- /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/model.js	(revision 20816)
+++ /issm/trunk-jpl/jenkins/javascript/karma/lib/bin/model.js	(revision 20817)
@@ -380,5 +380,5 @@
                 }
 
-                md.materials=extrude(md.materials,md);
+                md.materials=md.materials.extrude(md);
                 md.damage=extrude(md.damage,md);
                 md.mask=extrude(md.mask,md);
Index: /issm/trunk-jpl/src/m/classes/matice.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/matice.js	(revision 20816)
+++ /issm/trunk-jpl/src/m/classes/matice.js	(revision 20817)
@@ -84,4 +84,9 @@
 
 		} // }}}
+        this.extrude = function(md) {
+            this.rheology_B=project3d(md,'vector',this.rheology_B,'type','node');
+            this.rheology_n=project3d(md,'vector',this.rheology_n,'type','element');
+            return this;
+        }
 		this.classname = function() {// {{{
 			return "matice";
Index: /issm/trunk-jpl/src/m/classes/model.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.js	(revision 20816)
+++ /issm/trunk-jpl/src/m/classes/model.js	(revision 20817)
@@ -380,5 +380,5 @@
                 }
 
-                md.materials=extrude(md.materials,md);
+                md.materials=md.materials.extrude(md);
                 md.damage=extrude(md.damage,md);
                 md.mask=extrude(md.mask,md);