Skip to content
Snippets Groups Projects
Commit 846dc6e4 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

Correction on Z-pinch geometry

parent 665d0ec3
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 2 deletions
src/geometry_mod.F90
+ 2
− 2
View file @ 846dc6e4
...@@ -258,7 +258,7 @@ CONTAINS ...@@ -258,7 +258,7 @@ CONTAINS
hatB_NL(iz,eo) = 1._dp hatB_NL(iz,eo) = 1._dp
! Derivative of the magnetic field strenght ! Derivative of the magnetic field strenght
gradxB(iz,eo) = 0._dp ! Gene put a factor hatB^2 or 1/hatR^2 in this gradxB(iz,eo) = -1._dp ! Gene put a factor hatB^2 or 1/hatR^2 in this
gradyB(iz,eo) = 0._dp gradyB(iz,eo) = 0._dp
gradzB(iz,eo) = 0._dp ! Gene put a factor hatB or 1/hatR in this gradzB(iz,eo) = 0._dp ! Gene put a factor hatB or 1/hatR in this
...@@ -267,7 +267,7 @@ CONTAINS ...@@ -267,7 +267,7 @@ CONTAINS
ky = kyarray(iky) ky = kyarray(iky)
DO ikx= ikxs, ikxe DO ikx= ikxs, ikxe
kx = kxarray(ikx) kx = kxarray(ikx)
Ckxky(iky, ikx, iz,eo) = -ky * hatB(iz,eo) ! .. multiply by hatB to cancel the 1/ hatB factor in moments_eqs_rhs.f90 routine Ckxky(iky, ikx, iz,eo) = -ky
ENDDO ENDDO
ENDDO ENDDO
! coefficient in the front of parallel derivative ! coefficient in the front of parallel derivative
......
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