Dear Julien,
I glanced over your usrrnc.f. You open unit 81 in the first call, and
then write the output in unit 80 (which is the file where your
mgdraw_ju.f writes to as well).
Cheers,
Cesc
On Thu, Jul 25 2019, at 12:56 +0200, Julien Bonvalet wrote:
>
>Dear fluka users,
>
>I try to follow the production of radioisotopes in my target with mgdraw.f and usrrnc.f routines.
>But I have a problem with the usrrnc routine because my file declared in it is always empty at the end of the simualtion. I am interesting by all the nuclei produced and their position in my target.
>
>Someone can see the issue in my routine?
>
>Thank you very much
>
>Julien BONVALET
>CELIA - BORDEAUX
>*$ 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 '(EMFSTK)'
> 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. /
>*
> LOGICAL UDFIRST
> DATA UDFIRST / .TRUE. /
> SAVE UDFIRST
>*
> RETURN
>*
>*======================================================================*
>* *
> ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
>
> RETURN
>*
>*======================================================================*
>* *
> ENTRY EEDRAW ( ICODE )
> RETURN
>* *
>*======================================================================*
>* *
> ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
> RETURN
>*
>*======================================================================*
>* *
> ENTRY SODRAW
> RETURN
>*
>*======================================================================*
>* *
> ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
> IF ( UDFIRST ) THEN
> UDFIRST = .FALSE.
> OPEN ( UNIT = 80, FILE = "filejuju_usd", STATUS = 'UNKNOWN')
> ENDIF
> IF (NPHEAV/=0) THEN
> WRITE(80,'(1P,12G25.15)') 'c1 ',NCASE,JTRACK,NPFLKA,ICODE,
> & NP,NPHEAV,ETRACK-AM(JTRACK),XSCO,YSCO,ZSCO,WTRACK
> ENDIF
>*
>* Boucle sur particules secondaires
>* do 10 ip = 1, NP
>* IF (ZPART/=0) THEN
>* WRITE( 80,'(1P,12G25.15)')'c2 ',NCASE,JTRACK,
>* & ip,ICODE,ZPART,APART,TKI(ip)
>* ,CXR(ip),CYR(ip),
>* & CZR(ip),WEI(ip)
>* ENDIF
>*10 continue
>*
> do 20 ip=1,NPHEAV
> WRITE(80,'(1P,12G25.15)') 'c4', NCASE, JTRACK,
> & ip,ICODE,KHEAVY(ip),ICHEAV(KHEAVY(ip)),
> & IBHEAV(KHEAVY(ip)),TKHEAV(ip)
>20 continue
> DO 1 IZ = 1, MAXZ
> DO 2 NMZMK = MINNMZ-K, MAXNMZ-K
> WRITE(80,'(2I4,1P,G15.6)') 'cres',IZ, NMZMK+K+2*IZ,
> & ARRAY(IZ,NMZMK)
>2 CONTINUE
>1 CONTINUE
> RETURN
>*=== End of subrutine Mgdraw ==========================================*
> END
>*$ CREATE USRRNC.FOR
>*COPY USRRNC
>*
>*=== Usrrnc ===========================================================*
>*
> SUBROUTINE USRRNC ( IZ, IA, IS, X, Y, Z, MREG, WEE, ICALL )
>
> INCLUDE '(DBLPRC)'
> INCLUDE '(DIMPAR)'
> INCLUDE '(IOUNIT)'
>*
>*----------------------------------------------------------------------*
>* *
>* Copyright (C) 2005-2005 by Alfredo Ferrari & Paola Sala *
>* All Rights Reserved. *
>* *
>* *
>* USeR Residual NuClei: *
>* *
>* Created on 06 april 2005 by Alfredo Ferrari & Paola Sala *
>* Infn - Milan *
>* *
>* Last change on 06-apr-05 by Alfredo Ferrari *
>* *
>* *
>*----------------------------------------------------------------------*
>*
> LOGICAL LFIRST
> SAVE LFIRST
> DATA LFIRST /.TRUE./
>*
> IF (LFIRST) THEN
> WRITE(LUNOUT,*) 'Version juju usrrnc called'
> LFIRST = .FALSE.
> OPEN ( UNIT = 81, FILE = "recnucc", action='readwrite')
> ENDIF
>*
> IF (IZ .EQ. 9) THEN
> WRITE(80,'(1P,9G25.15)') 'c1 ',IZ, IA, IS, X, Y, Z, MREG, WEE
> ENDIF
> RETURN
>*=== End of subroutine Usrrnc =========================================*
> END
--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Jul 25 2019 - 16:05:21 CEST