Index: /issm/trunk-jpl/test/NightlyRun/test2084.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2084.m	(revision 22400)
+++ /issm/trunk-jpl/test/NightlyRun/test2084.m	(revision 22401)
@@ -1,5 +1,9 @@
 
 %Test Name: GiaCaron. 
-%Forward Love number solution for a viscoelastic earth, model M3-L70-V01 from Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z., Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A. and King, M. A. (2011), A benchmark study for glacial isostatic adjustment codes. Geophysical Journal International, 185: 106–132. doi:10.1111/j.1365-246X.2011.04952.x
+%Forward Love number solution for a viscoelastic earth, model M3-L70-V01 from
+%Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z.,
+%Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A. and King, M. A.
+%(2011), A benchmark study for glacial isostatic adjustment codes. Geophysical
+%Journal International, 185: 106--132. doi:10.1111/j.1365-246X.2011.04952.x
 
 md=model();
@@ -9,15 +13,15 @@
 md.miscellaneous.name='FourierLoveTest';
 
-md.verbose=verbose('111111101');
+md.verbose=verbose('all');
 cst=365.25*24*3600*1000;
 
-	md.materials.numlayers=6;
-        md.materials.radius =  [10 1222.5 3.4800e+03   5.7010e+03   5.9510e+03   6.3010e+03   6.3710e+03]'*1e3;
-	md.materials.density=  [1.0750e4 1.0750e+04   4.9780e+03   3.8710e+03   3.4380e+03   3.0370e+03]';
-	md.materials.lame_mu=  [1e-5         0   2.2834e+00   1.0549e+00   7.0363e-01   5.0605e-01]'*1e11;
-	md.materials.viscosity=[0             0   2.0000e+00   1.0000e+00   1.0000e+00   1.0000e+25]'*1e21;
-	md.materials.lame_lambda=md.materials.lame_mu*0+5e14;
-	md.materials.issolid=[1 0 1 1 1 1]';
-	md.materials.isburgers=zeros(md.materials.numlayers,1);
+md.materials.numlayers=6;
+md.materials.radius =  [10 1222.5 3.4800e+03   5.7010e+03   5.9510e+03   6.3010e+03   6.3710e+03]'*1e3;
+md.materials.density=  [1.0750e4 1.0750e+04   4.9780e+03   3.8710e+03   3.4380e+03   3.0370e+03]';
+md.materials.lame_mu=  [1e-5         0   2.2834e+00   1.0549e+00   7.0363e-01   5.0605e-01]'*1e11;
+md.materials.viscosity=[0            0   2.0000e+00   1.0000e+00   1.0000e+00   1.0000e+25]'*1e21;
+md.materials.lame_lambda=md.materials.lame_mu*0+5e14;
+md.materials.issolid=[1 0 1 1 1 1]';
+md.materials.isburgers=zeros(md.materials.numlayers,1);
 
 md.love.allow_layer_deletion=1;
@@ -30,6 +34,4 @@
 
 md=solve(md,'lv');
-
-
 
 %Fields and tolerances to track changes
@@ -51,6 +53,6 @@
 
 md=solve(md,'lv');
+
 %Fields and tolerances to track changes
-%loading love numbers
 field_names     ={field_names{:},'LoveH_loading_realpart','LoveK_loading_realpart','LoveL_loading_realpart','LoveH_loading_imagpart','LoveK_loading_imagpart','LoveL_loading_imagpart'};
 field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
@@ -65,8 +67,9 @@
 
 md.love.forcing_type=9;
-md.love.sh_nmin=2
+md.love.sh_nmin=2;
 md.love.frequencies=([0 1e-3 1e-2 1e-1 1 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;
 md.love.nfreq=length(md.love.frequencies);
 md=solve(md,'lv');
+
 %tidal love numbers
 field_names     ={field_names{:},'LoveH_tidal_elastic','LoveK_tidal_elastic','LoveL_tidal_elastic','LoveH_tidal_realpart','LoveK_tidal_realpart','LoveL_tidal_realpart','LoveH_tidal_imagpart','LoveK_tidal_imagpart','LoveL_tidal_imagpart'};
@@ -84,7 +87,5 @@
 	};
 
-
-
- %Many layers PREM-based model
+%Many layers PREM-based model
 %data=load('../Data/PREM_500layers');
 %md.love.sh_nmin=1;
@@ -126,5 +127,4 @@
 %	};
 
-
 %Model VSS96 from Vermeersen, L.L.A., Sabadini, R. & Spada, G., 1996a. Analytical visco-elastic relaxation models, Geophys. Res. Lett., 23, 697–700.
 md.materials.radius=[10, 1222.5, 3480., 3600., 3630.5, 3700., 3900., 4000., 4200., 4300., 4500., 4600., 4800., 4900., 5100., 5200., 5400., 5500., 5600.5, 5650., 5701., 5736., 5771.5, 5821., 5951., 5970.5, 6016., 6061., 6150.5, 6151.5, 6251., 6371.]'*1e3;
@@ -145,5 +145,4 @@
 md.love.nfreq=length(md.love.frequencies);
 
-
 field_names     ={field_names{:},'LoveH_loadingVSS96_elastic','LoveK_loadingVSS96_elastic','LoveL_loadingVSS96_elastic','LoveH_loadingVSS96_realpart','LoveK_loadingVSS96_realpart','LoveL_loadingVSS96_realpart','LoveH_loadingVSS96_imagpart','LoveK_loadingVSS96_imagpart','LoveL_loadingVSS96_imagpart'};
 field_tolerances={field_tolerances{:},1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
Index: /issm/trunk-jpl/test/NightlyRun/test2084.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2084.py	(revision 22400)
+++ /issm/trunk-jpl/test/NightlyRun/test2084.py	(revision 22401)
@@ -19,5 +19,5 @@
 md.miscellaneous.name='FourierLoveTest'
 
-md.verbose=verbose('111111101')
+md.verbose=verbose('all')
 cst=365.25*24*3600*1000
 
@@ -45,6 +45,4 @@
 md=solve(md,'lv')
 
-
-
 #Fields and tolerances to track changes
 #loading love numbers
@@ -65,6 +63,4 @@
 
 md=solve(md,'lv')
-
-
 
 #Fields and tolerances to track changes
@@ -89,6 +85,4 @@
 md=solve(md,'lv')
 
-
-
 #tidal love numbers, check
 field_names+=['LoveHtr','LoveKtr','LoveLtr','LoveHti','LoveKti','LoveLti']
@@ -105,6 +99,4 @@
 np.array(md.results.LoveSolution.LoveLi)[:,1:]
 ]
-
-
 
 #Many layers PREM-based model
@@ -189,6 +181,4 @@
 md=solve(md,'lv')
 
-
-
 md.love.frequencies=(np.array([0,1e-3,1e-2,1,-1e-3,-1e-2,
 				-1])*2*pi).reshape(-1,1)/cst