Changeset 7995
- Timestamp:
- Jul 27, 2006, 2:02:18 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/gastro2/src/gheader2.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/gastro2/src/gheader2.c
r7080 r7995 25 25 26 26 /*** use PutCoords to update header ***/ 27 if (Target[0].coords.Npolyterms > 1) { 28 gfits_modify (&header, "CTYPE1", "%s", 1, "RA---PLY"); 29 gfits_modify (&header, "CTYPE2", "%s", 1, "DEC--PLY"); 30 } else { 31 gfits_modify (&header, "CTYPE1", "%s", 1, "RA---TAN"); 32 gfits_modify (&header, "CTYPE2", "%s", 1, "DEC--TAN"); 33 } 34 gfits_modify (&header, "CDELT1", "%le", 1, Target[0].coords.cdelt1); 35 gfits_modify (&header, "CDELT2", "%le", 1, Target[0].coords.cdelt2); 36 gfits_modify (&header, "CRVAL1", "%lf", 1, Target[0].coords.crval1); 37 gfits_modify (&header, "CRVAL2", "%lf", 1, Target[0].coords.crval2); 38 gfits_modify (&header, "CRPIX1", "%lf", 1, Target[0].coords.crpix1); 39 gfits_modify (&header, "CRPIX2", "%lf", 1, Target[0].coords.crpix2); 40 gfits_modify (&header, "PC001001", "%le", 1, Target[0].coords.pc1_1); 41 gfits_modify (&header, "PC001002", "%le", 1, Target[0].coords.pc1_2); 42 gfits_modify (&header, "PC002001", "%le", 1, Target[0].coords.pc2_1); 43 gfits_modify (&header, "PC002002", "%le", 1, Target[0].coords.pc2_2); 44 gfits_modify (&header, "NPLYTERM", "%d", 1, Target[0].coords.Npolyterms); 45 if (Target[0].coords.Npolyterms > 1) { 46 /* RA Terms */ 47 gfits_modify (&header, "PCA1X2Y0", "%le", 1, Target[0].coords.polyterms[0][0]); /* polyterms[0]); */ 48 gfits_modify (&header, "PCA1X1Y1", "%le", 1, Target[0].coords.polyterms[1][0]); /* polyterms[1]); */ 49 gfits_modify (&header, "PCA1X0Y2", "%le", 1, Target[0].coords.polyterms[2][0]); /* polyterms[2]); */ 50 51 if (Target[0].coords.Npolyterms > 2) { 52 gfits_modify (&header, "PCA1X3Y0", "%le", 1, Target[0].coords.polyterms[3][0]); /* polyterms[3]); */ 53 gfits_modify (&header, "PCA1X2Y1", "%le", 1, Target[0].coords.polyterms[4][0]); /* polyterms[4]); */ 54 gfits_modify (&header, "PCA1X1Y2", "%le", 1, Target[0].coords.polyterms[5][0]); /* polyterms[5]); */ 55 gfits_modify (&header, "PCA1X0Y3", "%le", 1, Target[0].coords.polyterms[6][0]); /* polyterms[6]); */ 56 } 57 /* Dec Terms */ 58 gfits_modify (&header, "PCA2X2Y0", "%le", 1, Target[0].coords.polyterms[0][1]); /* polyterms[7]); */ 59 gfits_modify (&header, "PCA2X1Y1", "%le", 1, Target[0].coords.polyterms[1][1]); /* polyterms[8]); */ 60 gfits_modify (&header, "PCA2X0Y2", "%le", 1, Target[0].coords.polyterms[2][1]); /* polyterms[9]); */ 61 62 if (Target[0].coords.Npolyterms > 2) { 63 gfits_modify (&header, "PCA2X3Y0", "%le", 1, Target[0].coords.polyterms[3][1]); /* polyterms[10]); */ 64 gfits_modify (&header, "PCA2X2Y1", "%le", 1, Target[0].coords.polyterms[4][1]); /* polyterms[11]); */ 65 gfits_modify (&header, "PCA2X1Y2", "%le", 1, Target[0].coords.polyterms[5][1]); /* polyterms[12]); */ 66 gfits_modify (&header, "PCA2X0Y3", "%le", 1, Target[0].coords.polyterms[6][1]); /* polyterms[13]); */ 67 } 68 } 27 PutCoords (&Target[0].coords, &header); 69 28 70 29 dR = fabs (Target[0].answer.dR*Target[0].coords.cdelt1*3600.0);
Note:
See TracChangeset
for help on using the changeset viewer.
