RE: Rflctv.f not being called.

Thank you Joachim. Flair doesn't provide an explicit menu option for
this input card feature so I guess I kept missing it. Anyhow, after
making this change my RFLCTV function still doesn't get called.

Attached is the project I'm working on. As of now the photons being
traced are visible light photons however I eventually plan to include
infrared surface roughness/absorption/reflectivity curves for all the
materials with blackbody emission coming from a generic region (ie. the
room temperature portion of a cryostat). The included Source file has
optical primaries start in a disc region around the position 0,0,1982
with an isotropic angle distribution.

The attached RFLCTV function should just abort Fluka, however Fluka runs
without ever calling it.

Any advice?

--Nick

On Mon, 2012-02-20 at 17:18 +0100, Joachim Vollaire wrote:
> Dear Nicholas,
> It is hard to tell without your input but did you activate the call to the
> routine with OPT-PROP card ?
>
> RFLCTV: to specify the reflectivity of a material.
> This can be activated by card OPT-PROP with SDUM = METAL and
> WHAT(3) < -99.
>
> Cheers
> Joachim
>
> ________________________________________=
> From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Nicholas Bolibruch [nicholas.bolibruch_at_gmail.com]
> Sent: 08 February 2012 23:02
> To: fluka-discuss_at_fluka.org
> Subject: Rflctv.f not being called.
>
> Hi Fluka experts,
>
> After compiling a slightly modified rflctv.f into my fluka executable,
> I've noticed the routine is never called. Do I need a particular card=
> enabled to trigger the reflectivity user routine? Do I need to add
> another INCLUDE statement?
>
> I've added a FLABRT call to rflctv.f just to test if it is being
> accessed, and the abort request never gets triggered.
>
> I have optical photons generating no problem in a customized source
> routine.
>
> Any help would be greatly appreciated.
>
> Thanks,
>
> Nicholas Bolibruch
>
> Canadian Light Source Inc.
> Instrumentation Analyst
> Controls & Instrumentation Development
>

     0 0
* MlyrW
RCC MlyrW 0.000000 0.000000 1994.130005 0.000000 0.000000 0.010000 0.800000
* OuCuCnI
RCC OuCuCnI 0.000000 11.916800 2008.000000 0.000000 -16.920600 0.000000
               6.470000
* OuCuDnO
RCC OuCuDnO 0.000000 11.916800 2008.000000 0.000000 -16.920600 0.000000
               6.670000
* OuCuPlt
RCC OuCuPlt 0.000000 -5.000000 2008.000000 0.000000 -0.300000 0.000000
               6.670000
* OuStCnI
RCC OuStCnI 0.000000 10.929400 2008.000000 0.000000 -16.920600 0.000000
               8.255000
* OuStCnO
RCC OuStCnO 0.000000 10.929400 2008.000000 0.000000 -16.920600 0.000000
               8.400000
* OuStPlt
RCC OuStPlt 0.000000 -5.990600 2008.000000 0.000000 -0.300000 0.000000
               8.400000
* abConeI
TRC abConeI 0.000000 0.000000 2006.500000 0.000000 0.000000 2.950000 0.780000
               0.050000
* abConeO
TRC abConeO 0.000000 0.000000 2006.500000 0.000000 0.000000 3.000000 0.790000
               0.050000
* airsec
RPP airsec -50.000000 50.000000 -50.000000 50.000000 1845.000000 2145.000000
* fSntO2
RCC fSntO2 0.000000 0.000000 1983.800049 0.000000 0.000000 0.500000 0.700000
* fSntOut
RCC fSntOut 0.000000 0.000000 1983.849976 0.000000 0.000000 0.300000 1.800000
* inBlk
RPP inBlk -50000.000000 50000.000000 -50000.000000 50000.000000
               -50000.000000 50000.000000
* inCuCnI
RCC inCuCnI 0.000000 11.916800 2008.000000 0.000000 -15.925800 0.000000
               5.250000
