//
// $Id: models71to80.cpp,v 1.1 2010/09/01 19:49:05 thernis Exp $
//

#include "models71to80.h"
#include <cmath>


float CModel71::Density(const Cvec &v) {
  float r=v.norm();
  if (r <= 1.05) return 0.;

  // -- see article SolPhy 183: 165-180, 1998
/*  float a=3.3e5;
  float b=4.1e6;
  float c=8.0e7;*/
  
  float nel=3.3e5*pow(r,-2)+4.1e6*pow(r,-4)+8.0e7*pow(r,-6);

  return nel;
}
void CModel71::dumpDefaultParamForIDL(std::vector<moddefparam>& vp,int& flagcase)
{   
    flagcase=0x4;
    vp.push_back(moddefparam("","Leblanc, Dulk, Bougeret model","",""));
    return;
}