This topic has been migrated from the SCALE 5 ORIGEN notebook.
Date: Mon Aug 6 14:37:58 2007
Dear Sir
I have a question about using ORIGEN-S & OPUS of SCALE5.1.
Why are the result of ORIGEN and the result of OPUS different about curies of Se79. I attach my input file.
Y, Matsui
---ORIGEN---
charge discharge 1.0 y 2.0 y 3.0 y
se 75 0.000E+00 6.150E+00 7.430E-01 8.975E-02 1.084E-02
se 79 0.000E+00 1.541E-05 1.541E-05 1.541E-05 1.541E-05
kr 85 0.000E+00 6.888E-07 6.458E-07 6.053E-07 5.674E-07
kr 81 0.000E+00 2.914E-13 2.913E-13 2.913E-13 2.913E-13
h 3 0.000E+00 3.535E-15 3.342E-15 3.159E-15 2.986E-15
rb 86 0.000E+00 2.744E-08 3.446E-14 4.329E-20 0.000E+00
sr 89 0.000E+00 5.619E-18 3.755E-20 2.509E-22 1.677E-24
rb 87 0.000E+00 2.006E-22 2.006E-22 2.006E-22 2.006E-22
sr 90 0.000E+00 5.708E-24 5.569E-24 5.434E-24 5.302E-24
total 0.000E+00 2.009E+01 7.430E-01 8.977E-02 1.086E-02
--- Part of output ---
---OPUS---
0 nuclide time (years)
------- --------- --------- --------- ---------
0.000E+00 1.000E+00 2.000E+00 3.000E+00
se75 6.150E+00 7.430E-01 8.975E-02 1.084E-02
se79 1.377E-05 1.377E-05 1.377E-05 1.377E-05
kr85 6.888E-07 6.458E-07 6.053E-07 5.674E-07
kr81 2.914E-13 2.913E-13 2.913E-13 2.913E-13
h3 3.535E-15 3.342E-15 3.159E-15 2.986E-15
rb86 2.744E-08 3.446E-14 4.329E-20 0.000E+00
sr89 5.619E-18 3.755E-20 2.509E-22 0.000E+00
rb87 2.006E-22 2.006E-22 2.006E-22 2.006E-22
sr90 5.708E-24 5.569E-24 5.434E-24 5.302E-24
0 total 2.009E+01 7.430E-01 8.976E-02 1.086E-02
0 subtotal 6.150E+00 7.430E-01 8.976E-02 1.086E-02
--- Input ------------------------
=origens
0$$ a5 28 a11 71 e
1$$ 1
t
test
3$$ a10 0 0 0 a16 2 e
4** 0.65 3.0 0.8 e
t
35$$ 0
t
56$$ 10 10 1 a13 1 a15 3 a18 1 e
57** 0 a3 1e-05 1 e
t
test
test
59** 10r5.0e+13
60** 8i100 1000
66$$ a1 2 a5 2 a9 2 e
73$$ 340000
74** 1.0
75$$ 4
t
54$$ a8 1 e
56$$ a2 10 a10 10 a14 5 3 e
57** a3 1e-05 e
95$$ 0
t
test
test
60** 8i1 10
65$$ 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 e
t
56$$ 0 -10 a6 1 a10 0 e
t
56$$ f0
t
end
=opus
nrank=10
tmax=10
symnuc=se-79 end
end
Dear Y Matsui
OPUS is performing the activity calculation using the atom inventories calculated by ORIGEN and decay constants obtained from a binary library. In this case the decay data used in the OPUS conversion is different than that used by ORIGEN. There have been several recent corrections made to nuclear decay data, Se-79 being a nuclide that has an incorrect decay constant in ENDF/B-VI (see ORIGEN-ARP notebook entry 91 and SCALE notebook entry 110). The default binary library read by OPUS is on logical unit 21 (pwr33gwd), which is made from the card-image library on unit 28. These libraries should therefore have identical decay data (regardless of cross section updating). OPUS can not read the
card-image, but to maintain consistency it should be reading a binary library made from the same card image library on 28. If you have made any changes to the library or are using card-image and binary files of different origin, differences may be observed. To change the binary library used by OPUS set libu= to the correct unit.
Let me know if this does not correct the discrepancy you observed.
Regards
Ian Gauld
Dear Ian
My binary library (SCALE5.1/data/origen.rev02.pwr33gwd in UNIX) was too old. Its file date is 2004/03/31. Maybe it has not been updated.
Today, I tried to run script "makelibraries" for updating libraries. By using this new library, the difference between ORIGEN and OPUS became zero. Thanks.
My way to install SCALE5.1 to unix system was
1. make configure
2. installscale
3. runsmpls
according to README.
But this way did not update pwr33gwd library.
Add to this, it need "makelibraries".
On WINDOWS, same. Is it needed to run the script "scale5.1\work\makedata.bat" manually ?
Regards
Y, Matusi
Dear Yusuke,
This problem occurred because we failed to remake the pwr33gwd binary library after the last updates were made to the card-image decay library, origen.rev03.end6dec.data. We are grateful that you helped us identify this error.
On Windows, you should run the attached origen_data.inp case. The version distributed in \scale5.1\work has not been updated since SCALE 5.0.
We apologize for this mistake.
Regards
Steve Bowman