IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
May 5, 2006, 5:44:26 PM (20 years ago)
Author:
eugene
Message:

big move: fits_ to gfits_; minor updates to relphot

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/imregister/imphot/dumpfits.c

    r5014 r7080  
    2727
    2828  /* create primary header */
    29   fits_init_header (&header);   
     29  gfits_init_header (&header);   
    3030  header.extend = TRUE;
    31   fits_create_header (&header);
    32   fits_create_matrix (&header, &matrix);
    33   fits_print (&header, "NEXTEND", "%d", 1, 1);
     31  gfits_create_header (&header);
     32  gfits_create_matrix (&header, &matrix);
     33  gfits_print (&header, "NEXTEND", "%d", 1, 1);
    3434
    3535  ftable.header = &theader;
    36   fits_table_set_Image (&ftable, subset, Nmatch);
     36  gfits_table_set_Image (&ftable, subset, Nmatch);
    3737
    38   fits_fwrite_header  (f, &header);
    39   fits_fwrite_matrix  (f, &matrix);
    40   fits_fwrite_Theader (f, &theader);
    41   fits_fwrite_table   (f, &ftable);
     38  gfits_fwrite_header  (f, &header);
     39  gfits_fwrite_matrix  (f, &matrix);
     40  gfits_fwrite_Theader (f, &theader);
     41  gfits_fwrite_table   (f, &ftable);
    4242  fclose (f);
    4343  exit (0);
     
    5858
    5959  /* create primary header */
    60   fits_init_header (&header);   
     60  gfits_init_header (&header);   
    6161  header.extend = TRUE;
    62   fits_create_header (&header);
    63   fits_create_matrix (&header, &matrix);
    64   fits_print (&header, "NEXTEND", "%d", 1, 1);
     62  gfits_create_header (&header);
     63  gfits_create_matrix (&header, &matrix);
     64  gfits_print (&header, "NEXTEND", "%d", 1, 1);
    6565 
    6666  /* the ASCII table is always a little harder than the binary:
     
    6969
    7070  /* create table header */
    71   fits_create_table_header (&theader, "TABLE", "ZPTS");
     71  gfits_create_table_header (&theader, "TABLE", "ZPTS");
    7272
    7373  /* add current date/time to header */
    7474  str_to_time ("now", &tsecond);
    7575  datestr = sec_to_date (tsecond);
    76   fits_modify (&header,  "DATE", "%s", 1, datestr);
    77   fits_modify (&theader, "DATE", "%s", 1, datestr);
     76  gfits_modify (&header,  "DATE", "%s", 1, datestr);
     77  gfits_modify (&theader, "DATE", "%s", 1, datestr);
    7878 
    7979  /* define table layout */
    80   fits_define_table_column (&theader, "A20",   "START_TIME", "start time of measurement", "yyyy/mm/dd,hh:mm:ss");
    81   fits_define_table_column (&theader, "A10",   "FILTER",     "filter and camera name",    "");
    82   fits_define_table_column (&theader, "F8.4",  "ZP_OBS",     "measured zero point",       "mag");
    83   fits_define_table_column (&theader, "F7.4",  "ZP_ERR",     "error on zero point",       "mag");
    84   fits_define_table_column (&theader, "F11.6", "RA",         "RA (J2000)",                "dec. degrees");
    85   fits_define_table_column (&theader, "F11.6", "DEC",        "DEC (J2000)",               "dec. degrees");
    86   fits_define_table_column (&theader, "F7.3",  "C_AIRMASS",  "airmass coeff",             "mag per airmass");
    87   fits_define_table_column (&theader, "F7.1",  "SKY",        "median sky flux",           "counts");
    88   fits_define_table_column (&theader, "I6",    "NSTAR",      "Number of stars in image",  "stars");
     80  gfits_define_table_column (&theader, "A20",   "START_TIME", "start time of measurement", "yyyy/mm/dd,hh:mm:ss");
     81  gfits_define_table_column (&theader, "A10",   "FILTER",     "filter and camera name",    "");
     82  gfits_define_table_column (&theader, "F8.4",  "ZP_OBS",     "measured zero point",       "mag");
     83  gfits_define_table_column (&theader, "F7.4",  "ZP_ERR",     "error on zero point",       "mag");
     84  gfits_define_table_column (&theader, "F11.6", "RA",         "RA (J2000)",                "dec. degrees");
     85  gfits_define_table_column (&theader, "F11.6", "DEC",        "DEC (J2000)",               "dec. degrees");
     86  gfits_define_table_column (&theader, "F7.3",  "C_AIRMASS",  "airmass coeff",             "mag per airmass");
     87  gfits_define_table_column (&theader, "F7.1",  "SKY",        "median sky flux",           "counts");
     88  gfits_define_table_column (&theader, "I6",    "NSTAR",      "Number of stars in image",  "stars");
    8989
    9090  /* define TNULL, TNVAL values */
    91   fits_modify (&theader, "TNULL1",  "%s", 1, "NULL"); /* START_TIME */
    92   fits_modify (&theader, "TNULL2",  "%s", 1, "NULL"); /* FILTER     */
    93   fits_modify (&theader, "TNULL3",  "%s", 1, "NaN");  /* ZP_OBS     */
    94   fits_modify (&theader, "TNULL4",  "%s", 1, "NaN");  /* ZP_ERR     */
    95   fits_modify (&theader, "TNULL5",  "%s", 1, "NaN");  /* RA         */
    96   fits_modify (&theader, "TNULL6",  "%s", 1, "NaN");  /* DEC        */
    97   fits_modify (&theader, "TNULL7",  "%s", 1, "NaN");  /* C_AIRMASS  */
    98   fits_modify (&theader, "TNULL8",  "%s", 1, "NaN");  /* SKY        */
    99   fits_modify (&theader, "TNULL9",  "%s", 1,  "-1");  /* NSTAR      */
     91  gfits_modify (&theader, "TNULL1",  "%s", 1, "NULL"); /* START_TIME */
     92  gfits_modify (&theader, "TNULL2",  "%s", 1, "NULL"); /* FILTER     */
     93  gfits_modify (&theader, "TNULL3",  "%s", 1, "NaN");  /* ZP_OBS     */
     94  gfits_modify (&theader, "TNULL4",  "%s", 1, "NaN");  /* ZP_ERR     */
     95  gfits_modify (&theader, "TNULL5",  "%s", 1, "NaN");  /* RA         */
     96  gfits_modify (&theader, "TNULL6",  "%s", 1, "NaN");  /* DEC        */
     97  gfits_modify (&theader, "TNULL7",  "%s", 1, "NaN");  /* C_AIRMASS  */
     98  gfits_modify (&theader, "TNULL8",  "%s", 1, "NaN");  /* SKY        */
     99  gfits_modify (&theader, "TNULL9",  "%s", 1,  "-1");  /* NSTAR      */
    100100
    101   fits_modify (&theader, "TNVAL1",  "%s", 1, "NA");   /* START_TIME */
    102   fits_modify (&theader, "TNVAL2",  "%s", 1, "NA");   /* FILTER     */
    103   fits_modify (&theader, "TNVAL3",  "%s", 1, "Inf");  /* ZP_OBS     */
    104   fits_modify (&theader, "TNVAL4",  "%s", 1, "Inf");  /* ZP_ERR     */
    105   fits_modify (&theader, "TNVAL5",  "%s", 1, "Inf");  /* RA         */
    106   fits_modify (&theader, "TNVAL6",  "%s", 1, "Inf");  /* DEC        */
    107   fits_modify (&theader, "TNVAL7",  "%s", 1, "Inf");  /* C_AIRMASS  */
    108   fits_modify (&theader, "TNVAL8",  "%s", 1, "Inf");  /* SKY        */
    109   fits_modify (&theader, "TNVAL9",  "%s", 1,  "-2");  /* NSTAR      */
     101  gfits_modify (&theader, "TNVAL1",  "%s", 1, "NA");   /* START_TIME */
     102  gfits_modify (&theader, "TNVAL2",  "%s", 1, "NA");   /* FILTER     */
     103  gfits_modify (&theader, "TNVAL3",  "%s", 1, "Inf");  /* ZP_OBS     */
     104  gfits_modify (&theader, "TNVAL4",  "%s", 1, "Inf");  /* ZP_ERR     */
     105  gfits_modify (&theader, "TNVAL5",  "%s", 1, "Inf");  /* RA         */
     106  gfits_modify (&theader, "TNVAL6",  "%s", 1, "Inf");  /* DEC        */
     107  gfits_modify (&theader, "TNVAL7",  "%s", 1, "Inf");  /* C_AIRMASS  */
     108  gfits_modify (&theader, "TNVAL8",  "%s", 1, "Inf");  /* SKY        */
     109  gfits_modify (&theader, "TNVAL9",  "%s", 1,  "-2");  /* NSTAR      */
    110110
    111111  /* add data to table */
     
    121121
    122122    /* we should get an error here if we don't construct this line correctly */
    123     line = fits_table_print (&ftable, startstr, filtstr, zp, dzp, ra, dec, airmass, sky, subset[0].nstar);
    124     fits_add_rows (&ftable, line, 1, strlen(line));
     123    line = gfits_table_print (&ftable, startstr, filtstr, zp, dzp, ra, dec, airmass, sky, subset[0].nstar);
     124    gfits_add_rows (&ftable, line, 1, strlen(line));
    125125    free (line);
    126126    free (startstr);
     
    133133    exit (1);
    134134  }
    135   fits_fwrite_header  (f, &header);
    136   fits_fwrite_matrix  (f, &matrix);
    137   fits_fwrite_Theader (f, &theader);
    138   fits_fwrite_table   (f, &ftable);
     135  gfits_fwrite_header  (f, &header);
     136  gfits_fwrite_matrix  (f, &matrix);
     137  gfits_fwrite_Theader (f, &theader);
     138  gfits_fwrite_table   (f, &ftable);
    139139  fclose (f);
    140140  return (TRUE);
Note: See TracChangeset for help on using the changeset viewer.