Production of 39Ar by 39K(n,p) in a rock - RG03UCpub
c
c ******** BLOCK 1: CELL CARDS ************************************************
  1 10 -2.70 -1    $ rock of -rho[g/cm3] inside of sphere
  2  0       1     $ void outside of sphere

c ******** BLOCK 2: SURFACE CARDS *********************************************
  1 SO  300.     $ sphere or radius[cm] centered at origin

c ******** BLOCK 3: DATA CARDS ************************************************ 
c ======== MATERIAL CARDS ========
M0 nlib=70c    $ default data library
c 
READ file=zaid_atomfrac_RG03UCpub.mcnp6
c 
M91    19039 1.0
M92    12024 1.0
c 
c ======== PHYSICS CARDS ========
PHYS:N 9. 0 0 J J J 0 99.     $ neutron transport, max E 9 MeV, no models
c 
c ======== SOURCE CARDS ========
SDEF  pos=0 0 0  erg=d1       $ position, energy
READ file=sdef_si_232TH_18O.mcnp6
READ file=sdef_sp_232TH_18O.mcnp6
c 
c ======== TALLIES ========
F2:n 1
Fc2 >>>>>>>>>>>>> TALLY 2: FLUX OF NEUTRONS LEAVING CELL 1 <<<<<<<<<<<<<
Sd2 1
c 
F14:n 1
Fc14 >>>>>>>>>>>>> TALLY 14: # OF PROTONS FROM 39K(n,p) ONLY <<<<<<<<<<<<<
Fm14 8.9749508E-04 91 103     $ to account for atom density of 39K
Sd14 0.0231418576             $ to get 39Ar per wt.frac. K
c 
F24:n 1
Fc24 >>>>>>>>>>>>> TALLY 24: # OF ALPHAS FROM 24Mg(n,a) ONLY <<<<<<<<<<<<<
Fm24 7.8722718E-04 92 107     $ to account for atom density of 24Mg
Sd24 0.0148975700             $ to get 21Ne per wt.frac. Mg
c 
c ======== RELATED TO VARRED ========
IMP:n 1 0                     $ neutron importance in cells
c 
c ======== TERMINATION, OUTPUT, OTHER ========
STOP f14=0.005 nps=1000000    $ termination criteria
PRINT                         $ print everything about the calculation