Changeset 2250

Timestamp:

2008-07-23T21:38:40Z (15 years ago)

Author:

Gary J. Ferland

Message:

line pressure routine was moved to inline header but not removed from main source, so two copies existed.

File:

• branches/c08_branch/source/pressure_total.cpp (modified) (1 diff)

branches/c08_branch/source/pressure_total.cpp

@@ -1002 +1002 @@
         return;
 }
-
-double PressureRadiationLine( const transition *t, double densityDenorm )
-{
-        double PresssureReturned;
-        DEBUG_ENTRY( "PressureRadiationLine()" );
-
-        /* RT_LineWidth gets line width in terms of RT effects */
-        double width = RT_LineWidth(t);
-        long nelem = t->Hi->nelem-1;
-
-        double PopOpc = t->Emis->PopOpc/t->Lo->g;
-        /* return zero line radiation PresssureReturned if line mases or has
-         * zero opacity */
-        if( densityDenorm*PopOpc*t->Emis->opacity <= SMALLFLOAT )
-                return 0;
-
-        PresssureReturned = PI8 * HPLANCK / 3. * (powi(t->EnergyWN,4)) *
-                (t->Hi->Pop/t->Hi->g)/PopOpc * width;
-
-        /* this prevents line radiation PresssureReturned 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 = densityDenorm * t->Emis->PopOpc * t->Emis->opacity / DoppVel.doppler[nelem]/
-                        rfield.fine_opac_zone[ipLineCenter];
-                /* 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 );
-                PresssureReturned *= FractionThisLine;
-        }
-
-        return PresssureReturned;
-}