Changeset 25448

Timestamp:

08/23/20 22:29:05 (5 years ago)

Author:

adhikari

Message:

CHG: more accurate definition of eustatic sea level

File:

• issm/trunk-jpl/src/m/contrib/adhikari/sphericalharmonic_slr_solver/SHslr.m (modified) (1 diff)

issm/trunk-jpl/src/m/contrib/adhikari/sphericalharmonic_slr_solver/SHslr.m

@@ -94 +94 @@
         % Test whether mass is conserved (to ensure the solution is trustworthy!) {{{ 
         % ocean-average sea level and eustatic value must be equal! 
-        avg_eustatic = -force_lm(1)/oce_lm(1); % see equation B16 
-        barystatic = sum(sealevel.*ocean)/nnz(ocean); % we can do this because of equal-area pixelization 
-        sol_diff = abs(abs(barystatic) - abs(avg_eustatic))/abs(barystatic)*100; % per cent 
+        barystatic = -force_lm(1)/oce_lm(1); % see equation B16 
+        gmsl = sum(sealevel.*ocean)/nnz(ocean); % we can do this because of equal-area pixelization 
+        sol_diff = abs(abs(gmsl) - abs(barystatic))/abs(gmsl)*100; % per cent 
         if (sol_diff > 1) % per cent. 
                 error('Ocean-average (eustatic - sealevel) is NOT negligible. Returning.');