- Timestamp:
- May 3, 2010, 8:50:52 AM (16 years ago)
- Location:
- branches/simtest_nebulous_branches
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/simtest_nebulous_branches
- Property svn:mergeinfo changed
-
branches/simtest_nebulous_branches/Ohana
-
Property svn:mergeinfo
set to (toggle deleted branches)
/trunk/Ohana merged eligible /branches/eam_branches/Ohana.20100407 27635-27772 /branches/pap_delete/Ohana 27530-27595
-
Property svn:mergeinfo
set to (toggle deleted branches)
-
branches/simtest_nebulous_branches/Ohana/src/gastro2/src/gheader2.c
r8300 r27840 6 6 Header header; 7 7 FILE *f, *g; 8 int i, oldsize, nbytes, status;8 off_t i, oldsize, nbytes, status; 9 9 char line[1024]; 10 10 … … 13 13 exit(0); 14 14 } 15 oldsize = header. size;15 oldsize = header.datasize; 16 16 17 17 /* check for insufficient number of stars */ … … 21 21 if (Target[0].answer.N < 6) { 22 22 gfits_modify (&header, "NASTRO", "%d", 1, 0); 23 gfits_modify (&header, "NASTRO", "%C", 1, "number of stars used for astrometry");23 gfits_modify_alt (&header, "NASTRO", "%C", 1, "number of stars used for astrometry"); 24 24 goto skipstuff; 25 25 } … … 28 28 if (Target[0].answer.N < 12) { 29 29 gfits_modify (&header, "NASTRO", "%d", 1, 0); 30 gfits_modify (&header, "NASTRO", "%C", 1, "number of stars used for astrometry");30 gfits_modify_alt (&header, "NASTRO", "%C", 1, "number of stars used for astrometry"); 31 31 goto skipstuff; 32 32 } … … 35 35 if (Target[0].answer.N < 20) { 36 36 gfits_modify (&header, "NASTRO", "%d", 1, 0); 37 gfits_modify (&header, "NASTRO", "%C", 1, "number of stars used for astrometry");37 gfits_modify_alt (&header, "NASTRO", "%C", 1, "number of stars used for astrometry"); 38 38 goto skipstuff; 39 39 } … … 45 45 46 46 gfits_modify (&header, "NASTRO", "%d", 1, Target[0].answer.N); 47 gfits_modify (&header, "NASTRO", "%C", 1, "number of stars used for astrometry");47 gfits_modify_alt (&header, "NASTRO", "%C", 1, "number of stars used for astrometry"); 48 48 49 49 /*** use PutCoords to update header ***/ … … 52 52 dR = fabs (Target[0].answer.dR*Target[0].coords.cdelt1*3600.0); 53 53 gfits_modify (&header, "CERROR", "%lf", 1, dR); 54 gfits_modify (&header, "CERROR", "%C", 1, "scatter in astrometry soln (arcsec)");54 gfits_modify_alt (&header, "CERROR", "%C", 1, "scatter in astrometry soln (arcsec)"); 55 55 gfits_modify (&header, "CPRECISE", "%lf", 1, dR / sqrt(1.0*Target[0].answer.N)); 56 gfits_modify (&header, "CPRECISE", "%C", 1, "precision of astrometry soln (arcsec)");56 gfits_modify_alt (&header, "CPRECISE", "%C", 1, "precision of astrometry soln (arcsec)"); 57 57 gfits_modify (&header, "EQUINOX", "%lf", 1, 2000.0); 58 58 /* we force equinox to be 2000.0 for all images */ 59 59 60 60 skipstuff: 61 if (header. size > oldsize) {61 if (header.datasize > oldsize) { 62 62 if (VERBOSE) fprintf (stderr, "header expanded, creating new copy\n"); 63 63 sprintf (line, "mv %s %s~", file, file); … … 78 78 exit(0); 79 79 } 80 nbytes = fwrite (header.buffer, 1, header. size, g);81 fseek (f, oldsize, SEEK_SET);82 for (i = 0; (nbytes = fread (header.buffer, 1, header. size, f)) > 0; i++) {80 nbytes = fwrite (header.buffer, 1, header.datasize, g); 81 fseeko (f, oldsize, SEEK_SET); 82 for (i = 0; (nbytes = fread (header.buffer, 1, header.datasize, f)) > 0; i++) { 83 83 if (nbytes != fwrite (header.buffer, 1, nbytes, g)) { 84 84 fprintf (stderr, "ERROR: failure writing output data file\n"); … … 95 95 } 96 96 97 fseek (f, 0, SEEK_SET);98 nbytes = fwrite (header.buffer, 1, header. size, f);97 fseeko (f, 0, SEEK_SET); 98 nbytes = fwrite (header.buffer, 1, header.datasize, f); 99 99 100 100 fclose (f);
Note:
See TracChangeset
for help on using the changeset viewer.
