JSeaWater.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <cmath>
 
#include "JPhysics/JSeaWater.hh"
 
#include "JLang/JManip.hh"
 
#include "Jeep/JParser.hh"
 
 
namespace {
 
  /**
   * Write atomic parameters of sea water to output stream.
   *
   * \param  out                output stream
   * \param  object             atomic parameters
   * \return                    output stream
   */
  inline std::ostream& operator<<(std::ostream& out, const JPHYSICS::JSeaWater::atom_type& object)
  {
    using namespace std;
 
    return out << setw(2) << left << object.ps << right << ' '
               << FIXED(5,1) << object.Z    << ' '
               << FIXED(5,1) << object.A    << ' '
               << FIXED(7,4) << object.f;
  }
}
 
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
 
  double precision;
  int    debug;
 
  try {
 
    JParser<> zap;
 
    zap['p'] = make_field(precision)    = 1.0e-10;
    zap['d'] = make_field(debug)        = 3;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  {
    double F = 0.0;
 
#ifndef RADIATION
    const auto buffer = { JSeaWater::O, JSeaWater::H, JSeaWater::Na, JSeaWater::Cl };
#else
    const auto buffer = { JSeaWater::O, JSeaWater::H, JSeaWater::Na, JSeaWater::Cl, JSeaWater::Ca, JSeaWater::Mg, JSeaWater::K, JSeaWater::S };
#endif
    for (const auto& i : buffer) {
 
      cout << i << endl;
 
      F += i.f;
    }
 
    ASSERT(fabs(F - 1.0) <= precision, "Test total mass fraction " << SCIENTIFIC(12,3) << (F - 1.0));
  }
}