* inCuCnO
RCC inCuCnO 0.000000 11.916800 2008.000000 0.000000 -15.925800 0.000000
               5.397500
* inCuPlt
RCC inCuPlt 0.000000 -4.013200 2008.000000 0.000000 -0.300000 0.000000
               5.397500
* outPlO
RCC outPlO 0.000000 0.000000 1982.130005 0.000000 0.000000 0.700000 2.800000
* outPlOb
RCC outPlOb 0.000000 0.000000 1982.829956 0.000000 0.000000 0.700000 2.800000
* outPlin
RCC outPlin 0.000000 0.000000 1982.130005 0.000000 0.000000 0.700000 0.700000
* outPlinb
RCC outPlinb 0.000000 0.000000 1982.699951 0.000000 0.000000 0.900000 0.700000
* radBe
RCC radBe 0.000000 0.000000 1982.125000 0.000000 0.000000 0.010000 0.700000
* smTubeI
RCC smTubeI 0.000000 0.000000 1994.150024 0.000000 0.000000 9.600000 0.711200
* smTubeO
RCC smTubeO 0.000000 0.000000 1994.150024 0.000000 0.000000 8.602500 0.760000
* stTubeI
RCC stTubeI 0.000000 0.000000 1984.150024 0.000000 0.000000 17.830000
               1.700000
* stTubeO
RCC stTubeO 0.000000 0.000000 1984.150024 0.000000 0.000000 17.530000
               1.800000
* stelTuI
RCC stelTuI 0.000000 0.000000 1983.530029 0.000000 0.000000 16.720000
               2.600000
* stelTuO
RCC stelTuO 0.000000 0.000000 1983.530029 0.000000 0.000000 16.580000
               2.700000
END
* Reg # 1
* blkhole; assigned material: Blackhole; mat # (1)
blkhole 5 +inBlk -airsec
* Reg # 2
* airDiff; assigned material: Nitrogen; mat # (7)
airDiff 5 +airsec -radBe -( +outPlO -outPlin -outPlOb -outPlinb ) -( +( +inCuCnO
                -inCuCnI -smTubeI )| +inCuPlt | +( +smTubeO -smTubeI ))
                -( +( +OuCuDnO -OuCuCnI -stTubeI )| +OuCuPlt | +( +( +fSntOut -fSntO2 )|
                +( +stTubeO -stTubeI ))) -( +( +OuStCnO -OuStCnI -stelTuI )|+OuStPlt |
                +( +( +stelTuO -stelTuI )| +( +outPlOb -outPlinb ))) -MlyrW -( +abConeO
                -abConeI )
* Reg # 3
* radBe; assigned material: Beryllium; mat # (5)
radBe 5 +radBe
* Reg # 4
* outPlat; assigned material: Iron; mat # (11)
outPlat 5 +outPlO -outPlin -outPlOb -outPlinb
* Reg # 5
* inCuCan; assigned material: Copper; mat # (12)
inCuCan 5 +( +inCuCnO -inCuCnI -smTubeI )| +inCuPlt | +( +smTubeO -smTubeI )
* Reg # 6
* outCuCn; assigned material: Copper; mat # (12)
outCuCn 5 +( +OuCuDnO -OuCuCnI -stTubeI )| +OuCuPlt | +( +( +fSntOut -fSntO2 )|
                +( +stTubeO -stTubeI ))
* Reg # 7
* OStlCan; assigned material: Iron; mat # (11)
OStlCan 5 +( +OuStCnO -OuStCnI -stelTuI )| +OuStPlt | +( +( +stelTuO -stelTuI )|
                +( +outPlOb -outPlinb ))
* Reg # 8
* MlyrW; assigned material: Beryllium; mat # (5)
MlyrW 5 +MlyrW
* Reg # 9
* absorb; assigned material: Copper; mat # (12)
absorb 5 +abConeO -abConeI
END

