source: issm/oecreview/Archive/24684-25833/ISSM-25443-25444.diff

../trunk-jpl/src/m/contrib/adhikari/adhikari2016GMD_slr_solver/SESAWslr.m

@@ -110 +110 @@
 end 
 
 % compute bedrock motion aka VLM (vertical land motion) and geoid height change 
-rsl_no_mass_conservation = term1+term2; % eus-term3-term4 is excluded. See Tamisiea, 2011. 
+geoid = term1+term2; % eus-term3-term4 is excluded from rsl. See Tamisiea, 2011. 
 if strcmpi(para.solidearth,'rigid')
         vlm = 0.0*rsl;
 elseif strcmpi(para.solidearth,'elastic') 
@@ -121 +121 @@
         vlm = term1+term2; 
 end
 % now add VLM. 
-geoid = rsl_no_mass_conservation+vlm; 
+geoid = geoid+vlm;