27 vector<long>ipsave(LIMELM*LIMELM, INT_MIN);
28 vector<long>jpsave(LIMELM*LIMELM, INT_MIN);
29 vector<long>ipOrdered(LIMELM*LIMELM);
30 vector<string> chLabel(LIMELM*LIMELM);
61 for(
long i=0; i <
LIMELM; i++ )
63 for(
long j=0; j <
LIMELM; j++ )
87 fprintf( io,
" Negative coolant was %s %.2f %.2e\n",
94 fprintf( io,
"#>>>> Temperature not converged.\n" );
98 fprintf( io,
"#>>>> Electron density not converged.\n" );
102 fprintf( io,
"#>>>> Ionization not converged.\n" );
106 fprintf( io,
"#>>>> Pressure not converged.\n" );
114 for(
long k=0; k < ipnt; k++ )
125 fprintf(
ioQQQ,
"DISASTER assert about to be thrown - search for hit it\n");
135 else if( i == 0 && j == 1 )
141 else if( i == 0 && j == 3 )
147 else if( i == 0 && j == 7 )
152 else if( i == 0 && j == 8 )
158 else if( i == 0 && j == 17 )
165 else if( i == 0 && j == 9 )
170 else if( i == 0 && j == 20 )
175 else if( i == 0 && j == 21 )
180 else if( i == 0 && j == 11 )
185 else if( i == 0 && j == 12 )
190 else if( i == 0 && j == 13 )
195 else if( i == 0 && j == 14 )
200 else if( i == 0 && j == 15 )
205 else if( i == 0 && j == 16 )
210 else if( i == 0 && j == 18 )
215 else if( i == 0 && j == 19 )
220 else if( i == 0 && j == 22 )
225 else if( i == 0 && j == 23 )
231 else if( i == 0 && j == 24 )
236 else if( i == 1 && j == 3 )
241 else if( i == 1 && j == 5 )
246 else if( i == 1 && j == 6 )
251 else if( i == 25 && j == 27 )
259 oss <<
"[" << i <<
"][" << j <<
"]";
260 chLabel[k] = oss.str();
282 fprintf( io,
"%.5e\t%.4e\t%.4e\t%.4e",
291 for(
long k=0; k < ipnt; k++ )
293 int ip = ipOrdered[k];
296 ASSERT( i<LIMELM && j<LIMELM );
300 chLabel[ip].c_str(), SaveVal[ip] );
306 bool lgHeatLine =
false;
309 for( i=0; i < ipnt; i++ )
312 if( ipsave[ipOrdered[i]] == 0 && jpsave[ipOrdered[i]] == 22 )
331 fprintf( io,
" LHeat lv%2ld %s TIn%10.2e Pmp%9.1e EscP%9.1e CS%9.1e Hlin/tot%10.2e\n",
332 level,
chLineLbl(t).c_str(), TauIn, Pump, EscP, CS, ColHeat );
string chLineLbl(const TransitionProxy &t)
char chIonStage[LIMELM+1][CHARS_ION_STAGE]
t_elementnames elementnames
double xIonDense[LIMELM][LIMELM+1]
char chClntLab[NCOLNT][NCOLNT_LAB_LEN+1]
EmissionList::reference Emis() const
const TransitionProxy FndLineHt(long int *level)
sys_float safe_div(sys_float x, sys_float y, sys_float res_0by0)
double heating(long nelem, long ion)
char chElementSym[LIMELM][CHARS_ELEMENT_SYM]
realnum & col_str() const
CollisionProxy Coll() const
#define DEBUG_ENTRY(funcname)
int fprintf(const Output &stream, const char *format,...)
void spsort(realnum x[], long int n, long int iperm[], int kflag, int *ier)