32 if( strcmp( chType ,
"PRETTY")==0 )
37 fprintf( ioMEAN,
"\n" );
39 fprintf( ioMEAN,
"\n " );
40 for( i=1; i <= 17; i++ )
42 fprintf( ioMEAN,
"%7ld", i );
44 fprintf( ioMEAN,
"\n\n" );
47 for( nelem=0; nelem <
LIMELM; nelem++ )
49 if(
dense.lgElmtOn[nelem] )
53 fprintf( ioMEAN,
" %10.10s",
elementnames.chElementName[nelem] );
82 fprintf( ioMEAN,
"\n" );
84 fprintf( ioMEAN,
"%7.3f", aa );
89 if( i >
MAX2(3,nelem+2) )
94 fprintf(ioMEAN,
" (H2) Log10 Column density (cm^-2)");
97 fprintf( ioMEAN,
"\n" );
107 const int nEXCIT_COL = 12;
109 string strExcit_Col[nEXCIT_COL]={
110 "He1*",
"CII*",
"C11*",
"C12*",
"C13*",
"O11*",
"O12*",
"O13*",
"Si2*",
"C30*",
"C31*",
"C32*"};
113 fprintf(ioMEAN,
" Exc state ");
115 for(i=0; i<nEXCIT_COL; ++i )
117 ASSERT( strExcit_Col[i].size() <= 4 );
122 strExcit_Col[i].c_str(),
131 fprintf(ioMEAN,
"\n ");
134 fprintf(ioMEAN,
" %s%7.3f",
135 strExcit_Col[i].c_str(),
139 fprintf(ioMEAN,
"\n");
143 h2.H2_Prt_column_density(ioMEAN);
145 fprintf(ioMEAN,
"\n");
147 else if( strcmp( chType ,
"TABLE")==0 )
150 if(
save.lgPunHeader[ioPun] )
152 fprintf(ioMEAN,
"#Column densities %s",
156 fprintf(ioMEAN,
"\t%s",
159 fprintf(ioMEAN,
"\n");
160 save.lgPunHeader[ioPun] =
false;
162 fprintf(ioMEAN,
"%.3e",
163 mole.species[0].column );
166 fprintf(ioMEAN,
"\t%.3e",
167 mole.species[i].column );
169 fprintf(ioMEAN,
"\n");
NORETURN void TotalInsanity(void)
sys_float SDIV(sys_float x)
#define DEBUG_ENTRY(funcname)
int cdColm(const char *chLabel, long int ion, double *theocl)
t_elementnames elementnames
diatomics h2("h2", 4100., &hmi.H2_total, Yan_H2_CS)
void molcol(const char *chLabel, FILE *ioMEAN)
t_mole_global mole_global
void PrtColumns(FILE *ioMEAN, const char *chType, long int ioPun)