Index: /issm/trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 24284)
@@ -1,6 +1,6 @@
 %Test Name: EarthSlr
 
-%mesh earth: 
-md=model; 
+%mesh earth:
+md=model;
 md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',700.); %500 km resolution mesh
 
@@ -15,9 +15,9 @@
 pos=find(late <-80);
 md.slr.deltathickness(pos)=-100;
-%greenland 
+%greenland
 pos=find(late > 70 &  late < 80 & longe>-60 & longe<-30);
 md.slr.deltathickness(pos)=-100;
 
-%elastic loading from love numbers: 
+%elastic loading from love numbers:
 nlov=101;
 md.slr.love_h = love_numbers('h','CM'); md.slr.love_h(nlov+1:end)=[];
@@ -27,6 +27,6 @@
 %}}}
 %mask:  {{{
-md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset 
-mask=gmtmask(md.mesh.lat,md.mesh.long); 
+md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
+mask=gmtmask(md.mesh.lat,md.mesh.long);
 
 icemask=ones(md.mesh.numberofvertices,1);
@@ -39,5 +39,5 @@
 %make sure that the ice level set is all inclusive:
 md.mask.land_levelset=zeros(md.mesh.numberofvertices,1);
-md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
+md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1);
 
 %make sure that the elements that have loads are fully grounded:
@@ -45,10 +45,10 @@
 md.mask.groundedice_levelset(md.mesh.elements(pos,:))=1;
 
-%make sure wherever there is an ice load, that the mask is set to ice: 
-pos=find(md.slr.deltathickness); 
+%make sure wherever there is an ice load, that the mask is set to ice:
+pos=find(md.slr.deltathickness);
 md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
 % }}}
 
-md.slr.ocean_area_scaling=0; 
+md.slr.ocean_area_scaling=0;
 
 %geometry
@@ -68,4 +68,5 @@
 md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
 md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
+md.slr.hydro_rate = zeros(md.mesh.numberofvertices,1);
 
 %Miscellaneous
@@ -78,22 +79,22 @@
 
 % max number of iteration reverted back to 10 (i.e., the original default value)
-md.slr.maxiter=10; 
+md.slr.maxiter=10;
 
-%eustatic run: 
+%eustatic run:
 md.slr.rigid=0; md.slr.elastic=0;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Seustatic=md.results.SealevelriseSolution.Sealevel;
 
-%eustatic + rigid run: 
+%eustatic + rigid run:
 md.slr.rigid=1; md.slr.elastic=0;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Srigid=md.results.SealevelriseSolution.Sealevel;
 
-%eustatic + rigid + elastic run: 
+%eustatic + rigid + elastic run:
 md.slr.rigid=1; md.slr.elastic=1;md.slr.rotation=0;
 md=solve(md,'Sealevelrise');
 Selastic=md.results.SealevelriseSolution.Sealevel;
 
-%eustatic + rigid + elastic + rotation run: 
+%eustatic + rigid + elastic + rotation run:
 md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1;
 md=solve(md,'Sealevelrise');
Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24284)
@@ -79,4 +79,5 @@
 md.slr.Ngia = np.zeros((md.mesh.numberofvertices, ))
 md.slr.Ugia = np.zeros((md.mesh.numberofvertices, ))
+md.slr.hydro_rate = np.zeros((md.mesh.numberofvertices, ))
 
 #Miscellaneous
Index: /issm/trunk-jpl/test/NightlyRun/test2003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 24284)
@@ -1,6 +1,6 @@
 %Test Name: EarthSlr_rotationalFeedback
 
-%mesh earth: 
-md=model; 
+%mesh earth:
+md=model;
 md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',1000.); %500 km resolution mesh
 
@@ -16,5 +16,5 @@
 md.slr.deltathickness(pos)=-1;
 
-%elastic loading from love numbers: 
+%elastic loading from love numbers:
 nlov=1001;
 md.slr.love_h = love_numbers('h'); md.slr.love_h(nlov+1:end)=[];
@@ -24,6 +24,6 @@
 %}}}
 %mask:  {{{
-md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset 
-mask=gmtmask(md.mesh.lat,md.mesh.long); 
+md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
+mask=gmtmask(md.mesh.lat,md.mesh.long);
 
 icemask=ones(md.mesh.numberofvertices,1);
@@ -36,5 +36,5 @@
 %make sure that the ice level set is all inclusive:
 md.mask.land_levelset=zeros(md.mesh.numberofvertices,1);
-md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
+md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1);
 
 %make sure that the elements that have loads are fully grounded:
@@ -42,10 +42,10 @@
 md.mask.groundedice_levelset(md.mesh.elements(pos,:))=1;
 
-%make sure wherever there is an ice load, that the mask is set to ice: 
+%make sure wherever there is an ice load, that the mask is set to ice:
 pos=find(md.slr.deltathickness); md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
 % }}}
 
-% use model representation of ocea area (not the ture area) 
-md.slr.ocean_area_scaling = 0; 
+% use model representation of ocea area (not the ture area)
+md.slr.ocean_area_scaling = 0;
 
 %geometry
@@ -68,4 +68,5 @@
 md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
 md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
+md.slr.hydro_rate = zeros(md.mesh.numberofvertices,1);
 
 %Solution parameters
@@ -74,6 +75,6 @@
 md.slr.geodetic=1;
 
