nureg/cr-6700 - Oak Ridge National Laboratory
nureg/cr-6700 - Oak Ridge National Laboratory
nureg/cr-6700 - Oak Ridge National Laboratory
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Appendix A<br />
Computer Code Input and Output Listings<br />
Listing A.9 Example of OPUS input for total decay-heat calculations<br />
' $HOME/nrc/ranking/bwr/4.0_gd_rerun/opus_watts_all_bwr5e_4gd_70g.in<br />
'<br />
' this case begins with the concentrations in position 2 of ft71f001<br />
' which has 12 positions. position 1 has the end of irradiation<br />
' concentrations from the sas2 case. position 2 has the concentrations<br />
' from the end of the sas2 decay case (zero length). positions<br />
' 3 through 12 are from the origen-s case run separately.<br />
'<br />
' (1) first link after =shell is to link Ian's<br />
' latest opus executable<br />
'<br />
' (2) link below would use pwr 33 GWd/MTU<br />
' library but this does not use the<br />
' correct <strong>cr</strong>oss sections for absorption<br />
' ln -s $DATA/pwr33gwd<br />
ft33f001<br />
'<br />
' note that any library can be used to<br />
' get the other data needed for opus--<br />
' it takes the first library found in a<br />
' multiburnup library<br />
'<br />
' (3) second link after =shell uses last<br />
' temporary library from ft15f001 in<br />
' the sas2 run--ft15f001 is the unit<br />
' for temporary storage of binary<br />
' libraries before they are added to the<br />
' back of the multiburnup library on<br />
' ft33f001<br />
'<br />
' (3) third link after =shell uses ft71f001<br />
' saved from origen-s decay case run as<br />
' a separate job after sas2 run that<br />
' saved 12 positions on ft71f001<br />
'<br />
=shell<br />
ln -s /home/2ul/bin/opus opus<br />
ln -s $RTNDIR/bwr5e_4gd_70g.ft15f001 ft15f001<br />
ln -s $RTNDIR/bwr5e_4gd_70g.ft71f001 ft71f001<br />
end<br />
=opus<br />
default=no<br />
libunit=15<br />
units=watts<br />
time=year<br />
typarms=nucl<br />
tmax=100.1<br />
libtyp=all<br />
nposition= 2 3 4 5 6 7 8 9 10 11 12 end<br />
xlabel=cooling time (years)<br />
dbug=no<br />
nrank=25<br />
symnuc= ag-110m am-241 am-243 ba-137m cs-134 cs-137<br />
ce-144 cm-242 cm-243 cm-244 eu-154 kr-85<br />
np-239 pm-147 pr-144 pu-238 pu-239 pu-240<br />
pu-241 rh-106 ru-106 sb-125 sr-90 y-90<br />
co-60 end<br />
end<br />
' save plot data file<br />
=shell<br />
cp _plot000 $RTNDIR/_plot.watts_all_bwr5e_4gd_70g<br />
end<br />
79