diff --git a/experiments/long_runs/land_region.jl b/experiments/long_runs/land_region.jl
index 7dfe7962fa..a79aaddad0 100644
--- a/experiments/long_runs/land_region.jl
+++ b/experiments/long_runs/land_region.jl
@@ -660,7 +660,7 @@ setup_and_solve_problem(; greet = true);
 # read in diagnostics and make some plots!
 #### ClimaAnalysis ####
 simdir = ClimaAnalysis.SimDir(outdir)
-short_names = ["gpp", "ct", "lai", "swc", "si"]
+short_names = ["gpp", "ct", "swc", "si"]
 for short_name in short_names
     var = get(simdir; short_name)
     times = ClimaAnalysis.times(var)
@@ -668,7 +668,9 @@ for short_name in short_names
         fig = CairoMakie.Figure(size = (800, 600))
         kwargs = ClimaAnalysis.has_altitude(var) ? Dict(:z => 1) : Dict()
         tmp = ClimaAnalysis.slice(var, time = t; kwargs...)
-        viz.heatmap2D!(fig, tmp)
-        CairoMakie.save(joinpath(root_path, "$(short_name)_$t.png"), fig)
+        if ~all(isnan.(tmp.data))
+            viz.heatmap2D!(fig, tmp)
+            CairoMakie.save(joinpath(root_path, "$(short_name)_$t.png"), fig)
+        end
     end
 end