From 72b682229c85c9be6a48a90de407b77047e326d6 Mon Sep 17 00:00:00 2001
From: Daniela Farina <daniela.farina@istp.cnr.it>
Date: Wed, 1 Feb 2017 11:06:09 +0000
Subject: [PATCH] fixed check for N// > 1

---
 src/errcodes.f90 | 6 +++---
 src/graycore.f90 | 1 +
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/src/errcodes.f90 b/src/errcodes.f90
index 280399d..d257188 100644
--- a/src/errcodes.f90
+++ b/src/errcodes.f90
@@ -17,8 +17,8 @@ contains
     if(ibits(ierr,pnpl, lnpl )==2  .or.    & ! N// too large
        ibits(ierr,palph,lalph)==1) then      ! alpha < 0
       istop = 1
-    else
-      istop = 0
+!    else
+!      istop = 0
     end if
   end subroutine check_err
 
@@ -77,4 +77,4 @@ contains
       print*,i,'    IERR = ', ierrs*2**pcdfc,'  fcur integration error (2nd trapping region)'
   end subroutine print_errhcd
 
-end module errcodes
\ No newline at end of file
+end module errcodes
diff --git a/src/graycore.f90 b/src/graycore.f90
index 4771e91..7fcaec1 100644
--- a/src/graycore.f90
+++ b/src/graycore.f90
@@ -164,6 +164,7 @@ contains
 
       allout = .true. ! becomes false if at least part of the beam is inside the plsama
       ierr = 0
+      istop = 0
       do jk=1,nray
 ! compute derivatives with updated gradient and local plasma values
         xv = yw(1:3,jk)