Skip to content

Commit

Permalink
fix formula
Browse files Browse the repository at this point in the history
  • Loading branch information
bam241 committed Mar 17, 2017
1 parent 16ef3bd commit 4f00b59
Showing 1 changed file with 7 additions and 3 deletions.
10 changes: 7 additions & 3 deletions src/special_enrich.cc
Original file line number Diff line number Diff line change
Expand Up @@ -481,18 +481,22 @@ cyclus::Material::Ptr SEnrichment::Enrich_(cyclus::Material::Ptr mat,
// Re-Add the special nuc inside the product, and fix tails amount accordingly

double prod_mass = response->quantity();
double u5_raw_enrich = UraniumAssay(response);

cyclus::toolkit::MatQuery mq_resp(response);
double u5_raw_enrich = mq_resp.mass_frac(922350000);
cyclus::toolkit::MatQuery mq_flip(flip_mat.first);
double u5_flip_enrich = mq_flip.mass_frac(922350000);


std::map<cyclus::Nuc, double>::iterator it;
for (it = ux.begin(); it != ux.end(); it++) {

double nuc_i_enrich_factor =
u5_raw_enrich / UraniumAssay(flip_mat.first) * it->second;
u5_raw_enrich / u5_flip_enrich * it->second;

cyclus::toolkit::MatQuery mq_(natu_matl);
double nuc_i_feed_enrich = mq_.mass(it->first) / natu_matl->quantity();
double nuc_i_prod_enrich = nuc_i_feed_enrich * nuc_i_enrich_factor;

double nuc_i_prod_mass = nuc_i_prod_enrich * prod_mass;

// Remove come ux from the material pushed in the tails and add in into the
Expand Down

0 comments on commit 4f00b59

Please sign in to comment.