--=-Xxzvpp7KydpnuuIUTMcb
Content-Type: text/plain; charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
Dear Fluka experts
I am still having the problem described in a) crash in parse.f, using
fluka2011.2.5 with flair-0.9.4 and gfortran 4.5.2.
Output:
At line 6 of file parse.f
Fortran runtime error: Actual string length is shorter than the declared
one for dummy argument 'txline' (80/132)
Backtrace for this error:
+ function parse_ (0x41914F)
at line 55 of file parse.f
+ function geoskp_ (0x4D1B36)
at line 204 of file geoskp.f
+ function flukam_ (0x40C570)
at line 1328 of file flukam.f
+ function fluka (0x4081AC)
at line 289 of file fluka.f
+ /lib64/libc.so.6(__libc_start_main+0xfd) [0x7fbfc3a73cdd]
The problem only arises with certain input files, not with every one it
used to with the 2011.2.4 respin, when I try to plot the 2D Projection
and selecting "Geometry Use: Auto". I did recompile my user routine and
executable after the upgrade.
One of the inputs that does not work is attached with the needed files.
Is there anything I might have missed after the upgrade to the new
respin?
Thanks for your efforts,
Roger
On Wed, 2011-09-07 at 18:05 +0200, Alfredo Ferrari wrote:
> Dear Fluka users
>
> the fluka2011.2.5 respin has been pushed to the Fluka website.
> It is a bug fix respin which addresses the following issues:
>
> a) a problem when running Flair using the gfortran version (crash in
> parse.f) is fixed
> b) a problem when using the QUA geometric body to describe
> hyperboloid-like shapes is fixed
> c) the region-dependent group cut-off for low energy neutrons is now
> strictly enforced, eg at the boundaries in case of neighboring regions
> with different cutoffs
> d) the energy scored in the USRTRACK estimators for e+/e- is now the
> average one of the step (before it was the end step one), in order
> to make it consistent with what is scored for charged hadrons and
> muons
> e) a couple of rare random number non-reproducibility issues have
> solved in Fluka and rQMD
>
> All users are urged to update to this respin. Since there is no change in
> physics, only bug corrections, updating is harmless, or better it is
> critical if you were touched by one of these issues. Anyway, because of
> issues c) and e), the random number sequences could be different, hence
> debugging/help/fixing will be provided from now on only for this respin.
>
> The solution of b) required changes to one common block, hence recompiling
> of possible user routines is recommended.
>
> As usual, please report to fluka-discuss whichever problem you are
> experiencing
>
> Have fun
> The FLUKA developers
>
>
--=-Xxzvpp7KydpnuuIUTMcb
Content-Disposition: attachment;
filename="ambe_neutrons_water_phantom.flair"
Content-Type: text/plain; name="ambe_neutrons_water_phantom.flair";
charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
# FLUKA Project file
Version: None
Title: AmBe neutron source impinging on a water phantom
Input: ambe_neutrons_water_phantom.inp
Exec: myfluka
Submit: *Default
Debug: "Region #1" -15.1 -15.1 -5.1 15.1 15.1 15.1 100 100 100
Source: source.f
DefMain: True
LinkPrg: lfluka
F77bound: True
F77dline: False
# Run information
Run: <default>
Prev: 0
Last: 2
Status: 3
Pid: 0
StartRun: 1315482256
# USRxxx data file: water_phantom-neutron_fluence_usrbin
Data: water_phantom-neutron_fluence_usrbin
Unit: 31
Type: b
Recursive: 0
Rule: +,\I\d\d\d_fort\.\U
Rule: +,\I\d\d\d_ftn\.\U
Rule: +,_\I_\d\d\d\d\d_fort\.\U
End
End
# Geometry Information
Geometry:
End
# Geometry plot "Geometry AmBe neutron source impinging on a water phantom"
Plot: ambe_neutrons_water_phantom_geometry
Title: Geometry AmBe neutron source impinging on a water phantom
Format: .eps
Type: Geometry
LineType:
Keys: 1
origin: 0 0 0
showcenter: 0
vectscale: 0.1
labels: 1
coord: X-Y
boundaries: 0
subtype: Geometry
extends: 100 100
basisu: 1 0 0
basisv: 0 1 0
nu: 200
nv: 200
End
# USRBIN plot "Neutron fluence for AmBe neutron source and water phantom"
Plot: ambe_neutrons_water_phantom-neutron_fluence
Title: Neutron fluence for AmBe neutron source and water phantom
Format: .eps
Type: USRBIN
LineType:
Keys: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Neutron fluence
cbround: 1
cbcpd: 3
int: 3.7105618015068846
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 5.16555019E-06
xrebin: 1
axes: Auto
hist: 2D Projection
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 9.40733735E-05
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: water_phantom-neutron_fluence_usrbin
det: 1
cbpalette: FLUKA
pt.0: 1
End
--=-Xxzvpp7KydpnuuIUTMcb
Content-Disposition: attachment; filename="ambe_neutrons_water_phantom.inp"
Content-Type: text/plain; name="ambe_neutrons_water_phantom.inp";
charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7...
TITLE
AmBe neutron source impinging on a water phantom
*
* use names everywhere and free format for geometry
DEFAULTS
PRECISIO
RANDOMIZ 1. 0.0
START 1.E6
*
* beam definitions
BEAM -.0217993 0.0 0.0 1.0 1.0
1.0NEUTRON
BEAMPOS 0.0 0.0 -75.0 0.0 0.0
SOURCE
EMFRAY 4. BLKHOLE @LASTREG 1.0
DELTARAY 0.00001 100.0 1.03 AIR
NOPRINT
DELTARAY 0.00001 100.0 1.03 WATER
NOPRINT
DELTARAY 0.00001 100.0 1.03 PMMA
NOPRINT
DELTARAY 0.00001 100.0 1.03 Lexan
NOPRINT
DELTARAY 0.00001 100.0 1.03 Polyethy
NOPRINT
GEOBEGIN 21.
COMBNAME
water_phantom.geo
GEOEND
MATERIAL 1.0 WATER
COMPOUND 2.0 HYDROGEN 1.0 OXYGEN WATER
MATERIAL 0.1956 CO2
COMPOUND 1.0 CARBON 2.0 OXYGEN CO2
MATERIAL 0.001204 AIR
COMPOUND -.75518 NITROGEN -0.23135 OXYGEN -.01288 ARGONAIR
COMPOUND -0.00058 CO2 AIR
* Polymethyl Methacrylate PMMA Photoresist
* H
* Chemical H H-C-H
* Formula | |
* Polymethyl ----- C - C -----
* Methacrylate | | H
* H C - O - C-H
* C H O || H
* 5 8 2 O
MATERIAL 0.95 PMMA
COMPOUND 8.0 HYDROGEN 5.0 CARBON 2.0 OXYGENPMMA
* Polycarbonate Lexan,Makrofol
* H
* Chemical H-C - C-H H-C-H H-C - C-H O
* Formula // \\ | // \\ ||
* --- O - C C -- C -- C C - O - C ---
* C H 0 \ / | \ /
* 16 14 3 H-C = C-H H-C-H H-C = C-H
* H
* The density variation is very narrow, from 1.19 - 1.22 g/cm3.
MATERIAL 1.2 Lexan
COMPOUND 14.0 HYDROGEN 16.0 CARBON 3.0 OXYGENLexan
* Polyethylene Marlex
* Density variation of polyethylene is 0.91 - 1.05 g/cm3. "Low" density is
* 0.920, "medium" is .93, and "high" is .95. Special polyethelene is made
* for nuclear shielding, and this has loaded densities up to 1.08.
* Chemical H H
* Formula : | |
* |C H | ---- C -- C ----
* | 2 4|n | |
* H H
MATERIAL 0.93
Polyethy
COMPOUND 4.0 HYDROGEN 2.0 CARBON
Polyethy
STERNHEI 3.4354 .2211 2.7799 .08926 3.511
TISSUEIC
* Material assignments
ASSIGNMA BLCKHOLE BLKHOLE
ASSIGNMA VACUUM VACUUM
ASSIGNMA AIR AIR
EMFFIX AIR 0.1 WATER 0.05 PMMA 0.05
EMFFIX Lexan 0.05 Polyethy 0.05
EMFCUT -1.0E-06 1.0E-06 0.0 AIR
PROD-CUT
EMFCUT -1.0E-06 1.0E-06 0.0 WATER
PROD-CUT
EMFCUT -1.0E-06 1.0E-06 0.0 PMMA
PROD-CUT
EMFCUT -1.0E-06 1.0E-06 0.0 Lexan
PROD-CUT
EMFCUT -1.0E-06 1.0E-06 0.0 Polyethy
PROD-CUT
* Fluence
USRBIN 10. NEUTRON -31. 20.0 20.0 20.0FLUTOT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 110.0&
ASSIGNMA WATER WPHAN 1.0
ASSIGNMA PMMA PWALL
ASSIGNMA Lexan VSCOR0
ASSIGNMA Polyethy RAD1
ASSIGNMA Polyethy RAD2
STOP
--=-Xxzvpp7KydpnuuIUTMcb
Content-Disposition: attachment; filename="source.f"
Content-Type: text/x-fortran; name="source.f"; charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
SUBROUTINE SOURCE ( NOMORE )
INCLUDE '(DBLPRC)'
INCLUDE '(DIMPAR)'
INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2006 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* New source for FLUKA9x-FLUKA200x: *
* *
* Created on 07 january 1990 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 03-mar-06 by Alfredo Ferrari *
* *
* This is just an example of a possible user written source routine. *
* note that the beam card still has some meaning - in the scoring the *
* maximum momentum used in deciding the binning is taken from the *
* beam momentum. Other beam card parameters are obsolete. *
* *
*----------------------------------------------------------------------*
*
INCLUDE '(BEAMCM)'
INCLUDE '(FHEAVY)'
INCLUDE '(FLKSTK)'
INCLUDE '(IOIOCM)'
INCLUDE '(LTCLCM)'
INCLUDE '(PAPROP)'
INCLUDE '(SOURCM)'
INCLUDE '(SUMCOU)'
*
LOGICAL LFIRST
*
c defining and saving spectrum arrays
DIMENSION ENEPOI(0:1000),ENEPRO(0:1000),ENECUM(0:1000)
SAVE ENEPOI, ENEPRO, ENECUM
c saving spectrum dimension
SAVE IMAX
*
SAVE LFIRST
DATA LFIRST / .TRUE. /
*======================================================================*
* *
* BASIC VERSION *
* *
*======================================================================*
NOMORE = 0
* +-------------------------------------------------------------------*
* | First call initializations:
IF ( LFIRST ) THEN
* | *** The following 3 cards are mandatory ***
TKESUM = ZERZER
LFIRST = .FALSE.
LUSSRC = .TRUE.
* | *** User initialization ***
CALL OAUXFI('spectrum.dat',LUNRDB,'OLD',IERR)
c reading spectrum
DO I=0,1000
READ(LUNRDB,*,END=1972) ENEPOI(I),ENEPRO(I)
IMAX=I
ENDDO
STOP ' spectrum reading uncomplete!'
1972 CONTINUE
ENECUM(0)=ZERZER
c building cumulative spectrum
DO I=1,IMAX
ENECUM(I)=(ENEPRO(I-1)+ENEPRO(I))*(ENEPOI(I)-ENEPOI(I-1))
& +ENECUM(I-1)
ENDDO
c normalizing cumulative spectrum
DO I=1,IMAX
ENECUM(I)=ENECUM(I)/ENECUM(IMAX)
ENDDO
END IF
* |
* +-------------------------------------------------------------------*
* Push one source particle to the stack. Note that you could as well
* push many but this way we reserve a maximum amount of space in the
* stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
WTFLK (NPFLKA) = ONEONE
WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
* +-------------------------------------------------------------------*
* | (Radioactive) isotope:
IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
IARES = IPROA
IZRES = IPROZ
IISRES = IPROM
CALL STISBM ( IARES, IZRES, IISRES )
IJHION = IPROZ * 1000 + IPROA
IJHION = IJHION * 100 + KXHEAV
IONID = IJHION
CALL DCDION ( IONID )
CALL SETION ( IONID )
* |
* +-------------------------------------------------------------------*
* | Heavy ion:
ELSE IF ( IJBEAM .EQ. -2 ) THEN
IJHION = IPROZ * 1000 + IPROA
IJHION = IJHION * 100 + KXHEAV
IONID = IJHION
CALL DCDION ( IONID )
CALL SETION ( IONID )
ILOFLK (NPFLKA) = IJHION
* | Flag this is prompt radiation
LRADDC (NPFLKA) = .FALSE.
* |
* +-------------------------------------------------------------------*
* | Normal hadron:
ELSE
IONID = IJBEAM
ILOFLK (NPFLKA) = IJBEAM
* | Flag this is prompt radiation
LRADDC (NPFLKA) = .FALSE.
END IF
* |
* +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
LOFLK (NPFLKA) = 1
* User dependent flag:
LOUSE (NPFLKA) = 0
* User dependent spare variables:
DO 100 ISPR = 1, MKBMX1
SPAREK (ISPR,NPFLKA) = ZERZER
100 CONTINUE
* User dependent spare flags:
DO 200 ISPR = 1, MKBMX2
ISPARK (ISPR,NPFLKA) = 0
200 CONTINUE
* Save the track number of the stack particle:
ISPARK (MKBMX2,NPFLKA) = NPFLKA
NPARMA = NPARMA + 1
NUMPAR (NPFLKA) = NPARMA
NEVENT (NPFLKA) = 0
DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
AGESTK (NPFLKA) = +ZERZER
AKNSHR (NPFLKA) = -TWOTWO
* Group number for "low" energy neutrons, set to 0 anyway
IGROUP (NPFLKA) = 0
c
c sampling from the normalized cumulative spectrum
XYZ=FLRNDM(XYZ)
DO I=1,IMAX
IF(XYZ.LT.ENECUM(I)) GOTO 1973
ENDDO
STOP ' I did a big mistake'
1973 CONTINUE
c the sampled energy lies in the bin I (between ENEPOI(I-1) and ENEPOI(I))
c now determining the energy inside the bin I according to a linear
spectrum
XYZ=(XYZ-ENECUM(I-1))/(ENECUM(I)-ENECUM(I-1))
IF(ENEPRO(I).EQ.ENEPRO(I-1)) THEN
c indeed, inside the bin I the spectrum is flat!
ESAMPLE=XYZ
ELSE
ESAMPLE=XYZ*(ENEPRO(I)-ENEPRO(I-1))*(ENEPRO(I)+ENEPRO(I-1))
ESAMPLE=SQRT(ESAMPLE+ENEPRO(I-1)**2)-ENEPRO(I-1)
ESAMPLE=ESAMPLE/(ENEPRO(I)-ENEPRO(I-1))
ENDIF
ESAMPLE=ESAMPLE*(ENEPOI(I)-ENEPOI(I-1))+ENEPOI(I-1)
* Kinetic energy of the particle (GeV)
TKEFLK (NPFLKA) = ESAMPLE
* Particle momentum
PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
& + TWOTWO * AM (IONID) ) )
* Cosines (tx,ty,tz)
PI = 3.14159265
PHI = TWOTWO * PI * FLRNDM(XXX)
THETA = ACOS ( ONEONE - FLRNDM(XXX) )
TXFLK (NPFLKA) = SIN ( THETA ) * COS ( PHI )
TYFLK (NPFLKA) = SIN ( THETA ) * SIN ( PHI )
TZFLK (NPFLKA) = COS ( THETA )
* TZFLK (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
* & - TYFLK (NPFLKA)**2 )
* Polarization cosines:
TXPOL (NPFLKA) = -TWOTWO
TYPOL (NPFLKA) = +ZERZER
TZPOL (NPFLKA) = +ZERZER
* Particle coordinates
XFLK (NPFLKA) = XBEAM
YFLK (NPFLKA) = YBEAM
ZFLK (NPFLKA) = ZBEAM
* Calculate the total kinetic energy of the primaries: don't change
IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 )
& THEN
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
& * WTFLK (NPFLKA)
ELSE
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
END IF
RADDLY (NPFLKA) = ZERZER
* Here we ask for the region number of the hitting point.
* NREG (NPFLKA) = ...
* The following line makes the starting region search much more
* robust if particles are starting very close to a boundary:
CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
& NRGFLK(NPFLKA), IDISC )
* Do not change these cards:
CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
NLATTC (NPFLKA) = MLATTC
CMPATH (NPFLKA) = ZERZER
CALL SOEVSV
RETURN
*=== End of subroutine Source =========================================*
END
--=-Xxzvpp7KydpnuuIUTMcb
Content-Disposition: attachment; filename="spectrum.dat"
Content-Type: text/plain; name="spectrum.dat"; charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
1.03785E-011 0.0404233
1.13887E-011 0.0404233
1.24972E-011 0.0404233
1.37136E-011 0.0405672
1.50484E-011 0.0405672
1.65131E-011 0.0405672
1.81204E-011 0.0405672
1.98842E-011 0.0405672
2.18196E-011 0.0405672
2.39434E-011 0.0405672
2.62739E-011 0.0405672
2.88312E-011 0.0405672
3.16375E-011 0.0405672
3.47169E-011 0.0405672
3.8096E-011 0.0405672
4.1804E-011 0.0405672
4.5873E-011 0.0405672
5.0338E-011 0.0405672
5.52376E-011 0.0405672
6.06141E-011 0.0405672
6.65139E-011 0.0405672
7.2988E-011 0.0405672
8.00922E-011 0.0405672
8.78879E-011 0.0405672
9.64424E-011 0.0405672
1.05829E-010 0.0405672
1.1613E-010 0.0405672
1.27434E-010 0.0405672
1.39837E-010 0.0405672
1.53448E-010 0.0405672
1.68384E-010 0.0405672
1.84774E-010 0.0405672
2.02758E-010 0.0405672
2.22494E-010 0.0405672
2.4415E-010 0.0405672
2.67914E-010 0.0405672
2.93991E-010 0.0405672
3.22606E-010 0.0405672
3.54007E-010 0.0405672
3.88464E-010 0.0405672
4.18429E-010 0.00561036
4.3832E-010 0.00489108
4.94574E-010 0.0054665
5.42713E-010 0.0054665
5.95538E-010 0.0054665
6.53504E-010 0.0054665
7.17112E-010 0.0054665
7.86911E-010 0.0054665
8.63505E-010 0.0054665
9.47553E-010 0.0054665
0.000000001 0.0054665
1.14099E-009 0.0054665
1.25205E-009 0.0054665
1.37391E-009 0.0054665
1.50764E-009 0.00561036
1.65438E-009 0.00575421
1.81541E-009 0.00575421
0.000000002 0.00575421
2.18601E-009 0.00575421
2.39879E-009 0.00575421
2.63227E-009 0.00575421
2.88848E-009 0.00575421
3.14032E-009 0.00604192
3.41413E-009 0.00632963
3.74644E-009 0.00632963
4.11109E-009 0.00632963
4.51124E-009 0.00632963
0.000000005 0.00632963
5.43218E-009 0.00632963
0.000000006 0.00632963
6.42072E-009 0.0067612
6.85206E-009 0.00733662
7.51899E-009 0.00733662
8.25085E-009 0.00733662
9.05394E-009 0.00733662
9.93519E-009 0.00733662
1.09022E-008 0.00733662
0.000000012 0.00733662
1.31278E-008 0.00733662
0.000000014 0.00791204
0.000000015 0.00848746
1.64061E-008 0.00848746
0.000000018 0.00848746
1.97553E-008 0.00848746
2.16782E-008 0.00848746
2.37882E-008 0.00848746
2.61036E-008 0.00848746
2.83795E-008 0.00863132
2.97286E-008 0.00949445
0.000000032 0.00992602
3.51386E-008 0.00992602
3.85588E-008 0.00992602
4.23119E-008 0.00992602
4.64303E-008 0.00992602
5.09495E-008 0.00992602
5.59086E-008 0.00992602
6.13504E-008 0.0106453
6.42669E-008 0.0115084
7.05222E-008 0.0115084
7.73864E-008 0.0115084
8.49188E-008 0.0115084
9.31842E-008 0.0115084
1.02254E-007 0.0115084
1.12207E-007 0.0115084
1.23129E-007 0.0115084
1.30185E-007 0.0123716
1.36374E-007 0.0132347
1.49648E-007 0.0132347
1.64214E-007 0.0132347
1.80197E-007 0.0132347
1.97737E-007 0.0132347
0.000000217 0.0132347
2.38103E-007 0.0132347
2.61278E-007 0.0132347
2.76253E-007 0.013954
2.89385E-007 0.0148171
3.17552E-007 0.0148171
3.48461E-007 0.0148171
3.82378E-007 0.0148171
4.19596E-007 0.0148171
4.60437E-007 0.0148171
5.05254E-007 0.0148171
5.54432E-007 0.0148171
5.86208E-007 0.0155364
6.25589E-007 0.0162557
6.8648E-007 0.0162557
7.53297E-007 0.0162557
8.26619E-007 0.0162557
9.07077E-007 0.0162557
9.95367E-007 0.0162557
1.09225E-006 0.0162557
1.19856E-006 0.0162557
1.26726E-006 0.0169749
0.000001365 0.0174065
1.49787E-006 0.0174065
1.64366E-006 0.0174065
1.80365E-006 0.0174065
1.9792E-006 0.0174065
2.17185E-006 0.0174065
2.38324E-006 0.0174065
0.000002591 0.0178381
2.7651E-006 0.0182696
3.03423E-006 0.0182696
3.32957E-006 0.0182696
3.65365E-006 0.0182696
4.00927E-006 0.0182696
4.39951E-006 0.0182696
4.82774E-006 0.0182696
5.29764E-006 0.0182696
5.75953E-006 0.0187012
6.32013E-006 0.0187012
6.9353E-006 0.0187012
7.61034E-006 0.0187012
8.35108E-006 0.0187012
9.16392E-006 0.0187012
1.00559E-005 0.0187012
1.10347E-005 0.0187012
1.21087E-005 0.0187012
1.32873E-005 0.0187012
1.45806E-005 0.0187012
1.59998E-005 0.0187012
1.75571E-005 0.0187012
0.000019266 0.0187012
2.11413E-005 0.0187012
0.000023199 0.0187012
2.54571E-005 0.0187012
2.79349E-005 0.018845
0.000030654 0.018845
3.36376E-005 0.018845
3.69117E-005 0.018845
4.05045E-005 0.018845
4.44469E-005 0.018845
4.87731E-005 0.018845
5.25354E-005 0.0192766
5.71159E-005 0.0197082
6.26752E-005 0.0197082
6.87756E-005 0.0197082
7.54698E-005 0.0197082
8.28156E-005 0.0197082
9.08764E-005 0.0197082
9.97217E-005 0.0197082
0.000109428 0.0197082
0.000116779 0.0234484
0.000128146 0.0234484
0.000140619 0.0234484
0.000154306 0.0234484
0.000169325 0.0234484
0.000185806 0.0234484
0.000203891 0.0234484
0.000219619 0.0237361
0.000232206 0.0408549
0.000254807 0.0409988
0.000279609 0.0409988
0.000306824 0.0409988
0.000336689 0.0409988
0.00036946 0.0409988
0.000405421 0.0409988
0.000440769 0.0411426
0.000466031 0.0699137
0.000511392 0.0699137
0.000561167 0.0699137
0.000615788 0.0699137
0.000675725 0.0699137
0.000741496 0.0699137
0.000813669 0.0699137
0.000892867 0.0699137
0.000944039 0.143424
0.00103593 0.143424
0.00113676 0.143424
0.0012474 0.143424
0.00136882 0.143424
0.00150205 0.143424
0.00164825 0.143424
0.00180868 0.143424
0.00193019 0.246856
0.00209848 0.246856
0.00230273 0.246856
0.00252686 0.246856
0.00277281 0.246856
0.0030427 0.246856
0.00333886 0.247
0.00359641 0.246856
0.00383801 0.19291
0.00417265 0.192766
0.00457879 0.192766
0.00502446 0.192766
0.00551351 0.192766
0.00605016 0.192766
0.00663905 0.192766
0.0072179 0.192622
0.00770279 0.0201397
0.00829696 0.0202836
0.00910453 0.0202836
0.00999072 0.0202836
0.0109632 0.0202836
0.0120302 0.0202836
0.0132012 0.0202836
0.0144861 0.0202836
0.0154593 0.000143855
0.0164978 0.000143855
0.0181036 0.000143855
0.019682 0.000143855
0.0217993 0.000143855
--=-Xxzvpp7KydpnuuIUTMcb
Content-Disposition: attachment; filename="water_phantom.geo"
Content-Type: text/plain; name="water_phantom.geo"; charset="ISO-8859-15"
Content-Transfer-Encoding: 7bit
0 0 Rectangular phantom (30x30x15 cm^3) at zero with SSD75
*
* Bodies
* ------
*
* General definitions:
* blackhole to include geometry
SPH BLK 0.0 0.0 0.0 10000.0
* void
RPP VAC -1000.0 1000.0 -1000.0 1000.0 -1000.0 1000.0
* Treatment room air
RPP ROOM -100.0 100.0 -100.0 100.0 -100.0 100.0
* Body for phantom
RPP PHANTOM -15.0 15.0 -15.0 15.0 0.0 15.0
RPP PMMA -14.0 14.0 -14.0 14.0 0.25 14.0
* Scoring volumes
RPP SCOR0 -0.5 0.5 -0.5 0.5 -0.55001 -0.00001
XYP RAD1 -0.35001
XYP RAD2 -0.20001
END
*
* Regions
* -------
*
* Blackhole
BLKHOLE 5 +BLK -VAC
* Vacuum outside treatment room
VACUUM 5 +VAC -ROOM
* Air around target
AIR 5 +ROOM -PHANTOM -SCOR0
PWALL 5 +PHANTOM -PMMA
* Water phantom
WPHAN 5 +PMMA
* Radiator 1
RAD1 5 +SCOR0 +RAD1
* Radiator 2
RAD2 5 +SCOR0 -RAD2
* CR39
VSCOR0 5 +SCOR0 -RAD1 +RAD2
END
--=-Xxzvpp7KydpnuuIUTMcb--
Received on Thu Sep 08 2011 - 18:10:53 CEST
This archive was generated by hypermail 2.2.0 : Thu Sep 08 2011 - 18:10:54 CEST