Index: /issm/trunk-jpl/src/m/parameterization/sethydrostaticmask.m
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/sethydrostaticmask.m	(revision 28085)
+++ /issm/trunk-jpl/src/m/parameterization/sethydrostaticmask.m	(revision 28086)
@@ -1,3 +1,3 @@
-function md=sethydrostaticmask(md)
+function ocean_levelset = sethydrostaticmask(md)
 %SETHYDROSTATICMASK - establish ocean_levelset field
 %
@@ -6,8 +6,12 @@
 %
 %   Usage:
-%      md=sethydrostaticmask(md)
+%      ocean_levelset = sethydrostaticmask(md)
 %
 %   Examples:
-%      md=sethydrostaticmask(md);
+%      md.mask.ocean_levelset = sethydrostaticmask(md);
+
+%temporary warning for people that already use this function
+disp('Setting hydrostatic mask: WARNING - now returning ocean_levelset (not md)');
+disp('   -- you can blame Mathieu');
 
 if(length(md.geometry.bed)~=md.mesh.numberofvertices | length(md.geometry.thickness)~=md.mesh.numberofvertices | length(md.geometry.base)~=md.mesh.numberofvertices),
@@ -15,14 +19,14 @@
 end
 
-%grounded ice level set
-md.mask.ocean_levelset=md.geometry.thickness+md.geometry.bed*md.materials.rho_water/md.materials.rho_ice;
+%ocean level set based on height above floatation
+ocean_levelset=md.geometry.thickness+md.geometry.bed*md.materials.rho_water/md.materials.rho_ice;
 
 %Check consistency of geometry
-pos=find(md.mask.ocean_levelset>0);
+pos=find(ocean_levelset>0);
 if(any(md.geometry.base(pos)~=md.geometry.bed(pos))),
 	disp('WARNING: md.geometry.bed and md.geometry.base not equal on grounded ice');
 end
 
-pos=find(md.mask.ocean_levelset<=0);
+pos=find(ocean_levelset<=0);
 if(any(md.geometry.base(pos)<md.geometry.bed(pos))),
 	disp('WARNING: md.geometry.base < md.geometry.bed on floating ice');