64 fprintf( io,
" Cloudy punts, Te=%10.3e HTOT=%10.3e CTOT=%10.3e nzone=%4ld\n",
66 fprintf( io,
" COOLNG array is\n" );
70 coolpr(io,
"ZERO",1,0.,
"ZERO");
81 coolpr(io,
"DONE",1,0.,
"DONE");
82 fprintf( io,
" Line heating array follows\n" );
83 coolpr(io,
"ZERO",1,0.,
"ZERO");
95 coolpr(io,
"DONE",1,0.,
"DONE");
102 fprintf( io,
" map of heating, cooling, vs temp, follows.\n");
104 " Te Heat--------------------> Cool---------------------> dH/dT dC/DT Ne NH HII Helium \n" );
107 if( strcmp(chType,
"punt") == 0 )
122 else if( strcmp(chType,
" map") == 0 )
143 TeNew = (tlowst/factor);
151 fprintf(
ioQQQ,
" PUNT called with insane argument,=%4.4s\n",
166 printf(
" not enough memory to allocate hcmap.temap in map_do\n" );
172 printf(
" not enough memory to allocate hcmap.cmap in map_do\n" );
178 printf(
" not enough memory to allocate hcmap.hmap in map_do\n" );
188 fprintf(
ioQQQ,
" MAP called temp range %.4e %.4e in %li stops ===============================================\n",
225 fprintf(
ioQQQ,
" MAP new temp %.4e ===============================================\n",
263 for( k=0; k <
LIMELM; k++ )
265 for( j=0; j <
LIMELM; j++ )
281 fprintf( io,
" ----------------------------------------------- Probable thermal solution here. --------------------------------------------\n" );
302 fprintf( io,
"TRACE: te, htot, ctot%11.3e%11.3e%11.3e\n",
310 fprintf(io,
" [%2ld][%2ld]%6.3f",
314 fprintf(io,
" %s %.1f%c%6.3f",
325 fprintf(io,
"%6.2f%6.2f%6.2f%6.2f\n",
339 enum {DEBUG_LOC=
false};
342 static int kount = 0;
348 fprintf(
ioQQQ,
" exiting in map_do\n");
365 if( s1*s3 > 0. && s2*s3 < 0. )
374 " cooling curve had double inflection at T=%.2e. ",
376 fprintf( io,
" Slopes were %.2e %.2e %.2e", s1, s2, s3);
380 " error large, (rel slope of %.2e).\n",
387 " error is small, (rel slope of %.2e).\n",
395 if( s1*s3 > 0. && s2*s3 < 0. )
403 " heating curve had double inflection at T=%.2e.\n",