GROUP 1. Run title and other preliminaries                                  
TEXT(In-Form Represents Hot Moving Spheres                                      
TITLE                                                                           
libref=765                                                                      
  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                                                                        
  GROUP 2. Transience; time-step specification                                  
STEADY=F;GRDPWR(T,70,700.,1.0)                                                  
                                                                                
    GROUP 3. X-direction grid specification                                     
GRDPWR(X,80,4.,1.0)                                                             
                                                                                
    GROUP 4. Y-direction grid specification                                     
GRDPWR(Y,40,2.,1.0)                                                             
                                                                                
    GROUP 5. Z-direction grid specification                                     
GRDPWR(Z,1,0.1,1.0)                                                             
                                                                                
    GROUP 7. Variables stored, solved & named                                   
NAME(150)=TEM1                                                                  
SOLVE(TEM1)                                                                     
                                                                                
PRNDTL(TEM1)=0.1                                                                
    GROUP 9.Properties used
CP1=1000.                                                                    
    GROUP 11. Initialization of variable or porosity fields                     
FIINIT(TEM1)=20.                                                                
                                                                                
    GROUP 13. Boundary conditions and special sources                           
PATCH(WALLE,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)                                      
COVAL(WALLE,TEM1,1.0,20.0)                                                      
                                                                                
PATCH(WALLW,WWALL,1,1,1,NY,1,NZ,1,LSTEP)                                        
COVAL(WALLW,TEM1,1.0,20.0)                                                      
                                                                                
PATCH(WALLN,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)                                      
COVAL(WALLN,TEM1,1.0,20.0)                                                      
                                                                                
PATCH(WALLS,SWALL,1,NX,1,1,1,NZ,1,LSTEP)                                        
COVAL(WALLS,TEM1,1.0,20.0)                                                      
                                                                                
  INFORM13BEGIN                                                                 
  *** Hot sphere - 1                                                            
PATCH(FIRST,VOLUME,1,NX,1,NY/2,1,NZ,1,LSTEP)                                    
                                                                                
  *** Hot sphere - 2                                                            
PATCH(SECOND,VOLUME,1,NX,NY/2+1,NY,1,NZ,1,LSTEP)                                
  Note: The FIRST and SECOND patches cover the parts of the domain              
        in which the spheres move, namely the upper and the lower               
        halves.                                                                 
                                                                                
        Where the spheres actually are at any time is conveyed by               
        the INFOB statements, which also assign the (obligatory)                
        names INFOB_1 and INFOB_2.                                              
                                                                                
  Note: xce,yce and zce are the x,z and z coordinates of the sphere             
        centre, and radius has its obvious meaning.                             
        They are character variables which are evaluated by                     
        Satellite because they are enclosed within colons.                      
                                                                                
char(xce,yce,zce,radius,sour)                                                   
xce=.5+.5*(tim/100.-1); yce=.5; zce=.05; radius= .25                            
  The above setting of yce causes the first sphere to follow a                  
  horizontal path.                                                              
                                                                                
  Move the following statement one ot two spaces to the left if the             
  path is to have a gradient of 0.1/0.5; and also increase the                  
  PATCH(FIRST statement to NY, so as to ensure that the patch is                
  large enough.                                                                 
                                                                                
  yce=0.1 + 01*(tim/100.-1)                                                     
(INFOB at FIRST is SPHERE(:xce:,:yce:,:zce:,:radius:$                           
) with INFOB_1)                                                                 
                                                                                
xce=.5+.25*(tim/100.-1); yce=1.5; zce=.05; radius= .25                          
                                                                                
(INFOB at SECOND is SPHERE(:xce:,:yce:,:zce:,:radius:$                          
) with INFOB_2)                                                                 
                                                                                
  *** Setting of TEM1 sources into the two spheres                              
sour=100.                                                                       
(SOURCE of TEM1 at FIRST is :sour: with INFOB_1)                                
                                                                                
  The following statement, moved to the left so as to become active,            
  would cause the magnitude of the soulce to increase with                      
  temperature, a potentially explosive situation.                               
  sour=10.0 + 2.0*tem1                                                          
                                                                                
  Otherwise, the sour=100. setting will be used                                 
(SOURCE of TEM1 at SECOND is :sour: with INFOB_2)                               
  INFORM13END                                                                   
    GROUP 15. Termination of sweeps                                             
LSWEEP=30                                                                       
                                                                                
    GROUP 22. Spot-value print-out                                              
NPRMON=LSWEEP;IXMON=NX/4+1;IYMON=NY/4+1;IZMON=NZ/2+1                            
    GROUP 23. Print-out & plot control                                          
TSTSWP=-1;IDISPA=10; CSG1=P                                                     
SELREF=T; RESFAC=1.E-2                                                          
NXPRIN=1;NYPRIN=1                                                               
DISTIL=T;EX(TEM1)=2.295E+01                                                     
  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