test size h

Author: perazz

src/fitpack_core.f90

@@ -12232,7 +12232,7 @@ pure subroutine fpspgr(iopt,ider,u,mu,v,mv,r,mr,r0,r1,s, &
       end subroutine fpspgr
 
 
-      pure subroutine fpsphe(iopt,m,teta,phi,r,w,s,ntest,npest,eta,tol,maxit, &
+      subroutine fpsphe(iopt,m,teta,phi,r,w,s,ntest,npest,eta,tol,maxit, &
                              ib1,ib3,nc,ncc,intest,nrest,nt,tt,np,tp,c,fp,sup,fpint,coord,f, &
                              ff,row,coco,cosi,a,q,bt,bp,spt,spp,h,index,nummer,wrk,lwrk,ier)
 
@@ -12270,6 +12270,8 @@ pure subroutine fpsphe(iopt,m,teta,phi,r,w,s,ntest,npest,eta,tol,maxit, &
       lwest  = 0
       ntt    = 0
       iband1 = 0
+      
+      print *, 'ib3=',ib3,' size(h)=',size(h)
 
       bootstrap: if (iopt>=0) then
 
@@ -17233,7 +17235,7 @@ pure subroutine spgrid(iopt,ider,mu,u,mv,v,r,r0,r1,s, &
       end subroutine spgrid
 
 
-      pure subroutine sphere(iopt,m,teta,phi,r,w,s,ntest,npest, &
+      subroutine sphere(iopt,m,teta,phi,r,w,s,ntest,npest, &
                              eps,nt,tt,np,tp,c,fp,wrk1,lwrk1,wrk2,lwrk2,iwrk,kwrk,ier)
 
       !  subroutine sphere determines a smooth bicubic spherical spline