#  ^ Z
#  |                This is an example of two enclosures (boxes) joined
#  |---------       together to form a 1X1X1 cube. The walls which join
#  |         |      the two boxes together become virtual to the model
#  | Box 1   |      and do not interact in the radiation heat transfer.
#  |---------|      Note, Both boxes have a participating media with 
#  | Box 2   |      different properties.
#  |         |               
#  X-----------> Y Finally, all known temperatures and fluxes are based on
#                a value of 100 K. If the model runs correctly, all walls
#                and media should have a final temperature of 100 K and a
#                net flux of 0 Watts.

Seed(11516)     # The random number seed is 11516 this keyword is optional
Mode(Temperature) # The model will operate in the temperature mode  
Bundle(0.0001)    # Each photon bundle send out has 0.0001 Watts of energy
Model(2)          # '2` implies this model is built from two boxes

    # Enclosure 1 is a box where X = 1, Y = 1, and Z = 1
Box(1,1,1,0.5)
 
   # Wall 1 is a joining wall to box 2, wall 4. It is not a real wall, but
   # virtual
Wall(1,join,2,4)

    # Wall 2 is described in terms of flux sent into the system
    # 4.536 Watts/m/m => 100 K
Wall(2,flux,4.536,0.8,.2,0)

    # Wall 3 uses the initial temperature mode
    # Emissivity = 0.5 and diffuse scattering
Wall(3,temp,100,0.5,.5,0)

# Wall four has an unknown temperature, its value will be determined
# by iteration. It will be meshed into a 5X5 grid with an emis. of 0.6
Wall(4,unknown,5,5,.6,.4,0)

    # The keyword 'Iterate' allows for the following information
    # Max num. of iterations = 40 or stop at 5% error, Update factor = 1.0
    # Final flux value for the 5X5 mesh is 0.0
Iterate(40,0.5,1.0, 0,0,0,0,0 0,0,0,0,0 0,0,0,0,0 0,0,0,0,0 0,0,0,0,0)

Wall(5,temp,100,1.0,0,0)

# Wall 6 is broken down into a 3X3 mesh where each element property can be 
# described separately.
Wall(6,element,3,3)
Welement(temp, 100,1.0,0,0   100,0.7,.3,0   100,0.1,0,.9
               100,0.5,0,.5  100,1.0,0,0    100,1.0,0,0
               100,0.8,.2,0  100,0.1,0.9,0  100,0.3,.7,0)

# Box 1 has a Media which has a temp.= 100 K, Ke = 1.0 and ssal = 0.2 with
# isotropic scattering.
Media(1,temp,100,1.0,0.2,isotropic)

#---------- Second Enclosure ------------
Box(2,1,1,0.5)

    # Wall 1 is broken into a 2X2 grid with unknown temperatures
Wall(1,Element,2,2)
WELEMENT(Unknown  1.0,0,0,  0.5,.5,0,
                  0.1,0,.9, 0.9,.1,0 )

    #Iteration parameters are the same as the other box's wall
iterate(40,0.5,1.0, 0, 0,
                    0, 0)

    # Walls 2 - 6 use the temp keyword except wall 4 which joins box 2 to
    # box 1.
Wall(2,temp,100,0.3,.7,0)

   # Wall three has been removed from the model, the far field temperature
   # is 100 K
Wall(3,remove,100)

Wall(4,join,1,1)
Wall(5,temp,100,1.0,0,0)
Wall(6,temp,100,1.0,0,0)

 # Box 2 has an unknown Media temperature, notice only a 1 % tolerance is 
 # given in the iteration statement.
Media(2,unknown,2,2,4,0.7,0.9,isotropic)
Iterate(40,0.7,1.0,  0,0,
                     0,0,
           
                     0,0,
                     0,0,
           
                     0,0,
                     0,0,
           
                     0,0,
                     0,0 )