taulines.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 */
#include "cddefines.h"
#include "taulines.h"

INSTANTIATE_MULTI_ARR( qList, lgBOUNDSCHECKVAL );

vector<TransitionList> AllTransitions;

bool lgStatesAdded;
bool lgLinesAdded;
multi_arr<qList,2> StatesElemNEW;
char **chSpecies;
species *dBaseSpecies;
vector<qList> dBaseStates;
vector< multi_arr<int,2> > ipdBaseTrans;
vector<TransitionList> dBaseTrans;
multi_arr<CollRateCoeffArray,2> AtmolCollRateCoeff;
CollSplinesArray ****AtmolCollSplines;
StoutColls ****StoutCollData;
long int nSpecies;
qList AnonStates(1);
TransitionList TauLines("TauLines", &AnonStates);
multi_arr<int,3> ipExtraLymanLines;
vector<vector<TransitionList> > ExtraLymanLines;
long int nUTA;
TransitionList UTALines("UTALines", &AnonStates);
long int nLevel1;
TransitionList HFLines("HFLines", &AnonStates);
long int nHFLines;
//vector<vector<multi_arr<int,2> > > ipTransitions;
vector<vector<TransitionList> > Transitions;
multi_arr<int,2> ipFe2LevN;
static qList Fe2LevNStates;
TransitionList Fe2LevN("Fe2LevN", &Fe2LevNStates);
multi_arr<int,3> ipSatelliteLines; /* [ipISO][nelem][level] */
vector<vector<TransitionList> > SatelliteLines; /* [ipISO][nelem][level] */
TransitionList TauLine2("TauLine2", &AnonStates);
realnum *cs1_flag_lev2;

extern void checkTransitionListOfLists(vector<TransitionList>&list)
{ 
        for (vector<TransitionList>::iterator it=list.begin(); 
                  it != list.end(); ++it)
        {
                for (TransitionList::iterator tr = it->begin();
                          tr != it->end(); ++tr)
                {
                        (*tr).check();
                }
                for (EmissionList::iterator em = it->Emis().begin(); 
                          em != it->Emis().end(); ++em)
                {
                        (*em).check();
                }
        }
}

TransitionProxy::iterator TauDummy;
EmissionProxy DummyEmis;

namespace 
{
        class Init
        {
                EmissionList DummyEmisList;
                TransitionListImpl TauDummyTrans;
        public:
                Init(qList*states) : 
                        DummyEmisList(&TauDummyTrans, 1), TauDummyTrans("TauDummy",states, 1)
                {
                        DummyEmis = DummyEmisList[0];
                        TauDummy=TauDummyTrans.begin();
                };
        };
        qList TauDummyStates(1);
        Init TauDummyInit(&TauDummyStates);;
}