Skip to content

Commit

Permalink
fix a bug where i\omega\theta was missing from the A matrix in forced…
Browse files Browse the repository at this point in the history
… problems
  • Loading branch information
jrekier committed Jun 20, 2024
1 parent e7738b7 commit 007fb2d
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 2 deletions.
4 changes: 3 additions & 1 deletion bin/assemble.py
Original file line number Diff line number Diff line change
Expand Up @@ -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

Expand Down
2 changes: 1 addition & 1 deletion bin/operators_new.py
Original file line number Diff line number Diff line change
Expand Up @@ -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



Expand Down

0 comments on commit 007fb2d

Please sign in to comment.