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;