diff --git a/src/tracer/MOM_neutral_diffusion.F90 b/src/tracer/MOM_neutral_diffusion.F90 index 4e7b87886a..47f346a9d1 100644 --- a/src/tracer/MOM_neutral_diffusion.F90 +++ b/src/tracer/MOM_neutral_diffusion.F90 @@ -157,7 +157,7 @@ logical function neutral_diffusion_init(Time, G, param_file, diag, EOS, diabatic "That is, the algorithm will exclude the surface and bottom"//& "boundary layers.",default = .false.) - if (CS%continuous_reconstruction == .true. .and. CS%interior_only) then + if (CS%continuous_reconstruction .and. CS%interior_only) then call MOM_error(FATAL,"NDIFF_INTERIOR_ONLY=True only works with discontinuous" //& "reconstruction.") endif