no stop

src/epos/epos-fra.f

@@ -358,7 +358,12 @@ c          write(ifch,*)"nmax",nmax
           if(i.eq.nk1.and.nn.eq. 1  ) f=1.
           if(i.eq.nk2.and.nn.eq.nmax) f=1.
           nob=nob+1
-          if(nob.gt.mxnba-7)stop'gakfra: increase parameter mxnba '
+          if(nob.gt.mxnba-7) then !stop'gakfra: increase parameter mxnba '
+              write(*,*)'gaksbp: increase parameter mxnba'
+              write(*,*)'HACK! Original code stops here'
+              iok=1
+              return
+          endif
 c         if(nmax.ge.2) print *, nob,ff,f,i
 c         if(nmax.ge.2) write(ifch,*)'nob',nob,ff,f,i
           do k=1,4
@@ -371,14 +376,24 @@ c         if(nmax.ge.2) write(ifch,*)'nob',nob,ff,f,i
 
       do i1=nob-1,1,-1          ! copy again gluons
         nob=nob+1
-        if(nob.gt.mxnba-7)stop'gakfra: increase parameter mxnba '
+        if(nob.gt.mxnba-7) then !stop'gakfra: increase parameter mxnba '
+            write(*,*)'gaksbp: increase parameter mxnba'
+            write(*,*)'HACK! Original code stops here'
+            iok=1
+            return
+        endif
         do k=1,4
           pst(k,nob)=pst(k,i1)
           ipst(nob)=ipst(i1)
         enddo
       enddo
       nob=nob+1                 ! nob is number of kinks
-      if(nob.gt.mxnba-7)stop'gakfra: increase parameter mxnba '
+      if(nob.gt.mxnba-7) then !stop'gakfra: increase parameter mxnba '
+          write(*,*)'gaksbp: increase parameter mxnba'
+          write(*,*)'HACK! Original code stops here'
+          iok=1
+          return
+      endif
       if(ish.ge.6) then
         write (ifch,*) 'bands:'
         do i=0,nob-1
@@ -1161,7 +1176,12 @@ c......................................................................
 
 c.....breakpoint accepted......................
       nbr=nbr+1
-      if(nbr.gt.mxnbr-2) stop 'gaksbp: increase parameter mxnbr'
+      if(nbr.gt.mxnbr-2) then !stop 'gaksbp: increase parameter mxnbr'
+          write(*,*)'gaksbp: increase parameter mxnbr'
+          write(*,*)'HACK! Original code stops here'
+          iok=1
+          return
+      endif
       do i=nbr+1,ib1+1,-1
          do j=1,2
             ijb(j,i)=ijb(j,i-1)
@@ -1958,7 +1978,12 @@ c     &     -pst(2,i)*pst(2,j)-pst(1,i)*pst(1,j)
         endif
       enddo
       ni=ni+1
-      if(ni.gt.mxnin)stop'gakanp: increase parameter mxnin  '
+      if(ni.gt.mxnin) then !stop'gakanp: increase parameter mxnin  '
+          write(*,*)'gaksbp: increase parameter mxnin'
+          write(*,*)'HACK! Original code stops here'
+          iok=1
+          return
+      endif
       nin(1,ni)=ib
       nin(2,ni)=jb
       if(ish.ge.6)write(ifch,*)

src/epos/epos-ids.f

@@ -1938,7 +1938,7 @@ c wrapper for idtrafo to catch if ID not found
       idtrafo = idtrafostatus(code1,code2,idi,istatus)
       if (istatus.ne.0) then
          print *,'idtrafo: ',code1,' -> ', code2,idi,' not found.'
-         stop
+         ! stop
       endif
       return
       end

src/epos/epos-sem.f

@@ -5639,6 +5639,7 @@ c------------------------------------------------------------------------
       double precision amqt(4),yqm(4),yqm1(4),xlp(4),xlm(4),am23,sx,y2
      *,wp1,wp2,wm1,wm2,s,sq,psutz,yqmax,y,amjp,amjm,y1,s12,s34,x34,amqpt
       dimension samqt(4)
+      integer iteration
       include 'epos.inc'
 
       s=wp1*wm1
@@ -5649,7 +5650,14 @@ c------------------------------------------------------------------------
       enddo
       yqmax=max(yqm(1),yqm(2))
 
-1     y=yqmax*dble(rangen())
+      iteration=0
+1     iteration=iteration+1
+      if(iteration.gt.1000) then
+        write(*,*)'HACK! too many iterations in pslcsh, breaking out'
+        write(*,*)'HACK! this should not happen!'
+        goto 2
+      endif
+      y=yqmax*dble(rangen())
       j=int(1.5+rangen())
       if(y.gt.yqm(j))goto 1
 
@@ -5684,7 +5692,7 @@ c      write (*,*)'xlc',xlp(1),xlp(2),xlm(3),xlm(4)
       if(dble(rangen()).gt.(xlp(1)*xlp(2)*xlm(3)*xlm(4))**(-alpqua)*
      *(xlp(j)*(1.d0-xlp(j))*xlm(j1)*(1.d0-xlm(j1))))goto 1
 
-      wp2=xlp(2)*wp1
+2     wp2=xlp(2)*wp1
       wp1=xlp(1)*wp1
       wm2=xlm(4)*wm1
       wm1=xlm(3)*wm1

src/epos/epos-uti.f

@@ -3268,11 +3268,15 @@ c      parameter(itext=40)
      */1x,'***** current event number: ',i12
      */1x,'***** initial seed for current run:',i10
      */1x,'***** initial seed for current event:',d25.15
+     */1x,'HACK! Original code stops here.',
+     */1x,'HACK! MC should not stop at anomalous samples.',
+     */1x,'HACK! Be careful however.',
      */1x,72a1)
 1     continue
 c      c=0.
 c      b=a/c
-      stop
+      return ! modified
+      ! stop ! modified
       entry utmsg(txt)
       imsg=imsg+1
       write(ifch,'(1x,74a1)')('*',j=1,72)