33 for( nelem=0; nelem <
LIMELM; nelem++ )
57 for( ion=0; ion < (nelem + 2); ion++ )
125 for( ion=0; ion < (nelem + 2); ion++ )
218 static bool lgFirst=
true;
249 for( nelem=0; nelem<
LIMELM; ++nelem )
251 limit =
MAX2(3,nelem+2);
262 for( j=0; j < 2; j++ )
264 for( nelem=0; nelem <
LIMELM; nelem++ )
269 limit =
MAX2(3,nelem+2);
271 for( ion=0; ion < limit; ion++ )
330 double abund_radmean,
335 ASSERT( chType==
'i' || chType==
't' );
340 limit =
MAX2( 3, nelem+2 );
355 for( ion=0; ion < limit; ion++ )
372 if( chType==
'i' && lgDensity)
375 for( ion=0; ion < *n; ion++ )
388 else if( chType==
'i' )
391 for( ion=0; ion < *n; ion++ )
404 else if( chType==
't' && lgDensity )
407 for( ion=0; ion < *n; ion++ )
413 arlog[ion] = (
realnum)log10(
MAX2(1e-30,abund_radmean/normalize));
421 else if( chType==
't' )
424 for( ion=0; ion < *n; ion++ )
430 arlog[ion] = (
realnum)log10(
MAX2(1e-30,abund_radmean/normalize));
440 fprintf(
ioQQQ,
" MeanIonRadius called with insane job \n");
459 double abund_volmean,
465 ASSERT( chType==
'i' || chType==
't' );
469 limit =
MAX2( 3, nelem+2 );
483 for( ion=0; ion <= limit; ion++ )
502 if( chType==
'i' && lgDensity)
507 for( ion=0; ion < *n; ion++ )
520 else if( chType==
'i' )
525 for( ion=0; ion < *n; ion++ )
538 else if( chType==
't' && lgDensity )
543 for( ion=0; ion < *n; ion++ )
549 arlog[ion] = (
realnum)log10(
MAX2(1e-30,abund_volmean)/normalize);
557 else if( chType==
't' )
562 for( ion=0; ion < *n; ion++ )
568 arlog[ion] = (
realnum)log10(
MAX2(1e-30,abund_volmean)/normalize);
578 fprintf(
ioQQQ,
" MeanIonVolume called with insane job\n");
597 static char chLabavr[
NAVER][11];
602 double raver[
NAVER]={0.},
604 static double aversv[4][
NAVER];
608 if( strcmp(chWhat,
"zero") == 0 )
611 for( i=0; i <
NAVER; i++ )
613 for( j=0; j < 4; j++ )
620 else if( strcmp(chWhat,
"zone") == 0 )
624 else if( strcmp(chWhat,
"doit") == 0 )
629 fprintf(
ioQQQ,
" Too many values entered into AVER, increase NAVER\n" );
638 strcpy( chLabavr[n], chLabl );
641 else if( strcmp(chWhat,
"prin") == 0 )
647 fprintf(
ioQQQ,
"\n");
648 for( i=0; i < ioff; i++ )
651 fprintf(
ioQQQ,
" " );
655 fprintf(
ioQQQ,
"Averaged Quantities\n " );
657 fprintf(
ioQQQ,
" " );
658 for( i=0; i < n; i++ )
660 fprintf(
ioQQQ,
"%10.10s", chLabavr[i] );
662 fprintf(
ioQQQ,
" \n" );
664 for( i=0; i < n; i++ )
666 if( aversv[1][i] > 0. )
668 raver[i] = aversv[0][i]/aversv[1][i];
674 if( aversv[3][i] > 0. )
676 vaver[i] = aversv[2][i]/aversv[3][i];
684 fprintf(
ioQQQ,
" Radius:" );
685 for( i=0; i < n; i++ )
688 fprintf(
ioQQQ,
" ");
691 fprintf(
ioQQQ,
"\n" );
696 fprintf(
ioQQQ,
" Volume:" );
697 for( i=0; i < n; i++ )
700 fprintf(
ioQQQ,
" ");
703 fprintf(
ioQQQ,
"\n" );
708 fprintf(
ioQQQ,
" AVER does not understand chWhat=%4.4s\n",