Initial commit

1 files changed, 114 insertions, 0 deletions

diff --git a/src/analysis/modhead.c b/src/analysis/modhead.c
new file mode 100644
index 0000000..0ffc40b
--- /dev/null
+++ b/src/analysis/modhead.c
@@ -0,0 +1,114 @@
+/*****************************************************************************
+ *	      Johns Hopkins University
+ *	      Center For Astrophysical Sciences
+ *	      FUSE
+ *****************************************************************************
+ *
+ * Synopsis:	modhead file.fit[ext] KEYWORD NEW_VALUE
+ *
+ * Description: Returns or modifies the value of a header keyword.
+ *		If NEW_VALUE is not specified, current value of KEYWORD is
+ *		returned.  If KEYWORD is not found and NEW_VALUE is given,
+ *		KEYWORD is added to header.  If an extension beyond the
+ *		primary HDU is requested, file name must be given in
+ *		parentheses: "file.fit[3]" or an alternate construction
+ *		used: file.fit+3   Note that the primary HDU is extension 0.
+ *
+ * History:     05/15/05        tc	started work
+ *		08/30/05	wvd	Allow multi-word strings.
+ *
+ ****************************************************************************/
+
+#include <string.h>
+#include <stdio.h>
+#include "fitsio.h"
+
+int main(int argc, char *argv[])
+{
+    fitsfile *fptr;         /* FITS file pointer, defined in fitsio.h */
+    char card[FLEN_CARD], newcard[FLEN_CARD];
+    char oldvalue[FLEN_VALUE], comment[FLEN_COMMENT];
+    int status = 0;   /*  CFITSIO status value MUST be initialized to zero!  */
+    int i, iomode, keytype;
+
+    if (argc == 3) 
+      iomode = READONLY;
+    else if (argc >= 4)
+      iomode = READWRITE;
+    else if (argc < 3) {
+      printf("Usage:  modhead filename[ext] keyword new_value\n");
+      printf("\n");
+      printf("Write or modify the value of a header keyword.\n");
+      printf("If 'newvalue' is not specified, then just print \n");
+      printf("the current value. \n");
+      printf("\n");
+      printf("Examples: \n");
+      printf("  modhead file.fits dec      - list the DEC keyword \n");
+      printf("  modhead file.fits dec 30.0 - set DEC = 30.0 \n");
+      printf("\n");
+      printf("  modhead \"file.fits[2]\" OBSTIME  - list OBSTIME from extention 2 \n");
+      printf("  modhead file.fits+2 OBSTIME     - list OBSTIME from extention 2\n");
+      return(0);
+    }
+
+    if (!fits_open_file(&fptr, argv[1], iomode, &status))
+    {
+      if (fits_read_card(fptr,argv[2], card, &status))
+      {
+        printf("Keyword does not exist\n");
+        card[0] = '\0';
+        comment[0] = '\0';
+        status = 0;  /* reset status after error */
+      }
+      else
+        printf("%s\n",card);
+
+      if (argc >= 4)  /* write or overwrite the keyword */
+      {
+          /* check if this is a protected keyword that must not be changed */
+          if (*card && fits_get_keyclass(card) == TYP_STRUC_KEY)
+          {
+            printf("Protected keyword cannot be modified.\n");
+          }
+          else
+          {
+            /* get the comment string */
+            if (*card)fits_parse_value(card, oldvalue, comment, &status);
+
+            /* construct template for new keyword */
+            strcpy(newcard, argv[2]);     /* copy keyword name */
+            strcat(newcard, " = ");       /* '=' value delimiter */
+	    if (argc == 4) strcat(newcard, argv[3]);     /* new value */
+	    else {
+                strcat(newcard, "'");
+	        for (i = 3; i < argc-1; i++) {
+                    strcat(newcard, argv[i]);     /* new value */
+                    strcat(newcard, " ");
+		}
+                strcat(newcard, argv[i]);     /* new value */
+                strcat(newcard, "'");
+	    }
+
+            if (*comment) {
+              strcat(newcard, " / ");  /* comment delimiter */
+              strcat(newcard, comment);     /* append the comment */
+            }
+
+            /* reformat the keyword string to conform to FITS rules */
+            fits_parse_template(newcard, card, &keytype, &status);
+
+            /* overwrite the keyword with the new value */
+            fits_update_card(fptr, argv[2], card, &status);
+
+            printf("Keyword has been changed to:\n");
+            printf("%s\n",card);
+          }
+      }  /* if argc == 4 */
+      fits_close_file(fptr, &status);
+    }    /* open_file */
+
+    /* if error occured, print out error message */
+    if (status) fits_report_error(stderr, status);
+    return(status);
+}
+