Added minimum checks to see if maskfile matches the input raw data

scottransom · May 20, 2024 · a5286ef · a5286ef
1 parent 9a7c972
commit a5286ef
Show file tree

Hide file tree

Showing 4 changed files with 25 additions and 0 deletions.
diff --git a/src/mpiprepsubband.c b/src/mpiprepsubband.c
@@ -354,6 +354,11 @@ int main(int argc, char *argv[])
         if (myid == 0) {
             read_mask(cmd->maskfile, &obsmask);
             printf("Read mask information from '%s'\n\n", cmd->maskfile);
+            if ((obsmask.numchan != idata.num_chan) ||
+                (fabs(obsmask.mjd - (idata.mjd_i + idata.mjd_f)) > 1e-9)) {
+                    printf("WARNING!: maskfile has different number of channels or start MJD than raw data! Exiting.\n\n");
+                    exit(1);
+            }
             good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
         }
         broadcast_mask(&obsmask, myid);

diff --git a/src/prepdata.c b/src/prepdata.c
@@ -228,6 +228,11 @@ int main(int argc, char *argv[])
     if (cmd->maskfileP) {
         read_mask(cmd->maskfile, &obsmask);
         printf("Read mask information from '%s'\n\n", cmd->maskfile);
+        if ((obsmask.numchan != idata.num_chan) ||
+            (fabs(obsmask.mjd - (idata.mjd_i + idata.mjd_f)) > 1e-9)) {
+                printf("WARNING!: maskfile has different number of channels or start MJD than raw data! Exiting.\n\n");
+                exit(1);
+        }
         good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
     } else {
         obsmask.numchan = obsmask.numint = 0;

diff --git a/src/prepfold.c b/src/prepfold.c
@@ -285,6 +285,11 @@ int main(int argc, char *argv[])
         if (cmd->maskfileP) {
             read_mask(cmd->maskfile, &obsmask);
             printf("Read mask information from '%s'\n\n", cmd->maskfile);
+            if ((obsmask.numchan != idata.num_chan) ||
+                (fabs(obsmask.mjd - (idata.mjd_i + idata.mjd_f)) > 1e-9)) {
+                    printf("WARNING!: maskfile has different number of channels or start MJD than raw data! Exiting.\n\n");
+                    exit(1);
+            }
             good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
         } else {
             obsmask.numchan = obsmask.numint = 0;
@@ -328,6 +333,11 @@ int main(int argc, char *argv[])
             if (cmd->maskfileP) {
                 read_mask(cmd->maskfile, &obsmask);
                 printf("Read mask information from '%s'\n\n", cmd->maskfile);
+                if ((obsmask.numchan != idata.num_chan) ||
+                    (fabs(obsmask.mjd - (idata.mjd_i + idata.mjd_f)) > 1e-9)) {
+                        printf("WARNING!: maskfile has different number of channels or start MJD than raw data! Exiting.\n\n");
+                        exit(1);
+                }
                 good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
             } else {
                 obsmask.numchan = obsmask.numint = 0;

diff --git a/src/prepsubband.c b/src/prepsubband.c
@@ -209,6 +209,11 @@ int main(int argc, char *argv[])
         if (cmd->maskfileP) {
             read_mask(cmd->maskfile, &obsmask);
             printf("Read mask information from '%s'\n\n", cmd->maskfile);
+            if ((obsmask.numchan != idata.num_chan) ||
+                (fabs(obsmask.mjd - (idata.mjd_i + idata.mjd_f)) > 1e-9)) {
+                    printf("WARNING!: maskfile has different number of channels or start MJD than raw data! Exiting.\n\n");
+                    exit(1);
+            }
             good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
         } else {
             obsmask.numchan = obsmask.numint = 0;