TALK=T;RUN(1,1)
  display
     This Q1 demonstrates the use of 5 operations used
     for the creation of BFC grids, namely:
     * defining points,
     * defining lines,
     * defining frames,
     * setting grid dimensions, and
     * matching a grid to a frame.

     Those operations are used here in order to:
     1. create a circle,
     2. create a square, and
     3. to match them to the two sides of a constant-I plane,
        (namely I2 and I1 respectively).
     The relevant PIL commands will be displayed, then processed.
     The first sets the grid type to BFC


     Please enter the circle radius (default 0.25m)
     Please enter the side-length of the square (default 0.5m)
     The input data to be used in this example are..

     The centre of the circle will be placed at:
     x = 0.0
     y = - square root of 2 times the radius
     z = + square root of 2 times the radius


     Now define points A,B,C and D, positioned as follows:
       A at (0,-root,root)
       B at (0,root,root)
       C at (0,root,-root)
       D at (0,-root,-root)
     by the commands:
     for A
     for B
     for C
     and for D
     Now draw four arcs (AB,BC,CD,DA) from A to D
     to create a circle.
     Arcs AB and CD must have equal numbers of intervals;
     so must arcs BC and DA.
     The intervals will be supposed uniform on each arc.

     Please enter nab, the number of cells for AB and CD (default 5)
     The command to create arc AB, with nab uniform intervals,
     which also goes through the point (0,0,radius),
     is:

     Please enter nbc, the number of cells for BC and DA (default 5)
     The command to create arc BC, with nbc cells
     which also goes through the point (0,0,radius),
     is:
     The corresponding command for CD is:
     and for DA:

     Now define a frame ABCD which has corner points
     A,B,C & D by the command:


     Now set the BFC grid dimension as 1 x nab x nbc, with
     reference length 1m x side x side, by the command:


     Now match grid plane I2 on frame ABCD using the
     trans-finite method, by the command:

     Now create the square on plane I = 1
     Define points E,F,G and H, with positions:
       E at (-side,-side/2,side/2)
       F at (-side,side/2,side/2)
       G at (-side,side/2,-side/2)
       H at (-side,-side/2,-side/2)
     by the commands:
     for E:
     for F:
     for G:
     and for H:
     Now draw four lines (EF, FG, GH, HE) from E to H
     to create a square.
     Lines EF & GH each have nab uniformly-distributed intervals,
     and lines FG & HE each have nbc uniformly-distributed
     intervals, so as to match the numbers of intervals on the
     circle.

     The commands are:
     then
     then
     and finally

     Now define a frame EFGH (square) which has corner points
     E, F, G and H, by the command:

     Now match grid plane I1 on frame EFGH using the trans-finite
     by the command:

     This completes the exercise.

     However, you may now wish to type VIEW, so as to be able to
     inspect, or indeed alter, the grid by means of the graphical
     grid generator.
  enddis
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(GSET COMMAND TUTORIAL 1            :B538)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 1
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 5
 YVLAST =0.5
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 5
 ZWLAST =0.5
 ************************************************************
  Group 6. Body-Fitted Coordinates
 BFC = T ;NONORT = F
 NCRT = -1
 RSTGEO = F ;SAVGEO = F
 UUP = F ;VUP = F ;WUP = F
 NGEOM =CHAM
 NAMXYZ =CHAM
 ANGMIN =20.
 DOMAIN(   1,   1,   1,   6,   1,   6)
    * Set fixed sub-domain
 FIXDOM(   1,   0,   0,   0,   0,   0,   0)
 FIXDOM(   2,   0,   0,   0,   0,   0,   0)
 FIXDOM(   3,   0,   0,   0,   0,   0,   0)
 FIXDOM(   4,   0,   0,   0,   0,   0,   0)
 FIXDOM(   5,   0,   0,   0,   0,   0,   0)
 FIXDOM(   6,   0,   0,   0,   0,   0,   0)
 FIXDOM(   7,   0,   0,   0,   0,   0,   0)
 FIXDOM(   8,   0,   0,   0,   0,   0,   0)
 FIXDOM(   9,   0,   0,   0,   0,   0,   0)
 FIXDOM(  10,   0,   0,   0,   0,   0,   0)
 LIJ = F ;LJK = T ;LIK = F
 IMON = 1 ;JMON = 3 ;KMON = 3
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 ************************************************************
  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         
 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.0E-05 ;ENUT =0.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
   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
   No PATCHes yet used for this Group
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 ************************************************************
  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
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 1 ;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
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP