With reference to above discussion I am trying to sample the same neutron spectrum but in a cylindrical shape of Radius = 5 cm and Height = 10 cm.
I am able to compile the source file, but upon scoring the neutron flux using USRTRACK I get the following error.
“Abort called from DEDX reason Exceeded upper dp/dx abscissa. Run stopped!
STOP Exceeded upper dp/dx abscissa.”
Kindly guide me whether I have written the right code for cylindrical sampling in source_newgen.f or am I making mistake in scoring using USRTRACK.
Thank you so much for your kind reply and identifying my goof up. I was able to run the file. However the results from USRTRACK does not match the spectrum definition in the text file(kindly see the attached image). The spectrum definition in the text file is of the order of 1e14, while the results of USRTRACK scoring are in the order of 1E10.
neutron_spectrum.txt (2.9 KB) source_newgen.f (19.3 KB)
Also can you kindly proof-check the algorithm used for cylindrical sampling, I have modified it a little seeing the annular sampling algorithm:
after defining all necessary variables