wraparound in state2obs

ggebbie · Nov 3, 2021 · db0904a · db0904a
1 parent 8403614
commit db0904a
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/TMI.jl b/src/TMI.jl
@@ -977,7 +977,9 @@ function state2obs(cvec,wis,γ)
     c̃ = tracerinit(γ.wet)
     c̃[γ.wet] = cvec
     replace!(c̃,NaN=>0.0)
-    [ỹ[i] = c̃[wis[i]...]/sumwis[i] for i in 1:N]
+    cwrap = view(c̃,list,:,:)
+
+    [ỹ[i] = cwrap[wis[i]...]/sumwis[i] for i in eachindex(wis)]
     return ỹ
 end