Changeset 7054 for trunk/Ohana/src/libfits/doc/sample.c
- Timestamp:
- May 3, 2006, 12:01:09 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/libfits/doc/sample.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/libfits/doc/sample.c
r7039 r7054 25 25 26 26 /* create primary header */ 27 fits_init_header (&header); header.extend = TRUE;28 fits_create_header (&header);29 fits_create_matrix (&header, &matrix);30 fits_print (&header, "NEXTEND", "%d", 1, 2);27 gfits_init_header (&header); header.extend = TRUE; 28 gfits_create_header (&header); 29 gfits_create_matrix (&header, &matrix); 30 gfits_print (&header, "NEXTEND", "%d", 1, 2); 31 31 32 32 /* define bintable header & layout */ 33 fits_create_table_header (&theader1, "BINTABLE", "ZERO_POINTS");33 gfits_create_table_header (&theader1, "BINTABLE", "ZERO_POINTS"); 34 34 35 fits_define_bintable_column (&theader1, "E", "ZP_OBS", "measured zero point", "mag", 1.0, 0.0);36 fits_define_bintable_column (&theader1, "E", "ZP_REF", "measured zero point", "mag", 1.0, 0.0);37 fits_define_bintable_column (&theader1, "J", "TIME", "time of data", "seconds since Jan 1, 1970 UT", 1.0, 0.0);38 fits_define_bintable_column (&theader1, "8A", "TCHAR", "time of data", "seconds since Jan 1, 1970 UT", 1.0, 0.0);35 gfits_define_bintable_column (&theader1, "E", "ZP_OBS", "measured zero point", "mag", 1.0, 0.0); 36 gfits_define_bintable_column (&theader1, "E", "ZP_REF", "measured zero point", "mag", 1.0, 0.0); 37 gfits_define_bintable_column (&theader1, "J", "TIME", "time of data", "seconds since Jan 1, 1970 UT", 1.0, 0.0); 38 gfits_define_bintable_column (&theader1, "8A", "TCHAR", "time of data", "seconds since Jan 1, 1970 UT", 1.0, 0.0); 39 39 40 40 /* create table, add data values */ 41 fits_create_table (&theader1, &table1);41 gfits_create_table (&theader1, &table1); 42 42 43 43 /* set table column based on array, extend NAXIS2 as needed/appropriate */ 44 fits_set_bintable_column (&theader1, &table1, "ZP_OBS", zpobs, Nrow);45 fits_set_bintable_column (&theader1, &table1, "ZP_REF", zpref, Nrow);46 fits_set_bintable_column (&theader1, &table1, "TIME", time, Nrow);47 fits_set_bintable_column (&theader1, &table1, "TCHAR", tchar, Nrow);44 gfits_set_bintable_column (&theader1, &table1, "ZP_OBS", zpobs, Nrow); 45 gfits_set_bintable_column (&theader1, &table1, "ZP_REF", zpref, Nrow); 46 gfits_set_bintable_column (&theader1, &table1, "TIME", time, Nrow); 47 gfits_set_bintable_column (&theader1, &table1, "TCHAR", tchar, Nrow); 48 48 49 49 /* define ASCII table header */ 50 fits_create_table_header (&theader2, "TABLE", "ASCII_PTS");50 gfits_create_table_header (&theader2, "TABLE", "ASCII_PTS"); 51 51 52 fits_define_table_column (&theader2, "F5.2", "ZP_OBS", "measured zero point", "mag");53 fits_define_table_column (&theader2, "F5.2", "ZP_REF", "measured zero point", "mag");54 fits_define_table_column (&theader2, "I8", "TIME", "time of data", "seconds");55 fits_define_table_column (&theader2, "A8", "TCHAR", "time of data", "YYYYMMDD");52 gfits_define_table_column (&theader2, "F5.2", "ZP_OBS", "measured zero point", "mag"); 53 gfits_define_table_column (&theader2, "F5.2", "ZP_REF", "measured zero point", "mag"); 54 gfits_define_table_column (&theader2, "I8", "TIME", "time of data", "seconds"); 55 gfits_define_table_column (&theader2, "A8", "TCHAR", "time of data", "YYYYMMDD"); 56 56 57 57 /* create table, add data values */ 58 fits_create_table (&theader2, &table2);58 gfits_create_table (&theader2, &table2); 59 59 60 fits_set_table_column (&theader2, &table2, "ZP_OBS", zpobs, Nrow);61 fits_set_table_column (&theader2, &table2, "ZP_REF", zpref, Nrow);62 fits_set_table_column (&theader2, &table2, "TIME", time, Nrow);63 fits_set_table_column (&theader2, &table2, "TCHAR", tchar, Nrow);60 gfits_set_table_column (&theader2, &table2, "ZP_OBS", zpobs, Nrow); 61 gfits_set_table_column (&theader2, &table2, "ZP_REF", zpref, Nrow); 62 gfits_set_table_column (&theader2, &table2, "TIME", time, Nrow); 63 gfits_set_table_column (&theader2, &table2, "TCHAR", tchar, Nrow); 64 64 65 65 /* write header, matrix, table1, table2 */ 66 fits_write_header (filename, &header);67 fits_write_matrix (filename, &matrix);68 fits_write_Theader (filename, &theader1);69 fits_write_table (filename, &table1);70 fits_write_Theader (filename, &theader2);71 fits_write_table (filename, &table2);66 gfits_write_header (filename, &header); 67 gfits_write_matrix (filename, &matrix); 68 gfits_write_Theader (filename, &theader1); 69 gfits_write_table (filename, &table1); 70 gfits_write_Theader (filename, &theader2); 71 gfits_write_table (filename, &table2); 72 72 73 73 }
Note:
See TracChangeset
for help on using the changeset viewer.
