diff --git a/code/code_gen/post.py b/code/code_gen/post.py
index 6a0c4bf..0431106 100644
--- a/code/code_gen/post.py
+++ b/code/code_gen/post.py
@@ -472,13 +472,13 @@ class Stage3():
 
         pass1calc = set(map(repr, self.avg1))
         map (pass1calc.update, map(AveragedField.depClosure, self.avg1))
-        self.pass1 = self.sort_vars(self.dependency, pass1calc - self.primary)
+        self.pass1 = self.sort_vars_new(self.dependency, pass1calc - self.primary)
 
         # self.sort_vars(self.dependency, pass1calc - self.primary -set(map(repr, self.avg1)))
 
         pass2calc = set(map(repr, self.avg2))
         map (pass2calc.update, map(AveragedField.depClosure, self.avg2))
-        self.pass2 = self.sort_vars(self.dependency, pass2calc - self.primary)
+        self.pass2 = self.sort_vars_new(self.dependency, pass2calc - self.primary)
 
 
     def __repr__ (self):
@@ -498,6 +498,48 @@ class Stage3():
 
         return order
 
+    def calc_size (self, ordered, remaining):
+        count = 0
+        dep_union = set([])
+
+        g = self.dependency
+        
+        for v in remaining:
+            dep_union |= set(g[v])
+            
+        for v in ordered:
+            if v in dep_union:
+                count += 1
+            
+        return count
+
+
+    def sort_vars_new (self, dependency, group):
+        order = []
+        remain = list(group)
+        remain.sort()
+
+        while len(remain) > 0:
+            
+            candidate = []
+            for v in remain:
+                if set(dependency[v]).isdisjoint(remain):
+                    candidate.append(v)
+
+            impact = {}
+            
+            size0 = self.calc_size(set(order), set(remain))
+
+            for v in candidate:
+                impact[v] = self.calc_size(set(order) | set([v]), set(remain) - set([v])) - size0
+                
+            candidate.sort(key=impact.get)
+            
+            order.append(candidate[0])
+            remain.remove(candidate[0])
+
+        return order
+
 
     def print_program (self):
 
diff --git a/code/m_terms.f90 b/code/m_terms.f90
index 7d9b376..e31949b 100644
--- a/code/m_terms.f90
+++ b/code/m_terms.f90
@@ -174,23 +174,6 @@ real*8, allocatable, dimension(:,:,:) :: xyzbuffer16
 real*8, allocatable, dimension(:,:,:) :: xyzbuffer17
 real*8, allocatable, dimension(:,:,:) :: xyzbuffer18
 real*8, allocatable, dimension(:,:,:) :: xyzbuffer19
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer20
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer21
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer22
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer23
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer24
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer25
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer26
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer27
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer28
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer29
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer30
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer31
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer32
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer33
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer34
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer35
-real*8, allocatable, dimension(:,:,:) :: xyzbuffer36
 
 contains
 
@@ -364,23 +347,6 @@ allocate(xyzbuffer16(nxp,nyp,nzp), stat=ierr) ; xyzbuffer16 = 0.
 allocate(xyzbuffer17(nxp,nyp,nzp), stat=ierr) ; xyzbuffer17 = 0.
 allocate(xyzbuffer18(nxp,nyp,nzp), stat=ierr) ; xyzbuffer18 = 0.
 allocate(xyzbuffer19(nxp,nyp,nzp), stat=ierr) ; xyzbuffer19 = 0.
-allocate(xyzbuffer20(nxp,nyp,nzp), stat=ierr) ; xyzbuffer20 = 0.
-allocate(xyzbuffer21(nxp,nyp,nzp), stat=ierr) ; xyzbuffer21 = 0.
-allocate(xyzbuffer22(nxp,nyp,nzp), stat=ierr) ; xyzbuffer22 = 0.
-allocate(xyzbuffer23(nxp,nyp,nzp), stat=ierr) ; xyzbuffer23 = 0.
-allocate(xyzbuffer24(nxp,nyp,nzp), stat=ierr) ; xyzbuffer24 = 0.
-allocate(xyzbuffer25(nxp,nyp,nzp), stat=ierr) ; xyzbuffer25 = 0.
-allocate(xyzbuffer26(nxp,nyp,nzp), stat=ierr) ; xyzbuffer26 = 0.
-allocate(xyzbuffer27(nxp,nyp,nzp), stat=ierr) ; xyzbuffer27 = 0.
-allocate(xyzbuffer28(nxp,nyp,nzp), stat=ierr) ; xyzbuffer28 = 0.
-allocate(xyzbuffer29(nxp,nyp,nzp), stat=ierr) ; xyzbuffer29 = 0.
-allocate(xyzbuffer30(nxp,nyp,nzp), stat=ierr) ; xyzbuffer30 = 0.
-allocate(xyzbuffer31(nxp,nyp,nzp), stat=ierr) ; xyzbuffer31 = 0.
-allocate(xyzbuffer32(nxp,nyp,nzp), stat=ierr) ; xyzbuffer32 = 0.
-allocate(xyzbuffer33(nxp,nyp,nzp), stat=ierr) ; xyzbuffer33 = 0.
-allocate(xyzbuffer34(nxp,nyp,nzp), stat=ierr) ; xyzbuffer34 = 0.
-allocate(xyzbuffer35(nxp,nyp,nzp), stat=ierr) ; xyzbuffer35 = 0.
-allocate(xyzbuffer36(nxp,nyp,nzp), stat=ierr) ; xyzbuffer36 = 0.
 
 end subroutine m_terms_init
 
@@ -553,23 +519,6 @@ deallocate(xyzbuffer16)
 deallocate(xyzbuffer17)
 deallocate(xyzbuffer18)
 deallocate(xyzbuffer19)
-deallocate(xyzbuffer20)
-deallocate(xyzbuffer21)
-deallocate(xyzbuffer22)
-deallocate(xyzbuffer23)
-deallocate(xyzbuffer24)
-deallocate(xyzbuffer25)
-deallocate(xyzbuffer26)
-deallocate(xyzbuffer27)
-deallocate(xyzbuffer28)
-deallocate(xyzbuffer29)
-deallocate(xyzbuffer30)
-deallocate(xyzbuffer31)
-deallocate(xyzbuffer32)
-deallocate(xyzbuffer33)
-deallocate(xyzbuffer34)
-deallocate(xyzbuffer35)
-deallocate(xyzbuffer36)
 
 end subroutine m_terms_finalize
 
@@ -591,7 +540,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_w(i) = avg_w(i) + w(i,j,k)
+avg_v(i) = avg_v(i) + v(i,j,k)
 end do
 end do
 end do
@@ -600,7 +549,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_v(i) = avg_v(i) + v(i,j,k)
+avg_w(i) = avg_w(i) + w(i,j,k)
 end do
 end do
 end do
@@ -615,90 +564,6 @@ end do
 end do
 end do
 
