C+
C NAME:
C	Veleq3
C PURPOSE:
C	Determines the average value at the mid latitude of a three-dimensional array. 
C
C CALLING SEQUENCE:
	subroutine 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	FEB-2001, B. Jackson (UCSD)
C	NOV-2007, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
C	    Reversed order of do-loops
C-
	real 		Z(nLng,nLat,nT)

	NL = nLat/2+1

	Bad = BadR4()

	VELEQU = 0.0
	AN = 0.0

!$omp parallel do private(i,j,k) reduction(+:VELEQU,AN)
	do k=1,nT
	    do j=NL,NL
		do i=1,nLng
		    if (Z(i,j,k) .ne. Bad) then 
			VELEQU = VELEQU+Z(i,j,k)
			AN = AN+1.0
		    end if
		end do
	    end do
	end do
!$omp end parallel do

	if (AN .ne. 0.0) VELEQU = VELEQU/AN

	return
	end