TALK=T;RUN(1,1)
  photon use
  p
  parphi
  1 1 1000
 
  gr ou y 1
  msg velocity vectors of cooler gas; press return for hotter
  set vec comp
  1u v1 w1
  vec y 1 sh;pause;vec off;red
  msg velocity vectors of hotter gas; press return for pressures
  set vec comp 2u v1 w1
  vec y 1 sh;pause;vec off;red
  con p1 y 1 fi;0.01
  msg flow is upward at first; time is horizontal
  msg pressure contours; press return for hot-gas reactedness
  msg flow is upward at first; time is horizontal
  pause;con off;red;con rct2 y 1 fi;0.01
  msg hot-gas reactedness contours; press return for cold-gas
  msg flow is upward at first; time is horizontal
  pause;con off;red;con rct1 y 1 fi;0.01
  msg cold-gas reactedness contours
 
  enduse
  DISPLAY
  Intermingled fragments of burned and unburned combustible gas are
  set in relative motion by a shock wave. The increased entrainment
  rate leads to further chemical reaction, which, depending on the
  presumptions made for fragment size and reactivity, may lead to
  detonation. The following cases provide a parametric study of the
  influences of the interfluid friction factor (which increases as
  fragment size diminishes) and the inter-fluid entrainment factor
  (which is expected to be of the order of 1.0).
  PHOTON commands are provided to facilitate display.
 
    case 1  cfipc=0.0      ;cmdot=-1.0
    case 2  cfipc=1.0      ;cmdot=-1.0
    case 3  cfipc=1.0e1    ;cmdot=-1.0
    case 4  cfipc=2.0e1    ;cmdot=-1.0
    case 5  cfipc=5.0e1    ;cmdot=-1.0
    case 6  cfipc=1.0e2    ;cmdot=-1.0
    case 7  cfipc=1.0      ;cmdot=-1.0e1
    case 8  cfipc=5.0e1    ;cmdot=-1.0e1
    case 9  cfipc=1.0e1    ;cmdot=-1.0e1
    case 10 cfipc=1.0e2    ;cmdot=-1.0e1
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(1-D transient shock-induced prpgtn :W979)
 
 ************************************************************
 ************************************************************
 
 IRUNN   =         1 ;LIBREF =        14
 ************************************************************
  Group 2. Time dependence
 STEADY  =    F
    * Set overall time and no. of steps
 TFIRST  = 0.000000E+00 ;TLAST  = 1.000000E-03
 FSTEP   =         1 ;LSTEP  =       100
 TFRAC (   1) =  1.000000E-02 ;TFRAC (   2) =  2.000000E-02
 TFRAC (   3) =  3.000000E-02 ;TFRAC (   4) =  4.000000E-02
 TFRAC (   5) =  5.000000E-02 ;TFRAC (   6) =  6.000000E-02
 TFRAC (   7) =  7.000000E-02 ;TFRAC (   8) =  8.000000E-02
 TFRAC (   9) =  9.000000E-02 ;TFRAC (  10) =  1.000000E-01
 TFRAC (  11) =  1.100000E-01 ;TFRAC (  12) =  1.200000E-01
 TFRAC (  13) =  1.300000E-01 ;TFRAC (  14) =  1.400000E-01
 TFRAC (  15) =  1.500000E-01 ;TFRAC (  16) =  1.600000E-01
 TFRAC (  17) =  1.700000E-01 ;TFRAC (  18) =  1.800000E-01
 TFRAC (  19) =  1.900000E-01 ;TFRAC (  20) =  2.000000E-01
 TFRAC (  21) =  2.100000E-01 ;TFRAC (  22) =  2.200000E-01
 TFRAC (  23) =  2.300000E-01 ;TFRAC (  24) =  2.400000E-01
 TFRAC (  25) =  2.500000E-01 ;TFRAC (  26) =  2.600000E-01
 TFRAC (  27) =  2.700000E-01 ;TFRAC (  28) =  2.800000E-01
 TFRAC (  29) =  2.900000E-01 ;TFRAC (  30) =  3.000000E-01
 TFRAC (  31) =  3.100000E-01 ;TFRAC (  32) =  3.200000E-01
 TFRAC (  33) =  3.300000E-01 ;TFRAC (  34) =  3.400000E-01
 TFRAC (  35) =  3.500000E-01 ;TFRAC (  36) =  3.600000E-01
 TFRAC (  37) =  3.700000E-01 ;TFRAC (  38) =  3.800000E-01
 TFRAC (  39) =  3.900000E-01 ;TFRAC (  40) =  4.000000E-01
 TFRAC (  41) =  4.100000E-01 ;TFRAC (  42) =  4.200000E-01
 TFRAC (  43) =  4.300000E-01 ;TFRAC (  44) =  4.400000E-01
 TFRAC (  45) =  4.500000E-01 ;TFRAC (  46) =  4.600000E-01
 TFRAC (  47) =  4.700000E-01 ;TFRAC (  48) =  4.800000E-01
 TFRAC (  49) =  4.900000E-01 ;TFRAC (  50) =  5.000000E-01
 TFRAC (  51) =  5.100000E-01 ;TFRAC (  52) =  5.200000E-01
 TFRAC (  53) =  5.300000E-01 ;TFRAC (  54) =  5.400000E-01
 TFRAC (  55) =  5.500000E-01 ;TFRAC (  56) =  5.600000E-01
 TFRAC (  57) =  5.700000E-01 ;TFRAC (  58) =  5.800000E-01
 TFRAC (  59) =  5.900000E-01 ;TFRAC (  60) =  6.000000E-01
 TFRAC (  61) =  6.100000E-01 ;TFRAC (  62) =  6.199999E-01
 TFRAC (  63) =  6.300000E-01 ;TFRAC (  64) =  6.400000E-01
 TFRAC (  65) =  6.500000E-01 ;TFRAC (  66) =  6.600000E-01
 TFRAC (  67) =  6.700000E-01 ;TFRAC (  68) =  6.800000E-01
 TFRAC (  69) =  6.900000E-01 ;TFRAC (  70) =  7.000000E-01
 TFRAC (  71) =  7.100000E-01 ;TFRAC (  72) =  7.200000E-01
 TFRAC (  73) =  7.300000E-01 ;TFRAC (  74) =  7.400000E-01
 TFRAC (  75) =  7.499999E-01 ;TFRAC (  76) =  7.600000E-01
 TFRAC (  77) =  7.700000E-01 ;TFRAC (  78) =  7.800000E-01
 TFRAC (  79) =  7.900000E-01 ;TFRAC (  80) =  8.000000E-01
 TFRAC (  81) =  8.100000E-01 ;TFRAC (  82) =  8.200000E-01
 TFRAC (  83) =  8.300000E-01 ;TFRAC (  84) =  8.400000E-01
 TFRAC (  85) =  8.500000E-01 ;TFRAC (  86) =  8.600000E-01
 TFRAC (  87) =  8.700000E-01 ;TFRAC (  88) =  8.800000E-01
 TFRAC (  89) =  8.900000E-01 ;TFRAC (  90) =  9.000000E-01
 TFRAC (  91) =  9.100000E-01 ;TFRAC (  92) =  9.200000E-01
 TFRAC (  93) =  9.300000E-01 ;TFRAC (  94) =  9.400000E-01
 TFRAC (  95) =  9.499999E-01 ;TFRAC (  96) =  9.600000E-01
 TFRAC (  97) =  9.700000E-01 ;TFRAC (  98) =  9.800000E-01
 TFRAC (  99) =  9.900000E-01 ;TFRAC ( 100) =  1.000000E+00
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NX      =       100
 XULAST  = 1.000000E+00
 XFRAC (   1) =  1.000000E-02 ;XFRAC (   2) =  2.000000E-02
 XFRAC (   3) =  3.000000E-02 ;XFRAC (   4) =  4.000000E-02
 XFRAC (   5) =  5.000000E-02 ;XFRAC (   6) =  6.000000E-02
 XFRAC (   7) =  7.000000E-02 ;XFRAC (   8) =  8.000000E-02
 XFRAC (   9) =  9.000000E-02 ;XFRAC (  10) =  1.000000E-01
 XFRAC (  11) =  1.100000E-01 ;XFRAC (  12) =  1.200000E-01
 XFRAC (  13) =  1.300000E-01 ;XFRAC (  14) =  1.400000E-01
 XFRAC (  15) =  1.500000E-01 ;XFRAC (  16) =  1.600000E-01
 XFRAC (  17) =  1.700000E-01 ;XFRAC (  18) =  1.800000E-01
 XFRAC (  19) =  1.900000E-01 ;XFRAC (  20) =  2.000000E-01
 XFRAC (  21) =  2.100000E-01 ;XFRAC (  22) =  2.200000E-01
 XFRAC (  23) =  2.300000E-01 ;XFRAC (  24) =  2.400000E-01
 XFRAC (  25) =  2.500000E-01 ;XFRAC (  26) =  2.600000E-01
 XFRAC (  27) =  2.700000E-01 ;XFRAC (  28) =  2.800000E-01
 XFRAC (  29) =  2.900000E-01 ;XFRAC (  30) =  3.000000E-01
 XFRAC (  31) =  3.100000E-01 ;XFRAC (  32) =  3.200000E-01
 XFRAC (  33) =  3.300000E-01 ;XFRAC (  34) =  3.400000E-01
 XFRAC (  35) =  3.500000E-01 ;XFRAC (  36) =  3.600000E-01
 XFRAC (  37) =  3.700000E-01 ;XFRAC (  38) =  3.800000E-01
 XFRAC (  39) =  3.900000E-01 ;XFRAC (  40) =  4.000000E-01
 XFRAC (  41) =  4.100000E-01 ;XFRAC (  42) =  4.200000E-01
 XFRAC (  43) =  4.300000E-01 ;XFRAC (  44) =  4.400000E-01
 XFRAC (  45) =  4.500000E-01 ;XFRAC (  46) =  4.600000E-01
 XFRAC (  47) =  4.700000E-01 ;XFRAC (  48) =  4.800000E-01
 XFRAC (  49) =  4.900000E-01 ;XFRAC (  50) =  5.000000E-01
 XFRAC (  51) =  5.100000E-01 ;XFRAC (  52) =  5.200000E-01
 XFRAC (  53) =  5.300000E-01 ;XFRAC (  54) =  5.400000E-01
 XFRAC (  55) =  5.500000E-01 ;XFRAC (  56) =  5.600000E-01
 XFRAC (  57) =  5.700000E-01 ;XFRAC (  58) =  5.800000E-01
 XFRAC (  59) =  5.900000E-01 ;XFRAC (  60) =  6.000000E-01
 XFRAC (  61) =  6.100000E-01 ;XFRAC (  62) =  6.200000E-01
 XFRAC (  63) =  6.300000E-01 ;XFRAC (  64) =  6.400000E-01
 XFRAC (  65) =  6.500000E-01 ;XFRAC (  66) =  6.600000E-01
 XFRAC (  67) =  6.700000E-01 ;XFRAC (  68) =  6.800000E-01
 XFRAC (  69) =  6.900000E-01 ;XFRAC (  70) =  7.000000E-01
 XFRAC (  71) =  7.100000E-01 ;XFRAC (  72) =  7.200000E-01
 XFRAC (  73) =  7.300000E-01 ;XFRAC (  74) =  7.400000E-01
 XFRAC (  75) =  7.500000E-01 ;XFRAC (  76) =  7.600000E-01
 XFRAC (  77) =  7.700000E-01 ;XFRAC (  78) =  7.800000E-01
 XFRAC (  79) =  7.900000E-01 ;XFRAC (  80) =  8.000000E-01
 XFRAC (  81) =  8.100000E-01 ;XFRAC (  82) =  8.200000E-01
 XFRAC (  83) =  8.300000E-01 ;XFRAC (  84) =  8.400000E-01
 XFRAC (  85) =  8.500000E-01 ;XFRAC (  86) =  8.600000E-01
 XFRAC (  87) =  8.700000E-01 ;XFRAC (  88) =  8.800000E-01
 XFRAC (  89) =  8.900000E-01 ;XFRAC (  90) =  9.000000E-01
 XFRAC (  91) =  9.100000E-01 ;XFRAC (  92) =  9.200000E-01
 XFRAC (  93) =  9.300000E-01 ;XFRAC (  94) =  9.400000E-01
 XFRAC (  95) =  9.500000E-01 ;XFRAC (  96) =  9.600000E-01
 XFRAC (  97) =  9.700000E-01 ;XFRAC (  98) =  9.800000E-01
 XFRAC (  99) =  9.900000E-01 ;XFRAC ( 100) =  1.000000E+00
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY      =         1
 YVLAST  = 1.000000E+00
 YFRAC (   1) =  1.000000E+00
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB   =    F
 NZ      =         1
 ZWLAST  = 1.000000E+00
 ZFRAC (   1) =  1.000000E+00
 ************************************************************
  Group 6. Body-Fitted Coordinates
    * X-cyclic boundaries switched
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    F
 EQUVEL  =    F
 NAME(  1) =P1   ;NAME(  3) =1U
 NAME(  4) =2U   ;NAME(  9) =COLD
 NAME( 10) =HOT  ;NAME( 16) =RCT1
 NAME( 17) =RCT2 ;NAME(150) =MDOT
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging
 SOLUTN(P1  ,Y,Y,N,N,N,Y)
 SOLUTN(1U  ,Y,Y,N,N,N,Y)
 SOLUTN(2U  ,Y,Y,N,N,N,Y)
 SOLUTN(COLD,Y,Y,N,N,N,Y)
 SOLUTN(HOT ,Y,Y,N,N,N,Y)
 SOLUTN(RCT1,Y,Y,N,N,N,Y)
 SOLUTN(RCT2,Y,Y,N,N,N,Y)
 SOLUTN(MDOT,Y,N,N,N,N,Y)
 INTMDT  =       150
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport
 TERMS (P1  ,Y,Y,Y,N,Y,Y)
 TERMS (1U  ,Y,Y,Y,Y,Y,Y)
 TERMS (2U  ,Y,Y,Y,Y,N,Y)
 TERMS (COLD,Y,Y,Y,Y,Y,Y)
 TERMS (HOT ,Y,Y,Y,Y,N,Y)
 TERMS (RCT1,N,Y,Y,Y,Y,Y)
 TERMS (RCT2,N,Y,Y,Y,N,Y)
 DIFCUT  = 5.000000E-01 ;ZDIFAC = 1.000000E+00
 GALA    =    F  ;ADDDIF =    F
 NEWRH1  =    T
 NEWRH2  =    T
 EQDVDP  =    F
 ISOLX   =        -1 ;ISOLY  =        -1 ;ISOLZ  =        -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1    = GRND3        ;TMP1   = 0.000000E+00
 EL1     = 0.000000E+00
 TSURR   = 0.000000E+00 ;TEMP0  = 0.000000E+00
 PRESS0  = 1.000000E+05
 DVO1DT  = 0.000000E+00 ;DRH1DP = GRND3
 RHO1A   = 2.682695E-04 ;RHO1B  = 7.142857E-01
 RHO1C   = 0.000000E+00
 EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
 RADIA   = 0.000000E+00 ;RADIB  = 0.000000E+00
 RHO2    = GRND3        ;TMP2   = 0.000000E+00
 EL2     = 0.000000E+00
 DVO2DT  = 0.000000E+00 ;DRH2DP = GRND3
 RHO2A   = 5.365391E-05 ;RHO2B  = 7.142857E-01
 RHO2C   = 0.000000E+00
 ENUL    = 1.000000E-05 ;ENUT   = 0.000000E+00
 DRH2DP  = GRND3
 PHINT (1U  ) = -2.022000E+04 ;PHINT (2U  ) = -2.022000E+04
 PHINT (COLD) = -2.022000E+04 ;PHINT (HOT ) = -2.022000E+04
 PHINT (RCT1) = -2.022000E+04 ;PHINT (RCT2) = -2.022000E+04
 PRNDTL(1U  ) =  1.000000E+00 ;PRNDTL(2U  ) =  1.000000E+00
 PRNDTL(COLD) =  1.000000E+00 ;PRNDTL(HOT ) =  1.000000E+00
 PRNDTL(RCT1) =  1.000000E+00 ;PRNDTL(RCT2) =  1.000000E+00
 PRT   (1U  ) =  1.000000E+00 ;PRT   (2U  ) =  1.000000E+00
 PRT   (COLD) =  1.000000E+00 ;PRT   (HOT ) =  1.000000E+00
 PRT   (RCT1) =  1.000000E+00 ;PRT   (RCT2) =  1.000000E+00
 CP1     = 1.000000E+00 ;CP2    = 1.000000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS   = GRND2
 RLOLIM  = 0.000000E+00 ;CMDOT  =-1.000000E+00
 CFIPA   = 0.000000E+00 ;CFIPB  = 0.000000E+00
 CFIPC   = 0.000000E+00 ;CFIPD  = 0.000000E+00
 CINT  (P1  ) =  1.000000E+00 ;CINT  (1U  ) =  1.000000E+00
 CINT  (2U  ) =  1.000000E+00 ;CINT  (COLD) =  1.000000E+00
 CINT  (HOT ) =  1.000000E+00 ;CINT  (RCT1) =  1.000000E+00
 CINT  (RCT2) =  1.000000E+00
 PHINT (P1  ) = -2.022000E+04 ;PHINT (1U  ) = -2.022000E+04
 PHINT (2U  ) = -2.022000E+04 ;PHINT (COLD) = -2.022000E+04
 PHINT (HOT ) = -2.022000E+04 ;PHINT (RCT1) = -2.022000E+04
 PHINT (RCT2) = -2.022000E+04
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1  ) =  0.000000E+00 ;FIINIT(1U  ) =  1.000000E+01
 FIINIT(2U  ) =  1.000000E+01 ;FIINIT(COLD) =  9.000000E-01
 FIINIT(HOT ) =  1.000000E-01 ;FIINIT(RCT1) =  1.000000E-10
 FIINIT(RCT2) =  1.000000E+00 ;FIINIT(MDOT) =  1.000000E-10
   No PATCHes yet used for this Group
 INIADD  =    F
 FSWEEP  =         1
 NAMFI   =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(ENTRY   ,WEST  ,   1,   1,   1,   1,   1,   1,   1, 100)
 COVAL(ENTRY   ,1U  , 0.000000E+00, 1.000000E+01)
 COVAL(ENTRY   ,2U  , 0.000000E+00, 1.000000E+01)
 COVAL(ENTRY   ,COLD, FIXFLU      , 9.000000E+00)
 COVAL(ENTRY   ,HOT , FIXFLU      , 2.000000E-01)
 COVAL(ENTRY   ,RCT1, 0.000000E+00, 0.000000E+00)
 COVAL(ENTRY   ,RCT2, 0.000000E+00, 1.000000E+00)
 
 PATCH(CHSOTERM,VOLUME,   1, 100,   1,   1,   1,   1,   1, 100)
 COVAL(CHSOTERM,RCT2, GRND7       , 1.000000E+00)
   The following patch does not affect the calculation.
   It has been inserted as a signal to the Commander Search Engine
   that this is a two-fluid-model case
 PATCH(SHSO,CELL,1,1,1,1,1,1,1,1)
 XCYCLE  =    F
 EGWF    =    F
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =         4 ;ISWC1  =         1
 LITHYD  =      1 ;LITFLX =      1 ;LITC   =      1 ;ITHC1  =      1
 ISWR1   =         1 ;ISWR2  =     10000
 SELREF  =    T
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER (P1  ) =       20 ;LITER (1U  ) =       10
 LITER (2U  ) =       10 ;LITER (COLD) =        1
 LITER (HOT ) =        1 ;LITER (RCT1) =       20
 LITER (RCT2) =       20
 ENDIT (P1  ) =  1.000000E-03 ;ENDIT (1U  ) =  1.000000E-03
 ENDIT (2U  ) =  1.000000E-03 ;ENDIT (COLD) =  1.000000E-03
 ENDIT (HOT ) =  1.000000E-03 ;ENDIT (RCT1) =  1.000000E-03
 ENDIT (RCT2) =  1.000000E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(1U  ,FALSDT, 1.000000E-01)
 RELAX(2U  ,FALSDT, 1.000000E-01)
 RELAX(COLD,LINRLX, 4.000000E-01)
 RELAX(HOT ,LINRLX, 4.000000E-01)
 RELAX(RCT1,FALSDT, 1.000000E+09)
 RELAX(RCT2,FALSDT, 1.000000E+09)
 RELAX(MDOT,LINRLX, 1.000000E+00)
 OVRRLX  = 0.000000E+00
 EXPERT  =    F  ;NNORSL =    F
 ************************************************************
  Group 18. Limits
 VARMAX(P1  ) = 1.000000E+10 ;VARMIN(P1  ) =-1.000000E+10
 VARMAX(1U  ) = 1.000000E+06 ;VARMIN(1U  ) =-1.000000E+06
 VARMAX(2U  ) = 1.000000E+06 ;VARMIN(2U  ) =-1.000000E+06
 VARMAX(COLD) = 1.000000E+00 ;VARMIN(COLD) = 1.000000E-06
 VARMAX(HOT ) = 1.000000E+00 ;VARMIN(HOT ) = 1.000000E-06
 VARMAX(RCT1) = 1.000000E+10 ;VARMIN(RCT1) =-1.000000E+10
 VARMAX(RCT2) = 1.000000E+10 ;VARMIN(RCT2) =-1.000000E+10
 VARMAX(MDOT) = 1.000000E+10 ;VARMIN(MDOT) =-1.000000E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 USEGRD  =    T  ;USEGRX =    T
 CHSOA   = 1.000000E+05 ;CHSOB  = 6.000000E+00
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD  =    F  ;SUBWGR =    F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1  ,Y,N,Y,Y,Y,Y)
 OUTPUT(1U  ,Y,N,Y,Y,Y,Y)
 OUTPUT(2U  ,Y,N,Y,Y,Y,Y)
 OUTPUT(COLD,N,N,N,N,N,N)
 OUTPUT(HOT ,Y,N,Y,Y,Y,Y)
 OUTPUT(RCT1,Y,N,Y,Y,Y,Y)
 OUTPUT(RCT2,Y,N,Y,Y,Y,Y)
 OUTPUT(MDOT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        99 ;IYMON  =         1 ;IZMON  =         1
 NPRMON  =    100000 ;NPRMNT =         1 ;TSTSWP =        -1
 UWATCH  =    F  ;USTEER =    T
 HIGHLO  =    F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000 ;NUMCLS =         5
 NTPRIN  =        20 ;ISTPRF =         1 ;ISTPRL =     10000
 NXPRIN  =        10 ;IXPRF  =         1 ;IXPRL  =     10000
 IPLTF   =         2 ;IPLTL  =         4 ;NPLT   =        -1
 ISWPRF  =         1 ;ISWPRL =    100000
 ITABL   =         1 ;IPROF  =         1
 ABSIZ   = 5.000000E-01 ;ORSIZ  = 2.000000E-01
 NTZPRF  =         1 ;NCOLPF =        50
 ICHR    =         2 ;NCOLCO =        45 ;NROWCO =        20
 
 PATCH(PROFIL1 ,PROFIL,   1, 100,   1,   1,   1,   1,   1, 100)
 PLOT(PROFIL1 ,P1  , 0.000000E+00, 0.000000E+00)
 PLOT(PROFIL1 ,COLD, 0.000000E+00, 1.000000E+00)
 
 PATCH(PROFILE ,PROFIL,   1, 100,   1,   1,   1,   1,   1, 100)
 PLOT(PROFILE ,1U  ,-1.000000E+00,-1.000000E+00)
 PLOT(PROFILE ,2U  ,-1.000000E+00,-1.000000E+00)
 PLOT(PROFILE ,RCT2, 0.000000E+00, 1.000000E+00)
 PLOT(PROFILE ,MDOT, 0.000000E+00, 0.000000E+00)
 
 PATCH(TIMEPLOT,PROFIL, 100, 100,   1,   1,   1,   1,   1, 100)
 PLOT(TIMEPLOT,HOT , 0.000000E+00, 0.000000E+00)
   
  inform7begin
  In-Form statements: MAKE, STORE1 and PRINT
(make of ud75 is 0)
(store1 of ud75 is 1u[75]-2u[75])
(print of veldiff  ix=75 is ud75)
(make of pd75 is 0)
(store1 of pd75 is p1[75]-p1[76])
(print of p1_diff  ix=75 is pd75)
(print of veldif at timeplot  is ud75)
(print of p1_dif at timeplot  is pd75)
  inform7end
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE    =    T  ;AUTOPS =    F  ;NOWIPE =    F
 NSAVE   =CHAM
 IDISPA  =         5 ;IDISPB =         0 ;IDISPC =         0
STOP