Index: ../trunk-jpl/src/m/contrib/adhikari/adhikari2016GMD_slr_solver/SESAWslr.m =================================================================== --- ../trunk-jpl/src/m/contrib/adhikari/adhikari2016GMD_slr_solver/SESAWslr.m (revision 25443) +++ ../trunk-jpl/src/m/contrib/adhikari/adhikari2016GMD_slr_solver/SESAWslr.m (revision 25444) @@ -110,7 +110,7 @@ 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,5 +121,5 @@ vlm = term1+term2; end % now add VLM. -geoid = rsl_no_mass_conservation+vlm; +geoid = geoid+vlm;