From 56b8700af3994b0871a6a31b628685bad54f9e03 Mon Sep 17 00:00:00 2001
From: Daniele Micheletti <daniele.micheletti@istp.cnr.it>
Date: Thu, 14 Mar 2019 10:53:34 +0000
Subject: [PATCH] refl/gray/after_onestep: fixed test for search of ray that
 completed 1st pass in plasma

---
 src/gray.f | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/gray.f b/src/gray.f
index 637afd1..0b0302a 100644
--- 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+