| 1 |
cafagna |
1.1 |
* |
| 2 |
|
|
* $Id: gpustep.F,v 3.1.1.1 2002/07/11 16:02:00 cafagna Exp $ |
| 3 |
|
|
* |
| 4 |
|
|
* $Log: gpustep.F,v $ |
| 5 |
|
|
* Revision 3.1.1.1 2002/07/11 16:02:00 cafagna |
| 6 |
|
|
* First GPAMELA release on CVS |
| 7 |
|
|
* |
| 8 |
|
|
* |
| 9 |
|
|
*CMZ : 3.00/00 11/05/2001 13.44.09 by Marialuigia Ambriola |
| 10 |
|
|
*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola |
| 11 |
|
|
*CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna |
| 12 |
|
|
*CMZ : 1.01/00 30/04/96 18.28.33 by Francesco Cafagna |
| 13 |
|
|
*CMZ : 1.00/03 29/04/96 16.56.08 by Francesco Cafagna |
| 14 |
|
|
*-- Author : Francesco Cafagna 29/04/96 |
| 15 |
|
|
SUBROUTINE GPGIG |
| 16 |
|
|
************************************************************************ |
| 17 |
|
|
* * |
| 18 |
|
|
* User action in GUSTEP * |
| 19 |
|
|
* * |
| 20 |
|
|
* Variables definition: * |
| 21 |
|
|
* IN: * |
| 22 |
|
|
* IFLAG = 1, just entered GUSTEP * |
| 23 |
|
|
* IFLAG = 2, leaving GUSTEP * |
| 24 |
|
|
* * |
| 25 |
|
|
* Called by: GUSTEP * |
| 26 |
|
|
* Author: Francesco Cafagna, 29/04/96 16.56.08 * |
| 27 |
|
|
* * |
| 28 |
|
|
************************************************************************ |
| 29 |
|
|
#include "gpques.inc" |
| 30 |
|
|
#include "gctmed.inc" |
| 31 |
|
|
#include "gckine.inc" |
| 32 |
|
|
#include "gcking.inc" |
| 33 |
|
|
#include "gcflag.inc" |
| 34 |
|
|
#include "gctrak.inc" |
| 35 |
|
|
#include "gcvolu.inc" |
| 36 |
|
|
#include "gcsets.inc" |
| 37 |
|
|
#include "gpkey.inc" |
| 38 |
|
|
#include "gpsed.inc" |
| 39 |
|
|
|
| 40 |
|
|
* |
| 41 |
|
|
LOGICAL W,GPLOOK |
| 42 |
|
|
INTEGER IFLAG |
| 43 |
|
|
INTEGER IFL,IN,MECNAM(MAXMEC) |
| 44 |
|
|
IF(IFLAG.EQ.1) GO TO 999 |
| 45 |
|
|
IN = INWVOL |
| 46 |
|
|
IF(IN.NE.0) GO TO 999 |
| 47 |
|
|
C # open(20,file='GIG.dat',access='append') |
| 48 |
|
|
C # open(21,file='All.dat',access='append') |
| 49 |
|
|
* |
| 50 |
|
|
* Is it a gamma ? |
| 51 |
|
|
* |
| 52 |
|
|
IF ((IPART.NE.1).AND.(IPART.NE.13)) GO TO 999 |
| 53 |
|
|
IF (IPART.EQ.13) write(*,*)'particle=',IPART |
| 54 |
|
|
* |
| 55 |
|
|
* Store the mechanisms active for the current step |
| 56 |
|
|
* |
| 57 |
|
|
IF(NMEC.EQ.0)THEN |
| 58 |
|
|
MECNAM(1)=NAMEC(29) |
| 59 |
|
|
ELSE |
| 60 |
|
|
DO I=1,NMEC |
| 61 |
|
|
MEC=LMEC(I) |
| 62 |
|
|
IF(MEC.LE.MAXMEC) THEN |
| 63 |
|
|
MECNAM(I)=NAMEC(MEC) |
| 64 |
|
|
ELSEIF(MEC-100.LE.MAXME1.AND.MEC-100.GT.0) THEN |
| 65 |
|
|
MECNAM(I)=NAMEC1(MEC-100) |
| 66 |
|
|
ENDIF |
| 67 |
|
|
ENDDO |
| 68 |
|
|
ENDIF |
| 69 |
|
|
* |
| 70 |
|
|
W=GPLOOK('CAAB',NAMES,NLEVEL) |
| 71 |
|
|
IF(W) THEN |
| 72 |
|
|
c CALL GPCXYZ |
| 73 |
|
|
IF ((GETOT.GT.0.009097).AND.(GETOT.LT.0.028608).AND.(IPART.EQ.1)) |
| 74 |
|
|
+ then |
| 75 |
|
|
c PRINT *, ' GPUSTEP : WE FOUND A GAMMA INTO THE W ! Energy:' |
| 76 |
|
|
c + ,GETOT |
| 77 |
|
|
C # write(20,'(1x,f7.3,1x,f6.3,1x,i4,1x,a10,1x,i4,5(1x,e13.6),1x,i4)') |
| 78 |
|
|
C # + GETOT*1000., STEP, ISTAK, KCASE, NGKINE, |
| 79 |
|
|
C # + VECT(4)*VECT(7)*1000., VECT(5)*VECT(7)*1000., |
| 80 |
|
|
C # + VECT(6)*VECT(7)*1000., VECT(7)*1000., VECT(3), NUMBER(NLEVEL) |
| 81 |
|
|
ENDIF |
| 82 |
|
|
C # IF (IPART.EQ.13) then |
| 83 |
|
|
C # write(21,'(1x,f7.3,1x,f6.3,1x,i4,1x,a10,1x,i4,5(1x,e13.6),1x,i4)') |
| 84 |
|
|
C # + GETOT*1000., STEP, ISTAK, KCASE, NGKINE, |
| 85 |
|
|
C # + VECT(4)*VECT(7)*1000., VECT(5)*VECT(7)*1000., |
| 86 |
|
|
C # + VECT(6)*VECT(7)*1000., VECT(7)*1000., VECT(3), NUMBER(NLEVEL) |
| 87 |
|
|
C # ENDIF |
| 88 |
|
|
ENDIF |
| 89 |
|
|
999 CONTINUE |
| 90 |
|
|
C # close(20) |
| 91 |
|
|
C # close(21) |
| 92 |
|
|
|
| 93 |
|
|
RETURN |
| 94 |
|
|
END |
| 95 |
|
|
|
Parent Directory
| 
Revision Log