Index: /issm/trunk-jpl/src/m/classes/giacaron.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/giacaron.m	(revision 25055)
+++ /issm/trunk-jpl/src/m/classes/giacaron.m	(revision 25056)
@@ -185,4 +185,5 @@
 		end % }}}
 		function marshall(self,prefix,md,fid) % {{{
+			WriteData(fid,prefix,'name','md.gia.model','data',2,'format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','mantle_viscosity','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',self,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1,'scale',10^3); %from km to m
Index: /issm/trunk-jpl/src/m/classes/giaivins.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/giaivins.m	(revision 25055)
+++ /issm/trunk-jpl/src/m/classes/giaivins.m	(revision 25056)
@@ -58,4 +58,5 @@
 		end % }}}
 		function marshall(self,prefix,md,fid) % {{{
+			WriteData(fid,prefix,'name','md.gia.model','data',1,'format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','mantle_viscosity','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',self,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1,'scale',10^3); %from km to m
Index: /issm/trunk-jpl/src/m/classes/giamme.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/giamme.m	(revision 25056)
+++ /issm/trunk-jpl/src/m/classes/giamme.m	(revision 25056)
@@ -0,0 +1,72 @@
+%GIAMME class definition
+%
+%   Usage:
+%      gia=giamme(); %gia class based on a multi-model 
+%                     ensemble  (ex: Caron et al 2017 statistics)
+
+classdef giamme
+	properties (SetAccess=public) 
+
+		modelid; %index into the multi-model ensemble
+		Ngia;  %geoid ensemble
+		Ugia;  %uplift ensemble
+
+	end
+	methods
+		function self = extrude(self,md) % {{{
+			for i=1:length(self.Ngia),
+				self.Ngia{i}=project3d(md,'vector',self.Ngia{i},'type','node','layer',1);
+				self.Ugia{i}=project3d(md,'vector',self.Ugia{i},'type','node','layer',1);
+			end
+		end % }}}
+		function self = giamme(varargin) % {{{
+			switch nargin
+				case 0
+					self=setdefaultparameters(self);
+				case 1
+					self=structtoobj(gia(),varargin{1});
+				otherwise
+					error('constructor not supported');
+			end
+		end % }}}
+		function self = setdefaultparameters(self) % {{{
+
+			self.modelid=0;
+			self.Ngia={};
+			self.Ugia={};
+
+		end % }}}
+		function md = checkconsistency(self,md,solution,analyses) % {{{
+
+			%Early return
+			if ~ismember('SealevelriseAnalysis',analyses), return; end
+			if (strcmp(solution,'TransientSolution') & md.transient.isslr == 0), return; end
+			md = checkfield(md,'field',self.Ngia,'NaN',1,'Inf',1);
+			md = checkfield(md,'field',self.Ugia,'NaN',1,'Inf',1);
+			md = checkfield(md,'field',self.modelid,'NaN',1,'Inf',1,'>=',1,'<=',length(self.Ngia));
+		end % }}}
+		function disp(self) % {{{
+
+			disp(sprintf('   gia mme parameters:'));
+			fielddisplay(self,'modelid','index into the multi-model ensemble, determines which field will be used.');
+			fielddisplay(self,'Ngia','geoid (m/yr).');
+			fielddisplay(self,'Ugia','uplift (m/yr).');
+		end % }}}
+		function marshall(self,prefix,md,fid) % {{{
+
+			WriteData(fid,prefix,'name','md.gia.model','data',3,'format','Integer');
+			WriteData(fid,prefix,'name','md.gia.nummodels','data',length(self.Ngia),'format','Integer');
+			WriteData(fid,prefix,'object',self,'fieldname','modelid','format','Double');
+			WriteData(fid,prefix,'object',self,'fieldname','Ngia','format','DoubleMat');
+			WriteData(fid,prefix,'object',self,'fieldname','Ugia','format','DoubleMat');
+
+		end % }}}
+		function savemodeljs(self,fid,modelname) % {{{
+			
+			writejsdouble(fid,[modelname '.gia.modelid'],self.modelid);
+			writejscellarray(fid,[modelname '.gia.Ngia'],self.Ngia);
+			writejscellarray(fid,[modelname '.gia.Ugia'],self.Ugia);
+
+		end % }}}
+	end
+end