Changeset 42948 for trunk/ppSub/src/ppSubDefineOutput.c

Timestamp:

Dec 1, 2025, 3:27:47 PM (5 months ago)

Author:

tdeboer

Message:

adding changes to pswarp to propagate skycell and ZPT keywords

File:

• trunk/ppSub/src/ppSubDefineOutput.c (modified) (2 diffs)

trunk/ppSub/src/ppSubDefineOutput.c

@@ -30 +30 @@
 
     bool status = false;
+    float zptObs_in=0.;
+    float zptErr_in=0.;
+    float zptObs_ref=0.;
+    float zptErr_ref=0.;
 
     pmFPAview *view = ppSubViewReadout(); // View to readout
@@ -139 +143 @@
     psFree(inBase);
 
+    //read the zeropoint keywords
+    pmCell *cell_in = pmFPAfileThisCell(config->files, view, "PPSUB.INPUT");
+    pmHDU *hdu_in = pmHDUFromCell(cell_in);
+    zptObs_in = psMetadataLookupF32 (NULL, hdu_in->header, "ZPT_OBS");
+    zptErr_in = psMetadataLookupF32 (NULL, hdu_in->header, "ZPT_ERR");
+
+    pmCell *cell_ref = pmFPAfileThisCell(config->files, view, "PPSUB.REF");
+    pmHDU *hdu_ref = pmHDUFromCell(cell_ref);
+    zptObs_ref = psMetadataLookupF32 (NULL, hdu_ref->header, "ZPT_OBS");
+    zptErr_ref = psMetadataLookupF32 (NULL, hdu_ref->header, "ZPT_ERR");
+
+    if (zptObs_in != 0.) { 
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_OBS_IN", PS_META_REPLACE, "measured zero point",  zptObs_in);
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_ERR_IN", PS_META_REPLACE, "error on zero point",  zptErr_in);
+    }
+
+    if (zptObs_ref != 0.) { 
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_OBS_REF", PS_META_REPLACE, "measured zero point for template image",  zptObs_ref);
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_ERR_REF", PS_META_REPLACE, "error on zero point for template image",  zptErr_ref);
+    }
+
+    //also output the skycell id
+    char *skycell_id = psMetadataLookupStr(NULL, config->arguments, "SKYCELL_ID");
+    psLogMsg("ppSub", PS_LOG_INFO, "keykeykey o3: %s\n", skycell_id);
+
+    psMetadataAddStr (outHDU->header, PS_LIST_TAIL, "SKYCELL_ID", PS_META_REPLACE, "skycell_id",  skycell_id);
+
     ppSubVersionHeader(outHDU->header);