fix a bug where i\omega\theta was missing from the A matrix in forced…

… problems
repepo · Jun 20, 2024 · 007fb2d · 007fb2d
1 parent e7738b7
commit 007fb2d
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/bin/assemble.py b/bin/assemble.py
@@ -1042,7 +1042,9 @@ def main():
             basecol = (2*par.hydro+2*par.magnetic)*nb*ut.N1
 
             # Physics ----------------------------
-            mtx = op.thermal_diffusion(l,'h','',0)
+            iwtheta = 1j*ut.wf * op.theta(l,'h','', 0)
+            diffus = op.thermal_diffusion(l,'h','',0)
+            mtx = diffus - iwtheta
             # ------------------------------------
             col = basecol + col0
 

diff --git a/bin/operators_new.py b/bin/operators_new.py
@@ -754,7 +754,7 @@ def thermal_diffusion(l, section, component, offdiag):
 
             difus = - L*r0_kho0_D0_h + 2*r1_kho0_D1_h + r2_kho0_D2_h + r2_kho0_lnT1_D1_h + r2_kho1_D1_h
 
-    return difus * par.OmgTau * par.Ek / par.Prandtl
+    return difus * par.OmgTau * par.Etherm