C ++********************************************************************
C                                                                      *
C ROT3.F                                                                   *
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     ROT3(LUNI,LUNO,NSAM,NROW,NSLICE,AV1,IFANGL,ALPH,BETA,GAMA, 
C          IFSHIF,SHX,SHY,SHZ)
C
C     3-D ROTATION ACCORDING TO EULER ANGLES
C           
C     ANGLES ARE ACCORDING TO THE CONVENTION THAT THEY AGREE WITH WHAT 
C     YOU SEE IN THE SLICES
C     
C     IF IFANGL AND IFSHIF ARE EQUAL TO 1 THIS SOUBROUTINE DOES NOT ASK
C     FOR ANYTHING AND PROCREEDS WITH THE GIVEN DATA
C
C     (ALPH AND BETA ARE THE SAME THAN IN MICHAEL'S PROJECTION PROGRAM)
C     
C     10/27/86 J.M.CARAZO
C
C--*******************************************************************

      SUBROUTINE ROT3(LUNI,LUNO,NSAM,NROW,NSLICE,AV1,
     &            IFANGL,ALPH,BETA,GAMA, IFSHIF,SHX,SHY,SHZ)

      INCLUDE 'CMLIMIT.INC'
      COMMON /IOBUF/ RBUF(NBUFSIZ)
 
      INCLUDE 'CMBLOCK.INC' 

      CHARACTER NULL,ANS

      DATA PI/3.14159269/

      NULL = CHAR(0)

C     LUNI AND LUNO ARE MAIN INPUT/OUTPUT VOLUMES

C     ASK FOR TILTING ANGLES (IF IFANGL.NE.0)

      IF (IFANGL .NE. 1) THEN
         CALL RDPRM2(ALPH,BETA,NOT_USED,'ALPHA, BETA')
         CALL RDPRM (GAMA,NOT_USED,'GAMMA')
      ENDIF

      ALPH = ALPH*PI/180
      BETA = BETA*PI/180
      GAMA = GAMA*PI/180

      IF (IFSHIF .NE. 1) THEN
C        ASK FOR ROTATIONAL CENTER AND AVERAGE
         SHX = 0
         SHY = 0
         SHZ = 0
         CALL RDPRM2(SHX,SHY,NOT_USED,'SHIFT IN X-Y')

         CALL RDPRM (SHZ,NOT_USED,'SHIFT IN Z')
 
         CALL RDPRMC(ANS,NC,.TRUE.,
     &       'CHANGE DEFAULT MEAN VALUE? (N/Y)',NULL,IRT)

         IF (ANS .EQ. 'Y') THEN
            CALL RDPRM(AV1,NOT_USED,'VALUE OUTSIDE')
         ENDIF
      ENDIF

C     IT COPIES LUNI INTO LUNO (inefficient! al)

      IF (LUNI .NE. LUNO) THEN
        DO  I=1, NROW*NSLICE
           CALL REDLIN(LUNI,RBUF,NSAM,I)
           CALL WRTLIN(LUNO,RBUF,NSAM,I)  
	ENDDO
      ENDIF

C     MAKES THE FIRST ROTATION AROUND Z
      IF (ABS(ALPH) .GE. 0.0001) THEN
        DO  ISLICE=1, NSLICE  
            NROWS  = (ISLICE-1) * NROW + 1
            NROWE  = NROWS + NROW - 1
            NROWSK = -1
            CALL ROT32(LUNO,LUNO,NSAM,NROWS,NROWE,NROWSK,
     &                 ALPH,AV1,SHX,SHY)
         ENDDO
      ENDIF

C     could be speeded up by not writing to file each rotation!! al

C     MAKES THE ROTATION AROUND Y 
      IF (ABS(BETA) .GE. 0.0001) THEN  
         BETA=-BETA
         DO  IROW=1,NROW
            NROWS  = IROW         
            NROWE  = (NSLICE-1) * NROW + IROW
            NROWSK = -NROW
            CALL ROT32(LUNO,LUNO,NSAM,NROWS,NROWE,NROWSK,
     &                 BETA,AV1,SHX,SHZ)
         ENDDO
      ENDIF


C     THIS MAKES THE LAST ROTATION AROUND Z 
      IF (ABS(GAMA) .GE. 0.0001) THEN
         DO  ISLICE=1,NSLICE
            NROWS  = (ISLICE-1) * NROW + 1
            NROWE  = NROWS + NROW - 1
            NROWSK = -1
            CALL ROT32(LUNO,LUNO,NSAM,NROWS,NROWE,NROWSK,
     &                 GAMA,AV1,SHX,SHY)
         ENDDO
      ENDIF

      RETURN   
      END