INFN homepage
FLUKA: 13.2.3} dffcff.f: user-defined DiFFusion CoeFFicient Previous Index Next

13.2.3} dffcff.f: user-defined DiFFusion CoeFFicient


     Argument list (all variables are input only):
          WVLNGT : photon wavelength (in cm)
          OMGPHO : angular frequency (omega = 2pi nu) of the photon (in s-1)
          MMAT   : material index

 Function DFFCFF returns a user-defined diffusion coefficient for optical photons.
 It is activated by setting 
< -99 in command OPT-PROP, with
= blank. See option OPT-PROP and Chap. 12} for more information. 13.2.3b} elefld.f definition of an electric field ------------------------------------------------ Argument list: X,Y,Z,T : current position and time (input only) ETX,ETY,ETZ : direction cosines of the electric field vector (returned). E : electric field intensity in MV/m (returned) NREG : current region (input only) IDISC : if returned = 1, the particle will be discarded ELEFLD is activated by option ELCFIELD with WHAT(4-6) = 0.0 and is used to return intensity and direction of an electric field based on the current position, time and region. It is called only if the current region has been flagged as having a non-zero electric field by option ASSIGNMAt, with
= 2.0. The electric field spatial distribution is often read and interpolated from an external field map. Note that in any case the direction cosines MUST be properly normalised in double precision (e.g. ETX = SQRT(ONEONE - ETY**2 - ETZ**2)), even if E = 0.0. Please read carefully the notes on option ELCFIELD. If name-based input is being used, the name corresponding to NREG can be obtained via a call to routine GEOR2N: CALL GEOR2N (NUMREG, NAMREG, IERR) where NUMREG (input variable) is the region number, and NAMREG (returned variable) is the corresponding region name (to be declared as CHARACTER*8). IERR is a returned error code: if = 0, the conversion is successful. See example in the description of BXDRAW below.

Previous Index Next