* C:\Users\bolibrn\Documents\Fluka Optical Photons\Radiometer_with_SimpleAbsorber.dat
* Created: 13.1.2012
* At: 17:12:13
TITLE
MC-CAD Test
GLOBAL 1000.0 0.0 0.0 0.0 1.0 0.
GEOBEGIN
     0 0 MC-CAD
* MlyrW
   RCC MlyrW 0.000000 0.000000 1994.130005 0.000000 0.000000 0.010000
                   0.800000
* OuCuCnI
   RCC OuCuCnI 0.000000 11.916800 2008.000000 0.000000 -16.920600 0.000000
                   6.470000
* OuCuDnO
   RCC OuCuDnO 0.000000 11.916800 2008.000000 0.000000 -16.920600 0.000000
                   6.670000
* OuCuPlt
   RCC OuCuPlt 0.000000 -5.000000 2008.000000 0.000000 -0.300000 0.000000
                   6.670000
* OuStCnI
   RCC OuStCnI 0.000000 10.929400 2008.000000 0.000000 -16.920600 0.000000
                   8.255000
* OuStCnO
   RCC OuStCnO 0.000000 10.929400 2008.000000 0.000000 -16.920600 0.000000
                   8.400000
* OuStPlt
   RCC OuStPlt 0.000000 -5.990600 2008.000000 0.000000 -0.300000 0.000000
                   8.400000
* abConeI
   TRC abConeI 0.000000 0.000000 2006.500000 0.000000 0.000000 2.950000
                 0.780000 0.050000
* abConeO
   TRC abConeO 0.000000 0.000000 2006.500000 0.000000 0.000000 3.000000
                 0.790000 0.050000
* airsec
   RPP airsec -50.000000 50.000000 -50.000000 50.000000 1845.000000 2145.000000
* fSntO2
   RCC fSntO2 0.000000 0.000000 1983.800049 0.000000 0.000000 0.500000
                   0.700000
* fSntOut
   RCC fSntOut 0.000000 0.000000 1983.849976 0.000000 0.000000 0.300000
                   1.800000
* inBlk
   RPP inBlk -50000.000000 50000.000000 -50000.000000 50000.000000 -50000.000000 50000.000000
* inCuCnI
   RCC inCuCnI 0.000000 11.916800 2008.000000 0.000000 -15.925800 0.000000
                   5.250000
* inCuCnO
   RCC inCuCnO 0.000000 11.916800 2008.000000 0.000000 -15.925800 0.000000
                   5.397500
* inCuPlt
   RCC inCuPlt 0.000000 -4.013200 2008.000000 0.000000 -0.300000 0.000000
                   5.397500
* outPlO
   RCC outPlO 0.000000 0.000000 1982.130005 0.000000 0.000000 0.700000
                   2.800000
* outPlOb
   RCC outPlOb 0.000000 0.000000 1982.829956 0.000000 0.000000 0.700000
                   2.800000
* outPlin
   RCC outPlin 0.000000 0.000000 1982.130005 0.000000 0.000000 0.700000
                   0.700000
* outPlinb
   RCC outPlinb 0.000000 0.000000 1982.699951 0.000000 0.000000 0.900000
                   0.700000
* radBe
   RCC radBe 0.000000 0.000000 1982.125000 0.000000 0.000000 0.010000
                   0.700000
* smTubeI
   RCC smTubeI 0.000000 0.000000 1994.150024 0.000000 0.000000 9.600000
                   0.711200
* smTubeO
   RCC smTubeO 0.000000 0.000000 1994.150024 0.000000 0.000000 8.602500
                   0.760000
* stTubeI
   RCC stTubeI 0.000000 0.000000 1984.150024 0.000000 0.000000 17.830000
                   1.700000
* stTubeO
   RCC stTubeO 0.000000 0.000000 1984.150024 0.000000 0.000000 17.530000
                   1.800000
* stelTuI
   RCC stelTuI 0.000000 0.000000 1983.530029 0.000000 0.000000 16.720000
                   2.600000
