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

doubt on the NL geometric prefactor

parent 2c359b48
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 3 deletions
src/geometry_mod.F90
+ 2
− 3
View file @ f6c793cb
......@@ -139,9 +139,6 @@ CONTAINS
ENDDO
! coefficient in the front of parallel derivative
gradz_coeff(iz,eo) = 1._dp / jacobian(iz,eo) / hatB(iz,eo)
! Factor in front of the nonlinear term
hatB_NL(iz,eo) = Jacobian(iz,eo)&
*(gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)/hatB(iz,eo)
ENDDO
ENDDO
......@@ -213,6 +210,8 @@ CONTAINS
ENDDO
! coefficient in the front of parallel derivative
gradz_coeff(iz,eo) = 1._dp / Jacobian(iz,eo) / hatB(iz,eo)
! Nonlinear term prefactor
hatB_NL(iz,eo) = 1._dp !Jacobian(iz,eo)*(gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)/hatB(iz,eo)
ENDDO zloop
ENDDO parity
......
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