fixed check for N// > 1

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
 

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)