TALK=T;RUN(1,1)
DISPLAY
A rectangular plate with a centrally-placed CIRCULAR hole is
extended by uniform equal forces on its Top and Bottom
surfaces, which are allowed to distort.
By reason of symmetry, only one quarter of the plate is analysed.
For calculation is used POLAR coordinates and Boundary Flow
Cells.
Comparison - a analytical solution (Timoshenko )
Boundary Condition:
===================
V ( UR) :
---------
Y = RO : StrY = FX/2*(1-RI^2/RO^2) + FX/2*(1 + 3*RI^4/RO^4 -
4*RI^2/RO^2)*COS(2*X)
U ( UTeta) :
------------
X = 0, PI/2 : U = 0,
Y = RO : StXY = - FX/2 * (1 - 3*RI^4/RO^4 + 2*RI^2/RO^2)
*SIN(2*X)
Exact Solution:
===============
Displacements
-------------
V = FX/(2*E)* R* {(1-P) + RI^2/R^2*(1+P) + [(1 - RI^4/R^4)*(1+P)
+ 4*RI^2/R^2]*COS(2*X) }
U = -FX/(2*E)* R* [(1 + RI^4/R^4)*(1+P) + 2*RI^2/R^2*(1-P)]*
SIN(2*X)
Strains
-------
EpsY = FX/(2*E)* {(1-P) - RI^2/R^2*(1+P) + [(1+P)*(1 +
3*RI^4/R^4) - 4*RI^2/R^2]*COS(2*X)}
EpsX = FX/(2*E)* {(1-P) + RI^2/R^2*(1+P) - [(1+P)*(1 +
3*RI^4/R^4) - P * 4*RI^2/R^2]*COS(2*X)}
EpsZ = - P * FX/(2*E) * [ 2 - 4*RI^2/R^2*COS(2*X) ]
D = FX/E*(1-2*P) * [1 - 2*RI^2/R^2*COS(2*X) ]
Stress
------
StrY = FX/2*(1-RI^2/R^2) + FX/2*(1 + 3*RI^4/R^4 - 4*RI^2/R^2)*
COS(2*X)
StrX = FX/2*(1+RI^2/R^2) - FX/2*(1 + 3*RI^4/R^4)*COS(2*X)
StXY = - FX/2 * (1 - 3*RI^4/R^4 + 2*RI^2/R^2)*SIN(2*X)
StrZ = 0
ENDDIS
PHOTON USE
p;;;;
set prop off
msg
msg Grid
msg
gr Z 1
pause
cl
msg
msg DISY contours
msg
cont DISY z 1 fil;.0001
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 1 SH DASH 0
pause
cl
msg
msg V1TH contours
msg
cont V1TH z 1 fil;.0001
pause
cl
msg
msg V1/T contours
msg
cont V1/T z 1 fil;.0001
pause
cl
msg
msg DISX contours
msg
cont DISX z 1 fil;.0001
vec z 1 col 0
pause
cl
msg
msg U1TH contours
msg
cont U1TH z 1 fil;.0001
pause
cl
msg
msg U1/T contours
msg
cont U1/T z 1 fil;.0001
pause
cl
msg
msg STRX contours
msg
cont STRX z 1 fil;.0001
pause
cl
msg
msg SXTH contours
msg
cont SXTH z 1 fil;.0001
pause
cl
msg
msg SX/T contours
msg
cont SX/T z 1 fil;.0001
pause
cl
msg
msg STRY contours
msg
cont STRY z 1 fil;.0001
pause
cl
msg
msg SYTH contours
msg
cont SYTH z 1 fil;.0001
pause
cl
msg
msg SY/T contours
msg
cont SY/T z 1 fil;.0001
pause
cl
msg
msg STXY contours
msg
cont STXY z 1 fil;.0001
pause
cl
msg
msg XYTH contours
msg
cont XYTH z 1 fil;.0001
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Polar Plate with CIRCULAR hole )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 20
XULAST =1.5708
XFRAC(1)=0.05 ;XFRAC(2)=0.1
XFRAC(3)=0.15 ;XFRAC(4)=0.2
XFRAC(5)=0.25 ;XFRAC(6)=0.3
XFRAC(7)=0.35 ;XFRAC(8)=0.4
XFRAC(9)=0.45 ;XFRAC(10)=0.5
XFRAC(11)=0.55 ;XFRAC(12)=0.6
XFRAC(13)=0.65 ;XFRAC(14)=0.7
XFRAC(15)=0.75 ;XFRAC(16)=0.8
XFRAC(17)=0.85 ;XFRAC(18)=0.9
XFRAC(19)=0.95 ;XFRAC(20)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 30
YVLAST =1.5
RINNER =0.5 ;SNALFA =0.
YFRAC(1)=8.551264E-03 ;YFRAC(2)=0.022567
YFRAC(3)=0.039811 ;YFRAC(4)=0.059554
YFRAC(5)=0.081393 ;YFRAC(6)=0.105061
YFRAC(7)=0.130367 ;YFRAC(8)=0.157165
YFRAC(9)=0.18534 ;YFRAC(10)=0.214798
YFRAC(11)=0.24546 ;YFRAC(12)=0.277258
YFRAC(13)=0.310135 ;YFRAC(14)=0.34404
YFRAC(15)=0.378929 ;YFRAC(16)=0.414761
YFRAC(17)=0.451501 ;YFRAC(18)=0.489116
YFRAC(19)=0.527576 ;YFRAC(20)=0.566855
YFRAC(21)=0.606928 ;YFRAC(22)=0.647772
YFRAC(23)=0.689365 ;YFRAC(24)=0.731688
YFRAC(25)=0.774723 ;YFRAC(26)=0.818452
YFRAC(27)=0.862858 ;YFRAC(28)=0.907928
YFRAC(29)=0.953647 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(122)=EZ/T ;NAME(123)=EY/T
NAME(124)=EX/T ;NAME(125)=V1/T
NAME(126)=U1/T ;NAME(127)=V1TH
NAME(128)=U1TH ;NAME(129)=XY/T
NAME(130)=SY/T ;NAME(131)=SX/T
NAME(132)=XYTH ;NAME(133)=SYTH
NAME(134)=SXTH ;NAME(135)=STXY
NAME(136)=STRY ;NAME(137)=STRX
NAME(138)=EZTH ;NAME(139)=EYTH
NAME(140)=EXTH ;NAME(141)=EPSZ
NAME(142)=EPSY ;NAME(143)=EPSX
NAME(144)=DRH1 ;NAME(145)=DVO1
NAME(146)=ENUL ;NAME(147)=DEN1
NAME(148)=PRPS ;NAME(149)=DISY
NAME(150)=DISX
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(EZ/T,Y,N,N,N,N,Y)
SOLUTN(EY/T,Y,N,N,N,N,Y)
SOLUTN(EX/T,Y,N,N,N,N,Y)
SOLUTN(V1/T,Y,N,N,N,N,Y)
SOLUTN(U1/T,Y,N,N,N,N,Y)
SOLUTN(V1TH,Y,N,N,N,N,Y)
SOLUTN(U1TH,Y,N,N,N,N,Y)
SOLUTN(XY/T,Y,N,N,N,N,Y)
SOLUTN(SY/T,Y,N,N,N,N,Y)
SOLUTN(SX/T,Y,N,N,N,N,Y)
SOLUTN(XYTH,Y,N,N,N,N,Y)
SOLUTN(SYTH,Y,N,N,N,N,Y)
SOLUTN(SXTH,Y,N,N,N,N,Y)
SOLUTN(STXY,Y,N,N,N,N,Y)
SOLUTN(STRY,Y,N,N,N,N,Y)
SOLUTN(STRX,Y,N,N,N,N,Y)
SOLUTN(EZTH,Y,N,N,N,N,Y)
SOLUTN(EYTH,Y,N,N,N,N,Y)
SOLUTN(EXTH,Y,N,N,N,N,Y)
SOLUTN(EPSZ,Y,N,N,N,N,Y)
SOLUTN(EPSY,Y,N,N,N,N,Y)
SOLUTN(EPSX,Y,N,N,N,N,Y)
SOLUTN(DRH1,Y,N,N,N,N,Y)
SOLUTN(DVO1,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(DISY,Y,Y,N,N,N,Y)
SOLUTN(DISX,Y,Y,N,N,N,Y)
DEN1 = 147
VISL = 146
PRPS = 148
************************************************************
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(DISY,N,Y,Y,Y,N,Y)
TERMS(DISX,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 =7800. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =1.0E-05 ;DRH1DP =5.0E-12
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =0.3 ;ENUT =0.
PRNDTL(DISY)=1. ;PRNDTL(DISX)=1.
PRT(DISY)=1. ;PRT(DISX)=1.
CP1 =473. ;CP2 =1.
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=1
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(EZ/T)=1.001E-10 ;FIINIT(EY/T)=1.001E-10
FIINIT(EX/T)=1.001E-10 ;FIINIT(V1/T)=1.001E-10
FIINIT(U1/T)=1.001E-10 ;FIINIT(V1TH)=1.001E-10
FIINIT(U1TH)=1.001E-10 ;FIINIT(XY/T)=1.001E-10
FIINIT(SY/T)=1.001E-10 ;FIINIT(SX/T)=1.001E-10
FIINIT(XYTH)=1.001E-10 ;FIINIT(SYTH)=1.001E-10
FIINIT(SXTH)=1.001E-10 ;FIINIT(STXY)=1.001E-10
FIINIT(STRY)=1.001E-10 ;FIINIT(STRX)=1.001E-10
FIINIT(EZTH)=1.001E-10 ;FIINIT(EYTH)=1.001E-10
FIINIT(EXTH)=1.001E-10 ;FIINIT(EPSZ)=1.001E-10
FIINIT(EPSY)=1.001E-10 ;FIINIT(EPSX)=1.001E-10
FIINIT(DRH1)=1.001E-10 ;FIINIT(DVO1)=1.001E-10
FIINIT(ENUL)=1.001E-10 ;FIINIT(DEN1)=1.001E-10
FIINIT(PRPS)=160. ;FIINIT(DISY)=1.001E-10
FIINIT(DISX)=1.001E-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(OUTV ,NORTH , 1, 20, 30, 30, 1, 1, 1, 1)
COVAL(OUTV ,DISY,In-Form:source - see Grp 19)
PATCH(LOWX ,WWALL , 1, 1, 1, 30, 1, 1, 1, 1)
COVAL(LOWX ,DISX,1. ,0. )
PATCH(HIGHX ,EWALL , 20, 20, 1, 30, 1, 1, 1, 1)
COVAL(HIGHX ,DISX,1. ,0. )
PATCH(OUTU ,NORTH , 1, 20, 30, 30, 1, 1, 1, 1)
COVAL(OUTU ,DISX,In-Form:source - see Grp 19)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 400 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(DISY)=20 ;LITER(DISX)=20
ENDIT(DISY)=1.0E-03 ;ENDIT(DISX)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(EZ/T,LINRLX,1.)
RELAX(EY/T,LINRLX,1.)
RELAX(EX/T,LINRLX,1.)
RELAX(V1/T,LINRLX,1.)
RELAX(U1/T,LINRLX,1.)
RELAX(V1TH,LINRLX,1.)
RELAX(U1TH,LINRLX,1.)
RELAX(XY/T,LINRLX,1.)
RELAX(SY/T,LINRLX,1.)
RELAX(SX/T,LINRLX,1.)
RELAX(XYTH,LINRLX,1.)
RELAX(SYTH,LINRLX,1.)
RELAX(SXTH,LINRLX,1.)
RELAX(STXY,LINRLX,1.)
RELAX(STRY,LINRLX,1.)
RELAX(STRX,LINRLX,1.)
RELAX(EZTH,LINRLX,1.)
RELAX(EYTH,LINRLX,1.)
RELAX(EXTH,LINRLX,1.)
RELAX(EPSZ,LINRLX,1.)
RELAX(EPSY,LINRLX,1.)
RELAX(EPSX,LINRLX,1.)
RELAX(DRH1,LINRLX,1.)
RELAX(DVO1,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(DEN1,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
RELAX(DISY,FALSDT,1.0E+09)
RELAX(DISX,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(EZ/T)=1.0E+10 ;VARMIN(EZ/T)=-1.0E+10
VARMAX(EY/T)=1.0E+10 ;VARMIN(EY/T)=-1.0E+10
VARMAX(EX/T)=1.0E+10 ;VARMIN(EX/T)=-1.0E+10
VARMAX(V1/T)=1.0E+10 ;VARMIN(V1/T)=-1.0E+10
VARMAX(U1/T)=1.0E+10 ;VARMIN(U1/T)=-1.0E+10
VARMAX(V1TH)=1.0E+10 ;VARMIN(V1TH)=-1.0E+10
VARMAX(U1TH)=1.0E+10 ;VARMIN(U1TH)=-1.0E+10
VARMAX(XY/T)=1.0E+10 ;VARMIN(XY/T)=-1.0E+10
VARMAX(SY/T)=1.0E+10 ;VARMIN(SY/T)=-1.0E+10
VARMAX(SX/T)=1.0E+10 ;VARMIN(SX/T)=-1.0E+10
VARMAX(XYTH)=1.0E+10 ;VARMIN(XYTH)=-1.0E+10
VARMAX(SYTH)=1.0E+10 ;VARMIN(SYTH)=-1.0E+10
VARMAX(SXTH)=1.0E+10 ;VARMIN(SXTH)=-1.0E+10
VARMAX(STXY)=1.0E+10 ;VARMIN(STXY)=-1.0E+10
VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
VARMAX(EZTH)=1.0E+10 ;VARMIN(EZTH)=-1.0E+10
VARMAX(EYTH)=1.0E+10 ;VARMIN(EYTH)=-1.0E+10
VARMAX(EXTH)=1.0E+10 ;VARMIN(EXTH)=-1.0E+10
VARMAX(EPSZ)=1.0E+10 ;VARMIN(EPSZ)=-1.0E+10
VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
VARMAX(DRH1)=1.0E+10 ;VARMIN(DRH1)=-1.0E+10
VARMAX(DVO1)=1.0E+10 ;VARMIN(DVO1)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(DISY)=1.0E+10 ;VARMIN(DISY)=-1.0E+10
VARMAX(DISX)=1.0E+10 ;VARMIN(DISX)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
STRA = T
PARSOL = F
CONWIZ = T
ISG21 = 400
ISG52 = 3
ISG62 = 1
CSG10 ='q1'
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,160)
SPEDAT(SET,DOMAIN,PROPS_FILE_1,C,Q1)
SPEDAT(SET,STORED,SXTH,C,=5000.*((1+(0.5/RG)^2)-(1+3*(0.5/RG)^4)*$)
SPEDAT(SET,STORED,SXTH,C,COS(2*XG))!ZSLFIN)
SPEDAT(SET,STORED,SX/T,C,=STRX/SXTH!ZSLFIN)
SPEDAT(SET,STORED,SYTH,C,=5000.*((1-(0.5/RG)^2)+(1+3*(0.5/RG)^4-4$)
SPEDAT(SET,STORED,SYTH,C,*(0.5/RG)^2)*COS(2*XG))!ZSLFIN)
SPEDAT(SET,STORED,SY/T,C,=STRY/SYTH!ZSLFIN)
SPEDAT(SET,STORED,XYTH,C,=-5000.*(1-3*(0.5/RG)^4+2*(0.5/RG)^2)*SI$)
SPEDAT(SET,STORED,XYTH,C,N(2*XG)!ZSLFIN)
SPEDAT(SET,STORED,XY/T,C,=STXY/XYTH!ZSLFIN)
SPEDAT(SET,STORED,U1TH,C,=-2.5E-08*RG*(1.3*(1+(0.5/RG)^4)+2*(0.5/$)
SPEDAT(SET,STORED,U1TH,C,RG)^2*0.7)*SIN(2*XG)!ZSLFIN)
SPEDAT(SET,STORED,U1/T,C,=DISX/U1TH!ZSLFIN)
SPEDAT(SET,STORED,V1TH,C,=2.5E-08*RV*(0.7+(0.5/RG)^2*1.3+((1-(0.5$)
SPEDAT(SET,STORED,V1TH,C,/RG)^4)*1.3+4*(0.5/RG)^2)*COS(2*XG))!ZSL$)
SPEDAT(SET,STORED,V1TH,C,FIN)
SPEDAT(SET,STORED,V1/T,C,=DISY/V1TH!ZSLFIN)
SPEDAT(SET,STORED,EXTH,C,=2.5E-08*(0.7+(0.5/RG)^2*1.3-(1.3*(1+3*($)
SPEDAT(SET,STORED,EXTH,C,0.5/RG)^4)-4*(0.5/RG)^2*0.3)*COS(2*XG))!$)
SPEDAT(SET,STORED,EXTH,C,ZSLFIN)
SPEDAT(SET,STORED,EX/T,C,=EPSX/EXTH!ZSLFIN)
SPEDAT(SET,STORED,EYTH,C,=2.5E-08*(0.7-(0.5/RG)^2*1.3+(1.3*(1+3*($)
SPEDAT(SET,STORED,EYTH,C,0.5/RG)^4)-4*(0.5/RG)^2)*COS(2*XG))!ZSLF$)
SPEDAT(SET,STORED,EYTH,C,IN)
SPEDAT(SET,STORED,EY/T,C,=EPSY/EYTH!ZSLFIN)
SPEDAT(SET,STORED,EZTH,C,=-7.5E-09*(2-4*(0.5/RG)^2*COS(2*XG))!ZSL$)
SPEDAT(SET,STORED,EZTH,C,FIN)
SPEDAT(SET,STORED,EZ/T,C,=EPSZ/EZTH!ZSLFIN)
SPEDAT(SET,SOURCE,DISY!OUTV,C,=COVAL(2.0E-10&1.0/ANORTH*(5000.*(($)
SPEDAT(SET,SOURCE,DISY!OUTV,C,1-(0.5/RG)^2)+(1+3*(0.5/RG)^4-4*(0.$)
SPEDAT(SET,SOURCE,DISY!OUTV,C,5/RG)^2)*COS(2*XG)))))
SPEDAT(SET,SOURCE,DISX!OUTU,C,=COVAL(2.0E-10&-5000.*(1-3*(0.5/RG)$)
SPEDAT(SET,SOURCE,DISX!OUTU,C,^4+2*(0.5/RG)^2)*SIN(2*XG)))
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OUTV,C,OUTV)
SPEDAT(SET,OBJTYP,!OUTV,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!LOWX,C,LOWX)
SPEDAT(SET,OBJTYP,!LOWX,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!HIGHX,C,HIGHX)
SPEDAT(SET,OBJTYP,!HIGHX,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!OUTU,C,OUTV)
SPEDAT(SET,OBJTYP,!OUTU,C,USER_DEFINED)
SPEDAT(SET,FACETDAT,NUMOBJ,I,3)
SPEDAT(SET,MATERIAL,160,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
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(EZ/T,Y,N,Y,N,N,N)
OUTPUT(EY/T,Y,N,Y,N,N,N)
OUTPUT(EX/T,Y,N,Y,N,N,N)
OUTPUT(V1/T,Y,N,Y,N,N,N)
OUTPUT(U1/T,Y,N,Y,N,N,N)
OUTPUT(V1TH,Y,N,Y,N,N,N)
OUTPUT(U1TH,Y,N,Y,N,N,N)
OUTPUT(XY/T,Y,N,Y,N,N,N)
OUTPUT(SY/T,Y,N,Y,N,N,N)
OUTPUT(SX/T,Y,N,Y,N,N,N)
OUTPUT(XYTH,Y,N,Y,N,N,N)
OUTPUT(SYTH,Y,N,Y,N,N,N)
OUTPUT(SXTH,Y,N,Y,N,N,N)
OUTPUT(STXY,Y,N,Y,N,N,N)
OUTPUT(STRY,Y,N,Y,N,N,N)
OUTPUT(STRX,Y,N,Y,N,N,N)
OUTPUT(EZTH,Y,N,Y,N,N,N)
OUTPUT(EYTH,Y,N,Y,N,N,N)
OUTPUT(EXTH,Y,N,Y,N,N,N)
OUTPUT(EPSZ,Y,N,Y,N,N,N)
OUTPUT(EPSY,Y,N,Y,N,N,N)
OUTPUT(EPSX,Y,N,Y,N,N,N)
OUTPUT(DRH1,N,N,N,N,N,N)
OUTPUT(DVO1,N,N,N,N,N,N)
OUTPUT(ENUL,N,N,N,N,N,N)
OUTPUT(DEN1,N,N,N,N,N,N)
OUTPUT(PRPS,N,N,N,N,N,N)
OUTPUT(DISY,Y,N,Y,Y,Y,Y)
OUTPUT(DISX,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;IYMON = 29 ;IZMON = 1
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
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
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP