TALK=T;RUN(1,1)
  DISPLAY

  The two spheres move from left to right within the upper and
  lower halves of the domain, specified by PATCHes FIRST and
  SECOND.

  Their motion is specified by INFOB commands, wherein the
  formula takes the form:

    SPHERE(x-centre, y-centre, z-centre, radius) ,

  each of the arguments is itself a formula.

  SOURCE commands specify how much heat the spheres emit,
  also by way of formulae.

  Several formulae are provided from which the user can choose.

  WARNING:
  If the case is to be run by way of the VR-Editor, the patches
  should be extended to cover the whole domain; for otherwise
  the Editor will convert them to objects with different names.

  Then a single patch will suffice, e.g.

    PATCH(WHOLEDOM,1,NX,1,NY,1,LSTEP)

  to which the INFOB statements should refer.

  The Q1 contains PHOTON USE commands
  ENDDIS
  PHOTON USE
  p
  p10;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 100 sec.
  pause
  p
  p20;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 200 sec.
  pause
  p
  p30;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 300 sec.
  pause
  p
  p40;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 400 sec.
  pause
  p
  p50;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 500 sec.
  pause
  p
  p60;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 600 sec.
  pause
  p
  p70;

  gr ou z 1
  con tem1 z 1 fil;.00001
  msg at 700 sec.
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(In-Form Represents Hot Moving Spheres   )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 765
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =700.
 FSTEP = 1 ;LSTEP = 70
 TFRAC(1)=0.014286 ;TFRAC(2)=0.028571
 TFRAC(3)=0.042857 ;TFRAC(4)=0.057143
 TFRAC(5)=0.071429 ;TFRAC(6)=0.085714
 TFRAC(7)=0.1 ;TFRAC(8)=0.114286
 TFRAC(9)=0.128571 ;TFRAC(10)=0.142857
 TFRAC(11)=0.157143 ;TFRAC(12)=0.171429
 TFRAC(13)=0.185714 ;TFRAC(14)=0.2
 TFRAC(15)=0.214286 ;TFRAC(16)=0.228571
 TFRAC(17)=0.242857 ;TFRAC(18)=0.257143
 TFRAC(19)=0.271429 ;TFRAC(20)=0.285714
 TFRAC(21)=0.3 ;TFRAC(22)=0.314286
 TFRAC(23)=0.328571 ;TFRAC(24)=0.342857
 TFRAC(25)=0.357143 ;TFRAC(26)=0.371429
 TFRAC(27)=0.385714 ;TFRAC(28)=0.4
 TFRAC(29)=0.414286 ;TFRAC(30)=0.428571
 TFRAC(31)=0.442857 ;TFRAC(32)=0.457143
 TFRAC(33)=0.471429 ;TFRAC(34)=0.485714
 TFRAC(35)=0.5 ;TFRAC(36)=0.514286
 TFRAC(37)=0.528571 ;TFRAC(38)=0.542857
 TFRAC(39)=0.557143 ;TFRAC(40)=0.571429
 TFRAC(41)=0.585714 ;TFRAC(42)=0.6
 TFRAC(43)=0.614286 ;TFRAC(44)=0.628571
 TFRAC(45)=0.642857 ;TFRAC(46)=0.657143
 TFRAC(47)=0.671429 ;TFRAC(48)=0.685714
 TFRAC(49)=0.7 ;TFRAC(50)=0.714286
 TFRAC(51)=0.728571 ;TFRAC(52)=0.742857
 TFRAC(53)=0.757143 ;TFRAC(54)=0.771429
 TFRAC(55)=0.785714 ;TFRAC(56)=0.8
 TFRAC(57)=0.814286 ;TFRAC(58)=0.828571
 TFRAC(59)=0.842857 ;TFRAC(60)=0.857143
 TFRAC(61)=0.871429 ;TFRAC(62)=0.885714
 TFRAC(63)=0.9 ;TFRAC(64)=0.914286
 TFRAC(65)=0.928571 ;TFRAC(66)=0.942857
 TFRAC(67)=0.957143 ;TFRAC(68)=0.971429
 TFRAC(69)=0.985714 ;TFRAC(70)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 80
 XULAST =4.
 XFRAC(1)=0.0125 ;XFRAC(2)=0.025
 XFRAC(3)=0.0375 ;XFRAC(4)=0.05
 XFRAC(5)=0.0625 ;XFRAC(6)=0.075
 XFRAC(7)=0.0875 ;XFRAC(8)=0.1
 XFRAC(9)=0.1125 ;XFRAC(10)=0.125
 XFRAC(11)=0.1375 ;XFRAC(12)=0.15
 XFRAC(13)=0.1625 ;XFRAC(14)=0.175
 XFRAC(15)=0.1875 ;XFRAC(16)=0.2
 XFRAC(17)=0.2125 ;XFRAC(18)=0.225
 XFRAC(19)=0.2375 ;XFRAC(20)=0.25
 XFRAC(21)=0.2625 ;XFRAC(22)=0.275
 XFRAC(23)=0.2875 ;XFRAC(24)=0.3
 XFRAC(25)=0.3125 ;XFRAC(26)=0.325
 XFRAC(27)=0.3375 ;XFRAC(28)=0.35
 XFRAC(29)=0.3625 ;XFRAC(30)=0.375
 XFRAC(31)=0.3875 ;XFRAC(32)=0.4
 XFRAC(33)=0.4125 ;XFRAC(34)=0.425
 XFRAC(35)=0.4375 ;XFRAC(36)=0.45
 XFRAC(37)=0.4625 ;XFRAC(38)=0.475
 XFRAC(39)=0.4875 ;XFRAC(40)=0.5
 XFRAC(41)=0.5125 ;XFRAC(42)=0.525
 XFRAC(43)=0.5375 ;XFRAC(44)=0.55
 XFRAC(45)=0.5625 ;XFRAC(46)=0.575
 XFRAC(47)=0.5875 ;XFRAC(48)=0.6
 XFRAC(49)=0.6125 ;XFRAC(50)=0.625
 XFRAC(51)=0.6375 ;XFRAC(52)=0.65
 XFRAC(53)=0.6625 ;XFRAC(54)=0.675
 XFRAC(55)=0.6875 ;XFRAC(56)=0.7
 XFRAC(57)=0.7125 ;XFRAC(58)=0.725
 XFRAC(59)=0.7375 ;XFRAC(60)=0.75
 XFRAC(61)=0.7625 ;XFRAC(62)=0.775
 XFRAC(63)=0.7875 ;XFRAC(64)=0.8
 XFRAC(65)=0.8125 ;XFRAC(66)=0.825
 XFRAC(67)=0.8375 ;XFRAC(68)=0.85
 XFRAC(69)=0.8625 ;XFRAC(70)=0.875
 XFRAC(71)=0.8875 ;XFRAC(72)=0.9
 XFRAC(73)=0.9125 ;XFRAC(74)=0.925
 XFRAC(75)=0.9375 ;XFRAC(76)=0.95
 XFRAC(77)=0.9625 ;XFRAC(78)=0.975
 XFRAC(79)=0.9875 ;XFRAC(80)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 40
 YVLAST =2.
 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 = 1
 ZWLAST =0.1
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(149)=SPH1 ;NAME(150)=TEM1
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(SPH1,Y,N,N,N,N,N)
 SOLUTN(TEM1,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(TEM1,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 =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.
 PRNDTL(TEM1)=0.1
 PRT(TEM1)=1.
 CP1 =1000. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(SPH1)=1.0E-10 ;FIINIT(TEM1)=20.
   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(WALLE ,EWALL , 80, 80, 1, 40, 1, 1, 1, 70)
 COVAL(WALLE ,TEM1,1. ,20. )
 
 PATCH(WALLW ,WWALL , 1, 1, 1, 40, 1, 1, 1, 70)
 COVAL(WALLW ,TEM1,1. ,20. )
 
 PATCH(WALLN ,NWALL , 1, 80, 40, 40, 1, 1, 1, 70)
 COVAL(WALLN ,TEM1,1. ,20. )
 
 PATCH(WALLS ,SWALL , 1, 80, 1, 1, 1, 1, 1, 70)
 COVAL(WALLS ,TEM1,1. ,20. )
 
 PATCH(FIRST ,VOLUME, 1, 80, 1, 20, 1, 1, 1, 70)
 COVAL(FIRST ,TEM1,In-Form:source - see Grp 19)
 
 PATCH(SECOND ,VOLUME, 1, 80, 21, 40, 1, 1, 1, 70)
 COVAL(SECOND ,TEM1,In-Form:source - see Grp 19)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 30 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER(TEM1)=20
 ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(SPH1,LINRLX,1.)
 RELAX(TEM1,FALSDT,1.0E+09)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
 VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,INFOB,*!FIRST,C,=SPHERE(.5+.5*(TIM/100.-1)&.5&.05&.25&$)
 SPEDAT(SET,INFOB,*!FIRST,C,0&0&0)!INFOB_1)
 SPEDAT(SET,INFOB,*!SECOND,C,=SPHERE(.5+.25*(TIM/100.-1)&1.5&.05&.$)
 SPEDAT(SET,INFOB,*!SECOND,C,25&0&0&0)!INFOB_2)
 SPEDAT(SET,SOURCE,TEM1!FIRST,C,=100.!INFOB_1)
 SPEDAT(SET,SOURCE,TEM1!SECOND,C,=100.!INFOB_2)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(SPH1)=1000. ;EX(TEM1)=22.950001
 ************************************************************
  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(SPH1,N,N,Y,Y,N,N)
 OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 21 ;IYMON = 11 ;IZMON = 1
 NPRMON = 30 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 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
 IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='P'
STOP