Fix nh pressure index, as this is indexing into OffsetArray

CliMA · Jul 20, 2021 · 471fedd · 471fedd
1 parent ec66f9e
commit 471fedd
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Turbulence_PrognosticTKE.jl b/src/Turbulence_PrognosticTKE.jl
@@ -1135,7 +1135,7 @@ function compute_nh_pressure(self)
 
     @inbounds for i in xrange(self.n_updrafts)
         input.updraft_top = self.UpdVar.updraft_top[i]
-        alen = length(argwhere(self.UpdVar.Area.values[i, cinterior]))
+        alen = max(length(argwhere(self.UpdVar.Area.values[i, cinterior]))-1, 0)
         avals = off_arr(self.UpdVar.Area.values[i, cinterior])
         input.a_med = Statistics.median(avals[0:alen])
         @inbounds for k in xrange(grid(self).gw, grid(self).nzg - grid(self).gw)