46 fprintf(
ioQQQ,
" punching line data, then stopping\n" );
53 fprintf(
ioQQQ,
" WARNING - I am punching line data but element %s is turned off.\n",
60 fprintf(
ioQQQ,
"Some elements are turned off and punch line data requested.\n");
61 fprintf(
ioQQQ,
"Code is now designed to do punch line data only with all elements on.\n");
62 fprintf(
ioQQQ,
"Please try again with all elements on.\n");
81 fprintf( ioPUN,
"\n Recombination lines of C, N, O\n" );
82 fprintf( ioPUN,
" Ion WL(A) Coef Ion WL(A) Coef\n" );
83 for( i=0; i<471; i+=nskip)
86 limit =
MIN2(471,i+nskip);
87 fprintf( ioPUN,
" " );
88 for( j=i; j < limit; j++ )
90 fprintf( ioPUN,
"%2.2s%2ld%6ld%8.3f ",
96 fprintf( ioPUN,
" \n" );
98 fprintf( ioPUN,
"\n\n" );
118 for( i=0; i <
nUTA; i++ )
123 for( i=0; i <
LIMELM; i++ )
125 for( j=0; j < LIMELM+1; j++ )
134 fprintf( ioPUN,
" Level 1 transferred lines\n" );
136 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
147 fprintf( ioPUN,
"\n\n\n" );
148 fprintf( ioPUN,
" end level 1, start level 2\n" );
149 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
158 fprintf( ioPUN,
"\n\n\n" );
159 fprintf( ioPUN,
" end level 2, start inner shell\n" );
160 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
162 for( i=0; i <
nUTA; i++ )
168 fprintf( ioPUN,
"\n\n\n" );
169 fprintf( ioPUN,
" end inner shell, start h-like iso seq\n" );
170 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
174 for( nelem=0; nelem <
LIMELM; nelem++ )
191 fprintf( ioPUN,
"\n\n\n" );
192 fprintf( ioPUN,
" end h-like iso seq, start he-like iso seq\n" );
193 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
194 for( nelem=1; nelem <
LIMELM; nelem++ )
209 fprintf( ioPUN,
"\n\n\n" );
210 fprintf( ioPUN,
" end he-like iso seq, start hyperfine structure lines\n" );
211 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
218 fprintf( ioPUN,
"\n\n\n" );
219 fprintf( ioPUN,
" end hyperfine, start database lines\n" );
220 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
229 for(
long nelem = ipISO; nelem <
LIMELM; nelem++ )
252 fprintf( ioPUN,
"\n\n\n" );
253 fprintf( ioPUN,
" end database, start H2 lines\n" );
254 fprintf( ioPUN,
"Eu Vu Ju El Vl Jl WL gl gu gf A CS n(crt)\n" );
261 fprintf( ioPUN,
"\n\n\n" );
262 fprintf( ioPUN,
" end H2, start 12CO rotation lines\n" );
263 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
270 fprintf( ioPUN,
"\n\n\n" );
271 fprintf( ioPUN,
" end 12CO start 13CO rotation lines\n" );
272 fprintf( ioPUN,
"Ion WL gl gu gf A CS n(crt)\n" );
280 fprintf( ioPUN,
"\n\n\n" );
281 fprintf( ioPUN,
" end 13CO rotation lines, start FeII lines\n" );
282 fprintf( ioPUN,
" Lo Hi Ion label WL gl gu gf A CS n(crt)\n" );
318 fprintf(ioPUN,
"%s\t", chLbl );
325 fprintf( ioPUN,
" %3ld%3ld",
340 fprintf( ioPUN,
"%7.1f", t->
Coll.
cs );
342 else if( t->
Coll.
cs > 10. )
344 fprintf( ioPUN,
"%7.2f", t->
Coll.
cs );
346 else if( t->
Coll.
cs > 1. )
348 fprintf( ioPUN,
"%7.3f", t->
Coll.
cs );
350 else if( t->
Coll.
cs > .01 )
352 fprintf( ioPUN,
"%7.4f", t->
Coll.
cs );
354 else if( t->
Coll.
cs > 0.0 )
356 fprintf( ioPUN,
" %.3e", t->
Coll.
cs );
360 fprintf( ioPUN,
"%7.4f", 0. );
366 if( lgCS_2 && t->
Coll.
cs> 0. )
369 CritDen = log10(CritDen);
375 fprintf( ioPUN,
"%7.3f\n",CritDen );