SUBROUTINE sprsin_dp(a,thresh,sa)
	USE nrtype; USE nrutil, ONLY : arth,assert_eq
	IMPLICIT NONE
	REAL(dp), DIMENSION(:,:), INTENT(IN) :: a
	REAL(dp), INTENT(IN) :: thresh
	TYPE(sprs2_dp), INTENT(OUT) :: sa
	INTEGER(I4B) :: n,len
  integer :: arth_temp(size(a,1), size(a,2))  ! Added to get around ifort bug. DCS
	LOGICAL(LGT), DIMENSION(size(a,1),size(a,2)) :: mask
	n=assert_eq(size(a,1),size(a,2),'sprsin_dp')
	mask=abs(a)>thresh
	len=count(mask)
	allocate(sa%val(len),sa%irow(len),sa%jcol(len))
	sa%n=n
	sa%len=len
	sa%val=pack(a,mask)

  arth_temp = spread(arth(1,1,n),2,n)
	sa%irow=pack(arth_temp,mask)

  arth_temp = spread(arth(1,1,n),1,n)
	sa%jcol=pack(arth_temp,mask)

	END SUBROUTINE sprsin_dp