Changeset 26296 for issm/trunk-jpl/test/NightlyRun/test2004.m

Timestamp:

06/04/21 16:35:40 (4 years ago)

Author:

caronlam

Message:

BUG: BarystaticLoad was considering coastline elements that were fully covered by ice as also fully grounded; CHG: lots of renaming in sealevelchange_core & analysis; CHG: it is now possible to compute GRD with Barystatic contributions & ocean mass but without sea-level loading

File:

• issm/trunk-jpl/test/NightlyRun/test2004.m (modified) (5 diffs)

issm/trunk-jpl/test/NightlyRun/test2004.m

@@ -395 +395 @@
 md.solidearth.settings.reltol=NaN;
 md.solidearth.settings.abstol=1e-3;
-md.solidearth.settings.computesealevelchange=1;
+md.solidearth.settings.sealevelloading=1;
 md.solidearth.settings.isgrd=1;
 md.solidearth.settings.ocean_area_scaling=0;
@@ -423 +423 @@
 
 %eustatic run:
-md.solidearth.settings.rigid=0;
+md.solidearth.settings.selfattraction=0;
 md.solidearth.settings.elastic=0;
 md.solidearth.settings.rotation=0;
@@ -432 +432 @@
 Seustatic=md.results.TransientSolution.Sealevel;
 
-%eustatic + rigid run:
-md.solidearth.settings.rigid=1;
+%eustatic + selfattraction run:
+md.solidearth.settings.selfattraction=1;
 md.solidearth.settings.elastic=0;
 md.solidearth.settings.rotation=0;
 md=solve(md,'Transient');
-Srigid=md.results.TransientSolution.Sealevel;
-
-%eustatic + rigid + elastic run:
-md.solidearth.settings.rigid=1;
+Sselfattraction=md.results.TransientSolution.Sealevel;
+
+%eustatic + selfattraction + elastic run:
+md.solidearth.settings.selfattraction=1;
 md.solidearth.settings.elastic=1;
 md.solidearth.settings.rotation=0;
@@ -446 +446 @@
 Selastic=md.results.TransientSolution.Sealevel;
 
-%eustatic + rigid + elastic + rotation run:
-md.solidearth.settings.rigid=1;
+%eustatic + selfattraction + elastic + rotation run:
+md.solidearth.settings.selfattraction=1;
 md.solidearth.settings.elastic=1;
 md.solidearth.settings.rotation=1;
@@ -457 +457 @@
 field_names     ={'Eustatic','Rigid','Elastic','Rotation'};
 field_tolerances={1e-13,1e-13,1e-13,1e-13};
-field_values={Seustatic,Srigid,Selastic,Srotation};
+field_values={Seustatic,Sselfattraction,Selastic,Srotation};