TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 26/09/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
  A plate consisting of two materials has its temperature raised.
  Becuse the thermal expansion coefficients of the two materials 
  differ by a factor of 10, the plate distorts. its radius of 
  curvature being proportional to the width of the plate in the 
  direction normal to the line of separation between the materials, 
  divided by the differences in thermal strain.
  
  The opportunity is provided for the plate to be held either at its 
  high-y or low-y edges.
  
  ENDDIS
  PHOTON USE
  p;;;;
  
  

  set prop off
  cl
  msg prps distribution
  gr ou z 1
  cont prps z 1 fil;.0001
  pause
  cl
  msg x-displacement distribution
  gr ou z 1
  cont U1 z 1 fil;.0001
  vec z 1  col 1
  pause
  cl
  msg y-displacement distribution
  gr ou z 1
  cont V1 z 1 fil;.0001
  vec z 1  col 1
  pause
  cl
  msg reciprocal of radius of curvature
  gr ou z 1
  cont curv z 1 fil;.0001
  vec z 1  col 1
  pause
  cl
  msg epsy distribution
  gr ou z 1
  cont epsy z 1 fil;.0001
  vec z 1  col 1
  dump epsy
  pause
  cl
  msg epsx distribution
  gr ou z 1
  cont epsx z 1 fil;.0001
  vec z 1  col 1
  dump epsx
  ENDUSE

 ************************************************************
 IRUNN   =         1 ;LIBREF =       210
 ************************************************************
  Group 1. Run Title
 TEXT(2D xy Plate of two materials & TEM; s210 )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,18,18,1,1.000000E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(PRPS,STRX,STRY,STRZ,STXY,EPSY,EPSX,EPSZ,EPST)
 STORE(DEN1,ENUL,SPH1,KOND,DRH1,TEM1,CURV)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 SOLUTN(U1  ,Y,Y,Y,N,N,Y)
 SOLUTN(V1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Echo InForm settings for Group  7
  SAVE7BEGIN
 
  ***** SELECT CASENO ******
 integer(caseno)
     caseno = 1 : z free
     caseno = 2 : z fixed

 caseno=2
  ***** CASENO ******
 (stored var curv is (epsy-epst-epsy[-1,,]+epst[-1,,])/dxg$
 with if(ix.gt.1))
 
  SAVE7END
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
    * Domain material index is 160 signifying:
    * NOTSET
 SETPRPS(1,160,Q1)
 ENUT    = 0.000000E+00
 DRH1DP  = 5.000000E-12
 DVO1DT  = 1.000000E-05
   *  List of user-defined materials to be read by EARTH
    MATFLG=T;IMAT=2
   *  Name
   *Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
   *  
    160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
   *  
    161 7800.0 0.3 473.0 43.0 1.0E-4 0.5E-11
    
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1  ) =  0.000000E+00 ;FIINIT(U1  ) =  0.000000E+00
 FIINIT(V1  ) =  0.000000E+00 ;FIINIT(PRPS) =  1.600000E+02
 FIINIT(TEM1) =  1.000000E+01
   No PATCHes used for this Group
 
 
 INIADD  =    F
 ************************************************************
  Echo InForm settings for Group 11
  SAVE11BEGIN
(INITIAL PRPS at 2ndmat is 161)
SPEDAT(SET,MATERIAL,161,L,T)
  SAVE11END
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (UP      ,NWALL ,1,0,0,0,0,0,1,1)
 COVAL (UP      ,V1  , 1.000000E+00, 0.000000E+00)
 
 PATCH (POINTU  ,CELL  ,3,0,0,0,0,0,1,1)
 COVAL (POINTU  ,U1  , FIXVAL      , 0.000000E+00)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      200
 RESREF(U1  ) =  0.000000E+00 ;RESREF(V1  ) =  0.000000E+00
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 LITER (P1  ) =       50 ;LITER (U1  ) =       50
 LITER (V1  ) =       50
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 ISG21   =      1800
 ISG50   =         1
 ISG52   =         1
 CSG10   ='q1'
 SPEDAT(SET,MATERIAL,161,L,T)
 SPEDAT(SET,BOUNDARY,ZCONST,R,1.00000E+20)
 ************************************************************
  Echo InForm settings for Group 19
  SAVE19BEGIN
IF(CASENO.EQ.2) THEN
SPEDAT(BOUNDARY,ZCONST,R,1.e20)       ! z contraction prevented
ELSE
SPEDAT(BOUNDARY,ZCONST,R,0.0)         ! z contraction allowed
ENDIF
  SAVE19END
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =         9 ;IYMON  =         2 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NYPRIN  =         1
 NXPRIN  =         1
 NYPRIN = 1
 NXPRIN = 1
 ixprf=7
 ixprl=11

 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-1.810503E-02,5.071213E-02,9.985492E-01)
 GVIEW(UP,-6.446701E-02,9.965755E-01,-5.178077E-02)
 
> DOM,    SIZE,        9.000000E-02, 1.200000E-01, 1.000000E-03
> DOM,    MONIT,       7.650001E-02, 6.000001E-03, 5.000000E-04
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
  > GRID,   RSET_X_1,     18, 1.000000E+00
  > GRID,   RSET_Y_1,     18, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        UP
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 0.000000E+00, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        2NDMAT
> OBJ,    POSITION,    4.500000E-02, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        4.500000E-02, 1.200000E-01, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        POINTU
> OBJ,    POSITION,    4.010000E-02, 1.134000E-01, 0.000000E+00
> OBJ,    SIZE,        4.800000E-03, 5.666000E-03, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        CELLTYPE
 ************************************************************
STOP