Index: /issm/trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 24482)
@@ -10,6 +10,6 @@
 md.slr.sealevel=zeros(md.mesh.numberofvertices,1);
 md.dsl.global_average_thermosteric_sea_level_change=[0;0];
-md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices,1);
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices,1);
+md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
 %antarctica
 late=sum(md.mesh.lat(md.mesh.elements),2)/3;
Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24482)
@@ -19,6 +19,6 @@
 md.slr.sealevel = np.zeros((md.mesh.numberofvertices))
 md.dsl.global_average_thermosteric_sea_level_change=np.zeros((2, ))
-md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices, ))
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices, ))
+md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices+1, ))
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices+1, ))
 
 
Index: /issm/trunk-jpl/test/NightlyRun/test2003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 24482)
@@ -10,6 +10,6 @@
 md.slr.sealevel=zeros(md.mesh.numberofvertices,1);
 md.dsl.global_average_thermosteric_sea_level_change=[0;0];
-md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices,1);
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices,1);
+md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
 
 %antarctica
Index: /issm/trunk-jpl/test/NightlyRun/test2003.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24482)
@@ -19,6 +19,6 @@
 md.slr.sealevel = np.zeros((md.mesh.numberofvertices, ))
 md.dsl.global_average_thermosteric_sea_level_change=np.zeros((2, ))
-md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices, ))
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices, ))
+md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices+1, ))
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices+1, ))
 
 
Index: /issm/trunk-jpl/test/NightlyRun/test2010.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.m	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.m	(revision 24482)
@@ -17,6 +17,6 @@
 md.slr.sealevel=zeros(md.mesh.numberofvertices,1);
 md.dsl.global_average_thermosteric_sea_level_change=[0;0];
-md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices,1);
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices,1);
+md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
 
 md.slr.ocean_area_scaling = 1;
Index: /issm/trunk-jpl/test/NightlyRun/test2010.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24481)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24482)
@@ -27,6 +27,6 @@
 md.slr.sealevel = np.zeros((md.mesh.numberofvertices, ))
 md.dsl.global_average_thermosteric_sea_level_change=np.zeros((2, ))
-md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices, ))
-md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices, ))
+md.dsl.sea_surface_height_change_above_geoid=np.zeros((md.mesh.numberofvertices+1, ))
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros((md.mesh.numberofvertices+1, ))
 
 md.slr.ocean_area_scaling = 1