File:Puconcentrationsplot.jpg

Relative concentrations $C$ of various isotopes of Plutonium, measured at the unload of fuel from various reactors, versus burn-out $B$, measured in teraJoules per kilogram of the fuel loaded.

Data are by Plutonium Topic Strategy, 2008. .

In the original table, the burn-out parameter $B$ is measured in MW d/t that seems to be MegaWatt days per tonne. This parameter is converted to the International system, taking into account, that day is 24*3600 seconds, and tonne is 1000kg.

Table of the input
RType 86400 238  239  240  241  242 100 Magnox 3000 0.1 80   16.9  2.7  0.3 82.7 Magnox 5000 -1  68.5 25.0  5.3  1.2 73.8 CANDU 7500  -1  66.6 26.6  5.3  1.5 71.9 AGR  18000  0.6 53.7 30.8  9.9  5.0 63.6 BWR  27500  2.6 59.8 23.7 10.6  3.3 70.4 BWR  30400  -1  56.8 23.8 14.3  5.1 71.1 PWR  33000  1.3 56.6 23.2 13.9  4.7 70.5 PWR  43000  2.0 52.5 24.1 14.7  6.2 67.2 PWR  53000  2.7 50.4 24.1 15.2  7.1 65.6 Reactor Type & Mean fuel burn-up (MW d/t) & Percentage of & u Isotopes at Discharge & Fissile Content

C++ generator of curves

 * 1) include
 * 2) include
 * 3) include
 * 4) define DB double

void ado(FILE *O, int X, int Y) {      fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%'); fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y); fprintf(O,"/M {moveto} bind def\n"); fprintf(O,"/L {lineto} bind def\n"); fprintf(O,"/S {stroke} bind def\n"); fprintf(O,"/s {show newpath} bind def\n"); fprintf(O,"/C {closepath} bind def\n"); fprintf(O,"/F {fill} bind def\n"); fprintf(O,"/o {.01 0 360 arc C S} bind def\n"); fprintf(O,"/O {.03 0 360 arc C S} bind def\n"); fprintf(O,"/times-Roman findfont 20 scalefont setfont\n"); fprintf(O,"/W {setlinewidth} bind def\n"); fprintf(O,"/RGB {setrgbcolor} bind def\n");}

int main{int n,m; char name[10][10]; for(n=0;n<10;n++) sprintf(name[n],"00000\n"); DB E[10],e[10]; DB p238[10]; DB p239[10]; DB p240[10]; DB p241[10]; DB p242[10]; DB p243[10]; // (this is hust sum) DB MWt=1.e6; //Watt //DB Year=365.24*24*60*60; // second DB Day=24*60*60; //second printf("Day=%8.4lf %8.4le \n",Day, 53*Day); DB tonne=1000.; //kg

FILE *o; o=fopen("table.txt","r"); for(n=0;n<10;n++){ m=fscanf(o,"%s%lf%lf%lf%lf%lf%lf%lf",name[n], &E[n], &p238[n],&p239[n],&p240[n],&p241[n],&p242[n],&p243[n]); e[n]=E[n]*(MWt*Day/tonne); } fclose(o);

n=0; printf("%2d %6s %6.0lf %8.2le %4.0lf %4.0lf %4.0lf %4.0lf %4.0lf %4.0lf\n",n,name[n], E[n], e[n], p238[n],p239[n],p240[n],p241[n],p242[n],p243[n]);

for(n=1;n<10;n++){ printf("%2d %6s %6.0lf %8.2le %4.1lf %4.1lf %4.1lf %4.1lf %4.1lf %4.1lf\n",n,name[n], E[n], e[n],p238[n],p239[n],p240[n],p241[n],p242[n],p243[n]); } o=fopen("21.eps","w"); ado(o,504,204);
 * 1) define M(x,y) fprintf(o,"%6.4lf %6.4lf M\n",0.+x,0.+y);
 * 2) define L(x,y) fprintf(o,"%6.4lf %6.4lf L\n",0.+x,0.+y);
 * 3) define O(x,y) fprintf(o,"%6.4lf %6.4lf O\n",0.+x,0.+y);

