dense_fabden.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*dense_fabden called by dlaw command, returns density for any density law */
#include "cddefines.h"
#include "rfield.h"
#include "dense.h"
 
/*dense_fabden called by dlaw command, returns density for any density law */
double dense_fabden(double radius, 
  double depth)
{
        double fabden_v = pow(10.,dense.DensityLaw[0]);
        if( rfield.lgUSphON )
                fabden_v *= pow(radius/rfield.rstrom,dense.DensityLaw[1]);
        else
                fabden_v *= pow(depth/rfield.rstrom,dense.DensityLaw[1]);
 
        return fabden_v; 
}