#include <f77/f77_floatingpoint.h>
       subroutine fpfix(iret)
       integer iret
       external fp_handler
       iret=ieee_handler('set','underflow',fp_handler)
       iret=iret+ieee_handler('set','inexact',fp_handler)
       return
       end
       integer function fp_handler(sig,code,sigcontxt)
       integer sig,code,sigcontxt(5)
       return
       end