When using the source routine source.f are beam parameters defined with the BEAM card preserved?
By default, all settings defined with the beam card are ignored when using source, i.e., you must sample yourself any momentum spread, divergence or distribution of the beam spot. Only the beam momentum, defined with WHAT(1) of the BEAM card is available to the user as variable PBEAM (common block BEAMCM).
In case the particle position, direction and/or polarization should still be driven by the definitions in the input file the following calls to internal routines are available:
* set up the beam particle position LDIREC = .FALSE. CALL STBMPS( LDIREC ) * set up the beam particle direction CALL STBMDV( LDIREC ) * set up the beam particle polarization CALL STBMPL