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

Hyp coll coeff are now 1/(P+1)

to avoid /0 when P=0
parent 5c036d37
Branches
Tags
No related merge requests found
......@@ -363,7 +363,7 @@ CONTAINS
END DO
! Precomputations
pmax_xp = real(pmax,xp)
diff_p_coeff = pmax_xp*(1._xp/pmax_xp)**6
diff_p_coeff = pmax_xp*(1._xp/(pmax_xp+1._xp))**6
! Overwrite SOLVE_AMPERE flag if beta is zero
IF(.NOT. EM) THEN
SOLVE_AMPERE = .FALSE.
......@@ -386,7 +386,7 @@ CONTAINS
local_jmin = jarray(1+ngj/2)
! Precomputations
jmax_xp = real(jmax,xp)
diff_j_coeff = jmax_xp*(1._xp/jmax_xp)**6
diff_j_coeff = jmax_xp*(1._xp/(jmax_xp+1._xp))**6
! j=0 and j=1 indices
DO ij = 1,local_nj+ngj
IF(jarray(ij) .EQ. 0) ij0 = ij
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment