Index: /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/ismip6.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/ismip6.m	(revision 27091)
+++ /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/ismip6.m	(revision 27092)
@@ -22,4 +22,5 @@
 		deltathicknesshal = {};   %placeholder for delta thicknesses halosteric origins
 		deltathicknessbar = {};   %placeholder for delta thicknesses halosteric origins
+		vdeltathicknessbar = {};   %placeholder for delta thicknesses halosteric origins on vertices
 		thicknesscorrection={};  
 		icemask        = {};   %placeholder for ice masks
Index: /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/oceanmip.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/oceanmip.m	(revision 27091)
+++ /issm/branches/trunk-larour-SLPS2022/src/m/contrib/larour/oceanmip.m	(revision 27092)
@@ -258,4 +258,19 @@
 		end
 	end %}}}
+	function arrays= bottompressures(self,varargin) % {{{
+		options=pairoptions(varargin{:});
+		units=getfieldvalue(options,'units','m');
+		
+		arrays=cell(self.n,1);
+		for i=1:self.n,
+			pbo=self.mesh_pbo{i}; 
+			if strcmpi(units,'m'),
+				pbo=pbo*1e-3;
+			end
+			time=self.time{i};
+			array=[pbo;time'];
+			arrays{i}=array;
+		end
+	end %}}}
 	function arrays= dzoss(self,varargin) % {{{
 		options=pairoptions(varargin{:});
@@ -274,4 +289,20 @@
 		end
 	end	%}}}
+	function arrays= zoss(self,varargin) % {{{
+		
+		options=pairoptions(varargin{:});
+		units=getfieldvalue(options,'units','m');
+		
+		arrays=cell(self.n,1);
+		for i=1:self.n,
+			zos=self.mesh_zos{i}; 
+			if strcmpi(units,'m'),
+				zos=zos*1e-3;
+			end
+			time=self.time{i};
+			array=[zos;time'];
+			arrays{i}=array;
+		end
+	end	%}}}
 	function arrays= dzostogass(self,varargin) % {{{
 		options=pairoptions(varargin{:});
@@ -287,4 +318,19 @@
 			time=self.time{i};time=time(1:end-1);
 			array=[dzostoga';time'];
+			arrays{i}=array;
+		end
+	end	%}}}
+	function arrays= zostogas(self,varargin) % {{{
+		options=pairoptions(varargin{:});
+		units=getfieldvalue(options,'units','m');
+
+		arrays=cell(self.n,1);
+		for i=1:self.n,
+			z=self.zostoga{i}; 
+			if strcmpi(units,'m'),
+				z=z*1e-3;
+			end
+			time=self.time{i};
+			array=[z';time'];
 			arrays{i}=array;
 		end