C++*********************************************************************
C
C CRCSE3.FOR
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 IMAGE_PROCESSING_ROUTINE
C
C        1         2         3         4         5         6         7
C23456789012345678901234567890123456789012345678901234567890123456789012
C--*********************************************************************

         SUBROUTINE CRCSE3(LUN,LUN2,NSAM,NROW,NSLICE,IR)



         DIMENSION BUF(NSAM),SEC(IR),SNO(IR)
     
            SEC = 0.0
            SNO = 0.0
   

         DO K=1,NSLICE
            KK=K-NSLICE/2-1
            IF(IABS(KK).LE.IR-1) THEN
               DO J=1,NROW
                  KJ=J-NROW/2-1
                  NR1=J+(K-1)*NROW
                  IF(IABS(KJ).LE.IR-1) THEN
                     CALL REDLIN(LUN,BUF,NSAM,NR1)
                     DO I=1,NSAM
                        KI=I-NSAM/2-1
         R=SQRT(FLOAT(KJ*KJ)+FLOAT(KI*KI)+FLOAT(KK*KK))+1.0
                        L=R
                        IF(L.LE.IR-1) THEN
                           XD=R-L
                           SEC(L)=SEC(L)+BUF(I)*(1.0-XD)
                           SEC(L+1)=SEC(L+1)+BUF(I)*XD
                           SNO(L)=SNO(L)+1.0-XD
                           SNO(L+1)=SNO(L+1)+XD
                        ENDIF
                     ENDDO             
                  ENDIF
               ENDDO
            ENDIF
         ENDDO        

         DO    I=1,IR
            SEC(I)=SEC(I)/AMAX1(1.0,SNO(I))
         ENDDO

         CALL WRTLIN(LUN2,SEC,IR,1)

         END