fprintf(o,"2 2 translate 100 100 scale 2 setlinecap 1 setlinejoin\n"); DB sy=2; M(5,0) L(0,0) L(0,2) fprintf(o,".007 W S\n"); for(n=1;n<6;n++){ M(n,0)L(n,2) } for(n=1;n<11;n++){ M(0,.2*n)L(5,.2*n)} fprintf(o,".005 W S\n");

DB s=1.e12; // teraJoule DB p=.02; // 2 percent M(0,2) for(n=1;n<10;n++) L(e[n]/s,p*p239[n]) fprintf(o,".01 W S\n"); M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p240[n]) fprintf(o,"1 0 0 RGB .01 W S\n"); M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p241[n]) fprintf(o,"0 .7 0 RGB .01 W S\n"); M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p242[n]) fprintf(o,"0 0 .9 RGB .01 W S\n"); M(0,0) for(n=1;n<10;n++) {if(p238[n]>0) L(e[n]/s,p*p238[n])} fprintf(o,".8 0 .9 RGB .01 W S\n");

fprintf(o,"0 0 0 RGB .01 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p239[n]) fprintf(o,"1 0 0 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p240[n]) fprintf(o,"0 .6 0 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p241[n]) fprintf(o,"0 0 .9 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p242[n]) fprintf(o,".8 0 .8 RGB .02 W\n");for(n=1;n<10;n++){if(p238[n]>0){O(e[n]/s,p*p238[n])}}

fprintf(o,"showpage\n"); fclose(o); system("epstopdf 21.eps"); system( "open 21.pdf"); }

Latex generator of labels
\documentclass[12pt]{article} \usepackage{geometry} \topmargin -96pt \oddsidemargin -70pt \paperwidth 518pt \paperheight 216pt \textwidth 1800pt \textheight 1720pt \usepackage{amsmath} \usepackage{graphicx} \usepackage{rotating} \parindent 0pt \newcommand{\sx}{\scalebox} \newcommand \ing {\includegraphics} \newcommand \rot {\begin{rotate}} \newcommand \ero {\end{rotate}} \begin{document} \begin{picture}(510,201) \put(10,10){\ing{21}} %\put(0,924){\sx{6}{Bu, \%}} \put(-2,204){\sx{1.1}{$C,$}} \put(1,188){\sx{1.1}{$\%$}} \put(-1,168){\sx{1.}{$80$}} \put(-1,128){\sx{1}{$60$}} \put(-1,88){\sx{1}{$40$}} \put(-1,48){\sx{1}{$20$}} \put(2, 8){\sx{1}{$0$}} \put(9, 1){\sx{1.1}{$0$}} \put(109, 1){\sx{1.1}{$1$}} \put(209, 1){\sx{1.1}{$2$}} \put(309, 1){\sx{1.1}{$3$}} \put(409, 1){\sx{1.1}{$4$}} %\put(406,-6){\sx{1.3}{$4$}} \put(454,1){\sx{1.1}{$B$, TJ/kg}} % \put(40,76){\sx{1.2}{\rot{90}\bf Magnox \ero}} \put(60,76){\sx{1.2}{\rot{90}\bf Magnox \ero}} \put(81,76){\sx{1.2}{\rot{90}\bf CANDU\ero}} \put(172,78){\sx{1.2}{\rot{90}\bf AGR\ero}} \put(254,72){\sx{1.2}{\rot{90}\bf BWR\ero}} \put(280,72){\sx{1.2}{\rot{90}\bf BWR\ero}} \put(303,72){\sx{1.2}{\rot{90}\bf PWR\ero}} \put(391,68){\sx{1.2}{\rot{90}\bf PWR\ero}} \put(478,68){\sx{1.2}{\rot{90}\bf PWR\ero}} % \put(403,111){\sx{1.3}{Pu-239}} \put(403,56){\sx{1.3}{Pu-240}} \put(403,39){\sx{1.3}{Pu-241}} \put(403,22){\sx{1.3}{Pu-242}} \put(403,10){\sx{1.3}{Pu-238}} % \end{picture} \end{document}