C+
C NAME:
C	Veleq3
C PURPOSE:
C	Determines the average value at the mid latitude of a three-dimensional array.
C	Note - this program is only used to deconvolve velocities using mv^2 = constant
C	and will very seldom be used, since almost never will ipshtd be used to only
C	deconvolve velocities using mv^2 = constant.
C
C CALLING SEQUENCE:
C	call Veleq3(nLng,nLat,nT,Z,VELEQU)
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 OUTPUTS:
C	VELEQU		real		Average density at 1 AU at the solar equator
C PROCEDURE:
C MODIFICATION HISTORY:
C	February, 2001, B. Jackson (UCSD)
C-
	subroutine Veleq3(nLng,nLat,nT,Z,VELEQU)
	real 		Z(nLng,nLat,nT)
	NL = nLat/2 + 1
	Bad = BadR4()
	VELEQU = 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
			VELEQU = VELEQU + Z(i,j,k)
			AN = AN + 1.
		    end if
		end do
	    end do
	end do
	if(AN.ne.0.) VELEQU = VELEQU/AN
	return
	end