cool_nick.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 */
/*CoolNick compute nickel cooling */
#include "cddefines.h"
#include "taulines.h"
#include "phycon.h"
#include "lines_service.h"
#include "dense.h"
#include "atoms.h"
#include "cooling.h"

void CoolNick(void)
{
        double tused;
        realnum rate;

        DEBUG_ENTRY( "CoolNick()" );

        /*>>refer       Ni1     cs      Hollenbach, D. & McKee, C.F. 1989, ApJ, 342, 306 */
        /* rates are said to be ok over range 30 - 3000K */
        tused = MAX2( 30. , phycon.te );
        tused = MIN2( 3000. , phycon.te  );
        tused /= 100.;

        /* the 7.5 micron line */
        /* >>chng 03 nov 15, add these lines */
        rate = (realnum)(1.2e-7 * dense.eden + 
                /*8.0e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
                /* >>chng 05 jul 05, eden to cdsqte */
                8.0e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
        LineConvRate2CS( TauLines[ipNi1_7m]  , rate );

        /* the 11.3 micron line */
        rate = (realnum)(9.3e-8 * dense.eden + 
                /* >>chng 05 jul 05, eden to cdsqte */
                /*5.3e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
                5.3e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
        LineConvRate2CS( TauLines[ipNi1_11m]  , rate );

        rate = (realnum)(1.2e-7 * dense.eden + 
                /* >>chng 05 jul 05, eden to cdsqte */
                /*6.9e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
                6.9e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
        (*(*TauDummy).Hi()).g() = (*TauLines[ipNi1_11m].Hi()).g();
        LineConvRate2CS( *TauDummy  , rate );
        /* this says that line is a dummy, not real one */
        (*(*TauDummy).Hi()).g() = 0.;

        atom_level3(TauLines[ipNi1_7m],TauLines[ipNi1_11m],*TauDummy);

        return;
}