        subroutine STMESS(IUNIT,NMESS,NEQN,TINIT,TFINAL,TINCR
     &, RELERR,ABSERR)
        double precision T, TOUT, RELER, ABSER, DTSGN
        double precision TINIT,TFINAL,TINCR,RELERR,ABSERR
        integer NEQ, IFLAG, IOLD, INIT, KOP
        integer IUNIT,NMESS,NEQN
        common/SFDATA/T,TOUT,RELER,ABSER,DTSGN,NEQ,IFLAG,IOLD,INIT,KOP
      if (NMESS.gt.50) go to 100
      if (NMESS.eq.10) go to 110
      if (NMESS.eq.20) go to 120
      if (NMESS.eq.30) go to 130
      if (NMESS.eq.40) go to 140
      if (NMESS.eq.50) go to 150
100   continue
      if (NMESS.eq.60) go to 160
      if (NMESS.eq.70) go to 170
      if (NMESS.eq.80) go to 180
      if (NMESS.eq.90) go to 190
      if (NMESS.eq.92) go to 200
      if (NMESS.eq.95) go to 210
      write (IUNIT,5000)
      go to 220
5000  format (1X,'Malfunction in Subroutine "STMESS."')
110   continue
      write (IUNIT,5010) NEQN,TINIT,TFINAL,TINCR,RELERR,ABSERR
      go to 220
120   continue
      write (IUNIT,5020) T,IFLAG
      go to 220
130   continue
      write (IUNIT,5030) T,RELER,ABSER
      go to 220
140   continue
      write (IUNIT,5040) T,RELER,ABSER
      go to 220
150   continue
      write (IUNIT,5050) T
      write (IUNIT,5060) RELER,ABSER
      go to 220
160   continue
      write (IUNIT,5070) T,RELER,ABSER
      go to 220
170   continue
      write (IUNIT,5080) T
      go to 220
180   continue
      write (IUNIT,5090) T
      go to 220
190   continue
      write (IUNIT,5100) T,NEQ,RELER,ABSER,TINCR
      go to 220
200   continue
      write (IUNIT,5110) T
      go to 220
210   continue
      write (IUNIT,5120) T
      go to 220
5010  format (//' STFINT has been called to integrate',I5,' equations'/'
     & from initial time',1PG12.5/' to final time',G12.5/' with output a
     &t steps of',G12.5/' and initial error tolerances of'/' RELER =',  
     &G12.5,' and ABSER =',G12.5/)
5020  format (' ***',12('****')/' *** The program is terminated at T =',
     &1PG12.5,'***'/' *** because an unidentifiable error has caused  **
     &*'/' *** STIFF to return with the illegal flag value ***'/' *** IF
     &LAG =',I11,26X,'***'/' ***',12('****'))
5030  format (' *WARNING--at T =',1PG12.5,', the error tolerances'/' pas
     &sed to STIFF were found to be too small for the'/' machine precisi
     &on.  They have been increased to'/' RELER =',G12.5,' and ABSER =',
     &G12.5/' to continue integration.'/)
5040  format (' *WARNING--at T =',1PG12.5,', STIFF was unable to'/' meet
     & the requested error tolerances at the smallest'/' allowable steps
     &ize.  Integration is continuing with'/' the tolerances increased t
     &o RELER =',G12.5/' and ABSER =',G12.5,/)
5050  format (' *WARNING--at T =',1PG12.5,', the number of function'/' c
     &alls exceeded the maximum number, MAXNFE, allowed'/' by STIFF.  Th
     &e function evaluation counter has been'/' reduced by MAXNFE to all
     &ow integration to proceed.')
5060  format (' The error tolerances have been increased to'/' RELER =',
     &   1PG12.5,' and ABSER =',G12.5,/)
5070  format (' *WARNING--at T =',1PG12.5,', the convergence test'/' in 
     &GEAR could not be met at the smallest'/' allowable stepsize.  Inte
     &gration is continuing'/' with the error tolerances increased to RE
     &LER ='/' ',G12.5,' and ABSER =',G12.5/)
5080  format (' *WARNING--at T =',1PG12.5,', STIFF returned with an'/' i
     &ndication that this integration is unnecessarily'/' costly because
     & the requested output interval is'/' substantially smaller than th
     &e natural stepsize.'/' (Refer to program documentation for alterna
     &tive'/' approaches.)  The output frequency monitor has been'/' res
     &et to allow integration to proceed.',/)
5090  format (' ***',12('****')/' *** The program is terminated at T =',
     &1PG12.5,'***'/' *** because a component of the solution',9X,'***'/
     &' *** vanished while ABSER = 0.  it will be      ***'/' *** necess
     &ary to provide a nonzero absolute     ***'/' *** error tolerance t
     &o perform this integration.***'/' ***',12('****'))
5100  format (' ***',12('****')/' *** the program is terminated at T =',
     &1PG12.5,'***'/' *** because invalid input was received by STIFF.**
     &*'/' *** check the following parameters--',12X,'***'/' ***  NEQN =
     &',I11,26X,'***'/' ***  RELER =',G12.5,23X,'***'/' ***  ABSER =',  
     &G12.5,23X,'***'/' ***  TINCR =',G12.5,24X,'***'/' *** Also, check 
     &that the output routine does    ***'/' *** not reset IFLAG to an i
     &llegal value.',8X,'***'/' ***',12('****'))
5110  format (' ***',12('****')/' *** The integration has been truncated
     & at time  ***'/' *** ',1PG12.1,' (hours), because STIFF has    ***
     &'/' *** returned irregularly several times,         ***'/' *** ind
     &icating extreme difficulty in the        ***'/' *** integration.  
     &Please check the accuracy of  ***'/' *** the problem statement.   
     &                   ***'/' ***',12('****'))
5120  format (' ***',12('****')/' *** Integration complete at T =',     
     &   1PG12.5,5X,'***'/' ***',12('****')///)
220   continue
      return
      end
