source: issm/oecreview/Archive/20545-21336/ISSM-21043-21044.diff@ 21337

../trunk-jpl/test/NightlyRun/test2002.py

@@ -19 +19 @@
 
 #parameterize slr solution:
 #slr loading:  {{{
-md.slr.deltathickness=zeros((md.mesh.numberofvertices,1))
+md.slr.deltathickness=zeros((md.mesh.numberofelements,1))
 md.slr.sealevel=zeros((md.mesh.numberofvertices,1))
 #antarctica
-pos=numpy.nonzero(md.mesh.lat <-80)
+late=numpy.sum(md.mesh.lat[md.mesh.elements-1],axis=1)/3
+longe=numpy.sum(md.mesh.long[md.mesh.elements-1],axis=1)/3
+pos=numpy.nonzero(late <-80)
 md.slr.deltathickness[pos]=-100
-#greenlnd 
-pos=numpy.nonzero(logical_and_n(md.mesh.lat > 70,md.mesh.lat < 80,md.mesh.long>-60,md.mesh.long<-30))
+#greenland 
+pos=numpy.nonzero(logical_and_n(late > 70,late < 80,longe>-60,longe<-30))
 md.slr.deltathickness[pos]=-100
 
 #elastic loading from love numbers: 

../trunk-jpl/test/NightlyRun/test2002.m

@@ -14 +14 @@
         longe=sum(md.mesh.long(md.mesh.elements),2)/3;
         pos=find(late <-80);
         md.slr.deltathickness(pos)=-100;
-        %greenlnd 
+        %greenland 
         pos=find(late > 70 &  late < 80 & longe>-60 & longe<-30);
         md.slr.deltathickness(pos)=-100;