34 ASSERT( chType ==
'i' || chType ==
't' );
51 fprintf( ioMEAN,
"\n Hydrogen " );
52 for( i=0; i < 3; i++ )
54 fprintf( ioMEAN,
"%7.3f", aa[i] );
56 fprintf(ioMEAN,
" (H2)");
57 if( chType==
'i' && lgDensity )
60 " Log10 Mean Ionisation (over volume*electron density)\n" );
62 else if( chType==
'i' )
65 " Log10 Mean Ionisation (over volume)\n" );
67 else if( chType==
't' && lgDensity )
70 " Log10 Mean Temperature (over volume*electron density)\n" );
72 else if( chType==
't' )
75 " Log10 Mean Temperature (over volume)\n" );
79 fprintf(
ioQQQ,
"PrtMeanIon called with insane job\n");
91 for( i=0; i < limit; i++ )
93 fprintf( ioMEAN,
"%7.3f", aa[i] );
95 fprintf( ioMEAN,
"\n" );
100 fprintf( ioMEAN,
" " );
101 for( i=17; i < n; i++ )
103 fprintf( ioMEAN,
"%7.3f", aa[i] );
105 fprintf( ioMEAN,
"\n" );
111 fprintf( ioMEAN,
"\n " );
112 for( i=1; i <= 17; i++ )
114 fprintf( ioMEAN,
"%7ld", i );
116 fprintf( ioMEAN,
"\n" );
119 fprintf( ioMEAN,
" " );
120 for( i=18; i <=
LIMELM; i++ )
122 fprintf( ioMEAN,
"%7ld", i );
124 fprintf( ioMEAN,
"\n" );
133 fprintf( ioMEAN,
"\n Hydrogen " );
135 for( i=0; i < 3; i++ )
138 fprintf( ioMEAN,
"%7.3f", aa[i] );
140 fprintf(ioMEAN,
" (H2)");
141 if( chType==
'i' && lgDensity )
144 " Log10 Mean Ionisation (over radius*electron density)\n" );
146 else if( chType==
'i' )
149 " Log10 Mean Ionisation (over radius)\n" );
151 else if( chType==
't' && lgDensity )
154 " Log10 Mean Temperature (over radius*electron density)\n" );
156 else if( chType==
't' )
159 " Log10 Mean Temperature (over radius)\n" );
163 fprintf(
ioQQQ,
"PrtMeanIon called with insane job\n");
173 fprintf( ioMEAN,
" %10.10s",
175 for( i=0; i < limit; i++ )
177 fprintf( ioMEAN,
"%7.3f", aa[i] );
179 fprintf( ioMEAN,
"\n" );
184 fprintf( ioMEAN,
" " );
186 for( i=17; i < n; i++ )
188 fprintf( ioMEAN,
"%7.3f", aa[i] );
190 fprintf( ioMEAN,
"\n" );