122 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
135 for(nelem=0; nelem<
LIMELM; ++nelem )
139 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+2) );
141 (
double*)
MALLOC(
sizeof(
double)*(unsigned)(nelem+2) );
144 for( ion=0; ion<nelem+2; ++ion )
155 fprintf(
ioQQQ,
" IterStart called.\n" );
239 for( nelem=0; nelem <
LIMELM; nelem++ )
246 for( ion=0; ion < (nelem + 2); ion++ )
251 for( ion=0; ion <
LIMELM; ion++ )
262 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
316 for( ion=0; ion<nelem+2; ++ion )
325 for( ion2=0; ion2<nelem+2; ++ion2 )
426 for( i=0; i <
NCOLD; i++ )
446 for( i=0; i < 5; i++ )
454 for( i=0; i < 3; i++ )
465 for( i=0; i < 4; i++ )
499 for( ion=0; ion<nelem+1; ++ion )
509 for( ion=0; ion<nelem+1; ++ion )
529 aver(
"zero",0.,0.,
" ");
584 fprintf(
ioQQQ,
"IterStart has freed LindSv \n" );
614 fprintf(
ioQQQ,
"%7ld lines printed in main line array\n",
651 fprintf(
ioQQQ,
"PROBLEM could not find the normalisation line.\n");
653 "IterStart could not find a line with a wavelength of %f and a label of %s\n",
655 fprintf(
ioQQQ,
"Please check the emission line output to find the correct line identification.\n");
656 fprintf(
ioQQQ,
"Sorry.\n");
669 fprintf(
ioQQQ,
"PROBLEM could not find the default normalisation line.\n");
671 "IterStart could not find a line with a wavelength of 4861 and a label of TOTL\n" );
672 fprintf(
ioQQQ,
"Please check the emission line output to find the correct line identification.\n");
673 fprintf(
ioQQQ,
"Sorry.\n");
684 for(
long int nStopLine=0; nStopLine <
StopCalc.
nstpl; nStopLine++ )
686 double relint, absint ;
696 " IterStart could not find first line in STOP LINE command, number %ld with label *%s* and wl %.1f\n",
710 " IterStart could not find second line in STOP LINE command, line number %ld with label *%s* and wl %.1f\n",
720 " stop line 1 is number %5ld wavelength is %f label is %4.4s\n",
725 " stop line 2 is number %5ld wavelength is %f label is %4.4s\n",
755 fprintf(
ioQQQ,
" Outer optical depths defined.\n" );
759 fprintf(
ioQQQ,
" Outer optical depths NOT defined.\n" );
765 fprintf(
ioQQQ,
"%6f:%8.1e%8.1e%8.1e",
775 fprintf(
ioQQQ,
"%6f:%8.1e%8.1e%8.1e",
782 fprintf(
ioQQQ,
"\n" );
789 fprintf(
ioQQQ,
" IterStart does not change mid-zone optical depth since CASE B\n" );
829 fprintf(
ioQQQ,
" IterStart returns.\n" );
878 for( ion=0; ion<nelem+2; ++ion )
887 for( ion2=0; ion2<nelem+2; ++ion2 )
905 for( i=0; i <
NCOLD; i++ )
987 fprintf(
ioQQQ,
" IterRestart called.\n" );
997 for( ion=0; ion<nelem+1; ++ion )
1005 for( nelem=0; nelem<
LIMELM; ++nelem)
1007 for( ion=0; ion<nelem+1; ++ion )
1051 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
1066 fprintf(
ioQQQ,
" EDEN set to%12.4e by IterRestart.\n",
1087 for( ion=0; ion < (nelem + 2); ion++ )
1091 for( i=0; i <
LIMELM; i++ )
1186 fprintf(
ioQQQ,
"\f\n Start Iteration Number%2ld %75.75s\n\n\n",
1229 enum{DEBUG_LOC=
false};
1232 fprintf(
ioQQQ,
"i=%li opac %.2e \n", i,