diff --git a/src/input_xml.F90 b/src/input_xml.F90 index c7bf3ef17a..869f9fd63f 100644 --- a/src/input_xml.F90 +++ b/src/input_xml.F90 @@ -2647,6 +2647,12 @@ contains ! Get index of mesh filter k = t % find_filter(FILTER_MESH) + ! Check to make sure mesh filter was specified + if (k == 0) then + message = "Cannot tally surface current without a mesh filter." + call fatal_error() + end if + ! Get pointer to mesh i_mesh = t % filters(k) % int_bins(1) m => meshes(i_mesh)