diff --git a/testing/vcs/CMakeLists.txt b/testing/vcs/CMakeLists.txt
index 07ec860a23..6065f75135 100644
--- a/testing/vcs/CMakeLists.txt
+++ b/testing/vcs/CMakeLists.txt
@@ -981,6 +981,12 @@ cdat_add_test(test_vcs_settings_color_name_rgba
      "${UVCDAT_GIT_TESTDATA_DIR}/data/masked-file.nc"
      3
      )
+   cdat_add_test(test_vcs_isofill_masked_bounds
+     "${PYTHON_EXECUTABLE}"
+     ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_isofill_masked_bounds.py
+     "${BASELINE_DIR}/test_vcs_isofill_masked_bounds.png"
+     "${UVCDAT_GIT_TESTDATA_DIR}/data/coads_climatology.nc"
+     )
   cdat_add_test(test_vcs_issue_960_labels
     "${PYTHON_EXECUTABLE}"
     ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_issue_960_labels.py
diff --git a/testing/vcs/test_vcs_isofill_masked_bounds.py b/testing/vcs/test_vcs_isofill_masked_bounds.py
new file mode 100644
index 0000000000..6d165bcd8f
--- /dev/null
+++ b/testing/vcs/test_vcs_isofill_masked_bounds.py
@@ -0,0 +1,14 @@
+import cdms2
+import os
+import sys
+import vcs.testing.regression as regression
+import vcs
+import numpy
+
+data = sys.argv[2]
+x=regression.init(bg=1)
+f=cdms2.open(data)
+s=f["SST"]
+iso=x.getisofill('a_robinson_isofill')
+x.plot(s,iso)
+regression.run(x, "test_vcs_isofill_masked_bounds.png")