deletions | additions       diff --git a/begin_verbatim_real_dp_function__.tex b/begin_verbatim_real_dp_function__.tex  index 1d2bc61..1115ccb 100644  --- a/begin_verbatim_real_dp_function__.tex  +++ b/begin_verbatim_real_dp_function__.tex 
...
omega2 = pow2(1d-6*2*pi*s% nu_max*nu_factor)    ! k_inner and k_outer are bounds of evanescent region    ! k_outer is outermost k where Sl2 <= omega2 at k-1 and Sl2 > omega2 at k  ! 1st find outermost where Sl2 <= omega2  h1 = s% net_iso(ih1)  
...
n2 = s% brunt_N2(k)  sl2 = L2*cs2/r2  xh1 = s% xa(h1,k)  ! if (k == k_inner-1) then  ! ! interpolate to estimate where in cell k N2 == omega2  ! ! use that location in place of rmid(k)  ! if (n2 <= omega2 .and. omega2 <= s% brunt_N2(k+1)) then  ! r0_inner = pow_cr(&  ! find0(pow3(s% r(k)), n2-omega2, pow3(s% r(k+1)), s% brunt_N2(k+1)-omega2), 1d0/3d0)  ! if (r0_inner < s% r(k+1) .or. r0_inner > s% r(k)) then  ! if (dbg) write(*,1) 'r0_inner confusion', s% r(k+1), r0_inner, s% r(k)  ! end if  ! else  ! r0_inner = s% r(k+1)  ! if (dbg) write(*,2) 'N2 == omega2 failed for k_inner', k_inner  ! end if  ! dr = s% rmid(k-1) - r0_inner  ! if (k == k_outer) then ! find where in cell k-1 Sl2 == omega2  ! ! use that location in place of rmid(k-1)  ! sl2_next = L2*pow2(s% csound(k-1)/s% r(k-1))  ! if (sl2 >= omega2 .and. omega2 >= sl2_next) then  ! r0_outer = pow_cr(&  ! find0(pow3(s% r(k)), sl2-omega2, pow3(s% r(k-1)), sl2_next-omega2), 1d0/3d0)  ! if (r0_outer < s% r(k) .or. r0_outer > s% r(k-1)) then  ! if (dbg) write(*,1) 'r0_outer confusion', s% r(k), r0_outer, s% r(k-1)  ! end if  ! else  ! r0_outer = s% r(k-1)  ! if (dbg) write(*,2) 'Sl2 == omega2 failed for k_outer', k_outer  ! end if  ! dr = r0_outer - s% rmid(k)  ! else  ! dr = s% rmid(k-1) - s% rmid(k)  ! end if  dr = s% rmid(k-1) - s% rmid(k)   kr2 = (1 - n2/omega2)*(1 - Sl2/omega2)*omega2/cs2  if (kr2 < 0 .and. n2 < omega2 .and. omega2 < Sl2 .and. xh1 > 0.001) &