@@ -318,6 +318,14 @@ int main(int argc, char *argv[]) {
318
318
}
319
319
token = strtok (line , " \f\n\r\t\v" );
320
320
ncomponents = atoi (token );
321
+ sprintf (historyline ," %d hyperfine components" ,ncomponents );
322
+ checkfits (fits_write_history (fitsout ,historyline ,& status ));
323
+ checkfits (fits_write_history (tauout ,historyline ,& status ));
324
+ checkfits (fits_write_history (vlsrout ,historyline ,& status ));
325
+ checkfits (fits_write_history (sigmaout ,historyline ,& status ));
326
+ checkfits (fits_write_history (texout ,historyline ,& status ));
327
+ checkfits (fits_write_history (chisqout ,historyline ,& status ));
328
+
321
329
component_voffs = malloc (ncomponents * sizeof (double ));
322
330
component_relints = malloc (ncomponents * sizeof (double ));
323
331
if (component_voffs == NULL || component_relints == NULL ) {
@@ -333,6 +341,13 @@ int main(int argc, char *argv[]) {
333
341
component_voffs [i ] = atof (token );
334
342
token = strtok (NULL , " \f\n\r\t\v" );
335
343
component_relints [i ] = atof (token );
344
+ sprintf (historyline ," Component %d: Voff= %f, Rel. Intensity= %f" , i , component_voffs [i ], component_relints [i ]);
345
+ checkfits (fits_write_history (fitsout ,historyline ,& status ));
346
+ checkfits (fits_write_history (tauout ,historyline ,& status ));
347
+ checkfits (fits_write_history (vlsrout ,historyline ,& status ));
348
+ checkfits (fits_write_history (sigmaout ,historyline ,& status ));
349
+ checkfits (fits_write_history (texout ,historyline ,& status ));
350
+ checkfits (fits_write_history (chisqout ,historyline ,& status ));
336
351
}
337
352
fclose (fpin );
338
353
0 commit comments