TALK=T;RUN(1,1)
  DISPLAY


    This case  exemplifies  some  of  the  basic  techniques
    available in In-Form for field initialisations.

    Three examples are provided, each for a different xy-plane
    in a single 3D box.
    
    Example 1: Initialization of stagnation-point flow on IZ=1
               ===============================================

  PATCH(INI1,INIVAL,1,NX,1,NY,1,1,1,1)
  (INITIAL of U1 at INI1 is XU-10.)
  (INITIAL of V1 at INI1 is -(YV-10.))

    Example 2: Solid-body rotation flow on IZ=2
               ================================

  PATCH(INI2,INIVAL,1,NX,1,NY,2,2,1,1)
  (INITIAL of U1 at INI2 is YG-10.)
  (INITIAL of V1 at INI2 is -(XG-10.))

    Example 3:  The values on IZ=1 and IZ=2 are added together
                and placed on IZ=3
                ==============================================

  PATCH(INI3,INIVAL,1,NX,1,NY,3,3,1,1)
  (INITIAL of U1 at INI3 is -U1[&&-1]+U1[&&1])
  (INITIAL of V1 at INI3 is -V1[&&2]+V1[&&-2])

  The Q1 contains PHOTON USE commands
  ENDDIS
  PHOTON USE
  p



  msg
  msg
  msg   All examples which follow are assembled as different slabs
  msg   and number of slabs of one 3D domain called
  msg
  msg                INITBOX.
  msg
  msg   The series of PHOTON plots will follow.
  msg
  msg   They illustrate some of the existing initialisation
  msg   techniques available via In-Form.
  msg
  msg   The presentation covers:
  msg
  msg   * Flow field initialistions by formulae
  msg   * Manipulating with initial fields
  msg   * MARKing sub-domains to create arbitrary initial fields
  msg   * Geometry initialisations
  msg
  msg   Hit enter to proceed to 1st slab.
  msg
  pause
  clear
  vi z
  gr ou z 1
  set vec ref 25
  vec z 1 sh
  msg
  msg
  msg   The 1st slab of the box is initialised by stagnation
  msg   point flow
  msg   Hit enter to proceed to 2nd slab.
  pause
  vec off;red
  vec z 2 sh
  msg
  msg
  msg   The 2nd slab is initialised by solid body rotation flow
  msg
  msg   Hit enter to proceed to 3rd slab.
  pause
  vec off;red
  vec z 3 sh
  msg
  msg   The flow in the 3rd slab is constructed by adding the
  msg   velocity distributions of 1st and 2nd slabs.
  msg
  msg   Hit E to finish.
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Three examples of initialization :758   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 758
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 40
 XULAST =20.
 XFRAC(1)=0.025 ;XFRAC(2)=0.05
 XFRAC(3)=0.075 ;XFRAC(4)=0.1
 XFRAC(5)=0.125 ;XFRAC(6)=0.15
 XFRAC(7)=0.175 ;XFRAC(8)=0.2
 XFRAC(9)=0.225 ;XFRAC(10)=0.25
 XFRAC(11)=0.275 ;XFRAC(12)=0.3
 XFRAC(13)=0.325 ;XFRAC(14)=0.35
 XFRAC(15)=0.375 ;XFRAC(16)=0.4
 XFRAC(17)=0.425 ;XFRAC(18)=0.45
 XFRAC(19)=0.475 ;XFRAC(20)=0.5
 XFRAC(21)=0.525 ;XFRAC(22)=0.55
 XFRAC(23)=0.575 ;XFRAC(24)=0.6
 XFRAC(25)=0.625 ;XFRAC(26)=0.65
 XFRAC(27)=0.675 ;XFRAC(28)=0.7
 XFRAC(29)=0.725 ;XFRAC(30)=0.75
 XFRAC(31)=0.775 ;XFRAC(32)=0.8
 XFRAC(33)=0.825 ;XFRAC(34)=0.85
 XFRAC(35)=0.875 ;XFRAC(36)=0.9
 XFRAC(37)=0.925 ;XFRAC(38)=0.95
 XFRAC(39)=0.975 ;XFRAC(40)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 40
 YVLAST =20.
 YFRAC(1)=0.025 ;YFRAC(2)=0.05
 YFRAC(3)=0.075 ;YFRAC(4)=0.1
 YFRAC(5)=0.125 ;YFRAC(6)=0.15
 YFRAC(7)=0.175 ;YFRAC(8)=0.2
 YFRAC(9)=0.225 ;YFRAC(10)=0.25
 YFRAC(11)=0.275 ;YFRAC(12)=0.3
 YFRAC(13)=0.325 ;YFRAC(14)=0.35
 YFRAC(15)=0.375 ;YFRAC(16)=0.4
 YFRAC(17)=0.425 ;YFRAC(18)=0.45
 YFRAC(19)=0.475 ;YFRAC(20)=0.5
 YFRAC(21)=0.525 ;YFRAC(22)=0.55
 YFRAC(23)=0.575 ;YFRAC(24)=0.6
 YFRAC(25)=0.625 ;YFRAC(26)=0.65
 YFRAC(27)=0.675 ;YFRAC(28)=0.7
 YFRAC(29)=0.725 ;YFRAC(30)=0.75
 YFRAC(31)=0.775 ;YFRAC(32)=0.8
 YFRAC(33)=0.825 ;YFRAC(34)=0.85
 YFRAC(35)=0.875 ;YFRAC(36)=0.9
 YFRAC(37)=0.925 ;YFRAC(38)=0.95
 YFRAC(39)=0.975 ;YFRAC(40)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 3
 ZWLAST =0.5
 ZFRAC(1)=0.333333 ;ZFRAC(2)=0.666667
 ZFRAC(3)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(3)=U1 ;NAME(5)=V1
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(U1,Y,N,N,N,N,Y)
 SOLUTN(V1,Y,N,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         
 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)
 FIINIT(U1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 
 PATCH(INI1 ,INIVAL, 1, 40, 1, 40, 1, 1, 1, 1)
 
 PATCH(INI2 ,INIVAL, 1, 40, 1, 40, 2, 2, 1, 1)
 
 PATCH(INI3 ,INIVAL, 1, 40, 1, 40, 3, 3, 1, 1)
 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
 RELAX(U1,LINRLX,1.)
 RELAX(V1,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,INITIAL,U1!INI1,C,=XU-10.)
 SPEDAT(SET,INITIAL,V1!INI1,C,=-(YV-10.))
 SPEDAT(SET,INITIAL,U1!INI2,C,=YG-10.)
 SPEDAT(SET,INITIAL,V1!INI2,C,=-(XG-10.))
 SPEDAT(SET,INITIAL,U1!INI3,C,=-U1[&&-1]+U1[&&1])
 SPEDAT(SET,INITIAL,V1!INI3,C,=-V1[&&2]+V1[&&-2])
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(U1)=5.485 ;EX(V1)=5.485
 ************************************************************
  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(U1,Y,N,Y,N,N,N)
 OUTPUT(V1,Y,N,Y,N,N,N)
 ************************************************************
  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
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 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