-call d2dy ( xyzbuffer18, xyzbuffer19 )
-call ddx ( xyzbuffer13, xyzbuffer19 )
-call ddy ( xyzbuffer12, xyzbuffer19 )
-call ddz ( xyzbuffer11, xyzbuffer19 )
-
-! ( sqrt ( ( ( ((ddx_c(i,j,k))*(ddx_c(i,j,k))) + ((ddy_c(i,j,k))*(ddy_c(i,j,k))) ) + ((ddz_c(i,j,k))*(ddz_c(i,j,k))) ) ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer10(i,j,k) = ( sqrt ( ( ( ((xyzbuffer13(i,j,k))*(xyzbuffer13(i,j,k))) + ((xyzbuffer12(i,j,k))*(xyzbuffer12(i,j,k))) ) + ((xyzbuffer11(i,j,k))*(xyzbuffer11(i,j,k))) ) ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer17(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.05, 0.15 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.25, 0.35 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer15(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.45, 0.55 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer14(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.65, 0.75 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer9(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.85, 0.95 ) )
-end do
-end do
-end do
-
-
-! ( ( - ddy_c(i,j,k) ) / fsd_orig(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer8(i,j,k) = ( ( - xyzbuffer12(i,j,k) ) / xyzbuffer10(i,j,k) )
-end do
-end do
-end do
-
-
-! ( rxn_rate ( c(i,j,k) ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer3(i,j,k) = ( rxn_rate ( xyzbuffer19(i,j,k) ) )
-end do
-end do
-end do
-
 
 do k = 1, nzp
 do j = 1, nyp
@@ -708,146 +573,72 @@ end do
 end do
 end do
 
+call d2dx ( xyzbuffer18, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_d2dy_c(i) = avg_d2dy_c(i) + xyzbuffer18(i,j,k)
+avg_d2dx_c(i) = avg_d2dx_c(i) + xyzbuffer18(i,j,k)
 end do
 end do
 end do
 
+call d2dy ( xyzbuffer13, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_ddx_c(i) = avg_ddx_c(i) + xyzbuffer13(i,j,k)
+avg_d2dy_c(i) = avg_d2dy_c(i) + xyzbuffer13(i,j,k)
 end do
 end do
 end do
 
+call d2dz ( xyzbuffer12, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_ddz_c(i) = avg_ddz_c(i) + xyzbuffer11(i,j,k)
-end do
-end do
-end do
-
-call d2dx ( xyzbuffer2, xyzbuffer19 )
-call ddy ( xyzbuffer1, xyzbuffer8 )
-
-! ( fsd_orig(i,j,k) * mask1(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer0(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer17(i,j,k) )
+avg_d2dz_c(i) = avg_d2dz_c(i) + xyzbuffer12(i,j,k)
 end do
 end do
 end do
 
+call ddx ( xyzbuffer11, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_ny(i) = fsd1_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask3(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer7(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer16(i,j,k) )
+avg_ddx_c(i) = avg_ddx_c(i) + xyzbuffer11(i,j,k)
 end do
 end do
 end do
 
+call ddy ( xyzbuffer10, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd3_avg_ny(i) = fsd3_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask5(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer6(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer15(i,j,k) )
+avg_ddy_c(i) = avg_ddy_c(i) + xyzbuffer10(i,j,k)
 end do
 end do
 end do
 
+call ddz ( xyzbuffer17, xyzbuffer19 )
 
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd5_avg_ny(i) = fsd5_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer6(i,j,k)
+avg_ddz_c(i) = avg_ddz_c(i) + xyzbuffer17(i,j,k)
 end do
 end do
 end do
 
 
-! ( fsd_orig(i,j,k) * mask7(i,j,k) )
+! ( sqrt ( ( ( ((ddx_c(i,j,k))*(ddx_c(i,j,k))) + ((ddy_c(i,j,k))*(ddy_c(i,j,k))) ) + ((ddz_c(i,j,k))*(ddz_c(i,j,k))) ) ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer5(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer14(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_ny(i) = fsd7_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask9(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer4(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer9(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_ny(i) = fsd9_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer4(i,j,k)
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer17(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.15, 0.25 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.55, 0.65 ) )
+xyzbuffer16(i,j,k) = ( sqrt ( ( ( ((xyzbuffer11(i,j,k))*(xyzbuffer11(i,j,k))) + ((xyzbuffer10(i,j,k))*(xyzbuffer10(i,j,k))) ) + ((xyzbuffer17(i,j,k))*(xyzbuffer17(i,j,k))) ) ) )
 end do
 end do
 end do
@@ -857,385 +648,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer15(i,j,k) = ( ( - xyzbuffer13(i,j,k) ) / xyzbuffer10(i,j,k) )
+xyzbuffer15(i,j,k) = ( ( - xyzbuffer11(i,j,k) ) / xyzbuffer16(i,j,k) )
 end do
 end do
 end do
 
 
+! ( ( - ddy_c(i,j,k) ) / fsd_orig(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-avg_d2dx_c(i) = avg_d2dx_c(i) + xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_ddy_c(i) = avg_ddy_c(i) + xyzbuffer12(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd3(i) = avg_fsd3(i) + xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd5(i) = avg_fsd5(i) + xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd7(i) = avg_fsd7(i) + xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd9(i) = avg_fsd9(i) + xyzbuffer4(i,j,k)
-end do
-end do
-end do
-
-call d2dz ( xyzbuffer14, xyzbuffer19 )
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_nx(i) = fsd1_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask2(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer9(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer17(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_nx(i) = fsd2_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_nx(i) = fsd3_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_nx(i) = fsd5_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask6(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer13(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer16(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_nx(i) = fsd6_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_nx(i) = fsd7_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_nx(i) = fsd9_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k)
-end do
-end do
-end do
-
-
-! ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer12(i,j,k) = ( ( xyzbuffer2(i,j,k) + xyzbuffer18(i,j,k) ) + xyzbuffer14(i,j,k) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer17(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.75, 0.85 ) )
-end do
-end do
-end do
-
-
-! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer16(i,j,k) = ( ( ( ( ( xyzbuffer2(i,j,k) + xyzbuffer18(i,j,k) ) + xyzbuffer14(i,j,k) ) * rod ) + xyzbuffer3(i,j,k) ) / xyzbuffer10(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_d2dz_c(i) = avg_d2dz_c(i) + xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd2(i) = avg_fsd2(i) + xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd6(i) = avg_fsd6(i) + xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_lapc(i) = avg_lapc(i) + xyzbuffer12(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_sd(i) = fsd1_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_ny(i) = fsd2_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_sd(i) = fsd2_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_sd(i) = fsd3_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_sd(i) = fsd5_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_ny(i) = fsd6_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_sd(i) = fsd6_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_sd(i) = fsd7_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask8(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer14(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer17(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_nx(i) = fsd8_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_sd(i) = fsd8_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_sd(i) = fsd9_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k)
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer3(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.35, 0.45 ) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd1(i) = avg_fsd1(i) + xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd8(i) = avg_fsd8(i) + xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask4(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer2(i,j,k) = ( xyzbuffer10(i,j,k) * xyzbuffer3(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_nx(i) = fsd4_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_sd(i) = fsd4_avg_sd(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_ny(i) = fsd8_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer14(i,j,k)
+xyzbuffer14(i,j,k) = ( ( - xyzbuffer10(i,j,k) ) / xyzbuffer16(i,j,k) )
 end do
 end do
 end do
@@ -1245,291 +668,20 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer19(i,j,k) = ( ( - xyzbuffer11(i,j,k) ) / xyzbuffer10(i,j,k) )
-end do
-end do
-end do
-
-
-! ( ( ( u(i,j,k) * nx(i,j,k) ) + ( v(i,j,k) * ny(i,j,k) ) ) + ( w(i,j,k) * nz(i,j,k) ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer18(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer15(i,j,k) ) + ( v(i,j,k) * xyzbuffer8(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer19(i,j,k) ) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-avg_fsd4(i) = avg_fsd4(i) + xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-call ddz ( xyzbuffer12, xyzbuffer19 )
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_nz(i) = fsd1_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_nz(i) = fsd2_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_nz(i) = fsd3_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_ny(i) = fsd4_avg_ny(i) + xyzbuffer8(i,j,k) * xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_nz(i) = fsd5_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_nz(i) = fsd6_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_nz(i) = fsd7_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_nz(i) = fsd8_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_nz(i) = fsd9_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer4(i,j,k)
-end do
-end do
-end do
-
-
-! ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer11(i,j,k) = ( ( xyzbuffer18(i,j,k) + xyzbuffer16(i,j,k) ) * xyzbuffer15(i,j,k) )
-end do
-end do
-end do
-
-call ddx ( xyzbuffer10, xyzbuffer15 )
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_t1(i) = fsd1_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t1(i) = fsd2_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_t1(i) = fsd3_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_nz(i) = fsd4_avg_nz(i) + xyzbuffer19(i,j,k) * xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t1(i) = fsd5_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t1(i) = fsd6_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t1(i) = fsd7_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t1(i) = fsd8_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t1(i) = fsd9_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer4(i,j,k)
+xyzbuffer9(i,j,k) = ( ( - xyzbuffer17(i,j,k) ) / xyzbuffer16(i,j,k) )
 end do
 end do
 end do
 
+call ddx ( xyzbuffer8, xyzbuffer15 )
+call ddy ( xyzbuffer3, xyzbuffer14 )
+call ddz ( xyzbuffer2, xyzbuffer9 )
 
 ! ( ( ddx_nx(i,j,k) + ddy_ny(i,j,k) ) + ddz_nz(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer17(i,j,k) = ( ( xyzbuffer10(i,j,k) + xyzbuffer1(i,j,k) ) + xyzbuffer12(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_divn(i) = fsd1_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer0(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_divn(i) = fsd2_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_divn(i) = fsd3_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_divn(i) = fsd4_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer2(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_divn(i) = fsd5_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_divn(i) = fsd6_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer13(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_divn(i) = fsd7_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_divn(i) = fsd8_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_divn(i) = fsd9_avg_divn(i) + xyzbuffer17(i,j,k) * xyzbuffer4(i,j,k)
+xyzbuffer1(i,j,k) = ( ( xyzbuffer8(i,j,k) + xyzbuffer3(i,j,k) ) + xyzbuffer2(i,j,k) )
 end do
 end do
 end do
@@ -1539,7 +691,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer16(i,j,k) = ( dabs ( xyzbuffer17(i,j,k) ) )
+xyzbuffer0(i,j,k) = ( dabs ( xyzbuffer1(i,j,k) ) )
+end do
+end do
+end do
+
+
+! ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer7(i,j,k) = ( ( xyzbuffer18(i,j,k) + xyzbuffer13(i,j,k) ) + xyzbuffer12(i,j,k) )
 end do
 end do
 end do
@@ -1548,7 +710,27 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd2_avg_absk(i) = fsd2_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k)
+avg_lapc(i) = avg_lapc(i) + xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer6(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.05, 0.15 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask1(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer5(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k) )
 end do
 end do
 end do
@@ -1557,7 +739,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd4_avg_absk(i) = fsd4_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k)
+avg_fsd1(i) = avg_fsd1(i) + xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -1566,7 +748,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd5_avg_absk(i) = fsd5_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k)
+fsd1_avg_absk(i) = fsd1_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -1575,7 +757,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd7_avg_absk(i) = fsd7_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer5(i,j,k)
+fsd1_avg_divn(i) = fsd1_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -1584,7 +766,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd9_avg_absk(i) = fsd9_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k)
+fsd1_avg_nx(i) = fsd1_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -1593,7 +775,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_absk(i) = fsd1_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer0(i,j,k)
+fsd1_avg_ny(i) = fsd1_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -1602,7 +784,27 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd4_avg_t1(i) = fsd4_avg_t1(i) + xyzbuffer11(i,j,k) * xyzbuffer2(i,j,k)
+fsd1_avg_nz(i) = fsd1_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer4(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.15, 0.25 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask2(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k) )
 end do
 end do
 end do
@@ -1611,7 +813,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd8_avg_absk(i) = fsd8_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer14(i,j,k)
+avg_fsd2(i) = avg_fsd2(i) + xyzbuffer11(i,j,k)
 end do
 end do
 end do
@@ -1620,7 +822,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd3_avg_absk(i) = fsd3_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer7(i,j,k)
+fsd2_avg_absk(i) = fsd2_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer11(i,j,k)
 end do
 end do
 end do
@@ -1629,7 +831,754 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd6_avg_absk(i) = fsd6_avg_absk(i) + xyzbuffer16(i,j,k) * xyzbuffer13(i,j,k)
+fsd2_avg_divn(i) = fsd2_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_nx(i) = fsd2_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_ny(i) = fsd2_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_nz(i) = fsd2_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.25, 0.35 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask3(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer17(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd3(i) = avg_fsd3(i) + xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_absk(i) = fsd3_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_divn(i) = fsd3_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_nx(i) = fsd3_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_ny(i) = fsd3_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_nz(i) = fsd3_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer8(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.35, 0.45 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask4(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer3(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer8(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd4(i) = avg_fsd4(i) + xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_absk(i) = fsd4_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_divn(i) = fsd4_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_nx(i) = fsd4_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_ny(i) = fsd4_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_nz(i) = fsd4_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.45, 0.55 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask5(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer7(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd5(i) = avg_fsd5(i) + xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_absk(i) = fsd5_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_divn(i) = fsd5_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_nx(i) = fsd5_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_ny(i) = fsd5_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_nz(i) = fsd5_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer6(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.55, 0.65 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask6(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer4(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd6(i) = avg_fsd6(i) + xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_absk(i) = fsd6_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_divn(i) = fsd6_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_nx(i) = fsd6_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_ny(i) = fsd6_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_nz(i) = fsd6_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.65, 0.75 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask7(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer8(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd7(i) = avg_fsd7(i) + xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_absk(i) = fsd7_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_divn(i) = fsd7_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_nx(i) = fsd7_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_ny(i) = fsd7_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_nz(i) = fsd7_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.75, 0.85 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask8(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer6(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd8(i) = avg_fsd8(i) + xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_absk(i) = fsd8_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_divn(i) = fsd8_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_nx(i) = fsd8_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_ny(i) = fsd8_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_nz(i) = fsd8_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.85, 0.95 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask9(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_absk(i) = fsd9_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_divn(i) = fsd9_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+avg_fsd9(i) = avg_fsd9(i) + xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_nx(i) = fsd9_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_ny(i) = fsd9_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_nz(i) = fsd9_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! ( ( ( u(i,j,k) * nx(i,j,k) ) + ( v(i,j,k) * ny(i,j,k) ) ) + ( w(i,j,k) * nz(i,j,k) ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer15(i,j,k) ) + ( v(i,j,k) * xyzbuffer14(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer9(i,j,k) ) )
+end do
+end do
+end do
+
+
+! ( rxn_rate ( c(i,j,k) ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer0(i,j,k) = ( rxn_rate ( xyzbuffer19(i,j,k) ) )
+end do
+end do
+end do
+
+
+! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer19(i,j,k) = ( ( ( ( ( xyzbuffer18(i,j,k) + xyzbuffer13(i,j,k) ) + xyzbuffer12(i,j,k) ) * rod ) + xyzbuffer0(i,j,k) ) / xyzbuffer16(i,j,k) )
+end do
+end do
+end do
+
+
+! ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( ( xyzbuffer1(i,j,k) + xyzbuffer19(i,j,k) ) * xyzbuffer15(i,j,k) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_sd(i) = fsd1_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t1(i) = fsd1_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_sd(i) = fsd2_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_t1(i) = fsd2_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer11(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_sd(i) = fsd3_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t1(i) = fsd3_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_sd(i) = fsd4_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t1(i) = fsd4_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_sd(i) = fsd5_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t1(i) = fsd5_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer7(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_sd(i) = fsd6_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t1(i) = fsd6_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_sd(i) = fsd7_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t1(i) = fsd7_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_sd(i) = fsd8_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t1(i) = fsd8_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer6(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_sd(i) = fsd9_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t1(i) = fsd9_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k)
 end do
 end do
 end do
@@ -2073,61 +2022,23 @@ integer :: i, j, k
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer22(i,j,k) = ( 1.0 - y(i,j,k) )
+xyzbuffer18(i,j,k) = ( 1.0 - y(i,j,k) )
 end do
 end do
 end do
 
-call d2dy ( xyzbuffer23, xyzbuffer22 )
-call ddx ( xyzbuffer20, xyzbuffer22 )
-call ddy ( xyzbuffer21, xyzbuffer22 )
-call ddz ( xyzbuffer26, xyzbuffer22 )
+call d2dx ( xyzbuffer13, xyzbuffer18 )
+call d2dy ( xyzbuffer12, xyzbuffer18 )
+call d2dz ( xyzbuffer11, xyzbuffer18 )
+call ddx ( xyzbuffer10, xyzbuffer18 )
+call ddy ( xyzbuffer17, xyzbuffer18 )
+call ddz ( xyzbuffer16, xyzbuffer18 )
 
 ! ( sqrt ( ( ( ((ddx_c(i,j,k))*(ddx_c(i,j,k))) + ((ddy_c(i,j,k))*(ddy_c(i,j,k))) ) + ((ddz_c(i,j,k))*(ddz_c(i,j,k))) ) ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer27(i,j,k) = ( sqrt ( ( ( ((xyzbuffer20(i,j,k))*(xyzbuffer20(i,j,k))) + ((xyzbuffer21(i,j,k))*(xyzbuffer21(i,j,k))) ) + ((xyzbuffer26(i,j,k))*(xyzbuffer26(i,j,k))) ) ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.15, 0.25 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer25(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.35, 0.45 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer28(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.55, 0.65 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer29(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.75, 0.85 ) )
+xyzbuffer15(i,j,k) = ( sqrt ( ( ( ((xyzbuffer10(i,j,k))*(xyzbuffer10(i,j,k))) + ((xyzbuffer17(i,j,k))*(xyzbuffer17(i,j,k))) ) + ((xyzbuffer16(i,j,k))*(xyzbuffer16(i,j,k))) ) ) )
 end do
 end do
 end do
@@ -2137,256 +2048,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer31(i,j,k) = ( ( - xyzbuffer20(i,j,k) ) / xyzbuffer27(i,j,k) )
-end do
-end do
-end do
-
-
-! ( ( - ddz_c(i,j,k) ) / fsd_orig(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer35(i,j,k) = ( ( - xyzbuffer26(i,j,k) ) / xyzbuffer27(i,j,k) )
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd1_avg_nx(i)))*((nx(i,j,k) - fsd1_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer34(i,j,k) = (((xyzbuffer31(i,j,k) - fsd1_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd1_avg_nx(i))))
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd3_avg_nx(i)))*((nx(i,j,k) - fsd3_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer30(i,j,k) = (((xyzbuffer31(i,j,k) - fsd3_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd3_avg_nx(i))))
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd5_avg_nx(i)))*((nx(i,j,k) - fsd5_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer19(i,j,k) = (((xyzbuffer31(i,j,k) - fsd5_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd5_avg_nx(i))))
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd7_avg_nx(i)))*((nx(i,j,k) - fsd7_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer18(i,j,k) = (((xyzbuffer31(i,j,k) - fsd7_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd7_avg_nx(i))))
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd9_avg_nx(i)))*((nx(i,j,k) - fsd9_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer13(i,j,k) = (((xyzbuffer31(i,j,k) - fsd9_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd9_avg_nx(i))))
-end do
-end do
-end do
-
-
-! ( rxn_rate ( c(i,j,k) ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer12(i,j,k) = ( rxn_rate ( xyzbuffer22(i,j,k) ) )
-end do
-end do
-end do
-
-call d2dx ( xyzbuffer11, xyzbuffer22 )
-call ddx ( xyzbuffer10, xyzbuffer31 )
-call ddz ( xyzbuffer17, xyzbuffer35 )
-
-! ( fsd_orig(i,j,k) * mask2(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer16(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer24(i,j,k) )
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask4(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer15(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer25(i,j,k) )
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask6(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer14(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer28(i,j,k) )
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask8(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer9(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer29(i,j,k) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer8(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.05, 0.15 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer33(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.45, 0.55 ) )
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer3(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.85, 0.95 ) )
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd2_avg_nx(i)))*((nx(i,j,k) - fsd2_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer2(i,j,k) = (((xyzbuffer31(i,j,k) - fsd2_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd2_avg_nx(i))))
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd6_avg_nx(i)))*((nx(i,j,k) - fsd6_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer1(i,j,k) = (((xyzbuffer31(i,j,k) - fsd6_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd6_avg_nx(i))))
-end do
-end do
-end do
-
-call d2dz ( xyzbuffer0, xyzbuffer22 )
-
-! ( fsd_orig(i,j,k) * mask1(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer7(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer8(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_t6(i) = fsd1_avg_t6(i) + xyzbuffer34(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t6(i) = fsd2_avg_t6(i) + xyzbuffer2(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask5(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer6(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer33(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t6(i) = fsd5_avg_t6(i) + xyzbuffer19(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t6(i) = fsd6_avg_t6(i) + xyzbuffer1(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask9(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer5(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer3(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t6(i) = fsd9_avg_t6(i) + xyzbuffer13(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer4(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.25, 0.35 ) )
+xyzbuffer14(i,j,k) = ( ( - xyzbuffer10(i,j,k) ) / xyzbuffer15(i,j,k) )
 end do
 end do
 end do
@@ -2396,174 +2058,40 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer32(i,j,k) = ( ( - xyzbuffer21(i,j,k) ) / xyzbuffer27(i,j,k) )
+xyzbuffer9(i,j,k) = ( ( - xyzbuffer17(i,j,k) ) / xyzbuffer15(i,j,k) )
 end do
 end do
 end do
 
 
-! (((nx(i,j,k) - fsd4_avg_nx(i)))*((nx(i,j,k) - fsd4_avg_nx(i))))
+! ( ( - ddz_c(i,j,k) ) / fsd_orig(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer36(i,j,k) = (((xyzbuffer31(i,j,k) - fsd4_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd4_avg_nx(i))))
+xyzbuffer8(i,j,k) = ( ( - xyzbuffer16(i,j,k) ) / xyzbuffer15(i,j,k) )
 end do
 end do
 end do
 
+call ddx ( xyzbuffer3, xyzbuffer14 )
+call ddy ( xyzbuffer2, xyzbuffer9 )
 
 ! ( ( ( u(i,j,k) * nx(i,j,k) ) + ( v(i,j,k) * ny(i,j,k) ) ) + ( w(i,j,k) * nz(i,j,k) ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer20(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer31(i,j,k) ) + ( v(i,j,k) * xyzbuffer32(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer35(i,j,k) ) )
-end do
-end do
-end do
-
-call ddy ( xyzbuffer21, xyzbuffer32 )
-
-! ( fsd_orig(i,j,k) * mask3(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer26(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer4(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_t6(i) = fsd3_avg_t6(i) + xyzbuffer30(i,j,k) * xyzbuffer26(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t6(i) = fsd4_avg_t6(i) + xyzbuffer36(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = ( threshold_min_max ( xyzbuffer22(i,j,k), 0.65, 0.75 ) )
-end do
-end do
-end do
-
-
-! (((nx(i,j,k) - fsd8_avg_nx(i)))*((nx(i,j,k) - fsd8_avg_nx(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer22(i,j,k) = (((xyzbuffer31(i,j,k) - fsd8_avg_nx(i)))*((xyzbuffer31(i,j,k) - fsd8_avg_nx(i))))
+xyzbuffer1(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer14(i,j,k) ) + ( v(i,j,k) * xyzbuffer9(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer8(i,j,k) ) )
 end do
 end do
 end do
 
+call ddz ( xyzbuffer0, xyzbuffer8 )
 
 ! ( ( ddx_nx(i,j,k) + ddy_ny(i,j,k) ) + ddz_nz(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer25(i,j,k) = ( ( xyzbuffer10(i,j,k) + xyzbuffer21(i,j,k) ) + xyzbuffer17(i,j,k) )
-end do
-end do
-end do
-
-
-! ( fsd_orig(i,j,k) * mask7(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer28(i,j,k) = ( xyzbuffer27(i,j,k) * xyzbuffer24(i,j,k) )
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t6(i) = fsd7_avg_t6(i) + xyzbuffer18(i,j,k) * xyzbuffer28(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t6(i) = fsd8_avg_t6(i) + xyzbuffer22(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer22(i,j,k) = ( ( ( ( ( xyzbuffer11(i,j,k) + xyzbuffer23(i,j,k) ) + xyzbuffer0(i,j,k) ) * rod ) + xyzbuffer12(i,j,k) ) / xyzbuffer27(i,j,k) )
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd1_avg_divn(i)))*((divn(i,j,k) - fsd1_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer23(i,j,k) = (((xyzbuffer25(i,j,k) - fsd1_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd1_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd3_avg_divn(i)))*((divn(i,j,k) - fsd3_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer21(i,j,k) = (((xyzbuffer25(i,j,k) - fsd3_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd3_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd5_avg_divn(i)))*((divn(i,j,k) - fsd5_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer27(i,j,k) = (((xyzbuffer25(i,j,k) - fsd5_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd5_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd7_avg_divn(i)))*((divn(i,j,k) - fsd7_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = (((xyzbuffer25(i,j,k) - fsd7_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd7_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd9_avg_divn(i)))*((divn(i,j,k) - fsd9_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer29(i,j,k) = (((xyzbuffer25(i,j,k) - fsd9_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd9_avg_divn(i))))
+xyzbuffer7(i,j,k) = ( ( xyzbuffer3(i,j,k) + xyzbuffer2(i,j,k) ) + xyzbuffer0(i,j,k) )
 end do
 end do
 end do
@@ -2573,52 +2101,207 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer18(i,j,k) = ( dabs ( xyzbuffer25(i,j,k) ) )
+xyzbuffer6(i,j,k) = ( dabs ( xyzbuffer7(i,j,k) ) )
 end do
 end do
 end do
 
 
+! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_t7(i) = fsd1_avg_t7(i) + xyzbuffer23(i,j,k) * xyzbuffer7(i,j,k)
+xyzbuffer5(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.05, 0.15 ) )
 end do
 end do
 end do
 
 
+! ( fsd_orig(i,j,k) * mask1(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd3_avg_t7(i) = fsd3_avg_t7(i) + xyzbuffer21(i,j,k) * xyzbuffer26(i,j,k)
+xyzbuffer4(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k) )
 end do
 end do
 end do
 
 
+! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd5_avg_t7(i) = fsd5_avg_t7(i) + xyzbuffer27(i,j,k) * xyzbuffer6(i,j,k)
+xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.15, 0.25 ) )
 end do
 end do
 end do
 
 
+! ( fsd_orig(i,j,k) * mask2(i,j,k) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd7_avg_t7(i) = fsd7_avg_t7(i) + xyzbuffer24(i,j,k) * xyzbuffer28(i,j,k)
+xyzbuffer17(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer10(i,j,k) )
 end do
 end do
 end do
 
 
+! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) )
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd9_avg_t7(i) = fsd9_avg_t7(i) + xyzbuffer29(i,j,k) * xyzbuffer5(i,j,k)
+xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.25, 0.35 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask3(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer9(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer8(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.35, 0.45 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask4(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer3(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.45, 0.55 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask5(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer0(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer5(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.55, 0.65 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask6(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer10(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.65, 0.75 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask7(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer8(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.75, 0.85 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask8(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer5(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k) )
+end do
+end do
+end do
+
+
+! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.85, 0.95 ) )
+end do
+end do
+end do
+
+
+! ( fsd_orig(i,j,k) * mask9(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer2(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) )
+end do
+end do
+end do
+
+
+! ( rxn_rate ( c(i,j,k) ) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( rxn_rate ( xyzbuffer18(i,j,k) ) )
+end do
+end do
+end do
+
+
+! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( ( ( ( ( xyzbuffer13(i,j,k) + xyzbuffer12(i,j,k) ) + xyzbuffer11(i,j,k) ) * rod ) + xyzbuffer16(i,j,k) ) / xyzbuffer15(i,j,k) )
 end do
 end do
 end do
@@ -2628,445 +2311,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer35(i,j,k) = ( ( xyzbuffer20(i,j,k) + xyzbuffer22(i,j,k) ) * xyzbuffer31(i,j,k) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (nx(i,j,k) - fsd2_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer22(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd2_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (nx(i,j,k) - fsd4_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer23(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd4_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (nx(i,j,k) - fsd6_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer20(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd6_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (nx(i,j,k) - fsd8_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer21(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd8_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (divn(i,j,k) - fsd1_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer27(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd1_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (divn(i,j,k) - fsd3_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd3_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (divn(i,j,k) - fsd5_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer29(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd5_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (divn(i,j,k) - fsd7_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer12(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd7_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (divn(i,j,k) - fsd9_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer19(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd9_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (absk(i,j,k) - fsd2_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer10(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd2_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (absk(i,j,k) - fsd4_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer13(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd4_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (absk(i,j,k) - fsd6_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer34(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd6_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (absk(i,j,k) - fsd8_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer11(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd8_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd1_avg_t1(i)))*((t1(i,j,k) - fsd1_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer36(i,j,k) = (((xyzbuffer35(i,j,k) - fsd1_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd1_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd3_avg_t1(i)))*((t1(i,j,k) - fsd3_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer17(i,j,k) = (((xyzbuffer35(i,j,k) - fsd3_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd3_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd5_avg_t1(i)))*((t1(i,j,k) - fsd5_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer30(i,j,k) = (((xyzbuffer35(i,j,k) - fsd5_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd5_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd7_avg_t1(i)))*((t1(i,j,k) - fsd7_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer33(i,j,k) = (((xyzbuffer35(i,j,k) - fsd7_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd7_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd9_avg_t1(i)))*((t1(i,j,k) - fsd9_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer32(i,j,k) = (((xyzbuffer35(i,j,k) - fsd9_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd9_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd4_avg_divn(i)))*((divn(i,j,k) - fsd4_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer8(i,j,k) = (((xyzbuffer25(i,j,k) - fsd4_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd4_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd8_avg_divn(i)))*((divn(i,j,k) - fsd8_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer3(i,j,k) = (((xyzbuffer25(i,j,k) - fsd8_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd8_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd2_avg_absk(i)))*((absk(i,j,k) - fsd2_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer2(i,j,k) = (((xyzbuffer18(i,j,k) - fsd2_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd2_avg_absk(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd4_avg_absk(i)))*((absk(i,j,k) - fsd4_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer1(i,j,k) = (((xyzbuffer18(i,j,k) - fsd4_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd4_avg_absk(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd6_avg_absk(i)))*((absk(i,j,k) - fsd6_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer0(i,j,k) = (((xyzbuffer18(i,j,k) - fsd6_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd6_avg_absk(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd8_avg_absk(i)))*((absk(i,j,k) - fsd8_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer4(i,j,k) = (((xyzbuffer18(i,j,k) - fsd8_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd8_avg_absk(i))))
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_t3(i) = fsd1_avg_t3(i) + xyzbuffer27(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd1_avg_t5(i) = fsd1_avg_t5(i) + xyzbuffer36(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t2(i) = fsd2_avg_t2(i) + xyzbuffer22(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t4(i) = fsd2_avg_t4(i) + xyzbuffer10(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t8(i) = fsd2_avg_t8(i) + xyzbuffer2(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_t3(i) = fsd3_avg_t3(i) + xyzbuffer24(i,j,k) * xyzbuffer26(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_t5(i) = fsd3_avg_t5(i) + xyzbuffer17(i,j,k) * xyzbuffer26(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t2(i) = fsd4_avg_t2(i) + xyzbuffer23(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t4(i) = fsd4_avg_t4(i) + xyzbuffer13(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t7(i) = fsd4_avg_t7(i) + xyzbuffer8(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t3(i) = fsd5_avg_t3(i) + xyzbuffer29(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t5(i) = fsd5_avg_t5(i) + xyzbuffer30(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t2(i) = fsd6_avg_t2(i) + xyzbuffer20(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t4(i) = fsd6_avg_t4(i) + xyzbuffer34(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t8(i) = fsd6_avg_t8(i) + xyzbuffer0(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t3(i) = fsd7_avg_t3(i) + xyzbuffer12(i,j,k) * xyzbuffer28(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t5(i) = fsd7_avg_t5(i) + xyzbuffer33(i,j,k) * xyzbuffer28(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t2(i) = fsd8_avg_t2(i) + xyzbuffer21(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t4(i) = fsd8_avg_t4(i) + xyzbuffer11(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t7(i) = fsd8_avg_t7(i) + xyzbuffer3(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t3(i) = fsd9_avg_t3(i) + xyzbuffer19(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t5(i) = fsd9_avg_t5(i) + xyzbuffer32(i,j,k) * xyzbuffer5(i,j,k)
+xyzbuffer13(i,j,k) = ( ( xyzbuffer1(i,j,k) + xyzbuffer18(i,j,k) ) * xyzbuffer14(i,j,k) )
 end do
 end do
 end do
@@ -3076,127 +2321,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer22(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd1_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (nx(i,j,k) - fsd5_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer23(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd5_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (nx(i,j,k) - fsd9_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer20(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd9_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (divn(i,j,k) - fsd4_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer21(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd4_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (divn(i,j,k) - fsd8_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer27(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd8_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (absk(i,j,k) - fsd3_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd3_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (absk(i,j,k) - fsd7_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer29(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd7_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd2_avg_t1(i)))*((t1(i,j,k) - fsd2_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer12(i,j,k) = (((xyzbuffer35(i,j,k) - fsd2_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd2_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd6_avg_t1(i)))*((t1(i,j,k) - fsd6_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer19(i,j,k) = (((xyzbuffer35(i,j,k) - fsd6_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd6_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd2_avg_divn(i)))*((divn(i,j,k) - fsd2_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer10(i,j,k) = (((xyzbuffer25(i,j,k) - fsd2_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd2_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd1_avg_absk(i)))*((absk(i,j,k) - fsd1_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer13(i,j,k) = (((xyzbuffer18(i,j,k) - fsd1_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd1_avg_absk(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd5_avg_absk(i)))*((absk(i,j,k) - fsd5_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer34(i,j,k) = (((xyzbuffer18(i,j,k) - fsd5_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd5_avg_absk(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd9_avg_absk(i)))*((absk(i,j,k) - fsd9_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer11(i,j,k) = (((xyzbuffer18(i,j,k) - fsd9_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd9_avg_absk(i))))
+xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd1_avg_nx(i)) )
 end do
 end do
 end do
@@ -3205,7 +2330,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_t2(i) = fsd1_avg_t2(i) + xyzbuffer22(i,j,k) * xyzbuffer7(i,j,k)
+fsd1_avg_t2(i) = fsd1_avg_t2(i) + xyzbuffer12(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (nx(i,j,k) - fsd2_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd2_avg_nx(i)) )
 end do
 end do
 end do
@@ -3214,211 +2349,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_t8(i) = fsd1_avg_t8(i) + xyzbuffer13(i,j,k) * xyzbuffer7(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t5(i) = fsd2_avg_t5(i) + xyzbuffer12(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd3_avg_t4(i) = fsd3_avg_t4(i) + xyzbuffer24(i,j,k) * xyzbuffer26(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t3(i) = fsd4_avg_t3(i) + xyzbuffer21(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t8(i) = fsd4_avg_t8(i) + xyzbuffer1(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t8(i) = fsd5_avg_t8(i) + xyzbuffer34(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t5(i) = fsd6_avg_t5(i) + xyzbuffer19(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t4(i) = fsd7_avg_t4(i) + xyzbuffer29(i,j,k) * xyzbuffer28(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t3(i) = fsd8_avg_t3(i) + xyzbuffer27(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd8_avg_t8(i) = fsd8_avg_t8(i) + xyzbuffer4(i,j,k) * xyzbuffer9(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t8(i) = fsd9_avg_t8(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k)
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (nx(i,j,k) - fsd7_avg_nx(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer22(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd7_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (divn(i,j,k) - fsd6_avg_divn(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer21(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd6_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (absk(i,j,k) - fsd5_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer27(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd5_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd4_avg_t1(i)))*((t1(i,j,k) - fsd4_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer24(i,j,k) = (((xyzbuffer35(i,j,k) - fsd4_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd4_avg_t1(i))))
-end do
-end do
-end do
-
-
-! (((divn(i,j,k) - fsd6_avg_divn(i)))*((divn(i,j,k) - fsd6_avg_divn(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer29(i,j,k) = (((xyzbuffer25(i,j,k) - fsd6_avg_divn(i)))*((xyzbuffer25(i,j,k) - fsd6_avg_divn(i))))
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd7_avg_absk(i)))*((absk(i,j,k) - fsd7_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer34(i,j,k) = (((xyzbuffer18(i,j,k) - fsd7_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd7_avg_absk(i))))
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd2_avg_t7(i) = fsd2_avg_t7(i) + xyzbuffer10(i,j,k) * xyzbuffer16(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd4_avg_t5(i) = fsd4_avg_t5(i) + xyzbuffer24(i,j,k) * xyzbuffer15(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd5_avg_t4(i) = fsd5_avg_t4(i) + xyzbuffer27(i,j,k) * xyzbuffer6(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd6_avg_t7(i) = fsd6_avg_t7(i) + xyzbuffer29(i,j,k) * xyzbuffer14(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd7_avg_t8(i) = fsd7_avg_t8(i) + xyzbuffer34(i,j,k) * xyzbuffer28(i,j,k)
-end do
-end do
-end do
-
-
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-fsd9_avg_t2(i) = fsd9_avg_t2(i) + xyzbuffer20(i,j,k) * xyzbuffer5(i,j,k)
+fsd2_avg_t2(i) = fsd2_avg_t2(i) + xyzbuffer11(i,j,k) * xyzbuffer17(i,j,k)
 end do
 end do
 end do
@@ -3428,27 +2359,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer19(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer31(i,j,k) - fsd3_avg_nx(i)) )
-end do
-end do
-end do
-
-
-! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (absk(i,j,k) - fsd1_avg_absk(i)) )
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer10(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd1_avg_absk(i)) )
-end do
-end do
-end do
-
-
-! (((t1(i,j,k) - fsd8_avg_t1(i)))*((t1(i,j,k) - fsd8_avg_t1(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer13(i,j,k) = (((xyzbuffer35(i,j,k) - fsd8_avg_t1(i)))*((xyzbuffer35(i,j,k) - fsd8_avg_t1(i))))
+xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd3_avg_nx(i)) )
 end do
 end do
 end do
@@ -3457,7 +2368,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd1_avg_t4(i) = fsd1_avg_t4(i) + xyzbuffer10(i,j,k) * xyzbuffer7(i,j,k)
+fsd3_avg_t2(i) = fsd3_avg_t2(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (nx(i,j,k) - fsd4_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd4_avg_nx(i)) )
 end do
 end do
 end do
@@ -3466,7 +2387,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd3_avg_t2(i) = fsd3_avg_t2(i) + xyzbuffer19(i,j,k) * xyzbuffer26(i,j,k)
+fsd4_avg_t2(i) = fsd4_avg_t2(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (nx(i,j,k) - fsd5_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd5_avg_nx(i)) )
 end do
 end do
 end do
@@ -3475,7 +2406,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd5_avg_t2(i) = fsd5_avg_t2(i) + xyzbuffer23(i,j,k) * xyzbuffer6(i,j,k)
+fsd5_avg_t2(i) = fsd5_avg_t2(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (nx(i,j,k) - fsd6_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd6_avg_nx(i)) )
 end do
 end do
 end do
@@ -3484,7 +2425,83 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd7_avg_t2(i) = fsd7_avg_t2(i) + xyzbuffer22(i,j,k) * xyzbuffer28(i,j,k)
+fsd6_avg_t2(i) = fsd6_avg_t2(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (nx(i,j,k) - fsd7_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd7_avg_nx(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t2(i) = fsd7_avg_t2(i) + xyzbuffer12(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (nx(i,j,k) - fsd8_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd8_avg_nx(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t2(i) = fsd8_avg_t2(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (nx(i,j,k) - fsd9_avg_nx(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd9_avg_nx(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t2(i) = fsd9_avg_t2(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (divn(i,j,k) - fsd1_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd1_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t3(i) = fsd1_avg_t3(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k)
 end do
 end do
 end do
@@ -3494,17 +2511,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer22(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer25(i,j,k) - fsd2_avg_divn(i)) )
-end do
-end do
-end do
-
-
-! (((absk(i,j,k) - fsd3_avg_absk(i)))*((absk(i,j,k) - fsd3_avg_absk(i))))
-do k = 1, nzp
-do j = 1, nyp
-do i = 1, nxp
-xyzbuffer23(i,j,k) = (((xyzbuffer18(i,j,k) - fsd3_avg_absk(i)))*((xyzbuffer18(i,j,k) - fsd3_avg_absk(i))))
+xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd2_avg_divn(i)) )
 end do
 end do
 end do
@@ -3513,7 +2520,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd2_avg_t3(i) = fsd2_avg_t3(i) + xyzbuffer22(i,j,k) * xyzbuffer16(i,j,k)
+fsd2_avg_t3(i) = fsd2_avg_t3(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (divn(i,j,k) - fsd3_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd3_avg_divn(i)) )
 end do
 end do
 end do
@@ -3522,7 +2539,273 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd6_avg_t3(i) = fsd6_avg_t3(i) + xyzbuffer21(i,j,k) * xyzbuffer14(i,j,k)
+fsd3_avg_t3(i) = fsd3_avg_t3(i) + xyzbuffer18(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (divn(i,j,k) - fsd4_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd4_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t3(i) = fsd4_avg_t3(i) + xyzbuffer12(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (divn(i,j,k) - fsd5_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd5_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t3(i) = fsd5_avg_t3(i) + xyzbuffer11(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (divn(i,j,k) - fsd6_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd6_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t3(i) = fsd6_avg_t3(i) + xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (divn(i,j,k) - fsd7_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd7_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t3(i) = fsd7_avg_t3(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (divn(i,j,k) - fsd8_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd8_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t3(i) = fsd8_avg_t3(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (divn(i,j,k) - fsd9_avg_divn(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd9_avg_divn(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t3(i) = fsd9_avg_t3(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (absk(i,j,k) - fsd1_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd1_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t4(i) = fsd1_avg_t4(i) + xyzbuffer12(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (absk(i,j,k) - fsd2_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd2_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_t4(i) = fsd2_avg_t4(i) + xyzbuffer11(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (absk(i,j,k) - fsd3_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd3_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t4(i) = fsd3_avg_t4(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (absk(i,j,k) - fsd4_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd4_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t4(i) = fsd4_avg_t4(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (absk(i,j,k) - fsd5_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd5_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t4(i) = fsd5_avg_t4(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (absk(i,j,k) - fsd6_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd6_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t4(i) = fsd6_avg_t4(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (absk(i,j,k) - fsd7_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd7_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t4(i) = fsd7_avg_t4(i) + xyzbuffer12(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (absk(i,j,k) - fsd8_avg_absk(i)) )
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd8_avg_absk(i)) )
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t4(i) = fsd8_avg_t4(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k)
 end do
 end do
 end do
@@ -3532,7 +2815,7 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-xyzbuffer22(i,j,k) = ( (xyzbuffer35(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer18(i,j,k) - fsd9_avg_absk(i)) )
+xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd9_avg_absk(i)) )
 end do
 end do
 end do
@@ -3541,7 +2824,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd3_avg_t8(i) = fsd3_avg_t8(i) + xyzbuffer23(i,j,k) * xyzbuffer26(i,j,k)
+fsd9_avg_t4(i) = fsd9_avg_t4(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd1_avg_t1(i)))*((t1(i,j,k) - fsd1_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = (((xyzbuffer13(i,j,k) - fsd1_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd1_avg_t1(i))))
 end do
 end do
 end do
@@ -3550,7 +2843,17 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd9_avg_t4(i) = fsd9_avg_t4(i) + xyzbuffer22(i,j,k) * xyzbuffer5(i,j,k)
+fsd1_avg_t5(i) = fsd1_avg_t5(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd2_avg_t1(i)))*((t1(i,j,k) - fsd2_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = (((xyzbuffer13(i,j,k) - fsd2_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd2_avg_t1(i))))
 end do
 end do
 end do
@@ -3559,7 +2862,653 @@ end do
 do k = 1, nzp
 do j = 1, nyp
 do i = 1, nxp
-fsd8_avg_t5(i) = fsd8_avg_t5(i) + xyzbuffer13(i,j,k) * xyzbuffer9(i,j,k)
+fsd2_avg_t5(i) = fsd2_avg_t5(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd3_avg_t1(i)))*((t1(i,j,k) - fsd3_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = (((xyzbuffer13(i,j,k) - fsd3_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd3_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t5(i) = fsd3_avg_t5(i) + xyzbuffer18(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd4_avg_t1(i)))*((t1(i,j,k) - fsd4_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = (((xyzbuffer13(i,j,k) - fsd4_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd4_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t5(i) = fsd4_avg_t5(i) + xyzbuffer12(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd5_avg_t1(i)))*((t1(i,j,k) - fsd5_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = (((xyzbuffer13(i,j,k) - fsd5_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd5_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t5(i) = fsd5_avg_t5(i) + xyzbuffer11(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd6_avg_t1(i)))*((t1(i,j,k) - fsd6_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = (((xyzbuffer13(i,j,k) - fsd6_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd6_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t5(i) = fsd6_avg_t5(i) + xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd7_avg_t1(i)))*((t1(i,j,k) - fsd7_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = (((xyzbuffer13(i,j,k) - fsd7_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd7_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t5(i) = fsd7_avg_t5(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd8_avg_t1(i)))*((t1(i,j,k) - fsd8_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = (((xyzbuffer13(i,j,k) - fsd8_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd8_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t5(i) = fsd8_avg_t5(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! (((t1(i,j,k) - fsd9_avg_t1(i)))*((t1(i,j,k) - fsd9_avg_t1(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = (((xyzbuffer13(i,j,k) - fsd9_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd9_avg_t1(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t5(i) = fsd9_avg_t5(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd1_avg_nx(i)))*((nx(i,j,k) - fsd1_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer13(i,j,k) = (((xyzbuffer14(i,j,k) - fsd1_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd1_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t6(i) = fsd1_avg_t6(i) + xyzbuffer13(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd2_avg_nx(i)))*((nx(i,j,k) - fsd2_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = (((xyzbuffer14(i,j,k) - fsd2_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd2_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_t6(i) = fsd2_avg_t6(i) + xyzbuffer12(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd3_avg_nx(i)))*((nx(i,j,k) - fsd3_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = (((xyzbuffer14(i,j,k) - fsd3_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd3_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t6(i) = fsd3_avg_t6(i) + xyzbuffer11(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd4_avg_nx(i)))*((nx(i,j,k) - fsd4_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = (((xyzbuffer14(i,j,k) - fsd4_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd4_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t6(i) = fsd4_avg_t6(i) + xyzbuffer16(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd5_avg_nx(i)))*((nx(i,j,k) - fsd5_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = (((xyzbuffer14(i,j,k) - fsd5_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd5_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t6(i) = fsd5_avg_t6(i) + xyzbuffer15(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd6_avg_nx(i)))*((nx(i,j,k) - fsd6_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = (((xyzbuffer14(i,j,k) - fsd6_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd6_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t6(i) = fsd6_avg_t6(i) + xyzbuffer1(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd7_avg_nx(i)))*((nx(i,j,k) - fsd7_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = (((xyzbuffer14(i,j,k) - fsd7_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd7_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t6(i) = fsd7_avg_t6(i) + xyzbuffer18(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd8_avg_nx(i)))*((nx(i,j,k) - fsd8_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer13(i,j,k) = (((xyzbuffer14(i,j,k) - fsd8_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd8_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t6(i) = fsd8_avg_t6(i) + xyzbuffer13(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! (((nx(i,j,k) - fsd9_avg_nx(i)))*((nx(i,j,k) - fsd9_avg_nx(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = (((xyzbuffer14(i,j,k) - fsd9_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd9_avg_nx(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t6(i) = fsd9_avg_t6(i) + xyzbuffer12(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd1_avg_divn(i)))*((divn(i,j,k) - fsd1_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = (((xyzbuffer7(i,j,k) - fsd1_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd1_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t7(i) = fsd1_avg_t7(i) + xyzbuffer11(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd2_avg_divn(i)))*((divn(i,j,k) - fsd2_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = (((xyzbuffer7(i,j,k) - fsd2_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd2_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_t7(i) = fsd2_avg_t7(i) + xyzbuffer16(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd3_avg_divn(i)))*((divn(i,j,k) - fsd3_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = (((xyzbuffer7(i,j,k) - fsd3_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd3_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t7(i) = fsd3_avg_t7(i) + xyzbuffer15(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd4_avg_divn(i)))*((divn(i,j,k) - fsd4_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer14(i,j,k) = (((xyzbuffer7(i,j,k) - fsd4_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd4_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t7(i) = fsd4_avg_t7(i) + xyzbuffer14(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd5_avg_divn(i)))*((divn(i,j,k) - fsd5_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = (((xyzbuffer7(i,j,k) - fsd5_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd5_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t7(i) = fsd5_avg_t7(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd6_avg_divn(i)))*((divn(i,j,k) - fsd6_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer18(i,j,k) = (((xyzbuffer7(i,j,k) - fsd6_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd6_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t7(i) = fsd6_avg_t7(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd7_avg_divn(i)))*((divn(i,j,k) - fsd7_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer13(i,j,k) = (((xyzbuffer7(i,j,k) - fsd7_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd7_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t7(i) = fsd7_avg_t7(i) + xyzbuffer13(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd8_avg_divn(i)))*((divn(i,j,k) - fsd8_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer12(i,j,k) = (((xyzbuffer7(i,j,k) - fsd8_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd8_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t7(i) = fsd8_avg_t7(i) + xyzbuffer12(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! (((divn(i,j,k) - fsd9_avg_divn(i)))*((divn(i,j,k) - fsd9_avg_divn(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer11(i,j,k) = (((xyzbuffer7(i,j,k) - fsd9_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd9_avg_divn(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t7(i) = fsd9_avg_t7(i) + xyzbuffer11(i,j,k) * xyzbuffer2(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd1_avg_absk(i)))*((absk(i,j,k) - fsd1_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer16(i,j,k) = (((xyzbuffer6(i,j,k) - fsd1_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd1_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd1_avg_t8(i) = fsd1_avg_t8(i) + xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd2_avg_absk(i)))*((absk(i,j,k) - fsd2_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer15(i,j,k) = (((xyzbuffer6(i,j,k) - fsd2_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd2_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd2_avg_t8(i) = fsd2_avg_t8(i) + xyzbuffer15(i,j,k) * xyzbuffer17(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd3_avg_absk(i)))*((absk(i,j,k) - fsd3_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer14(i,j,k) = (((xyzbuffer6(i,j,k) - fsd3_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd3_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd3_avg_t8(i) = fsd3_avg_t8(i) + xyzbuffer14(i,j,k) * xyzbuffer9(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd4_avg_absk(i)))*((absk(i,j,k) - fsd4_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer9(i,j,k) = (((xyzbuffer6(i,j,k) - fsd4_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd4_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd4_avg_t8(i) = fsd4_avg_t8(i) + xyzbuffer9(i,j,k) * xyzbuffer3(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd5_avg_absk(i)))*((absk(i,j,k) - fsd5_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer3(i,j,k) = (((xyzbuffer6(i,j,k) - fsd5_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd5_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd5_avg_t8(i) = fsd5_avg_t8(i) + xyzbuffer3(i,j,k) * xyzbuffer0(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd6_avg_absk(i)))*((absk(i,j,k) - fsd6_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer1(i,j,k) = (((xyzbuffer6(i,j,k) - fsd6_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd6_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd6_avg_t8(i) = fsd6_avg_t8(i) + xyzbuffer1(i,j,k) * xyzbuffer10(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd7_avg_absk(i)))*((absk(i,j,k) - fsd7_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer0(i,j,k) = (((xyzbuffer6(i,j,k) - fsd7_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd7_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd7_avg_t8(i) = fsd7_avg_t8(i) + xyzbuffer0(i,j,k) * xyzbuffer8(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd8_avg_absk(i)))*((absk(i,j,k) - fsd8_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer7(i,j,k) = (((xyzbuffer6(i,j,k) - fsd8_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd8_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd8_avg_t8(i) = fsd8_avg_t8(i) + xyzbuffer7(i,j,k) * xyzbuffer5(i,j,k)
+end do
+end do
+end do
+
+
+! (((absk(i,j,k) - fsd9_avg_absk(i)))*((absk(i,j,k) - fsd9_avg_absk(i))))
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+xyzbuffer5(i,j,k) = (((xyzbuffer6(i,j,k) - fsd9_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd9_avg_absk(i))))
+end do
+end do
+end do
+
+
+do k = 1, nzp
+do j = 1, nyp
+do i = 1, nxp
+fsd9_avg_t8(i) = fsd9_avg_t8(i) + xyzbuffer5(i,j,k) * xyzbuffer2(i,j,k)
 end do
 end do
 end do