#cls
text(In-Form Equivalent Of PLANT Case Z111
title
libref=711
DISPLAY
This In-Form case does what Z111 does, but without fortran coding
or re-compilation.
First, case z111 is loaded. Then NAMSAT = NONE is set, to
deactivate.
The Q1 uses the In-Form functions:
SOURCE and STORED.
ENDDIS
load(z111)
text(In-Form equivalent of PLANT case Z111
libref=711
INFORM13BEGIN
** The next line de-activates PLANT **
namsat=none
COVAL(SOURCE,C1,0.,0.)
INFORM13END
** The following settings have been chosen as being
preferable to those of z111, in respect of convergence
and printout. **
LSWEEP=7
RELAX(C1,FALSDT,10.)
** In-Form formulae which are equivalent to the
PLANT formulae of z111 **
INFORM13BEGIN
PATCH(ISOURCE,CELL,1,NX,1,NY,1,1,1,lstep)
(SOURCE of C1 at ISOURCE is EXP(TIM)-5.*C1)
(STORED of EXCT is 1./6*(EXP(TIM)+5.*EXP(-5.*TIM)))
INFORM13END
DISTIL=T
EX(C1)=4.695E-01; EX(EXCT)=4.587E-01
STOP