*$ CREATE MGDRAW.FOR *COPY MGDRAW * * *=== mgdraw ===========================================================* * * SUBROUTINE MGDRAW ( ICODE, MREG ) INCLUDE '(DBLPRC)' INCLUDE '(DIMPAR)' INCLUDE '(IOUNIT)' * *----------------------------------------------------------------------* * * * Copyright (C) 1990-2006 by Alfredo Ferrari * * All Rights Reserved. * * * * * * MaGnetic field trajectory DRAWing: actually this entry manages * * all trajectory dumping for * * drawing * * * * Created on 01 march 1990 by Alfredo Ferrari * * INFN - Milan * * Last change 05-may-06 by Alfredo Ferrari * * INFN - Milan * * * *----------------------------------------------------------------------* * INCLUDE '(CASLIM)' INCLUDE '(COMPUT)' INCLUDE '(SOURCM)' INCLUDE '(FHEAVY)' INCLUDE '(FLKSTK)' INCLUDE '(GENSTK)' INCLUDE '(MGDDCM)' INCLUDE '(PAPROP)' INCLUDE '(QUEMGD)' INCLUDE '(SUMCOU)' INCLUDE '(TRACKR)' * DIMENSION DTQUEN ( MXTRCK, MAXQMG ) * CHARACTER*20 FILNAM CHARACTER*8 MRGNAM, NRGNAM LOGICAL LFCOPE SAVE LFCOPE DATA LFCOPE / .FALSE. / * RETURN * *======================================================================* * * ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO ) IF ( .NOT. LFCOPE ) THEN LFCOPE = .TRUE. OPEN ( UNIT = 88, FILE = "collfile", STATUS = 'UNKNOWN') END IF CALL GEOR2N ( MREG, MRGNAM, IERR1 ) CALL GEOR2N ( NEWREG, NRGNAM, IERR2 ) IF(IERR1 .NE. 0 .OR. IERR2 .NE. 0) STOP "Error in name conversion" IF(MRGNAM .EQ. "VpNEHaVa" .AND. nrgnam .EQ. "DummyC4d") THEN IF(JTRACK .EQ. 7) WRITE(88,'(1P,4G25.15)') & XSCO, YSCO, ETRACK, WTRACK END IF RETURN * *======================================================================* * * ENTRY EEDRAW ( ICODE ) RETURN * * *======================================================================* * * ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO ) RETURN * *======================================================================* * * ENTRY SODRAW RETURN * *======================================================================* * * ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO ) RETURN *=== End of subrutine Mgdraw ==========================================* END