29 size_t sizeof_sanity =
sizeof( sanity );
40 if( (n=fread( pnt , 1 , sizeof_pnt ,
ioSTATE )) - sizeof_pnt )
42 fprintf(
ioQQQ,
" state_do failed reading state file, wanted %lu got %lu\n",
43 (
unsigned long)sizeof_pnt ,
48 if( (n=fread( &chk_sanity , 1 , sizeof_sanity ,
ioSTATE )) - sizeof_sanity )
50 fprintf(
ioQQQ,
" state_do failed reading sanity par of state file, wanted %lu got %lu\n",
51 (
unsigned long)sizeof_sanity ,
55 if( !
fp_equal( sanity, chk_sanity ) )
57 fprintf(
ioQQQ,
" state_do sanity fails in state file, wanted %g got %g\n",
66 fwrite( pnt , 1 , sizeof_pnt ,
ioSTATE );
68 fwrite( &sanity , 1 , sizeof_sanity ,
ioSTATE );
77 long int ipISO , nelem , ipHi, i ,
82 if( (strcmp( chJob ,
"get" ) == 0) )
87 else if( (strcmp( chJob ,
"put" ) == 0) )
103 sprintf( chIteration ,
"_%li",
iteration );
104 strcat( chFilename , chIteration );
112 fprintf(
ioQQQ,
" Print state quantities, start iso seq \n");
121 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
130 fprintf(
ioQQQ,
" start ISO ipISO= %li, nelem= %li, ipHi %li \n",
134 for( n=0; n< ipHi; ++n )
136 fprintf(
ioQQQ,
" ISO %li %li %li %li %.4e %.4e \n",
137 ipISO , nelem , ipHi , n ,
141 fprintf(
ioQQQ,
" end ISO ipISO\n");
147 fprintf(
ioQQQ,
" start Ext ipISO= %li, nelem= %li, got %li \n",
156 fprintf(
ioQQQ,
" Ext %li %li %li %.4e %.4e \n",
161 fprintf(
ioQQQ,
" end Ext ipISO\n");
170 for( n=0; n< (
nLevel1+1); ++n )
172 fprintf(
ioQQQ,
" Taulines %li %.4e %.4e \n",
184 fprintf(
ioQQQ,
" TauLine2 %li %.4e %.4e \n",
194 for( n=0; n<
nUTA; ++n )
196 fprintf(
ioQQQ,
" UTALines %li %.4e %.4e \n",
208 fprintf(
ioQQQ,
" HFLines %li %.4e %.4e \n",
220 fprintf(
ioQQQ,
" C12O16Rotate %li %.4e %.4e \n",
232 fprintf(
ioQQQ,
" C13O16Rotate %li %.4e %.4e \n",
244 for( n=0; n< ipHi; ++n )
246 fprintf(
ioQQQ,
" Fe2LevN %li %li %.4e %.4e \n",
260 fprintf(
ioQQQ,
" TauAbsGeo %li %li %.4e \n",
271 fprintf(
ioQQQ,
" TauScatGeo %li %li %.4e \n",
282 fprintf(
ioQQQ,
" TauTotalGeo %li %li %.4e \n",
321 for( ion=0; ion<(
LIMELM+1); ++ion )