dynamics.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2023 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */

#ifndef DYNAMICS_H_
#define DYNAMICS_H_

#include "module.h"
#include "container_classes.h"

void DynaIterStart(void);

void DynaIterEnd(void);

void DynaStartZone(void);

void DynaEndZone(void);

void DynaIonize(void);

void DynaCreateArrays( void );

class Parser;
void ParseDynaWind( Parser &p );

void ParseDynaTime( Parser &p );

void DynaPrtZone( void );

void DynaSave(FILE* ipPnunit , char chJob );

void DynaPunchTimeDep( FILE* ipPnunit , const char *chJob );

realnum DynaFlux(double depth);

struct t_dynamics : public module
{
        const char *chName() const
        {
                return "dynamics";
        }
        void zero();
        void comment(t_warnings&) {}

        bool isInitialRelaxIteration( long int iteration );

        bool doNonEquilibriumSolve( long int iteration );


        bool lgAdvection;

        double Cool_r, Heat_v, dHeatdT;

        double Cool(), Heat(), dCooldT();

        double CoolMax, HeatMax;

        double Rate;

        multi_arr<double,2> Source;

        multi_arr<double,3> StatesElem;

        vector<double> molecules;

        bool lgISO[NISO];

        bool lgMETALS;

        bool lgCoolHeat;

        bool lg_coronal_time_init;

        bool lgTimeDependentStatic;

        double time_elapsed;

        bool lgRecom;

        bool lgStatic_completed;

        double AdvecLengthInit;

        double FluxCenter;

        char chPresMode[20];

        double ShockDepth;

        double ShockMach;

        long int n_initial_relax;

        double FluxScale;

        bool lgFluxDScale;

        double FluxIndex;

        double dRad;

        double oldFullDepth;

        /* the error from comparing this iteration with the previous one */
        double convergence_error;

        double convergence_tolerance;

        double discretization_error;

        double error_scale1 , error_scale2;

        bool lgSetPresMode;

        realnum Upstream_density;

        realnum DivergePresInteg;

        bool lgEquilibrium;

        /* set true with trace option on time command */
        bool lgTracePrint;

         /* initial timestep (seconds) read off time command,
          * each iteration accounts for this much time */
        double timestep_init,
                 timestep,
                 timestep_stop,
                 timestep_factor;
};
extern t_dynamics dynamics;

#endif /* DYNAMICS_H_ */