Ticket #67: dense_fabden.cpp
| File dense_fabden.cpp, 511 bytes (added by peter, 4 months ago) |
|---|
| Line | |
|---|---|
| 1 | /* This file is part of Cloudy and is copyright (C)1978-2007 by Gary J. Ferland |
| 2 | * For conditions of distribution and use see copyright notice in license.txt */ |
| 3 | /*dense_fabden called by dlaw command, returns density for any density law */ |
| 4 | #include "cddefines.h" |
| 5 | #include "rfield.h" |
| 6 | #include "dense.h" |
| 7 | |
| 8 | double dense_fabden(double radius, |
| 9 | double depth) |
| 10 | { |
| 11 | double fabden_v = pow(10.,dense.DensityLaw[0]); |
| 12 | if( rfield.lgUSphON ) |
| 13 | fabden_v *= pow(radius/rfield.rstrom,dense.DensityLaw[1]); |
| 14 | return fabden_v; |
| 15 | } |
