Changeset 26288

Timestamp:

05/24/21 11:37:16 (4 years ago)

Author:

adhikari

Message:

CHG: grace function requires info about whether we want to map data on the vertices or elemental centroids

File:

• issm/trunk-jpl/examples/EsaGRACE/runme.m (modified) (2 diffs)

issm/trunk-jpl/examples/EsaGRACE/runme.m

@@ -2 +2 @@
 addpath('../Data','../Functions');
 
-steps=[1];
+steps=[1:5];
 
 if any(steps==1)
@@ -27 +27 @@
         time_range = [year_month year_month];
 
-        water_load = grace(md.mesh.elements,md.mesh.lat,md.mesh.long,time_range(1),time_range(2));
+        onvertex = 0; % map data on vertex. If 0, it maps on the elemental centroid. 
+        water_load = grace(md.mesh.elements,md.mesh.lat,md.mesh.long,time_range(1),time_range(2),onvertex);
 
         md.esa.deltathickness = water_load*md.materials.rho_freshwater/md.materials.rho_ice; % ice height equivalent