09.05.2014 Views

nureg/cr-6700 - Oak Ridge National Laboratory

nureg/cr-6700 - Oak Ridge National Laboratory

nureg/cr-6700 - Oak Ridge National Laboratory

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!