Index: /issm/trunk-jpl/examples/SlrGRACE/runme.m
===================================================================
--- /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 26879)
+++ /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 26880)
@@ -2,5 +2,5 @@
 addpath('../Data','../Functions');
 
-steps=[1:7];
+steps=[7:7];
 
 if any(steps==1) % {{{ 
Index: /issm/trunk-jpl/src/c/modules/FourierLoveCorex/love_numbers.f90
===================================================================
--- /issm/trunk-jpl/src/c/modules/FourierLoveCorex/love_numbers.f90	(revision 26879)
+++ /issm/trunk-jpl/src/c/modules/FourierLoveCorex/love_numbers.f90	(revision 26880)
@@ -308,5 +308,5 @@
    if (allow_layer_del.eqv..true.) then
 
-      print*, layerrap, epsdb, nbc 
+      !print*, layerrap, epsdb, nbc 
       !do while ((layerrap<=epsdb).and.(nbc>12*6).or.(isnan(layerrap)) )
       do while ((layerrap<=epsdb).and.(nbc>12).or.(isnan(layerrap)) )
Index: /issm/trunk-jpl/src/c/modules/FourierLoveCorex/lovenb_sub.f90
===================================================================
--- /issm/trunk-jpl/src/c/modules/FourierLoveCorex/lovenb_sub.f90	(revision 26879)
+++ /issm/trunk-jpl/src/c/modules/FourierLoveCorex/lovenb_sub.f90	(revision 26880)
@@ -48,4 +48,5 @@
    endif
 
+	valini=1.d0;
   ! Which among the model interfaces is the first to display boundary conditions
   ifirst = 0
@@ -65,5 +66,5 @@
    xmin= radius(i)/ra
    xmax= radius(i+1)/ra - epsdb
-   nstep=100
+   nstep=1
   
    if (soliddim(i)) then
@@ -193,5 +194,4 @@
 	!write(89,*),xsen,r,ra,cst,mu0,mu,n	
 
-
 	bc(1,1)=xsen*r/ra
 	bc(1,2)=xsen/(r*ra)
@@ -334,9 +334,12 @@
 	if (radbc(i)) ifirst = i
   end do
-! do i = 1,nbc 
-!	do j = 1,nbc
-!	if (bc(i,j)/=0.d0) print*, i,j,bc(i,j)
-!	enddo
-! enddo
+
+ open(unit=200, file="/home/caronlam/Codes/trunk-exp/trunk-jpl/test/NightlyRun/BC_fourier", status='unknown')
+ do i = 1,nbc 
+ 	do j = 1,nbc
+	!if (bc(i,j)/=0.d0) 
+ 	write(200,*), deg,nbc,i,j, dble(f(11,j)) ,dble(bc(i,j))
+ 	enddo
+ enddo
 ! write(*,*)
   
@@ -375,4 +378,5 @@
 	lovel = sc(nbc-1,1)*ra*go_surf
 	lovek = sc(nbc,1)*ra*go0
+	!print*,loveh,lovel, lovek
 	delta = (1.d0-dble(n+1)/dble(n)*(lovek-1.d0)+2.d0/dble(n)*loveh)
 
@@ -382,5 +386,7 @@
       if (j<nbc-2) then ! no need for valini restoration for surface. 
          sc(j,1) = sc(j,1)*valini**((nbc-j-3)/6+1)
+	
       endif
+	!print*,ipiv(j),dble(sc(j,1))
 	   love_kernels(j,1)=sc(j,1) !*valini
 	end do
Index: /issm/trunk-jpl/test/NightlyRun/test2071.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2071.m	(revision 26879)
+++ /issm/trunk-jpl/test/NightlyRun/test2071.m	(revision 26880)
@@ -104,4 +104,26 @@
 field_values={h,k,l};
 
+return
+subplot(3,1,1)
+pcolor(0:256,t,log10(abs((hts-h(:,1:257))./hts)))       
+set(gca,'yscale', 'log')
+shading flat
+colorbar 
+title('h')
+
+subplot(3,1,2)
+pcolor(0:256,t,log10(abs((kts-k(:,1:257))./kts)))       
+set(gca,'yscale', 'log')
+shading flat
+colorbar 
+title('k')
+
+subplot(3,1,3)
+pcolor(0:256,t,log10(abs((lts-l(:,1:257))./lts)))       
+set(gca,'yscale', 'log')
+shading flat
+colorbar 
+title('l')
+
 
 return
Index: /issm/trunk-jpl/test/NightlyRun/test2092.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2092.m	(revision 26879)
+++ /issm/trunk-jpl/test/NightlyRun/test2092.m	(revision 26880)
@@ -22,5 +22,5 @@
 
 defaultoptions={'KeepVertices',0,'MaxCornerAngle',0.0000000001,'NoBoundaryRefinement',1}; 
-md2d=bamg(model,'domain',dom,'subdomains',co,'hmin',50e3,'hmax',10000e3,defaultoptions{:});
+md2d=bamg(model,'domain',dom,'subdomains',co,'hmin',100e3,'hmax',10000e3,defaultoptions{:});
 
 colat=sqrt(md2d.mesh.x.^2+md2d.mesh.y.^2)/re;