* stelTuO
   RCC stelTuO 0.000000 0.000000 1983.530029 0.000000 0.000000 16.580000
                   2.700000
   END
* Reg # 1
* blkhole; assigned material: Blackhole; mat # (1)
blkhole 5 +inBlk -airsec
* Reg # 2
* airDiff; assigned material: Nitrogen; mat # (7)
airDiff 5 +airsec -radBe -( +outPlO -outPlin -outPlOb -outPlinb ) -( +( +inCuCnO
            -inCuCnI -smTubeI )| +inCuPlt | +( +smTubeO -smTubeI )) -( +( +OuCuDnO -OuCuCnI -stTubeI )|
            +OuCuPlt | +( +( +fSntOut -fSntO2 )| +( +stTubeO -stTubeI ))) -( +( +OuStCnO -OuStCnI -stelTuI )|
            +OuStPlt | +( +( +stelTuO -stelTuI )| +( +outPlOb -outPlinb ))) -MlyrW -( +abConeO -abConeI )
* Reg # 3
* radBe; assigned material: Beryllium; mat # (5)
radBe 5 +radBe
* Reg # 4
* outPlat; assigned material: Iron; mat # (11)
outPlat 5 +outPlO -outPlin -outPlOb -outPlinb
* Reg # 5
* inCuCan; assigned material: Copper; mat # (12)
inCuCan 5 +( +inCuCnO -inCuCnI -smTubeI )| +inCuPlt | +( +smTubeO -smTubeI )
* Reg # 6
* outCuCn; assigned material: Copper; mat # (12)
outCuCn 5 +( +OuCuDnO -OuCuCnI -stTubeI )| +OuCuPlt | +( +( +fSntOut -fSntO2 )| +( +stTubeO -stTubeI ))
* Reg # 7
* OStlCan; assigned material: Iron; mat # (11)
OStlCan 5 +( +OuStCnO -OuStCnI -stelTuI )| +OuStPlt | +( +( +stelTuO -stelTuI )| +( +outPlOb -outPlinb ))
* Reg # 8
* MlyrW; assigned material: Beryllium; mat # (5)
MlyrW 5 +MlyrW
* Reg # 9
* absorb; assigned material: Copper; mat # (12)
absorb 5 +abConeO -abConeI
   END
GEOEND
ASSIGNMAT BLCKHOLE blkhole
ASSIGNMAT NITROGEN airDiff
ASSIGNMAT BERYLLIU radBe
ASSIGNMAT IRON outPlat
ASSIGNMAT COPPER inCuCan
ASSIGNMAT COPPER outCuCn
ASSIGNMAT IRON OStlCan
ASSIGNMAT BERYLLIU MlyrW
ASSIGNMAT COPPER absorb

TITLE
Inner Snout Test
* Set the defaults for precision simulations
DEFAULTS NEW-DEFA
* Define the beam characteristics
BEAM 3.3E-10 PHOTON
SOURCE
* Define the beam position
BEAMPOS 0.0 0.0 1980.
GEOBEGIN 21. COMBNAME
IR_absorber3.par.geo
GEOEND
OPT-PROP -110. METAL
*OPT-PROP 300E-9 750E-9 1000E-9 COPPER GOLD 1.WV-LIMIT
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA BLCKHOLE blkhole
ASSIGNMA VACUUM airDiff
ASSIGNMA VACUUM radBe
ASSIGNMA COPPER outPlat OStlCan
ASSIGNMA VACUUM MlyrW
ASSIGNMA GOLD absorb
PART-THR -1.E-11 PHOTON 0.0
SCORE ENERGY PHOTON
USRBIN 10. OPTIPHOT -35. 6. 6. 2010.df
USRBIN -6. -6. 1975. 100. 100. 200.&
* Set the random number seed
RANDOMIZ 1.0
* Set the number of primary histories to be simulated in the run
START 1E6 0.0
STOP