Skip to content
Snippets Groups Projects
Commit 5c727b90 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann
Browse files

Debugging 1D linear

parent c3598dd5
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 4 additions and 0 deletions
src/moments_eq_rhs.F90
+ 4
− 0
View file @ 5c727b90
...@@ -127,6 +127,8 @@ SUBROUTINE moments_eq_rhs ...@@ -127,6 +127,8 @@ SUBROUTINE moments_eq_rhs
kzloope : DO ikz = ikzs,ikze kzloope : DO ikz = ikzs,ikze
kr = krarray(ikr) ! Poloidal wavevector kr = krarray(ikr) ! Poloidal wavevector
kz = kzarray(ikz) ! Toroidal wavevector kz = kzarray(ikz) ! Toroidal wavevector
IF (Nkz .EQ. 1) kz = krarray(ikr) ! If 1D simulation we put kr as kz
kperp2 = kr**2 + kz**2 ! perpendicular wavevector kperp2 = kr**2 + kz**2 ! perpendicular wavevector
be_2 = kperp2 * sigmae2_taue_o2 ! Kernel argument be_2 = kperp2 * sigmae2_taue_o2 ! Kernel argument
...@@ -338,6 +340,8 @@ SUBROUTINE moments_eq_rhs ...@@ -338,6 +340,8 @@ SUBROUTINE moments_eq_rhs
kzloopi : DO ikz = ikzs,ikze kzloopi : DO ikz = ikzs,ikze
kr = krarray(ikr) ! Poloidal wavevector kr = krarray(ikr) ! Poloidal wavevector
kz = kzarray(ikz) ! Toroidal wavevector kz = kzarray(ikz) ! Toroidal wavevector
IF (Nkz .EQ. 1) kz = krarray(ikr) ! If 1D simulation we put kr as kz
kperp2 = kr**2 + kz**2 ! perpendicular wavevector kperp2 = kr**2 + kz**2 ! perpendicular wavevector
bi_2 = kperp2 * sigmai2_taui_o2 ! Kernel argument bi_2 = kperp2 * sigmai2_taui_o2 ! Kernel argument
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment