Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

the new suffix entered for checking ACER run is ignored for Thermal Neutron Scattering Data (NSUB = 12) #163

Closed
MAS-OUD opened this issue May 19, 2020 · 2 comments
Assignees

Comments

@MAS-OUD
Copy link

MAS-OUD commented May 19, 2020

running the following, with the Thermal Neutron Scattering Data sublibrary in tape 21,

acer
-21 -28 0 29 30 /
2 0 1 .00 /
'Processed by NJOY' /
125 500.0 'h-h2o' /
1001 /
222 64 0 0 10.0 /

writes h-h2o.00t in the first line of tape 29. Running

acer
0 29 0 80 90 /
7 0 1 .80 /
'Processed by NJOY' /

also writes h-h2o.00t in the first line of tape 80, but it must be h-h2o.80t.

@whaeck
Copy link
Member

whaeck commented May 30, 2020

FYI: I'm not ignoring this issue, I'm still working on the other one you submitted.

I may pick this one up anyway since the binary atomic relaxation fix is giving me unexpected problems.

@whaeck whaeck self-assigned this May 30, 2020
@MAS-OUD
Copy link
Author

MAS-OUD commented May 31, 2020

replacing the following lines of code in aceth module (starting from line 599, based on 2016.57 version of NJOY)

!--adjust zaid
if (mcnpx.gt.0) then
    read(hz,'(f10.0,a3)') zaid,ht
else
    read(hz,'(a9,a1)') str,ht
    if (ht(1:1).ne.'t') then
        read(hz,'(f9.0,a1)') zaid,ht
    endif
endif
if (suff.ge.zero.and.ht(1:1).ne.'t') then
    iza=nint(zaid)
    zaid=iza+suff
    if (mcnpx.gt.0) then
        write(hz,'(f10.3,a3)') zaid,ht
    else
        write(hz,'(f9.2,a1)') zaid,ht
    endif
endif

with these ones

!--adjust zaid
if (mcnpx.gt.0) then
    read(hz,'(f10.0,a3)') zaid,ht
else
    read(hz,'(a9,a1)') str,ht
endif
if (suff.ge.zero) then
    iza=nint(zaid)
    zaid=iza+suff
    if (mcnpx.gt.0) then
        write(hz,'(f10.3,a3)') zaid,ht
    else
        write(hz,'(a6,f3.2,a1)') str(1:6),suff,ht
    endif
endif

solved this issue for me.

@whaeck whaeck mentioned this issue Jul 27, 2020
@whaeck whaeck closed this as completed Aug 5, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants