TALK=T;RUN(1,1)
  PHOTON USE
  p



  gr ou z 1
  con thex z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Linear thermal expansion coefficient
  pause
  con cl;red
  con stif z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Stiffnes ( Young's module ) contours
  pause
  con cl;red
  vec z 1 sh
  msg Velocity and displacement vectors
  pause
  vec cl;red
  con u1 z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg X-component displacement contours
  pause
  con cl;red
  con v1 z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Y-component displacement contours
  pause
  con cl;red
  con dil z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Dilatation contours
  pause
  con cl;red
  con epsx z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg X-strain contours
  pause
  con cl;red
  con epsy z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Y-strain contours
  pause
  con cl;red
  con epst z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Thermal strain contours
  pause
  con cl;red
  con strx z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg X-normal stress contours
  pause
  con cl;red
  con stry z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Y-normal stress contours
  pause
  con cl;red
  con taxy z 1 fil;.001
  gr ou z 1 x 1 3 y 6 m
  gr ou z 1 x 8 m y 6 m
  msg Shear stress contours
  pause
  msg Press E to end
  enduse
  DISPLAY
    This case, deals     with  the calculations of
    the displacements and stresses
    in uniformly heated  block  composed  the  materials  of
    different mechanical properties.

    It is  used  here  as an example of extensive processing
    the  calculation  results  with  the  aim  to  get   the
    distribution   of  normal  and  shear  stresses  by  the
    numerical differentiation of displacement field.

    The outcome is the stress fields which may be viewed  by
    PHOTON. The use file for the latter is also supplied.
  ENDDIS
      PLANTBEGIN
PATCH(STELC,INIVAL,1,NX,1,NY,1,1,1,1)
    VAL=(1.-RG(4))/((1.+RG(4))*(1-2.*RG(4)))
INIT (STELC,CON1,0.0,GRND)
    VAL=RG(4)/((1.+RG(4))*(1-2.*RG(4)))
INIT (STELC,CON2,0.0,GRND)
    VAL=1./(1.-2.*RG(4))
INIT (STELC,CON3,0.0,GRND)
    VAL=1./(2.*(1+RG(4)))
INIT (STELC,CON4,0.0,GRND)
              Stress - in -Solid BC
    ---------------- WEST ----------------------------
   ** Zero displacement at the west boundary
PATCH(WESZD,WWALL,1,1,1,NY,1,1,1,1)
COVAL(WESZD,U1,1.,0.0)
    ---------------- EAST ----------------------------
   ** Normal stress at the east solid boundary
PATCH(EASZN,EAST,NX-1,NX-1,1,NY,1,1,1,1)
    VAL=U1[-1,,]+0.333*(P1+4.*(THEX*TEM+RG(3)/STIF))*DXU2D
COVAL(EASZN,U1,FIXVAL,GRND)
   ** Zero shear stress at the east boundary
PATCH(EASZS,EAST,NX,NX,1,NY-1,1,1,1,1)
    VAL=-(U1[-1,+1,]-U1[-1,,])/DYG2D
COVAL(EASZS,V1,FIXFLU,GRND)
   ** Dilatation at the east solid boundary
PATCH(FXDEAS,EAST,NX,NX,2,NY-1,1,1,1,1)
    VAL=2.*(THEX*TEM-RG(3)/STIF)-1.5*(V1-V1[,-1,])/DYV2D
COVAL(FXDEAS,P1,FIXVAL, grnd)

PATCH(FXDEAC,EAST,NX,NX,1,1,1,1,1,1)
    VAL=2.*(THEX*TEM-RG(3)/STIF)-            1.5*(V1[,+1,]-$
V1)/DYV2D[,+1,]
COVAL(FXDEAC,P1,FIXVAL, grnd)

    ---------------- NORTH ----------------------------
   ** Normal stress at the free north boundary
PATCH(NORZN,CELL,1,NX,NY-1,NY-1,1,1,1,1)
    VAL=V1[,-1,]+0.333*(P1+4.*(THEX*TEM+RG(2)/STIF))*DYV2D
COVAL(NORZN,V1,FIXVAL,grnd)
   ** Zero shear stress at the free north boundary
PATCH(NORZS,NORTH,1,NX-1,NY,NY,1,1,1,1)
    VAL=-(V1[+1,-1,]-V1[,-1,])/DXG2D
COVAL(NORZS,U1,FIXFLU,GRND)
   ** Dilatation at the free north boundary
PATCH(FDNW,NORTH,1,1,NY,NY,1,1,1,1)
    VAL=2.*(THEX*TEM-RG(2)/STIF)-            1.5*(U1[+1,,]-$
U1)/DXU2D[+1,,]
COVAL(FDNW,P1,FIXVAL, grnd)
PATCH(FDNOR,NORTH,2,NX-1,NY,NY,1,1,1,1)
    VAL=2.*(THEX*TEM-RG(2)/STIF)-1.5*(U1-U1[-1,,])/DXU2D
COVAL(FDNOR,P1,FIXVAL, grnd)
PATCH(FDNEC,NORTH,NX,NX,NY,NY,1,1,1,1)
    VAL=2.*(THEX*TEM-RG(2)/STIF)-       1.5*(U1[-1,,]-U1[-2$
,,])/DXU2D[-1,,]
COVAL(FDNEC,P1,FIXVAL, grnd)
    ---------------- SOUTH ----------------------------
   ** Zero displacement at the south symmetry
PATCH(SOUZD,SWALL,1,NX,1,1,1,1,1,1)
COVAL(SOUZD,V1,1.,0.0)
    ==================================================
   ** Volumetric mass source for dilatation
PATCH(DILAT,VOLUME,1,nx,1,ny,1,1,1,1)
    VAL=-P1/CON3+THEX*TEM/CON4
COVAL(DILAT,P1,FIXFLU,GRND)
    ===================================================
            OUTPUT INFORMATION
   ** Dilatation
STORE(DIL)
    DIL=-P1/CON3+THEX*TEM/CON4
   ** Stress and strain calculations
STORE(STRX,STRY,TAXY)
STORE(EPSX,EPSY,EPST)
      * Thermal strains
   EPST=THEX*TEM
  REGION(1,NX,1,NY,1,1)
      * Normal strains
   EPSX=(U1-WEST(U1))/DXU2D
  REGION(2,NX-1,1,NY,1,1)
   EPSX=(U1[+1,,]-U1)/DXU2D[+1,,]
  REGION(1,1,1,NY,1,1)
   EPSX=(U1[-1,,]-U1[-2,,])/DXU2D[-1,,]
  REGION(NX,NX,1,NY,1,1)
   EPSY=(V1-SOUTH(V1))/DYV2D
  REGION(1,NX,2,NY-1,1,1)
   EPSY=(V1[,+1,]-V1)/DYV2D[,+1,]
  REGION(1,NX,1,1,1,1)
   EPSY=(V1[,-1,]-V1[,-2,])/DYV2D[,-1,]
  REGION(1,NX,NY,NY,1,1)
      * Normal stresses
   STRX=STIF*(CON1*EPSX+CON2*EPSY-CON3*EPST)
  REGION(1,NX,1,NY,1,1)
   STRY=STIF*(CON1*EPSY+CON2*EPSX-CON3*EPST)
  REGION(1,NX,1,NY,1,1)
      * Shear stresses
STORE(DDYU,DDXV)
              == ddy.u ==
   DDYU=0.25*((U1-U1[,-1,])/DYG2D[,-1,]+$
 (U1[,+1,]-U1)/DYG2D+(U1[-1,,]-                         U1[-1,-1,])$
/DYG2D[,-1,]+                    (U1[-1,+1,]-U1[-1,,])/DYG2D)
  REGION(2,NX-1,2,NY-1,1,1)

   DDYU=0.5*((U1-U1[,-1,])/DYG2D[,-1,]+                (U1[$
,+1,]-U1)/DYG2D)
  REGION(1,1,2,NY-1,1,1)

   DDYU=(U1[,+1,]-U1)/DYG2D
  REGION(1,1,1,1,1,1)

   DDYU=(U1-U1[,-1,])/DYG2D[,-1,]
  REGION(1,1,NY,NY,1,1)

   DDYU=0.5*((U1[-1,,]-U1[-1,-1,])/DYG2D[,-1,]+$
              (U1[-1,+1,]-U1[-1,,])/DYG2D)
  REGION(NX,NX,2,NY-1,1,1)

   DDYU=(U1[-1,+1,]-U1[-1,,])/DYG2D
  REGION(NX,NX,1,1,1,1)

   DDYU=(U1[-1,,]-U1[-1,-1,])/DYG2D[,-1,]
  REGION(NX,NX,NY,NY,1,1)

   DDYU=0.5*((U1-U1[,-1,])/DYG2D[,-1,]+             (U1[-1,$
,]-U1[-1,-1,])/DYG2D[,-1,])
  REGION(2,NX-1,NY,NY,1,1)

   DDYU=0.5*((U1[,+1,]-U1)/DYG2D+             (U1[-1,+1,]-U$
1[-1,,])/DYG2D)
  REGION(2,NX-1,1,1,1,1)

              == ddx.v ==
   DDXV=0.25*((V1[+1,-1,]-V1[,-1,]  )/DXG2D+$
       (V1[+1,,]-V1)/DXG2D+(V1[,-1,]-                            V1$
[-1,-1,])/DXG2D[-1,,]+                          (V1-V1[-1,,])/DXG2D$
[-1,,])
  REGION(2,NX-1,2,NY-1,1,1)

   DDXV=(V1[+1,,] -V1)/DXG2D
  REGION(1,1,1,1,1,1)

   DDXV=0.5*((V1[+1,-1,]-V1[,-1,])/DXG2D+$
         (V1[+1,,]-V1)/DXG2D)
  REGION(1,1,2,NY-1,1,1)

   DDXV=(V1[+1,-1,]-V1[,-1,]  )/DXG2D
  REGION(1,1,NY,NY,1,1)

   DDXV=(V1-V1[-1,,])/DXG2D[-1,,]
  REGION(NX,NX,1,1,1,1)

   DDXV=0.5*((V1[,-1,]-V1[-1,-1,])/DXG2D[-1,,]+$
                (V1-V1[-1,,])/DXG2D[-1,,])
  REGION(NX,NX,2,NY-1,1,1)

   DDXV=(V1[,-1,] -V1[-1,-1,])/DXG2D[-1,,]
  REGION(NX,NX,NY,NY,1,1)

   DDXV=0.5*((V1[+1,,]-V1)/DXG2D+               (V1-V1[-1,,$
])/DXG2D[-1,,])
  REGION(2,NX-1,1,1,1,1)

   DDXV=0.5*((V1[+1,-1,]-V1[,-1,])/DXG2D+               (V1$
[,-1,]-V1[-1,-1,])/DXG2D[-1,,])
  REGION(2,NX-1,NY,NY,1,1)

   TAXY=STIF*CON4*(DDYU+DDXV)
  REGION(1,NX,1,NY,1,1)
    PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2D Stress in solid, XY plane            )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 611
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 10
 XULAST =1.
 XFRAC(1)=0.1 ;XFRAC(2)=0.2
 XFRAC(3)=0.3 ;XFRAC(4)=0.4
 XFRAC(5)=0.5 ;XFRAC(6)=0.6
 XFRAC(7)=0.7 ;XFRAC(8)=0.8
 XFRAC(9)=0.9 ;XFRAC(10)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 10
 YVLAST =1.
 YFRAC(1)=0.1 ;YFRAC(2)=0.2
 YFRAC(3)=0.3 ;YFRAC(4)=0.4
 YFRAC(5)=0.5 ;YFRAC(6)=0.6
 YFRAC(7)=0.7 ;YFRAC(8)=0.8
 YFRAC(9)=0.9 ;YFRAC(10)=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(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(134)=DDXV
 NAME(135)=DDYU ;NAME(136)=EPST
 NAME(137)=EPSY ;NAME(138)=EPSX
 NAME(139)=TAXY ;NAME(140)=STRY
 NAME(141)=STRX ;NAME(142)=DIL
 NAME(143)=CON4 ;NAME(144)=CON3
 NAME(145)=CON2 ;NAME(146)=CON1
 NAME(147)=POIS ;NAME(148)=STIF
 NAME(149)=THEX ;NAME(150)=TEM
    * 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(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(DDXV,Y,N,N,N,N,Y)
 SOLUTN(DDYU,Y,N,N,N,N,Y)
 SOLUTN(EPST,Y,N,N,N,N,Y)
 SOLUTN(EPSY,Y,N,N,N,N,Y)
 SOLUTN(EPSX,Y,N,N,N,N,Y)
 SOLUTN(TAXY,Y,N,N,N,N,Y)
 SOLUTN(STRY,Y,N,N,N,N,Y)
 SOLUTN(STRX,Y,N,N,N,N,Y)
 SOLUTN(DIL,Y,N,N,N,N,Y)
 SOLUTN(CON4,Y,N,N,N,N,Y)
 SOLUTN(CON3,Y,N,N,N,N,Y)
 SOLUTN(CON2,Y,N,N,N,N,Y)
 SOLUTN(CON1,Y,N,N,N,N,Y)
 SOLUTN(POIS,Y,N,N,N,N,Y)
 SOLUTN(STIF,Y,N,N,N,N,Y)
 SOLUTN(THEX,Y,N,N,N,N,Y)
 SOLUTN(TEM,Y,Y,N,N,N,Y)
 ************************************************************
  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,N,Y,Y,Y,Y)
 TERMS(V1,Y,N,Y,Y,Y,Y)
 TERMS(TEM,N,N,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 =1. ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(TEM)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(TEM)=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(DDXV)=1.0E-10
 FIINIT(DDYU)=1.0E-10 ;FIINIT(EPST)=1.0E-10
 FIINIT(EPSY)=1.0E-10 ;FIINIT(EPSX)=1.0E-10
 FIINIT(TAXY)=1.0E-10 ;FIINIT(STRY)=1.0E-10
 FIINIT(STRX)=1.0E-10 ;FIINIT(DIL)=1.0E-10
 FIINIT(CON4)=1.0E-10 ;FIINIT(CON3)=1.0E-10
 FIINIT(CON2)=1.0E-10 ;FIINIT(CON1)=1.0E-10
 FIINIT(POIS)=1.0E-10 ;FIINIT(STIF)=1.0E-10
 FIINIT(THEX)=1.0E-10 ;FIINIT(TEM)=100.
 
 PATCH(STELP ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
 INIT(STELP ,POIS,0. ,0.3333 )
 
 PATCH(STELT ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
 INIT(STELT ,THEX,0. ,1.0E-05 )
 
 PATCH(STELS ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
 INIT(STELS ,STIF,0. ,2.0E+11 )
 
 PATCH(ALMT1 ,INIVAL, 1, 3, 6, 10, 1, 1, 1, 1)
 INIT(ALMT1 ,STIF,0. , FIXVAL )
 INIT(ALMT1 ,THEX,0. ,1.0E-04 )
 
 PATCH(ALMT2 ,INIVAL, 8, 10, 6, 10, 1, 1, 1, 1)
 INIT(ALMT2 ,STIF,0. , FIXVAL )
 INIT(ALMT2 ,THEX,0. ,1.0E-04 )
 
 PATCH(STELC ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
 INIT(STELC ,CON4,0. , GRND )
 INIT(STELC ,CON3,0. , GRND )
 INIT(STELC ,CON2,0. , GRND )
 INIT(STELC ,CON1,0. , GRND )
 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(WESZD ,WWALL , 1, 1, 1, 10, 1, 1, 1, 1)
 COVAL(WESZD ,U1 ,1. ,0. )
 
 PATCH(EASZN ,EAST , 9, 9, 1, 10, 1, 1, 1, 1)
 COVAL(EASZN ,U1 , FIXVAL , GRND )
 
 PATCH(EASZS ,EAST , 10, 10, 1, 9, 1, 1, 1, 1)
 COVAL(EASZS ,V1 , FIXFLU , GRND )
 
 PATCH(FXDEAS ,EAST , 10, 10, 2, 9, 1, 1, 1, 1)
 COVAL(FXDEAS ,P1 , FIXVAL , GRND )
 
 PATCH(FXDEAC ,EAST , 10, 10, 1, 1, 1, 1, 1, 1)
 COVAL(FXDEAC ,P1 , FIXVAL , GRND )
 
 PATCH(NORZN ,CELL , 1, 10, 9, 9, 1, 1, 1, 1)
 COVAL(NORZN ,V1 , FIXVAL , GRND )
 
 PATCH(NORZS ,NORTH , 1, 9, 10, 10, 1, 1, 1, 1)
 COVAL(NORZS ,U1 , FIXFLU , GRND )
 
 PATCH(FDNW ,NORTH , 1, 1, 10, 10, 1, 1, 1, 1)
 COVAL(FDNW ,P1 , FIXVAL , GRND )
 
 PATCH(FDNOR ,NORTH , 2, 9, 10, 10, 1, 1, 1, 1)
 COVAL(FDNOR ,P1 , FIXVAL , GRND )
 
 PATCH(FDNEC ,NORTH , 10, 10, 10, 10, 1, 1, 1, 1)
 COVAL(FDNEC ,P1 , FIXVAL , GRND )
 
 PATCH(SOUZD ,SWALL , 1, 10, 1, 1, 1, 1, 1, 1)
 COVAL(SOUZD ,V1 ,1. ,0. )
 
 PATCH(DILAT ,VOLUME, 1, 10, 1, 10, 1, 1, 1, 1)
 COVAL(DILAT ,P1 , FIXFLU , GRND )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 3000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = F
 RESREF(P1)=1.0E-08 ;RESREF(U1)=1.0E-08
 RESREF(V1)=1.0E-08 ;RESREF(TEM)=1.0E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(TEM)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(TEM)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,5.0E-03)
 RELAX(U1,FALSDT,5.0E-04)
 RELAX(V1,FALSDT,5.0E-04)
 RELAX(DDXV,LINRLX,1.)
 RELAX(DDYU,LINRLX,1.)
 RELAX(EPST,LINRLX,1.)
 RELAX(EPSY,LINRLX,1.)
 RELAX(EPSX,LINRLX,1.)
 RELAX(TAXY,LINRLX,1.)
 RELAX(STRY,LINRLX,1.)
 RELAX(STRX,LINRLX,1.)
 RELAX(DIL,LINRLX,1.)
 RELAX(CON4,LINRLX,1.)
 RELAX(CON3,LINRLX,1.)
 RELAX(CON2,LINRLX,1.)
 RELAX(CON1,LINRLX,1.)
 RELAX(POIS,LINRLX,1.)
 RELAX(STIF,LINRLX,1.)
 RELAX(THEX,LINRLX,1.)
 RELAX(TEM,FALSDT,1.0E+09)
 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(DDXV)=1.0E+10 ;VARMIN(DDXV)=-1.0E+10
 VARMAX(DDYU)=1.0E+10 ;VARMIN(DDYU)=-1.0E+10
 VARMAX(EPST)=1.0E+10 ;VARMIN(EPST)=-1.0E+10
 VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
 VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
 VARMAX(TAXY)=1.0E+10 ;VARMIN(TAXY)=-1.0E+10
 VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
 VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
 VARMAX(DIL)=1.0E+10 ;VARMIN(DIL)=-1.0E+10
 VARMAX(CON4)=1.0E+10 ;VARMIN(CON4)=-1.0E+10
 VARMAX(CON3)=1.0E+10 ;VARMIN(CON3)=-1.0E+10
 VARMAX(CON2)=1.0E+10 ;VARMIN(CON2)=-1.0E+10
 VARMAX(CON1)=1.0E+10 ;VARMIN(CON1)=-1.0E+10
 VARMAX(POIS)=1.0E+10 ;VARMIN(POIS)=-1.0E+10
 VARMAX(STIF)=1.0E+10 ;VARMIN(STIF)=-1.0E+10
 VARMAX(THEX)=1.0E+10 ;VARMIN(THEX)=-1.0E+10
 VARMAX(TEM)=1.0E+10 ;VARMIN(TEM)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 POISSN =0.3333
 RG( 1) =1.0E-05
 RG( 2) =-1.0E+08
 RG( 4) =0.3333
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=4.671E-03 ;EX(U1)=3.079E-03
 EX(V1)=6.59E-04 ;EX(DDXV)=2.347E-03
 EX(DDYU)=3.255E-03 ;EX(EPST)=3.7E-03
 EX(EPSY)=3.955E-03 ;EX(EPSX)=6.382E-03
 EX(TAXY)=1.362E+08 ;EX(STRY)=3.664E+08
 EX(STRX)=1.83E+08 ;EX(DIL)=8.618E-03
 EX(CON4)=0.375 ;EX(CON3)=2.999
 EX(CON2)=0.7498 ;EX(CON1)=1.5
 EX(POIS)=0.3333 ;EX(STIF)=1.46E+11
 EX(THEX)=3.7E-05 ;EX(TEM)=100.
 ************************************************************
  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(DDXV,Y,N,Y,N,N,N)
 OUTPUT(DDYU,Y,N,Y,N,N,N)
 OUTPUT(EPST,Y,N,Y,N,N,N)
 OUTPUT(EPSY,Y,N,Y,N,N,N)
 OUTPUT(EPSX,Y,N,Y,N,N,N)
 OUTPUT(TAXY,Y,N,Y,N,N,N)
 OUTPUT(STRY,Y,N,Y,N,N,N)
 OUTPUT(STRX,Y,N,Y,N,N,N)
 OUTPUT(DIL,Y,N,Y,N,N,N)
 OUTPUT(CON4,Y,N,Y,N,N,N)
 OUTPUT(CON3,Y,N,Y,N,N,N)
 OUTPUT(CON2,Y,N,Y,N,N,N)
 OUTPUT(CON1,Y,N,Y,N,N,N)
 OUTPUT(POIS,Y,N,Y,N,N,N)
 OUTPUT(STIF,Y,N,Y,N,N,N)
 OUTPUT(THEX,Y,N,Y,N,N,N)
 OUTPUT(TEM,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 9 ;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