Skip to content
Permalink
Browse files

improve (or at least make more robust) characteristic speeds

  • Loading branch information...
evanoconnor committed Feb 15, 2019
1 parent 1677b0b commit 7ca205e64f8954359466e8804477979bd9e41637
Showing with 8 additions and 8 deletions.
  1. +8 −8 src/M1/M1_explicitterms.F90
@@ -211,11 +211,11 @@ subroutine M1_explicitterms(dts,implicit_factor)
endif

!actual speed
l_min(1) = (3.0d0*M1chi_space_minus(k+1)-1.0d0)*0.5d0*l_min_thin(1) + &
(1.0d0 - M1chi_space_minus(k+1))*1.5d0*l_min_thick(1)
l_min(1) = (3.0d0*M1chi_space(k+1)-1.0d0)*0.5d0*l_min_thin(1) + &
(1.0d0 - M1chi_space(k+1))*1.5d0*l_min_thick(1)

l_max(1) = (3.0d0*M1chi_space_minus(k+1)-1.0d0)*0.5d0*l_max_thin(1) + &
(1.0d0 - M1chi_space_minus(k+1))*1.5d0*l_max_thick(1)
l_max(1) = (3.0d0*M1chi_space(k+1)-1.0d0)*0.5d0*l_max_thin(1) + &
(1.0d0 - M1chi_space(k+1))*1.5d0*l_max_thick(1)

!plus interface (k zone)
!thin limit:
@@ -274,10 +274,10 @@ subroutine M1_explicitterms(dts,implicit_factor)
endif

!actualy speed
l_min(2) = (3.0d0*M1chi_space_plus(k)-1.0d0)*0.5d0*l_min_thin(2) + &
(1.0d0 - M1chi_space_plus(k))*1.5d0*l_min_thick(2)
l_max(2) = (3.0d0*M1chi_space_plus(k)-1.0d0)*0.5d0*l_max_thin(2) + &
(1.0d0 - M1chi_space_plus(k))*1.5d0*l_max_thick(2)
l_min(2) = (3.0d0*M1chi_space(k)-1.0d0)*0.5d0*l_min_thin(2) + &
(1.0d0 - M1chi_space(k))*1.5d0*l_min_thick(2)
l_max(2) = (3.0d0*M1chi_space(k)-1.0d0)*0.5d0*l_max_thin(2) + &
(1.0d0 - M1chi_space(k))*1.5d0*l_max_thick(2)


!check for NaNs

0 comments on commit 7ca205e

Please sign in to comment.
You can’t perform that action at this time.