🚀 go-pugleaf

The following short program compiled with MS-FORTRAN 5.0 with 
coprocessor crashes OS2 Ver. 2 on my machine:

c DBZ.FOR - Division by zero test
      double precision x, y
      x = 1.0d0
      y = 0.0d0
      call devide( x, y )
      end

      subroutine devide( x, y )
      double precision x, y
      double precision z
      z = x / y
      return
      end

However, on OS2 1.3 and plain DOS the program exits gracefully with
the following message:

run-time error M6103: MATH
- floating-point error: divide by 0

Does anybody know a workaround other then masking the numerical 
cooprocessor ?

My machine is 386/33 with 387 copprocessor.

Thanks,