    0                                     IPRT1  gureghian 3a (1983, ex 3)
    1        1        .F.                 IFLOW    ITRANS  OSC_FLAG
    0        0.01                         KSLOPE   TOLKSL  nt = 9000
    -3.0     -1.0      -3.0      -1.0     PKRL     PKRR    PSEL PSER
    -3.0     -2.5      -1.5      -1.0     PDSE1L   PDSE1R  PDSE2L PDSE2R
    0        1                            ISFONE   ISFCVG 
    0.5      0         0.5      0         TETAF    LUMP    TETAC    LUMPC
  100        12        20       1.e-4     ITMX     ITMX1   ITMX2    TOLNL
  100        7         12       1.e-4     ITMXC    ITMXC1  ITMXC2   TOLNLC
    0        1.0e+30                      L2NORM   ERNLMX
    0        0         0.8                ISOLV    IRELAX  OMEGA
    0        0         1.0                ISOLVC   IRELAXC OMEGAC
   100    1.0e-10   500  1.0e-10          ITMXCGSY TOLCGSY ITMXCGNS TOLCGNS
   0.0  1.0E+2 1.0e+2 1.0E+7 6.30E+4      TIMEP  DELTAT  DTMIN    DTMAX  TMAX
   0.0       2.0       0.0       0.6      DTMAGA   DTMAGM  DTREDS   DTREDM
    0        0         0      500.0       IPRT     ISEC    NPRT  TIMPR(1:NPRT)
    4   607   654   623   659             NUMVP    NODVP(1:NUMVP)
    0                                     NR (# of output nodes)
    6       534      1238      1766
  138       666      1370      1898
  171       699      1403      1931       CONTR(I),I=1,NR

   0.0  1.0E+6 1.0e+6 1.0E+7 6.30E+8      TIMEP  DELTAT  DTMIN    DTMAX  TMAX
   0.0  1.0E+36 1.0e+6 1.0E+7 6.30E+8      TIMEP  DELTAT  DTMIN    DTMAX  TMAX
   0.0  1.0E+2 1.0e+2 1.0E+7 6.30E+8      TIMEP  DELTAT  DTMIN    DTMAX  TMAX
