L($100)
TEXT(Case 100; Various Solver Options 
lit=-liter(temp)
OVRRLX=1.0;isolbk=0.; isolx=0.; isoly=0.; isolz=0.
do ii=1,5
+ mesg(
enddo
mesg( Initial data that can be changed:
mesga( Number of iterations in the solver is :lit:
mesg( Over-relaxation factor is :ovrrlx:
mesg( Frequency of block correction is :isolbk:
mesg( Frequency of x-wise correction is :isolx:
mesg( Frequency of y-wise correction is :isoly:
mesg( Frequency of z-wise correction is :isolz:
mesga( Do you want to change settings (y/n)? (Default n)
readvdu(ans,char,n)
if(:ans:.eq.y) then
mesga(Change NX, NY or NZ ? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y) then
mesg(change NX, currently :nx: to ...
readvdu(nx,int,nx)
mesg(change NY, currently :ny: to ...
readvdu(ny,int,ny)
mesg(change NZ, currently :nZ: to ...
readvdu(nZ,int,nZ)
goto top
endif
mesga(Number of iterations in the solver is :lit:. OK?
mesg(If not, insert new value
readvdu(lit,int,:lit:)
lit
liter(temp)=-lit
  ** Set the over-relaxation factor for the linear-equation
     solver.
OVRRLX=1.0
mesga(Over-relaxation factor is :ovrrlx:. OK?
mesg(If not, insert new value
readvdu(ovrrlx,real,:ovrrlx:)
ovrrlx
isolbk=0
mesga(Frequency of block correction is :isolbk:. OK?
mesg(If not, insert new value
readvdu(isolbk,int,:isolbk:)
isolbk
isolx=0
mesga(Frequency of x-wise correction is :isolx:. OK?
mesg(If not, insert new value
readvdu(isolx,int,:isolx:)
isolx
isoly=0
mesga(Frequency of y-wise correction is :isoly:. OK?
mesg(If not, insert new value
readvdu(isoly,int,:isoly:)
isoly
mesga(Frequency of z-wise correction is :isolz:. OK?
mesg(If not, insert new value
readvdu(isolz,int,:isolz:)
isolz
endif
 
selref=t; resfac=1.e-3;TSTSWP=1
mesga(Are settings now OK? Enter n to start again
ans=y
readvdu(ans,char,y)
if(:ans:.eq.n) then
 goto top
endif
 **END OF LIBRARY CASE 100