pressure.h

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 */

#ifndef PRESSURE_H_
#define PRESSURE_H_

#include "rt.h"
#include "rfield.h"
#include "doppvel.h"
#include "physconst.h"
#include "transition.h"

void PresTotCurrent(void);

inline double PressureRadiationLine( const TransitionProxy &t, realnum DopplerWidth )
{
        DEBUG_ENTRY( "PressureRadiationLine()" );

        /* return zero if below plasma frequency */     
        if( t.EnergyErg() / EN1RYD <= rfield.plsfrq )
                return 0.;

        /* RT_LineWidth gets line width in terms of RT effects */
        double width = RT_LineWidth(t, DopplerWidth);

        double PopOpc = t.Emis().PopOpc()/(*t.Lo()).g();
        /* return zero line radiation PressureReturned if line mases or has
         * zero opacity */
        /* \todo 2 1e-22 is arbtrary but roughly 1/kpc.  Replace with a cloud width if available? */
        if( PopOpc*t.Emis().opacity()/ DopplerWidth <= 1.e-22 || width<=0. )
                return 0.;

        double PressureReturned = PI8 * HPLANCK / 3. * POW4(t.EnergyWN()) *
                ((*t.Hi()).Pop()/(*t.Hi()).g())/PopOpc * width;

        /* this prevents line radiation PressureReturned from being very large when line 
         * is not optically thick but total opacity at that energy is large 
         * due to overlapping transitions */
        long int ipLineCenter = t.Emis().ipFine() + rfield.ipFineConVelShift;
        if( ipLineCenter > 0 && ipLineCenter < rfield.nfine && rfield.lgOpacityFine &&
                rfield.fine_opac_zone[ipLineCenter] > SMALLFLOAT )
        {
                double FractionThisLine = t.Emis().PopOpc() * t.Emis().opacity() / DopplerWidth/
                        rfield.fine_opac_zone[ipLineCenter];
                if( FractionThisLine<1e-5 )
                        FractionThisLine = 0.;
                /* fine opacities are only reevaluated one time per zone due
                 * to the expense - PopOpc is for the current solution - but the two
                 * may be out of step by a few percent, due to the variation in
                 * abundance from zone to zone.  This prevents the change
                 * in solution from increasing the radiation pressure.
                 * This correction is mainly an order of magnitude scaler to prevent
                 * optically thin lines from appearing to be optically thick due to
                 * overlapping lines */
                FractionThisLine = MIN2(1., FractionThisLine);
                ASSERT( FractionThisLine >= 0. && FractionThisLine <= 1.0 );
                PressureReturned *= FractionThisLine;
        }

        return PressureReturned;
}

struct t_pressure {

        realnum PresLow, 
          PresHigh;

        realnum PresPowerlaw;

        double PresRamCurr;

        double PresTurbCurr;

        double PresTotlCurr, 
          PresTotlError,
          PresGasCurr;

        double PresTotlInit;

        bool lgPressureInitialSpecified;
        double PressureInitialSpecified;

        double pres_radiation_lines_curr;

        bool lgContRadPresOn;

        realnum PresInteg, 
          pinzon;

        realnum PresIntegElecThin,
                pinzon_PresIntegElecThin;

        double RhoGravity_dark;
        double RhoGravity_self;
        double RhoGravity_external;
        double RhoGravity;
        double IntegRhoGravity;
        int gravity_symmetry;
        double self_mass_factor;
        vector<double> external_mass[3];

        bool lgPres_radiation_ON;
        bool lgPres_magnetic_ON;
        bool lgPres_ram_ON;

        realnum 
          RadBetaMax, 
          pbeta, 
          PresMax;

        long int ipPradMax_line;

        long int ipPradMax_nzone;

        char chLineRadPres[101];

        bool lgPradCap, 
          lgPradDen;

        bool lgLineRadPresOn;

        bool lgRadPresAbortOK;

        bool lgSonicPointAbortOK; 

        bool lgSonicPoint;

        bool lgStrongDLimbo;

        };
extern t_pressure pressure;


#endif /* PRESSURE_H_ */