C ********************************************************************** C C GOTBRANCH.FOR -- CREATED OCT 90 C ********************************************************************** C * AUTHOR: ArDean Leith 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 GOTBRANCH(NCLUS,NOW,IEQUIV,NEQUIV,NEQMAX,IRTFLG) C C PURPOSE: SETS BRANCH CONNECTIVITY ARRAY C C PARAMETERS: C C CALLED BY: CCONECT C C CALLS: C C23456789012345678901234567890123456789012345678901234567890123456789012 C--******************************************************************** SUBROUTINE GOTBRANCH(NCLUS,NOW,IEQUIV,NEQUIV,NEQMAX,IRTFLG) COMMON /UNITS/LUN,NIN,NOUT INTEGER IEQUIV(2,*) NOWMIN = MIN(NCLUS,NOW) NOWMAX = MAX(NCLUS,NOW) DO I = 1,NEQUIV C CHECK IF ALREADY HAVE RECORDED THIS Y IN NEQUIV IF (IEQUIV(2,I) .EQ. NOWMIN .AND. & IEQUIV(1,I) .EQ. NOWMAX) RETURN END DO c********************************** C** WRITE(NOUT,*) ' NOW,NCLUS: ',NOW,NCLUS c********************************** NEQUIV = NEQUIV + 1 IERR = 0 IF (NEQUIV .GT. NEQMAX) THEN C RAN OUT OF SPACE IN IEQUIV ARRAY IF (IERR .LE. 1) THEN WRITE(NOUT,*) '*** BRANCH LIMIT (NEQMAX) IS:', NEQMAX WRITE(NOUT,*) ' OUTPUT IS JUNK !!!!!!!!!!!!!!!' IERR = IERR + 1 ENDIF IRTFLG = 1 RETURN ENDIF IEQUIV(1,NEQUIV) = NOWMAX IEQUIV(2,NEQUIV) = NOWMIN IRTFLG = 0 RETURN END