c----------------------------------------------------------------------
	function DGINTGRL1(n,r0,f,x)
c integrates a function from 0 to x using two successive grids: (y=1/(1+r/r0) 
c to scale into [1,r0/(r0+x)] and a uniform grid to scale to [0,1]
c
c author: Kieron Burke
c     in: n  - number of gridpoints
c         r0 - parameter for y grid
c         f  - array with function values
c         x  - upper bound for integral
c    out: return value of function: integral value of f from 0 to x
c  calls: STDNTGRL
	implicit real*8(a-h,o-z)
	parameter(m=1000000)
	dimension f(0:n)
	dimension dum(0:m)