IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 7995


Ignore:
Timestamp:
Jul 27, 2006, 2:02:18 PM (20 years ago)
Author:
eugene
Message:

replaced header updates with PutCoords

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/gastro2/src/gheader2.c

    r7080 r7995  
    2525
    2626  /*** 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);
    6928
    7029  dR = fabs (Target[0].answer.dR*Target[0].coords.cdelt1*3600.0);
Note: See TracChangeset for help on using the changeset viewer.