TALK=T;RUN(1,1) PHOTON USE p red msg u1 contours. Press return for w1 con u1 z 1 fi;0.001;con u1 z 5 fi;0.001;con u1 z 10 fi;0.001 gr ou z m;gr ou y m;gr ou x m pause con off;red msg w1 contours, and velocity vectors on the top. msg Press return for h1 con w1 y 1 fi;0.001;con w1 y 5 fi;0.001;con w1 y m fi;0.001 vec y m pause vec off;con off;red msg temperature contours. Press e to end con h1 z 1 fi;0.001;con h1 z 5 fi;0.001;con h1 z 10 fi;0.001 enduse DISPLAY This simulation is similar to that of case 249; but the cavity is limited in the third direction by adiabatic end walls, at rest. Economy could be affected by noting that symmetry must prevail about the NZ/2 plane; but this is not done here. The symmetry of the solution can therefore be checked.. _ _ _ _ _ _ _ / /| / north wall / | / at H1=1. --/--|--->moving /_ _ _ _ _ _ _/ | y | | | ^ ^ | |east wall | /z | low wall | /H1=0. | / | adiabatic | / |/----> x |_ _ _ _ _ _ _|/ EXPERT is in use The Q1 contains PHOTON USE commands ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Cubical Cavity With Moving Lid ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 788 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 5 XULAST =0.1 XFRAC(1)=0.08 ;XFRAC(2)=0.32 XFRAC(3)=0.68 ;XFRAC(4)=0.92 XFRAC(5)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 5 YVLAST =0.1 YFRAC(1)=0.08 ;YFRAC(2)=0.32 YFRAC(3)=0.68 ;YFRAC(4)=0.92 YFRAC(5)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 5 ZWLAST =0.1 ZFRAC(1)=0.08 ;ZFRAC(2)=0.32 ZFRAC(3)=0.68 ;ZFRAC(4)=0.92 ZFRAC(5)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(7)=W1 NAME(14)=H1 ;NAME(107)=GAH1 NAME(108)=ALH1 ;NAME(109)=AHH1 NAME(110)=ASH1 ;NAME(111)=ANH1 NAME(112)=AWH1 ;NAME(113)=AEH1 NAME(114)=APH1 ;NAME(115)=SUH1 NAME(116)=GAW1 ;NAME(117)=ALW1 NAME(118)=AHW1 ;NAME(119)=ASW1 NAME(120)=ANW1 ;NAME(121)=AWW1 NAME(122)=AEW1 ;NAME(123)=APW1 NAME(124)=SUW1 ;NAME(125)=GAV1 NAME(126)=ALV1 ;NAME(127)=AHV1 NAME(128)=ASV1 ;NAME(129)=ANV1 NAME(130)=AWV1 ;NAME(131)=AEV1 NAME(132)=APV1 ;NAME(133)=SUV1 NAME(134)=GAU1 ;NAME(135)=ALU1 NAME(136)=AHU1 ;NAME(137)=ASU1 NAME(138)=ANU1 ;NAME(139)=AWU1 NAME(140)=AEU1 ;NAME(141)=APU1 NAME(142)=SUU1 ;NAME(143)=ALP1 NAME(144)=AHP1 ;NAME(145)=ASP1 NAME(146)=ANP1 ;NAME(147)=AWP1 NAME(148)=AEP1 ;NAME(149)=APP1 NAME(150)=SUP1 * 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,Y,N,N,N) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(H1,Y,Y,N,N,N,Y) SOLUTN(GAH1,Y,N,N,N,N,N) SOLUTN(ALH1,Y,N,N,N,N,N) SOLUTN(AHH1,Y,N,N,N,N,N) SOLUTN(ASH1,Y,N,N,N,N,N) SOLUTN(ANH1,Y,N,N,N,N,N) SOLUTN(AWH1,Y,N,N,N,N,N) SOLUTN(AEH1,Y,N,N,N,N,N) SOLUTN(APH1,Y,N,N,N,N,N) SOLUTN(SUH1,Y,N,N,N,N,N) SOLUTN(GAW1,Y,N,N,N,N,N) SOLUTN(ALW1,Y,N,N,N,N,N) SOLUTN(AHW1,Y,N,N,N,N,N) SOLUTN(ASW1,Y,N,N,N,N,N) SOLUTN(ANW1,Y,N,N,N,N,N) SOLUTN(AWW1,Y,N,N,N,N,N) SOLUTN(AEW1,Y,N,N,N,N,N) SOLUTN(APW1,Y,N,N,N,N,N) SOLUTN(SUW1,Y,N,N,N,N,N) SOLUTN(GAV1,Y,N,N,N,N,N) SOLUTN(ALV1,Y,N,N,N,N,N) SOLUTN(AHV1,Y,N,N,N,N,N) SOLUTN(ASV1,Y,N,N,N,N,N) SOLUTN(ANV1,Y,N,N,N,N,N) SOLUTN(AWV1,Y,N,N,N,N,N) SOLUTN(AEV1,Y,N,N,N,N,N) SOLUTN(APV1,Y,N,N,N,N,N) SOLUTN(SUV1,Y,N,N,N,N,N) SOLUTN(GAU1,Y,N,N,N,N,N) SOLUTN(ALU1,Y,N,N,N,N,N) SOLUTN(AHU1,Y,N,N,N,N,N) SOLUTN(ASU1,Y,N,N,N,N,N) SOLUTN(ANU1,Y,N,N,N,N,N) SOLUTN(AWU1,Y,N,N,N,N,N) SOLUTN(AEU1,Y,N,N,N,N,N) SOLUTN(APU1,Y,N,N,N,N,N) SOLUTN(SUU1,Y,N,N,N,N,N) SOLUTN(ALP1,Y,N,N,N,N,N) SOLUTN(AHP1,Y,N,N,N,N,N) SOLUTN(ASP1,Y,N,N,N,N,N) SOLUTN(ANP1,Y,N,N,N,N,N) SOLUTN(AWP1,Y,N,N,N,N,N) SOLUTN(AEP1,Y,N,N,N,N,N) SOLUTN(APP1,Y,N,N,N,N,N) SOLUTN(SUP1,Y,N,N,N,N,N) ************************************************************ 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(U1,Y,Y,Y,Y,Y,Y) TERMS(V1,Y,Y,Y,Y,Y,Y) TERMS(W1,Y,Y,Y,Y,Y,Y) TERMS(H1,N,Y,Y,Y,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = 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 =1. ;TMP1 =0. EL1 =0. TSURR =0. ;TEMP0 =0. PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =9.999999E-04 ;ENUT =0. PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRNDTL(W1)=1. ;PRNDTL(H1)=0.7 PRT(U1)=1. ;PRT(V1)=1. PRT(W1)=1. ;PRT(H1)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10 FIINIT(H1)=1.0E-10 ;FIINIT(GAH1)=1.0E-10 FIINIT(ALH1)=1.0E-10 ;FIINIT(AHH1)=1.0E-10 FIINIT(ASH1)=1.0E-10 ;FIINIT(ANH1)=1.0E-10 FIINIT(AWH1)=1.0E-10 ;FIINIT(AEH1)=1.0E-10 FIINIT(APH1)=1.0E-10 ;FIINIT(SUH1)=1.0E-10 FIINIT(GAW1)=1.0E-10 ;FIINIT(ALW1)=1.0E-10 FIINIT(AHW1)=1.0E-10 ;FIINIT(ASW1)=1.0E-10 FIINIT(ANW1)=1.0E-10 ;FIINIT(AWW1)=1.0E-10 FIINIT(AEW1)=1.0E-10 ;FIINIT(APW1)=1.0E-10 FIINIT(SUW1)=1.0E-10 ;FIINIT(GAV1)=1.0E-10 FIINIT(ALV1)=1.0E-10 ;FIINIT(AHV1)=1.0E-10 FIINIT(ASV1)=1.0E-10 ;FIINIT(ANV1)=1.0E-10 FIINIT(AWV1)=1.0E-10 ;FIINIT(AEV1)=1.0E-10 FIINIT(APV1)=1.0E-10 ;FIINIT(SUV1)=1.0E-10 FIINIT(GAU1)=1.0E-10 ;FIINIT(ALU1)=1.0E-10 FIINIT(AHU1)=1.0E-10 ;FIINIT(ASU1)=1.0E-10 FIINIT(ANU1)=1.0E-10 ;FIINIT(AWU1)=1.0E-10 FIINIT(AEU1)=1.0E-10 ;FIINIT(APU1)=1000. FIINIT(SUU1)=1.0E-10 ;FIINIT(ALP1)=1.0E-10 FIINIT(AHP1)=1.0E-10 ;FIINIT(ASP1)=1.0E-10 FIINIT(ANP1)=1.0E-10 ;FIINIT(AWP1)=1.0E-10 FIINIT(AEP1)=1.0E-10 ;FIINIT(APP1)=1.0E-10 FIINIT(SUP1)=1.0E-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(SOUTH ,SWALL , 1, 10, 1, 1, 1, 10, 1, 1) COVAL(SOUTH ,U1 ,1. ,0. ) COVAL(SOUTH ,W1 ,1. ,0. ) COVAL(SOUTH ,H1 ,1. ,0. ) PATCH(MOVING ,NWALL , 1, 10, 10, 10, 1, 10, 1, 1) COVAL(MOVING ,U1 ,1. ,1. ) COVAL(MOVING ,W1 ,1. ,0. ) COVAL(MOVING ,H1 ,1. ,1. ) PATCH(WEST ,WWALL , 1, 1, 1, 10, 1, 10, 1, 1) COVAL(WEST ,V1 ,1. ,0. ) COVAL(WEST ,W1 ,1. ,0. ) COVAL(WEST ,H1 ,1. ,0. ) PATCH(EAST ,EWALL , 10, 10, 1, 10, 1, 10, 1, 1) COVAL(EAST ,V1 ,1. ,0. ) COVAL(EAST ,W1 ,1. ,0. ) COVAL(EAST ,H1 ,1.428571 ,0. ) PATCH(HIGH ,HWALL , 1, 10, 1, 10, 10, 10, 1, 1) COVAL(HIGH ,U1 ,1. ,0. ) COVAL(HIGH ,V1 ,1. ,0. ) PATCH(LOW ,LWALL , 1, 10, 1, 10, 1, 1, 1, 1) COVAL(LOW ,U1 ,1. ,0. ) COVAL(LOW ,V1 ,1. ,0. ) PATCH(RELIEF ,CELL , 5, 5, 5, 5, 1, 1, 1, 1) COVAL(RELIEF ,P1 ,1. ,0. ) COVAL(RELIEF ,H1 ,0. , SAME ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 10 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1. ************************************************************ Group 16. Terminate Iterations LITER(P1)=100 ;LITER(U1)=20 LITER(V1)=20 ;LITER(W1)=20 LITER(H1)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03 ENDIT(H1)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,1.0E-03) RELAX(V1,FALSDT,1.0E-03) RELAX(W1,FALSDT,1.0E-03) RELAX(H1,LINRLX,0.1) RELAX(GAH1,LINRLX,1.) RELAX(ALH1,LINRLX,1.) RELAX(AHH1,LINRLX,1.) RELAX(ASH1,LINRLX,1.) RELAX(ANH1,LINRLX,1.) RELAX(AWH1,LINRLX,1.) RELAX(AEH1,LINRLX,1.) RELAX(APH1,LINRLX,1.) RELAX(SUH1,LINRLX,1.) RELAX(GAW1,LINRLX,1.) RELAX(ALW1,LINRLX,1.) RELAX(AHW1,LINRLX,1.) RELAX(ASW1,LINRLX,1.) RELAX(ANW1,LINRLX,1.) RELAX(AWW1,LINRLX,1.) RELAX(AEW1,LINRLX,1.) RELAX(APW1,LINRLX,1.) RELAX(SUW1,LINRLX,1.) RELAX(GAV1,LINRLX,1.) RELAX(ALV1,LINRLX,1.) RELAX(AHV1,LINRLX,1.) RELAX(ASV1,LINRLX,1.) RELAX(ANV1,LINRLX,1.) RELAX(AWV1,LINRLX,1.) RELAX(AEV1,LINRLX,1.) RELAX(APV1,LINRLX,1.) RELAX(SUV1,LINRLX,1.) RELAX(GAU1,LINRLX,1.) RELAX(ALU1,LINRLX,1.) RELAX(AHU1,LINRLX,1.) RELAX(ASU1,LINRLX,1.) RELAX(ANU1,LINRLX,1.) RELAX(AWU1,LINRLX,1.) RELAX(AEU1,LINRLX,1.) RELAX(APU1,LINRLX,1.0E-03) RELAX(SUU1,LINRLX,1.) RELAX(ALP1,LINRLX,1.) RELAX(AHP1,LINRLX,1.) RELAX(ASP1,LINRLX,1.) RELAX(ANP1,LINRLX,1.) RELAX(AWP1,LINRLX,1.) RELAX(AEP1,LINRLX,1.) RELAX(APP1,LINRLX,1.) RELAX(SUP1,LINRLX,1.) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(GAH1)=1.0E+10 ;VARMIN(GAH1)=-1.0E+10 VARMAX(ALH1)=1.0E+10 ;VARMIN(ALH1)=-1.0E+10 VARMAX(AHH1)=1.0E+10 ;VARMIN(AHH1)=-1.0E+10 VARMAX(ASH1)=1.0E+10 ;VARMIN(ASH1)=-1.0E+10 VARMAX(ANH1)=1.0E+10 ;VARMIN(ANH1)=-1.0E+10 VARMAX(AWH1)=1.0E+10 ;VARMIN(AWH1)=-1.0E+10 VARMAX(AEH1)=1.0E+10 ;VARMIN(AEH1)=-1.0E+10 VARMAX(APH1)=1.0E+10 ;VARMIN(APH1)=-1.0E+10 VARMAX(SUH1)=1.0E+10 ;VARMIN(SUH1)=-1.0E+10 VARMAX(GAW1)=1.0E+10 ;VARMIN(GAW1)=-1.0E+10 VARMAX(ALW1)=1.0E+10 ;VARMIN(ALW1)=-1.0E+10 VARMAX(AHW1)=1.0E+10 ;VARMIN(AHW1)=-1.0E+10 VARMAX(ASW1)=1.0E+10 ;VARMIN(ASW1)=-1.0E+10 VARMAX(ANW1)=1.0E+10 ;VARMIN(ANW1)=-1.0E+10 VARMAX(AWW1)=1.0E+10 ;VARMIN(AWW1)=-1.0E+10 VARMAX(AEW1)=1.0E+10 ;VARMIN(AEW1)=-1.0E+10 VARMAX(APW1)=1.0E+10 ;VARMIN(APW1)=-1.0E+10 VARMAX(SUW1)=1.0E+10 ;VARMIN(SUW1)=-1.0E+10 VARMAX(GAV1)=1.0E+10 ;VARMIN(GAV1)=-1.0E+10 VARMAX(ALV1)=1.0E+10 ;VARMIN(ALV1)=-1.0E+10 VARMAX(AHV1)=1.0E+10 ;VARMIN(AHV1)=-1.0E+10 VARMAX(ASV1)=1.0E+10 ;VARMIN(ASV1)=-1.0E+10 VARMAX(ANV1)=1.0E+10 ;VARMIN(ANV1)=-1.0E+10 VARMAX(AWV1)=1.0E+10 ;VARMIN(AWV1)=-1.0E+10 VARMAX(AEV1)=1.0E+10 ;VARMIN(AEV1)=-1.0E+10 VARMAX(APV1)=1.0E+10 ;VARMIN(APV1)=-1.0E+10 VARMAX(SUV1)=1.0E+10 ;VARMIN(SUV1)=-1.0E+10 VARMAX(GAU1)=1.0E+10 ;VARMIN(GAU1)=-1.0E+10 VARMAX(ALU1)=1.0E+10 ;VARMIN(ALU1)=-1.0E+10 VARMAX(AHU1)=1.0E+10 ;VARMIN(AHU1)=-1.0E+10 VARMAX(ASU1)=1.0E+10 ;VARMIN(ASU1)=-1.0E+10 VARMAX(ANU1)=1.0E+10 ;VARMIN(ANU1)=-1.0E+10 VARMAX(AWU1)=1.0E+10 ;VARMIN(AWU1)=-1.0E+10 VARMAX(AEU1)=1.0E+10 ;VARMIN(AEU1)=-1.0E+10 VARMAX(APU1)=1.0E+10 ;VARMIN(APU1)=-1.0E+10 VARMAX(SUU1)=1.0E+10 ;VARMIN(SUU1)=-1.0E+10 VARMAX(ALP1)=1.0E+10 ;VARMIN(ALP1)=-1.0E+10 VARMAX(AHP1)=1.0E+10 ;VARMIN(AHP1)=-1.0E+10 VARMAX(ASP1)=1.0E+10 ;VARMIN(ASP1)=-1.0E+10 VARMAX(ANP1)=1.0E+10 ;VARMIN(ANP1)=-1.0E+10 VARMAX(AWP1)=1.0E+10 ;VARMIN(AWP1)=-1.0E+10 VARMAX(AEP1)=1.0E+10 ;VARMIN(AEP1)=-1.0E+10 VARMAX(APP1)=1.0E+10 ;VARMIN(APP1)=-1.0E+10 VARMAX(SUP1)=1.0E-03 ;VARMIN(SUP1)=-1.0E-03 ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F CONWIZ = T ISG50 = 1 ISG52 = 2 ISG62 = 1 SPEDAT(SET,STORED,SUP1,C,=RESI(P1)!RESIDU) SPEDAT(SET,STORED,APP1,C,=APCO(P1)!EQCOEF) SPEDAT(SET,STORED,AEP1,C,=AECO(P1)!EQCOEF) SPEDAT(SET,STORED,AWP1,C,=AWCO(P1)!EQCOEF) SPEDAT(SET,STORED,ANP1,C,=ANCO(P1)!EQCOEF) SPEDAT(SET,STORED,ASP1,C,=ASCO(P1)!EQCOEF) SPEDAT(SET,STORED,AHP1,C,=AHCO(P1)!EQCOEF) SPEDAT(SET,STORED,ALP1,C,=ALCO(P1)!EQCOEF) SPEDAT(SET,STORED,SUU1,C,=RESI(U1)!RESIDU) SPEDAT(SET,STORED,APU1,C,=APCO(U1)!EQCOEF) SPEDAT(SET,STORED,AEU1,C,=AECO(U1)!EQCOEF) SPEDAT(SET,STORED,AWU1,C,=AWCO(U1)!EQCOEF) SPEDAT(SET,STORED,ANU1,C,=ANCO(U1)!EQCOEF) SPEDAT(SET,STORED,ASU1,C,=ASCO(U1)!EQCOEF) SPEDAT(SET,STORED,AHU1,C,=AHCO(U1)!EQCOEF) SPEDAT(SET,STORED,ALU1,C,=ALCO(U1)!EQCOEF) SPEDAT(SET,STORED,GAU1,C,=GAMM(U1)!GAMCOF) SPEDAT(SET,STORED,SUV1,C,=RESI(V1)!RESIDU) SPEDAT(SET,STORED,APV1,C,=APCO(V1)!EQCOEF) SPEDAT(SET,STORED,AEV1,C,=AECO(V1)!EQCOEF) SPEDAT(SET,STORED,AWV1,C,=AWCO(V1)!EQCOEF) SPEDAT(SET,STORED,ANV1,C,=ANCO(V1)!EQCOEF) SPEDAT(SET,STORED,ASV1,C,=ASCO(V1)!EQCOEF) SPEDAT(SET,STORED,AHV1,C,=AHCO(V1)!EQCOEF) SPEDAT(SET,STORED,ALV1,C,=ALCO(V1)!EQCOEF) SPEDAT(SET,STORED,GAV1,C,=GAMM(V1)!GAMCOF) SPEDAT(SET,STORED,SUW1,C,=RESI(W1)!RESIDU) SPEDAT(SET,STORED,APW1,C,=APCO(W1)!EQCOEF) SPEDAT(SET,STORED,AEW1,C,=AECO(W1)!EQCOEF) SPEDAT(SET,STORED,AWW1,C,=AWCO(W1)!EQCOEF) SPEDAT(SET,STORED,ANW1,C,=ANCO(W1)!EQCOEF) SPEDAT(SET,STORED,ASW1,C,=ASCO(W1)!EQCOEF) SPEDAT(SET,STORED,AHW1,C,=AHCO(W1)!EQCOEF) SPEDAT(SET,STORED,ALW1,C,=ALCO(W1)!EQCOEF) SPEDAT(SET,STORED,GAW1,C,=GAMM(W1)!GAMCOF) SPEDAT(SET,STORED,SUH1,C,=RESI(H1)!RESIDU) SPEDAT(SET,STORED,APH1,C,=APCO(H1)!EQCOEF) SPEDAT(SET,STORED,AEH1,C,=AECO(H1)!EQCOEF) SPEDAT(SET,STORED,AWH1,C,=AWCO(H1)!EQCOEF) SPEDAT(SET,STORED,ANH1,C,=ANCO(H1)!EQCOEF) SPEDAT(SET,STORED,ASH1,C,=ASCO(H1)!EQCOEF) SPEDAT(SET,STORED,AHH1,C,=AHCO(H1)!EQCOEF) SPEDAT(SET,STORED,ALH1,C,=ALCO(H1)!EQCOEF) SPEDAT(SET,STORED,GAH1,C,=GAMM(H1)!GAMCOF) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=0.4606 ;EX(U1)=0.03798 EX(V1)=0.01968 ;EX(W1)=6.408E-05 EX(H1)=0.1072 ;EX(GAH1)=1.429E-03 EX(ALH1)=2.189E-05 ;EX(AHH1)=2.189E-05 EX(ASH1)=2.19E-05 ;EX(ANH1)=2.19E-05 EX(AWH1)=2.237E-05 ;EX(AEH1)=2.237E-05 EX(APH1)=1.581E-04 ;EX(SUH1)=1.462E-05 EX(GAW1)=8.0E-04 ;EX(ALW1)=8.882E-06 EX(AHW1)=8.882E-06 ;EX(ASW1)=1.244E-05 EX(ANW1)=1.244E-05 ;EX(AWW1)=1.15E-05 EX(AEW1)=1.15E-05 ;EX(APW1)=7.461E-03 EX(SUW1)=2.179E-07 ;EX(GAV1)=1.0E-03 EX(ALV1)=1.532E-05 ;EX(AHV1)=1.532E-05 EX(ASV1)=7.645E-06 ;EX(ANV1)=7.645E-06 EX(AWV1)=1.319E-05 ;EX(AEV1)=1.319E-05 EX(APV1)=0.01312 ;EX(SUV1)=8.753E-06 EX(GAU1)=1.0E-03 ;EX(ALU1)=1.532E-05 EX(AHU1)=1.532E-05 ;EX(ASU1)=1.397E-05 EX(ANU1)=1.397E-05 ;EX(AWU1)=6.775E-06 EX(AEU1)=6.775E-06 ;EX(APU1)=8.119E-03 EX(SUU1)=1.853E-05 ;EX(ALP1)=1.531E-05 EX(AHP1)=1.531E-05 ;EX(ASP1)=1.532E-05 EX(ANP1)=1.532E-05 ;EX(AWP1)=1.531E-05 EX(AEP1)=1.531E-05 ;EX(APP1)=8.0E-03 EX(SUP1)=3.272E-08 ************************************************************ 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(U1,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(GAH1,Y,N,Y,N,N,N) OUTPUT(ALH1,Y,N,Y,N,N,N) OUTPUT(AHH1,Y,N,Y,N,N,N) OUTPUT(ASH1,Y,N,Y,N,N,N) OUTPUT(ANH1,Y,N,Y,N,N,N) OUTPUT(AWH1,Y,N,Y,N,N,N) OUTPUT(AEH1,Y,N,Y,N,N,N) OUTPUT(APH1,Y,N,Y,N,N,N) OUTPUT(SUH1,Y,N,Y,N,N,N) OUTPUT(GAW1,Y,N,Y,N,N,N) OUTPUT(ALW1,Y,N,Y,N,N,N) OUTPUT(AHW1,Y,N,Y,N,N,N) OUTPUT(ASW1,Y,N,Y,N,N,N) OUTPUT(ANW1,Y,N,Y,N,N,N) OUTPUT(AWW1,Y,N,Y,N,N,N) OUTPUT(AEW1,Y,N,Y,N,N,N) OUTPUT(APW1,Y,N,Y,N,N,N) OUTPUT(SUW1,Y,N,Y,N,N,N) OUTPUT(GAV1,Y,N,Y,N,N,N) OUTPUT(ALV1,Y,N,Y,N,N,N) OUTPUT(AHV1,Y,N,Y,N,N,N) OUTPUT(ASV1,Y,N,Y,N,N,N) OUTPUT(ANV1,Y,N,Y,N,N,N) OUTPUT(AWV1,Y,N,Y,N,N,N) OUTPUT(AEV1,Y,N,Y,N,N,N) OUTPUT(APV1,Y,N,Y,N,N,N) OUTPUT(SUV1,Y,N,Y,N,N,N) OUTPUT(GAU1,Y,N,Y,N,N,N) OUTPUT(ALU1,Y,N,Y,N,N,N) OUTPUT(AHU1,Y,N,Y,N,N,N) OUTPUT(ASU1,Y,N,Y,N,N,N) OUTPUT(ANU1,Y,N,Y,N,N,N) OUTPUT(AWU1,Y,N,Y,N,N,N) OUTPUT(AEU1,Y,N,Y,N,N,N) OUTPUT(APU1,Y,N,Y,N,N,N) OUTPUT(SUU1,Y,N,Y,N,N,N) OUTPUT(ALP1,Y,N,Y,N,N,N) OUTPUT(AHP1,Y,N,Y,N,N,N) OUTPUT(ASP1,Y,N,Y,N,N,N) OUTPUT(ANP1,Y,N,Y,N,N,N) OUTPUT(AWP1,Y,N,Y,N,N,N) OUTPUT(AEP1,Y,N,Y,N,N,N) OUTPUT(APP1,Y,N,Y,N,N,N) OUTPUT(SUP1,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 5 ;IYMON = 9 ;IZMON = 5 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = T ;YZPR = F IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 PATCH(MIDX ,PROFIL, 5, 5, 1, 10, 5, 5, 1, 1) PLOT(MIDX ,U1 ,0. ,0. ) PLOT(MIDX ,H1 ,0. ,0. ) PATCH(MAP ,CONTUR, 1, 10, 1, 10, 5, 5, 1, 1) PLOT(MAP ,U1 ,0. ,10. ) PLOT(MAP ,V1 ,0. ,10. ) PLOT(MAP ,W1 ,0. ,10. ) PLOT(MAP ,H1 ,0. ,10. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP