496 for( ion=0; ion<nelem+1; ++ion )
513 (
long**)
MALLOC(
sizeof(
long*)*(unsigned)LIMELM );
515 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
517 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
519 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
521 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
523 (
double****)
MALLOC(
sizeof(
double***)*(unsigned)LIMELM );
525 (
double***)
MALLOC(
sizeof(
double**)*(unsigned)LIMELM );
527 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
529 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
534 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
536 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)LIMELM );
559 for(nelem=0; nelem<
LIMELM; ++nelem )
573 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
575 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
577 (
long*)
MALLOC(
sizeof(
long)*(unsigned)(nelem+1) );
579 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
581 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
583 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
585 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+1) );
587 (
double***)
MALLOC(
sizeof(
double**)*(unsigned)(nelem+1) );
589 (
double**)
MALLOC(
sizeof(
double*)*(unsigned)(nelem+1) );
592 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+2) );
594 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+2) );
597 for( ion=0; ion<nelem+2; ++ion )
603 for( ion=0; ion<nelem+1; ++ion )
620 (
double*)
MALLOC(
sizeof(
double)*(unsigned)2 );
624 (
double*)
MALLOC(
sizeof(
double)*(unsigned)3 );
640 for(ion=0; ion<
LIMELM; ++ion )
646 for( nelem=0; nelem<
LIMELM; ++nelem )
648 for( ion=0; ion<nelem+1; ++ion )
658 for( ns=0; ns <
NSHELLS; ++ns )
668 for( ion=0; ion<nelem+2; ++ion )
674 for( ion2=0; ion2<nelem+2; ++ion2 )
706 for( i=0; i <
LIMELM; i++ )
769 for( ion=0; ion < LIMELM+1; ion++ )
987 int hmi_protons[
N_H_MOLEC] = {1,1,1,2,2,3,2,1};
991 int hmi_electrons[
N_H_MOLEC] = {0,0,-1,0,1,1,0,1};
1138 for( i=0; i <
NCOLD; i++ )
1152 for( i=0; i < 5; i++ )
1160 for( i=0; i < 3; i++ )
1171 for( i=0; i < 4; i++ )
1201 for( i=0; i <
LIMSPC; ++i )
1232 for( i=0; i <
LIMSPC; i++ )
1408 for( nelem=0; nelem<
LIMELM; ++nelem )
1410 for( ion=0; ion<
LIMELM; ++ion )
1421 for( nelem=0; nelem <
LIMELM; nelem++ )
1428 fprintf(
ioQQQ,
" ZERO finds insane abundance or depletion.\n" );
1429 fprintf(
ioQQQ,
" atomic number=%6ld abundance=%10.2e depletion=%10.2e\n",
1543 lgFirstCall =
false;
1561 unsigned long n=(
unsigned long)(ihi-lo+1);
1594 memset(&
opac.
albedo[lo] , 0 , n*
sizeof(
double) );
1602 for( i=lo; i <= ihi; i++ )