njoy / NJOY2016

Nuclear data processing with legacy NJOY
https://www.njoy21.io/NJOY2016
Other
95 stars 83 forks source link

Read THERMR::Card2::natom parameter from file instead of from user input #110

Open jlconlin opened 5 years ago

jlconlin commented 5 years ago

The natom parameter in Card2 of THERMR is present on the evaluation file. It should be read from the evaluation file instead of taken from the user. This should be an easy thing to fix.

jlconlin commented 5 years ago

@kahlerac Has the natom always been in ENDF? That is, do evaluations before ENDF/B-VI have that parameter in it? If so, we should check the format number to determine if we can ignore the user input or not.

jlconlin commented 5 years ago

This is the same as the first point in #80

kahlerac commented 5 years ago

natom is not defined in legacy (ENDF/B-III) files but B(1) which is sig(el)*natom is present which is why natom is a user input

jlconlin commented 5 years ago

So if we just check the ENDF version and if it is greater than 3 then we can ignore the natom value, right?

kahlerac commented 5 years ago

e3 was carried forward to e4, e5 and first couple of e6 releases. natom definition appeared when MacFarlane updated the thermal files in e63 (although his report says e62 since that was the release he hoped to have things done by). Skip

On Mon, Nov 5, 2018, 5:09 PM Jeremy Lloyd Conlin <notifications@github.com wrote:

So if we just check the ENDF version and if it is greater than 3 then we can ignore the natom value, right?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/njoy/NJOY2016/issues/110#issuecomment-436054113, or mute the thread https://github.com/notifications/unsubscribe-auth/AQ6UScyihWP4GnMewf_SV26dHByrHqPEks5usLcQgaJpZM4YPQ0o .

jlconlin commented 5 years ago

And now we begin to see why this hasn't been taken out previously—it isn't backwards compatible.

What was B(6) before it was defined to be natom? If it was always 0, then we can check against that.

kahlerac commented 5 years ago

I believe B(6) is unused ( and therefore set to zero) in e3. That's what's in the e3 h-h2o file but I'd have to check an old manual to be sure. Skip

On Mon, Nov 5, 2018, 6:04 PM Jeremy Lloyd Conlin <notifications@github.com wrote:

And now we begin to see why this hasn't been taken out previously—it isn't backwards compatible.

What was B(6) before it was defined to be natom? If it was always 0, then we can check against that.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/njoy/NJOY2016/issues/110#issuecomment-436069000, or mute the thread https://github.com/notifications/unsubscribe-auth/AQ6USYDSGt3Gmk-XPPSHJl7EG5KtUFmgks5usMPggaJpZM4YPQ0o .