-%eustatic + rigid + elastic run: 
-md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=0; 
+%eustatic + rigid + elastic run:
+md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=0;
 md.cluster=generic('name',oshostname(),'np',3);
 %md.verbose=verbose('111111111');
@@ -81,5 +82,5 @@
 SnoRotation=md.results.SealevelriseSolution.Sealevel;
 
-%eustatic + rigid + elastic + rotation run: 
+%eustatic + rigid + elastic + rotation run:
 md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk-jpl/test/NightlyRun/test2003.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24284)
@@ -88,4 +88,5 @@
 md.slr.Ngia = np.zeros((md.mesh.numberofvertices, ))
 md.slr.Ugia = np.zeros((md.mesh.numberofvertices, ))
+md.slr.hydro_rate = np.zeros((md.mesh.numberofvertices, ))
 
 #Solution parameters
Index: /issm/trunk-jpl/test/NightlyRun/test2010.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.m	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.m	(revision 24284)
@@ -1,8 +1,8 @@
 %Test Name: MomentOfIntertia
 
-%mesh earth: 
-md=model; 
+%mesh earth:
+md=model;
 rad_e = 6.371012*10^3; % mean radius of Earth, km
-md.mesh=gmshplanet('radius',rad_e,'resolution',1000.0);  % km resolution 
+md.mesh=gmshplanet('radius',rad_e,'resolution',1000.0);  % km resolution
 
 %parameterize slr solution:
@@ -17,7 +17,7 @@
 md.slr.sealevel=zeros(md.mesh.numberofvertices,1);
 md.slr.steric_rate=zeros(md.mesh.numberofvertices,1);
-md.slr.ocean_area_scaling = 1; 
+md.slr.ocean_area_scaling = 1;
 
-%elastic loading from love numbers: 
+%elastic loading from love numbers:
 nlov=1001;
 md.slr.love_h = love_numbers('h'); md.slr.love_h(nlov+1:end)=[];
@@ -27,6 +27,6 @@
 %}}}
 %mask:  {{{
-md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset 
-mask=gmtmask(md.mesh.lat,md.mesh.long); 
+md.mask=maskpsl(); % use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
+mask=gmtmask(md.mesh.lat,md.mesh.long);
 
 icemask=ones(md.mesh.numberofvertices,1);
@@ -39,5 +39,5 @@
 %make sure that the ice level set is all inclusive:
 md.mask.land_levelset=zeros(md.mesh.numberofvertices,1);
-md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
+md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1);
 
 %make sure that the elements that have loads are fully grounded:
@@ -45,7 +45,7 @@
 md.mask.groundedice_levelset(md.mesh.elements(pos,:))=1;
 
-%make sure wherever there is an ice load, that the mask is set to ice: 
+%make sure wherever there is an ice load, that the mask is set to ice:
 md.mask.ice_levelset=ones(md.mesh.numberofvertices,1);
-pos=find(md.slr.deltathickness); 
+pos=find(md.slr.deltathickness);
 md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
 % }}}
@@ -75,22 +75,23 @@
 md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
 md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
+md.slr.hydro_rate = zeros(md.mesh.numberofvertices,1);
 
-%eustatic + rigid + elastic run: 
-md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1; 
+%eustatic + rigid + elastic run:
+md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1;
 md.cluster=generic('name',oshostname(),'np',3);
 
-% uncomment following 2 lines for 
+% uncomment following 2 lines for
 md=solve(md,'Sealevelrise');
 eus=md.results.SealevelriseSolution.SealevelRSLEustatic;
 slr=md.results.SealevelriseSolution.Sealevel;
-moixz=md.results.SealevelriseSolution.SealevelInertiaTensorXZ; 
-moiyz=md.results.SealevelriseSolution.SealevelInertiaTensorYZ; 
-moizz=md.results.SealevelriseSolution.SealevelInertiaTensorZZ; 
+moixz=md.results.SealevelriseSolution.SealevelInertiaTensorXZ;
+moiyz=md.results.SealevelriseSolution.SealevelInertiaTensorYZ;
+moizz=md.results.SealevelriseSolution.SealevelInertiaTensorZZ;
 
 % analytical moi => just checking FOR ICE only!!! {{{
-% ...have to mute ** slr induced MOI in Tria.cpp ** prior to the comparison 
-%rad_e = rad_e*1e3; % now in meters 
+% ...have to mute ** slr induced MOI in Tria.cpp ** prior to the comparison
+%rad_e = rad_e*1e3; % now in meters
 %areas=GetAreasSphericalTria(md.mesh.elements,md.mesh.lat,md.mesh.long,rad_e);
-%lat=late*pi/180; lon=longe*pi/180; 
+%lat=late*pi/180; lon=longe*pi/180;
 %moi_xz = sum(-md.materials.rho_freshwater.*md.slr.deltathickness.*areas.*rad_e^2.*sin(lat).*cos(lat).*cos(lon));
 %moi_yz = sum(-md.materials.rho_freshwater.*md.slr.deltathickness.*areas.*rad_e^2.*sin(lat).*cos(lat).*sin(lon));
Index: /issm/trunk-jpl/test/NightlyRun/test2010.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24283)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24284)
@@ -90,4 +90,5 @@
 md.slr.Ngia = np.zeros((md.mesh.numberofvertices, ))
 md.slr.Ugia = np.zeros((md.mesh.numberofvertices, ))
+md.slr.hydro_rate = np.zeros((md.mesh.numberofvertices, ))
 
 #eustatic + rigid + elastic run: