13.2.5} 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 
WHAT(5)
 = 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