GCR extended layers and scoring through mgdraw file

Advanced Features and User Routines
3

Dear @tuneu ,

In reply to your questions:

  1. Indeed the normalization factor is applied by default only to the (100) standard atmospheric layers. If you want to add more layers, you can bypass this step and apply the normalization factor in post-process.
    Firstly, make sure you select the option NO-NORM in the GCR-SPE card. Please note that you may need to rename the *.spc files as 01NO-NORM.spc, 02NO-NORM.spc, etc, in order to run the simulation.
    Secondly, in the .out you can find the global normalization factor FLUXST (integral over energy and angle) and the equivalent flux per primary particle FLUX. The normalization factor is then given by the ratio FLUXST/FLUX. Please note that your USRTRACK/USRBDX scoring must be properly normalized. For example, if you don’t provide the volume in the USRTRACK card you must add this factor in the final post-process normalization.
    Below an example of what you can find the *.out:
  *** Fluxst:  0.651134312 ions/cm^2/s
  *** for full earth shadowing ***
  *** <Z_ion>:  1.16797233 ***
  *** Emission radius:  644899968. cm ***
  *** Flux:  3.82679652E-19 Part/pr/cmq ***

  1. Please make sure that the regions that you have specified in the mgdraw are correct. Without the *.inp is difficult to cross-check. Please note also the large curvature radius of the model.

  2. Please note that the numbers provided in the example shouldn’t be considered as universal/all-purposes values. You may change these parameters depending on your specific problem/application.

Cheers,

Angelo