/home66/gary/public_html/cloudy/c08_branch/source/rfield.h

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

#ifndef _RFIELD_H_
#define _RFIELD_H_

/* rfield.h */

#define WL_V_FILT       5500.

#define WL_B_FILT       4400.

void rfield_opac_zero( long lo , long ihi );

extern bool lgRfieldMalloced;

EXTERN struct t_rfield {

        long int nflux;

        long int nupper;

        realnum emm;

        realnum egamry;

        realnum FluxFaint;

        realnum *anu;

        realnum *AnuOrg;

        realnum *widflx;

        long int *line_count;

        realnum *OccNumbBremsCont;

        realnum *OccNumbContEmitOut;

        realnum *anulog, 
          *anusqr, 
          *anu2, 
          *anu3;

        realnum *flux;

        realnum *flux_isotropic;

        realnum *flux_beam_time , *flux_beam_const;

        realnum *flux_accum;

        bool lgMustBlockHIon;

        /* this is set true if H-ionizing radiation is blocked with extinguish
         * command */
        bool lgBlockHIon;

        bool lgDoLineTrans;

        bool lgOpacityReevaluate;

        bool lgCMB_set;

        bool lgIonizReevaluate;

        realnum *convoc;

        realnum *OccNumbIncidCont;

        realnum *OccNumbDiffCont;

        double *ContBoltz;

        realnum **ConEmitLocal/* [depth][energy]*/;

        realnum *ConEmitReflec;

        realnum *ConEmitOut;

        realnum *ConInterOut;

        realnum *ConRefIncid;

        realnum *SummedCon;
        realnum *SummedDif;
        realnum *SummedOcc;
        realnum *SummedDifSave;

        realnum *ConOTS_local_photons,
                *ConOTS_local_OTS_rate;

        realnum *TotDiff2Pht;

        realnum 
                *otslin/*[NC_ELL]*/, 
                *otscon/*[NC_ELL]*/, 

                 *otslinNew/*[NC_ELL]*/,
                 *otsconNew/*[NC_ELL]*/,
                **otssav/*[NC_ELL][2]*/;

        realnum *outlin,
                *outlin_noplot;

        realnum *reflin/*[NC_ELL]*/;

        realnum *flux_total_incident;
        realnum *flux_beam_const_save , *flux_time_beam_save , *flux_isotropic_save;

        realnum time_continuum_scale;

        realnum DiffPumpOn;

        char **chLineLabel/*[NC_ELL][5]*/;

        char **chContLabel/*[NC_ELL][5]*/;

        realnum **gff/*[LIMELM][NC_ELL]*/;

        bool lgCompileGauntFF;

        char chDffTrns[4];

        bool lgOutOnly;

        long int ipEnergyBremsThin;
        realnum EnergyBremsThin;

        long int ipMaxBolt;

        bool lgInducProcess;

        double *comup, 
          *comdn;

        long int ipB_filter , ipV_filter;

        long int ipG0_TH85_lo , ipG0_TH85_hi;

        long int ipG0_DB96_lo , ipG0_DB96_hi;

        long int ipG0_spec_lo , ipG0_spec_hi;

        long int ip1000A;

        double extin_mag_B_point , extin_mag_V_point;

        double extin_mag_B_extended , extin_mag_V_extended;

        double opac_mag_B_point, opac_mag_V_point, opac_mag_B_extended , opac_mag_V_extended;

        realnum *csigh, 
          *csigc;

        double comtot, 
          cmheat, 
          cmcool, 
          comoff, 
          cinrat;

        bool lgComUndr;

#       define  LIMSPC  100

#       define  NCELL   130000

        double totpow[LIMSPC], 
          slope[LIMSPC], 
          cutoff[LIMSPC][3], 
          spfac[LIMSPC];

        bool lgTimeVary[LIMSPC];

        /* beamed or isotropic continuum?  if isotropic then does not vary
         * with time */
        bool lgBeamed[LIMSPC];

        long int nspec, 
          ipspec;

        realnum **tNuRyd  , 
          **tslop , 
          **tFluxLog ;

        int *lgContMalloc;

        double range[LIMSPC][2];

        realnum *ConTabRead;

        char chSpNorm[LIMSPC][5], 
          chRSpec[LIMSPC][5], 
          chSpType[LIMSPC][6];

        string ioTableRead[LIMSPC];

        realnum qhtot, 
          qhe, 
          qheii, 
          qbal, 
          qrad, 
          qtot;

        realnum uh;

        realnum uheii;

        bool lgUSphON;
        realnum rstrom;

        bool lgHabing;

        long int fine_opac_nelem;

        long int fine_opac_nresolv;

        realnum fine_opac_velocity_width;

        realnum *trans_coef_zone;
        realnum *trans_coef_total;

        long int *ipnt_coarse_2_fine;

        realnum fine_ener_lo, fine_ener_hi;
        long nfine;
        long nfine_malloc;
        double fine_resol;
        realnum *fine_opac_zone;
        realnum *fine_opt_depth;
        realnum *fine_anu;

        long int ipFineConVelShift;

        bool lgOpacityFine;
        bool lgPunchOpacityFine;

        bool lgFine_opac_update;

        bool lgLyaOTS;

        bool lgHeIIOTS;

        bool lgKillOutLine;

        bool lgKillOutCont;

        bool lgKillOTSLine;

        bool lgTableRead;

        bool lgPlasNu;
        realnum plsfrq, 
                plsfrqmax;
        long int ipPlasma, 
          ipPlasmax;

        bool lgMMok, 
          lgHPhtOK, 
          lgXRayOK, 
          lgGamrOK;

        realnum EnerGammaRay;
        long int ipEnerGammaRay;

        bool lgHionRad;

        realnum occmax, 
          occmnu, 
          tbrmax, 
          tbrmnu, 
          tbr4nu, 
          occ1nu;

        bool lgOcc1Hi;

}       rfield;


#endif /* _RFIELD_H_ */

---