/home66/gary/public_html/cloudy/c08_branch/source/rt_line_one_tau_reset.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2008 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 "lines_service.h"
#include "prt.h"
#include "rt.h"

void RT_line_one_tau_reset(transition * t, 
  double WeightNew)
{
        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 );
        }

        /* iteration is 1 when this routine is called at end of first iteration
         * lgStatic is false by default, set true with command sphere SPHERE STATIC */
        if( iteration == 1 && (!geometry.lgStatic) )
        {
                /* end of first iteration, first 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 */
                t->Emis->TauTot = t->Emis->TauIn*rt.DoubleTau;
                t->Emis->TauIn = MIN2(opac.taumin,t->Emis->TauTot/2.f);
        }
        else if( geometry.lgSphere && geometry.lgStatic )
        {
                /* static sphere, both sides interact */
                if( iteration == 1 )
                {
                        /* on first iteration TauIn was 0 for first zone */
                        t->Emis->TauTot = 2.f*t->Emis->TauIn;
                }
                else
                {
                        t->Emis->TauTot = (realnum)(t->Emis->TauIn*WeightNew + t->Emis->TauTot*
                          (1. - WeightNew));
                }
                t->Emis->TauIn = t->Emis->TauTot/2.f;
        }
        else
        {
                /* take some sort of mean of old and new limiting optical depths */
                if( t->Emis->TauIn > 0. && t->Emis->TauTot > 0. )
                {
                        /* this is case where no maser occurs */
                        /* >>chng 03 nov 06, try taking linear average */
                        t->Emis->TauTot = (realnum)(t->Emis->TauTot*(1. - WeightNew) + 
                                t->Emis->TauIn*rt.DoubleTau*WeightNew);
                        t->Emis->TauIn = MIN2(opac.taumin,t->Emis->TauTot/ 2.f);
                }
                else
                {
                        /* case where maser did occr */
                        t->Emis->TauTot = (t->Emis->TauIn*rt.DoubleTau + t->Emis->TauTot)/2.f;
                        t->Emis->TauIn = MIN2(opac.taumin,t->Emis->TauTot/2.f);
                }
        }

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

        /* this is destruction probability
         * >>chng 96 sep 4, was not here, needed to add since now taking
         * mean of old and new dest prob */
        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;
}

---