diff --git a/parm/jcb-gdas b/parm/jcb-gdas
index 0b41d937b..a74480b57 160000
--- a/parm/jcb-gdas
+++ b/parm/jcb-gdas
@@ -1 +1 @@
-Subproject commit 0b41d937b3697aab37a6088f02ba28668f989abe
+Subproject commit a74480b57fc5dc04ea77eb21812081599f7a547e
diff --git a/utils/fv3jedi/fv3jedi_analcalc.h b/utils/fv3jedi/fv3jedi_analcalc.h
index dfa5aacab..1181b221e 100644
--- a/utils/fv3jedi/fv3jedi_analcalc.h
+++ b/utils/fv3jedi/fv3jedi_analcalc.h
@@ -97,8 +97,8 @@ namespace gdasapp {
         dx.read(incrConfig);
 
         // Initialize ensemble mean analysis
-        fv3jedi::State xxAnlEnsMean(anlEnsMeanGeom, stateVars, currentCycle);
-        dx.read(anlEnsMeanConfig);
+        fv3jedi::State xxAnlEnsMean(anlEnsMeanGeom, incrVars, currentCycle);
+        xxAnlEnsMean.read(anlEnsMeanConfig);
 
         // Perform calcuations
         // -------------------
@@ -107,8 +107,8 @@ namespace gdasapp {
         fv3jedi::State xxAnlDet(anlDetGeom, xxBkg);
         xxAnlDet += dx;
 
-        // Interpolate full resolution analysis to ensemble resolution
-        fv3jedi::State xxAnlDetEnsRes(incrCorGeom, xxAnlDet);
+        // Interpolate full resolution analysis to ensemble resolution and then change variables
+        fv3jedi::State xxAnlDetEnsRes(incrCorGeom, fv3jedi::State(incrVars, xxAnlDet));
 
         // Compute correction increment
         fv3jedi::Increment dxCor(incrCorGeom, incrVars, xxBkg.validTime());