Dear Fluka experts,
I want to split the transport problem into two or more sequential phases.
Frist,I want to use USERDUMP to record all particles crossing a given boundary,so I write the
mgdraw.f like this:
ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
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. "TO" .AND. NRGNAM .EQ. "TW") THEN
IF ( .NOT. LFCOPE ) THEN
LFCOPE = .TRUE.
IF ( KOMPUT .EQ. 2 ) THEN
FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
ELSE
FILNAM = CFDRAW
END IF
OPEN ( UNIT = 66, FILE = FILNAM, STATUS = 'NEW', FORM =
& 'UNFORMATTED' )
END IF
IF (ICODE .EQ. 29) THEN
WRITE (66) 0, ICODE, JTRACK, SNGL (ETRACK), SNGL (WTRACK)
WRITE (66) SNGL (XSCO), SNGL (YSCO), SNGL (ZSCO)
END IF
END IF
and I use USERDUMP in the inp file like this£º
USERDUMP 200 49.0 0.0 TRAKFILE
I can get the TRAKFILE,but it's a binary file,I don't know whether it true or not.
When I meet this, I think I have two chooses:
1¡¢I should read the binary file at next phase
2¡¢I should change a way to write a file which is not a binary file
I choose the first one,and I write the source.f like this:
Open(60,File='TRACKFILE',Access='Direct',Form='Unformatted')
READ(60,*) IPART,ICODE,JTRACK,ETRACK,WTRACK,XSCO,YSCO,ZSCO
The source.f is ok, but it can't run ,It always keep running like this:
======================= Running FLUKA for cycle # 1 =======================
I don't know how to deal with it
If I choose the second one,how should I do?
Please help me! Thanks!
Thanks in advance
xuchao
University Chinese Academy of Sciences
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Nov 08 2016 - 11:31:30 CET