kolmo_law: good S2 equation and raise

Clovis Lambert · Clovis Lambert · commit 062a1e7d7c7f · 2026-05-13T10:09:12.000+02:00
diff --git a/fluidsim/base/output/kolmo_law.py b/fluidsim/base/output/kolmo_law.py
@@ -217,7 +217,7 @@ def compute(self):
 
         # Compute second-order structure function
         val = sum(fft_vi[i] * fft_vi[i].conj() for i in range(3))
-        S2_k_r = 2 * E_k_mean - 2 * self.sim.oper.ifft(val)
+        S2_k_r = 4 * E_k_mean - 2 * self.sim.oper.ifft(val)
 
         # If buoyancy field exists, compute J_p
         if "b" in keys_state_phys:
@@ -264,6 +264,8 @@ def compute(self):
             Jk_r_array[0], Jk_r_array[1], Jk_r_array[2]
         )
 
+        raise
+
         Jh_k, Jt_k, Jv_k = self.coord_conv.compute_cylindrical_components(
             Jk_r_array[0], Jk_r_array[1], Jk_r_array[2]
         )

Original file line number	Diff line number	Diff line change
`@@ -217,7 +217,7 @@ def compute(self):`
`217`	`217`
`218`	`218`	`# Compute second-order structure function`
`219`	`219`	`val = sum(fft_vi[i] * fft_vi[i].conj() for i in range(3))`
`220`		`- S2_k_r = 2 * E_k_mean - 2 * self.sim.oper.ifft(val)`
	`220`	`+ S2_k_r = 4 * E_k_mean - 2 * self.sim.oper.ifft(val)`
`221`	`221`
`222`	`222`	`# If buoyancy field exists, compute J_p`
`223`	`223`	`if "b" in keys_state_phys:`
`@@ -264,6 +264,8 @@ def compute(self):`
`264`	`264`	`Jk_r_array[0], Jk_r_array[1], Jk_r_array[2]`
`265`	`265`	`)`
`266`	`266`
	`267`	`+ raise`
	`268`	`+`
`267`	`269`	`Jh_k, Jt_k, Jv_k = self.coord_conv.compute_cylindrical_components(`
`268`	`270`	`Jk_r_array[0], Jk_r_array[1], Jk_r_array[2]`
`269`	`271`	`)`