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