TALK=F;RUN( 1, 1)
Group 1. Run Title
TEXT(sub-grid objects POINT, LINE and PLANE
DISPLAY
The use if illustrated of the commands for introducing, placing
and sizing POINT, LINE and PLANE objects, the parameters of which
are supplied via PIL DO loops.
It is recommended that the VR-Viewer should be used for displaying
them visually.
ENDDIS
LIBREF=385
************************************************************
Group 2. Transience
STEADY=T
************************************************************
Groups 3, 4, 5 Grid Information
NX=50;NY=50;NZ=50
XULAST=1.0; YVLAST=1.0; ZWLAST=1.0
#UNIGRID
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Solved variables list
* Stored variables list
STORE(PRPS,MARK) ! MARK is set to non-zero values inside In-Form
objects
STORE(VPOR,EPOR,NPOR,HPOR)
************************************************************
Group 9 Properties
SPEDAT(SET,MATERIAL,100,L,T) ! needed to warn earth to load the
! properties of material to be set
! in in-form objects
************************************************************
Group 11.Initialise Var/Porosity Fields
patch(patch1,volume,1,nx,1,ny,1,nz,1,1) ! the space within which
! infob_1 will be placed
FIINIT(PRPS) = 0.0
FIINIT(PRPS) = 0.0 !
inform11begin
real(x0,y0,z0,nomrad) ! declarations
Point objects at fixed y and z but varying x
do ii = 1,10
(stored of mark at patch1 is 1.0 with infob_:ii:) ! MARK object
(initial of prps is 100 with infob_:ii:) ! initialise PRPS
x0=:ii:*xulast/10 ! x/y/z position of point object
y0=:ii:*yvlast/10 !
z0=zwlast/10 !
nomrad=:ii:*xulast/100 ! x/y/z radius of POINT
(infob at patch1 is point(x0,y0,z0,nomrad) with infob_:ii:!poro)
enddo
line objects
real(x1,y1,z1)
x0=xulast/10
z0=zwlast/10 !
z0=zwlast-z0
z1=z0
do ii = 1,10
(stored of mark at patch1 is 1.0 with infob_:ii+10:) ! MARK object
(initial of prps is 100 with infob_:ii+10:) ! initialise PRPS
x1=xulast-:ii:*xulast/10
y0=:ii:*yvlast/10 !
y1=y0
nomrad=:ii:*xulast/100 ! x/y/z radius of LINE
(infob at patch1 is LINE(x0,y0,z0,x1,y1,z1,nomrad) with infob_:ii+$
10:!poro)
enddo
plane quadrilateral objects
real(x2,y2,z2,x3,y3,z3)
do ii=1,3
(stored of mark at patch1 is 1.0 with infob_:ii+20:) ! MARK object
(initial of prps is 100 with infob_21) ! initialise PRPS
x0=xulast/2
x0=x0+0.2*:ii-1:
x1=x0
x2=x0
x3=x0
y0=3*yvlast/10
y1=yvlast-y0*:ii:
y2=yvlast-y1
y3=y1
z0=3*zwlast/10
z1=zwlast-z0*:ii:
z2=z1
z3=z1
nomrad=xulast/100 ! x/y/z radius of LINE
(infob at patch1 is PLANE(x0,y0,z0,x1,y1,z1,x2,y2,z2,x3,y3,z3,nomrad$
) with infob_:ii+20:!poro)
enddo
inform11end
************************************************************
Group 13. Boundary & Special Sources
************************************************************
Group 15. Terminate Sweeps
yzpr=t
tstswp=-1
#endpause
#maxabs
stop