C++*********************************************************************
C
C VOEA.F              USED REG_SET_NSEL  AUG 00 ARDEAN LEITH
C
C **********************************************************************
C=* FROM: SPIDER - MODULAR IMAGE PROCESSING SYSTEM.   AUTHOR: J.FRANK  *
C=* Copyright (C) 1985-2005  Health Research Inc.                      *
C=*                                                                    *
C=* HEALTH RESEARCH INCORPORATED (HRI),                                *   
C=* ONE UNIVERSITY PLACE, RENSSELAER, NY 12144-3455.                   *
C=*                                                                    *
C=* Email:  spider@wadsworth.org                                       *
C=*                                                                    *
C=* This program is free software; you can redistribute it and/or      *
C=* modify it under the terms of the GNU General Public License as     *
C=* published by the Free Software Foundation; either version 2 of the *
C=* License, or (at your option) any later version.                    *
C=*                                                                    *
C=* This program is distributed in the hope that it will be useful,    *
C=* but WITHOUT ANY WARRANTY; without even the implied warranty of     *
C=* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU  *
C=* General Public License for more details.                           *
C=*                                                                    *
C=* You should have received a copy of the GNU General Public License  *
C=* along with this program; if not, write to the                      *
C=* Free Software Foundation, Inc.,                                    *
C=* 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.      *
C=*                                                                    *
C **********************************************************************
C
C  PREPARE ANG DOC FILE FOR PJ 3Q (QUASI-EQUAL ANGULAR SPACING) 04/23/92
C  MODIFIED ON 10/23/96 TO GET EVEN SPACING IN PHI.
C
C
C23456789012345678901234567890123456789012345678901234567890123456789012
C--*********************************************************************

         SUBROUTINE VOEA(MAXMEM)

         INCLUDE 'CMBLOCK.INC'
         COMMON        DUMMY(1024),FIO(4)
         CHARACTER*80  FILNAM,FINPAT
         EQUIVALENCE   (FIO(2),PSI),(FIO(3),THETA),(FIO(4),PHI)
	 LOGICAL       SKIP
	 PARAMETER (QUADPI = 3.141592653589793238462643383279502884197)
	 PARAMETER (DGR_TO_RAD = (QUADPI/180))

         DATA  NDOUT/56/
#ifdef USE_MPI
         INCLUDE 'mpif.h'
         INTEGER MYPID, COMM, IERR
         COMM = MPI_COMM_WORLD
         CALL MPI_COMM_RANK(COMM, MYPID, IERR)
#endif
         PSI = 0.0

         CALL  RDPRM(DELTA,NOT_USED,'DELTA THETA')
         CALL  RDPRM2(T1,T2,NOT_USED,'RANGE OF THETA (0,90)')
         CALL  RDPRM2(P1,P2,NOT_USED,'RANGE OF PHI (0,359.9)')
         IF(T1.EQ.0.0.AND.T2.EQ.0.0.OR.T1.GE.T2) THEN
            T1=0.0
            T2=90.0
         ENDIF
         IF(P1.EQ.0.0.AND.P2.EQ.0.0.OR.P1.GE.P2) THEN
            P1=0.0
            P2=359.9
         ENDIF

         SKIP=T1.LT.90.0.AND.T2.EQ.90.0.AND.P1.EQ.0.0.AND.P2.GT.180.0
         LITER=0
         NLIST=4
C        do  1  theta=0.,90.0,delta
         DO  1  THETA=T1,T2,DELTA
            IF (THETA.EQ.0.0.OR.THETA.EQ.180.0)  THEN
		DETPHI=360.0
		LT=1
            ELSE
		DETPHI=DELTA/SIN(THETA*DGR_TO_RAD)
		LT=MAX0(INT((P2-P1)/DETPHI)-1,1)
		DETPHI=(P2-P1)/LT
            ENDIF
C           do  1  phi=0.,359.9,detphi
C           DO  1  PHI=P1,P2,DETPHI
	       DO  1  I=1,LT
	         PHI=P1+(I-1)*DETPHI
	         IF (SKIP.AND.THETA.EQ.90.0.AND.PHI.GE.180.0)  GOTO  1
                 LITER=LITER+1
                 FIO(1)=LITER
                 CALL  SAVD(NDOUT,FIO,NLIST,IRTFLG)
                 IF (IRTFLG.EQ.-1) GOTO  5
1        CONTINUE
#ifdef USE_MPI
         IF (MYPID .EQ. 0) THEN
            WRITE(NOUT,33)  LITER
         ENDIF
#else
         WRITE(NOUT,33)  LITER
#endif
33       FORMAT('   Total number of directions =',i10)
5        CALL  SAVDC
#ifdef USE_MPI
         IF (MYPID .EQ. 0) THEN
            CLOSE(NDOUT)
            CALL FLUSHRESULTS
         ENDIF
c
#else
         CLOSE(NDOUT)
#endif

C        IF (NSEL(1).NE.0)  PARAM(NSEL(1))=LITER
         CALL REG_SET_NSEL(1,1,FLOAT(LITER),0.0,0.0,0.0,0.0,IRTFLG)

         RETURN
         END