C+
C NAME:
C	Deneq3
C PURPOSE:
C	Determines the average value at the mid latitude of a three-dimensional array
C	from the density at the source surface and assuming the radial fall-off given
C	in the density fall-off model.  Note - this is only used when mv^2 is assumed 
C	constant in order to obtain a velocity map.
C
C CALLING SEQUENCE:
C	call Deneq3(nLng,nLat,nT,Z,RR,FALLOFF,DENEQU) 
C INPUTS:
C	nLng		integer		# longitudes
C	nLat		integer		# latitudes
C	nT		integer		# times
C	Z(nLng,nLat,nT)	real		3D array of function values
C	RR		real		Height of source surface in AU
C	FALLOFF		real		Radial falloff of density
C OUTPUTS:
C	DENEQU		real		Average density at 1 AU at the solar equator
C PROCEDURE:
C MODIFICATION HISTORY:
C	February, 2001, B. Jackson (UCSD)
C-
	subroutine Deneq3(nLng,nLat,nT,Z,RR,FALLOFF,DENEQU)
	real 		Z(nLng,nLat,nT)
	NL = nLat/2 + 1
	Bad = BadR4()
	R1AU = 1.0
	DENEQU = 0.
	AN = 0.
	do i=1,nLng
	    do j=NL-1,NL+1
		do k=1,nT
		    if(Z(i,j,k).ne.Bad) then 
			DENEQU = DENEQU + Z(i,j,k)
			AN = AN + 1.
		    end if
		end do
	    end do
	end do
	if(AN.ne.0.) DENEQU = (DENEQU/AN)*((RR/R1AU)**FALLOFF)		! Average equatorial density at 1AU
	return
	end