#cls
text(In-Form Equivalent Of PLANT Case Z108
title
libref=708
DISPLAY
This In-Form case does what Z108 does, but without fortran coding
or re-compilation.
First, case z108 is loaded.
Then NAMSAT = NONE is set, so as to deactivate PLANT.
Then, the following In-Form statement is provided:
(SOURCE of H1 at iRAD is radco*(TMPX^4-H1^4) with LINE)
This ensures that the patch iRAD (at the surface of a block) loses
heat by radiation at a rate dependent on the fourth power of its
temperature (H1) .
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
p
msg Temperature contours
con h1 z m fi;0.001
con h1 y m fi;0.001
con h1 x m fi;0.001
msg Press e to END
ENDUSE
load(z108)
text(In-Form equivalent of PLANT case Z108
libref=708
** The next line de-activates PLANT **
namsat=none
RAD=SKIP
** The following settings have been chosen as being
preferable to those of z108, in respect of convergence
and printout. **
FIINIT(H1)=TMPX
RELAX(H1,LINRLX,0.35)
** In-Form formulae which are equivalent to the
PLANT formulae of z108 **
INFORM13BEGIN
PATCH(iRAD,high,1,NX,1,NY,NZ,NZ,1,1)
(SOURCE of H1 at iRAD is radco*(TMPX^4-H1^4) with LINE)
INFORM13END
DISTIL=T
EX(H1)=6.898E+01
STOP