56 SaveVal[i] = -FLT_MAX;
57 chLabel[i] = (
char *)
CALLOC(10,
sizeof(
char));
88 for( i=0; i <
LIMELM; i++ )
90 for( j=0; j <
LIMELM; j++ )
105 for( i=0; i <
thermal.ncltot; i++ )
114 fprintf( io,
" Negative coolant was %s %.2f %.2e\n",
119 if( !
conv.lgConvTemp )
121 fprintf( io,
"#>>>> Temperature not converged.\n" );
123 else if( !
conv.lgConvEden )
125 fprintf( io,
"#>>>> Electron density not converged.\n" );
127 else if( !
conv.lgConvIoniz() )
129 fprintf( io,
"#>>>> Ionization not converged.\n" );
131 else if( !
conv.lgConvPres )
133 fprintf( io,
"#>>>> Pressure not converged.\n" );
141 for( k=0; k < ipnt; k++ )
152 fprintf(
ioQQQ,
"DISASTER assert about to be thrown - search for hit it\n");
162 else if( i == 0 && j == 1 )
166 strcpy( chLabel[k],
"Hn=2" );
168 else if( i == 0 && j == 3 )
172 strcpy( chLabel[k],
"Hion" );
174 else if( i == 0 && j == 7 )
177 strcpy( chLabel[k],
" UTA" );
179 else if( i == 0 && j == 8 )
183 strcpy( chLabel[k],
"H2vH" );
185 else if( i == 0 && j == 17 )
190 strcpy( chLabel[k],
"H2dH" );
192 else if( i == 0 && j == 9 )
195 strcpy( chLabel[k],
"COds" );
197 else if( i == 0 && j == 20 )
200 strcpy( chLabel[k],
"extH" );
202 else if( i == 0 && j == 21 )
205 strcpy( chLabel[k],
"pair" );
207 else if( i == 0 && j == 11 )
210 strcpy( chLabel[k],
"H FF" );
212 else if( i == 0 && j == 12 )
215 strcpy( chLabel[k],
"Hcol" );
217 else if( i == 0 && j == 13 )
220 strcpy( chLabel[k],
"GrnP" );
222 else if( i == 0 && j == 14 )
225 strcpy( chLabel[k],
"GrnC" );
227 else if( i == 0 && j == 15 )
230 strcpy( chLabel[k],
"H- " );
232 else if( i == 0 && j == 16 )
235 strcpy( chLabel[k],
"H2+ " );
237 else if( i == 0 && j == 18 )
240 strcpy( chLabel[k],
"H2ph" );
242 else if( i == 0 && j == 19 )
245 strcpy( chLabel[k],
"Comp" );
247 else if( i == 0 && j == 22 )
250 strcpy( chLabel[k],
"line" );
252 else if( i == 0 && j == 23 )
256 strcpy( chLabel[k],
"Hlin" );
258 else if( i == 0 && j == 24 )
261 strcpy( chLabel[k],
"ChaT" );
263 else if( i == 1 && j == 3 )
266 strcpy( chLabel[k],
"He3l" );
268 else if( i == 1 && j == 5 )
271 strcpy( chLabel[k],
"adve" );
273 else if( i == 1 && j == 6 )
276 strcpy( chLabel[k],
"CR H" );
278 else if( i == 25 && j == 27 )
281 strcpy( chLabel[k],
"Fe 2" );
285 sprintf( chLabel[k],
"[%ld][%ld]" , i , j );
307 fprintf( io,
"%.5e\t%.4e\t%.4e\t%.4e",
316 for( k=0; k < ipnt; k++ )
318 int ip = ipOrdered[k];
322 if(k > 4 &&
thermal.heating[i][j]/
SDIV(heat_total) <
save.WeakHeatCool )
324 fprintf( io,
"\t%s\t%.7f ",
325 chLabel[ip], SaveVal[ip] );
327 fprintf( io,
" \n" );
334 for( i=0; i < ipnt; i++ )
337 if( ipsave[i] == 0 && jpsave[i] == 22 )
357 fprintf( io,
" LHeat lv%2ld %10.10s TIn%10.2e Pmp%9.1e EscP%9.1e CS%9.1e Hlin/tot%10.2e\n",
358 level, chLbl, TauIn, Pump, EscP, CS, ColHeat );
void spsort(realnum x[], long int n, long int iperm[], int kflag, int *ier)
sys_float SDIV(sys_float x)
#define DEBUG_ENTRY(funcname)
realnum & col_str() const
CollisionProxy Coll() const
EmissionList::reference Emis() const
t_elementnames elementnames
const TransitionProxy FndLineHt(long int *level)
char * chLineLbl(const TransitionProxy &t)