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

Index: ../trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test2002.py    (revision 21043)
+++ ../trunk-jpl/test/NightlyRun/test2002.py    (revision 21044)
@@ -19,13 +19,15 @@
 
 #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: 
Index: ../trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test2002.m     (revision 21043)
+++ ../trunk-jpl/test/NightlyRun/test2002.m     (revision 21044)
@@ -14,7 +14,7 @@
        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;