refl/gray/after_onestep: fixed test for search of ray that completed 1st pass in plasma

2019-03-14 10:53:34 +00:00 · 2019-03-14 10:53:34 +00:00 · 56b8700af3
commit 56b8700af3
parent 49303ac2e3
1 changed files with 1 additions and 1 deletions
--- a/src/gray.f
+++ b/src/gray.f
@ -468,7 +468,7 @@ c     if found, use its polarization state to compute mode coupling
                end do powloop
 c     if no ray completed a first pass in the plasma, use central ray
 c     initial polarization (possibly reflected)
-                if (qqout.le.0.0_wp_) then
+                if (qqout**2+uuout**2+vvout**2.le.0.0_wp_) then
                  call stokes(ext(1,1,0),eyt(1,1,0),qqout,uuout,vvout)
                end if
                powrfl=0.5_wp_*(1.0_wp_+vvout*vvin2+