JSeaWater.hh

Go to the documentation of this file.

#ifndef __JPHYSICS__SEAWATER__
#define __JPHYSICS__SEAWATER__
 
/**
 * \author mdejong
 */
 
namespace JPHYSICS {}
namespace JPP { using namespace JPHYSICS; }
 
namespace JPHYSICS {
 
  /**
   * Sea water composition.
   */
  struct JSeaWater 
  {
    struct atom_type {
 
      constexpr double operator()() const { return f; }
 
      const char* const ps;    //!< short name
      double            Z;     //!< electric charge
      double            A;     //!< mass number
      double            f;     //!< mass fraction  
    };
 
#ifndef RADIATION
    static constexpr atom_type H   =  { "H",   1.0,  1.0, 0.10870 };
    static constexpr atom_type O   =  { "O",   8.0, 16.0, 0.85890 };
    static constexpr atom_type Na  =  { "Na", 11.0, 23.0, 0.01180 };
    static constexpr atom_type Cl  =  { "Cl", 17.0, 35.5, 0.02060 };
#else
    static constexpr atom_type H   =  { "H",   1.0,  1.0, 0.10790 };
    static constexpr atom_type O   =  { "O",   8.0, 16.0, 0.85630 };
    static constexpr atom_type Na  =  { "Na", 11.0, 23.0, 0.01160 };
    static constexpr atom_type Mg  =  { "Mg", 12.0, 24.3, 0.00140 };
    static constexpr atom_type S   =  { "S",  16.0, 32.0, 0.00100 };
    static constexpr atom_type Cl  =  { "Cl", 17.0, 35.5, 0.02090 };
    static constexpr atom_type K   =  { "K",  19.0, 39.0, 0.00045 };
    static constexpr atom_type Ca  =  { "Ca", 20.0, 40.0, 0.00045 };
#endif
  };
}
 
#endif