rt_line_one_tau_reset.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 */
/*RT_line_one_tau_reset computes average of old and new optical depths for new scale at end of iter,
 * called by update, also FeIILevelPops */
#include "cddefines.h"
#include "opacity.h"
#include "geometry.h"
#include "transition.h"
#include "prt.h"
#include "rt.h"

void RT_line_one_tau_reset(const TransitionProxy& t)
{
        char chSave[11];

        DEBUG_ENTRY( "RT_line_one_tau_reset()" );

        if( t.ipCont() <= 0 )
        {
                return;
        }

        /* option to print masing lines, set with print maser */
        if( prt.lgPrtMaser && 
            ( t.Emis().TauTot() < -0.01 || t.Emis().TauIn() < -0.01 ) )
        {
                strcpy( chSave, chLineLbl(t) );
                fprintf( ioQQQ, " Masing line:%10.10s t(in, out)=%10.2e%10.2e\n", 
                  chSave, t.Emis().TauIn(), t.Emis().TauTot() );
        }

        realnum TauNext;

        /* lgStatic is false by default, set true with command sphere SPHERE STATIC */
        if( geometry.lgSphere && geometry.lgStatic )
        {
                /* static sphere, both sides interact */
                TauNext = 2.f*t.Emis().TauIn();
        }
        else
        {
                /* end of iteration, the next estimate of total optical depth 
                 * is now the inward optical depth - set T(ipLnTauTot) to T(1)
                 * DoubleTau normally 1, set to 2 by DoubleTau command in order
                 * to simulate two-sided photoionization */
                TauNext = rt.DoubleTau*t.Emis().TauIn();
        }

        /* iteration is 1 when this routine is called at end of first iteration
         * estimate is bad when starting first iteration, so don't track that one */
        if( iteration > 1 )
                TauNext = t.Emis().TauTrack().next_val( t.Emis().TauTot(), TauNext );

        if( geometry.lgSphere && geometry.lgStatic )
                t.Emis().TauIn() = TauNext/2.f;
        else
                t.Emis().TauIn() = min(opac.taumin,TauNext/2.f);

        t.Emis().TauTot() = TauNext;

        /* this is escape prob */
        t.Emis().Pesc() = 0.5f*(1.f + 1.f/MAX2(1.f,t.Emis().TauTot()));

        /* this is fraction inward */
        t.Emis().FracInwd() = MIN2(1.f,1.5f-t.Emis().Pesc());

        /*  destruction probability */
        t.Emis().Pdest() = 0.;
        t.Emis().Pelec_esc() = 0.;

        /* optical depth to the continuum source */
        t.Emis().TauCon() = opac.taumin;

        /* >>chng 01 sep 01, zero out some pops and energies */
        (*t.Lo()).Pop() = 0.;
        /* >>chng 97 jul 21, added following zero
         * population of upper level */
        (*t.Hi()).Pop() = 0.;
        /* population of lower level with correction for stim emission */
        t.Emis().PopOpc() = 0.;
        /* following two heat exchange excitation, deexcitation */
        t.Coll().cool() = 0.;
        t.Coll().heat() = 0.;
        /* intensity of line */
        t.Emis().xIntensity() = 0.;
        /* number of photons emitted in line */
        t.Emis().phots() = 0.;
        return;
}