TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2014, Date 26/09/14
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
up z
gr ou z 1;gr ou z m;gr ou x 1;gr ou x m
gr ou z 5 x 2 4 y 2 4; gr ou z 5 x 2 4 y 7 9
gr ou z 5 x 7 9 y 7 9; gr ou z 5 x 7 9 y 2 4
gr ou z m x 2 4 y 2 4; gr ou z m x 2 4 y 7 9
gr ou z m x 7 9 y 7 9; gr ou z m x 7 9 y 2 4
gr ou x 2 z 5 m y 2 4; gr ou x 2 z 5 m y 7 9
gr ou x 5 z 5 m y 2 4; gr ou x 5 z 5 m y 7 9
gr ou x 7 z 5 m y 2 4; gr ou x 7 z 5 m y 7 9
gr ou x 10 z 5 m y 2 4; gr ou x 10 z 5 m y 7 9
gr ou z 3 col 12; gr ou z 6 col 13
block;4
2 4 2 4 5 7
2 4 7 9 5 7
7 9 2 4 5 7
7 9 7 9 5 7
msg The geometry
pause
gr del; gr del
vi z
set con sca pl
set por on
con hi z 2 sh; int 14
msg The metal-eletrolyte interface height
pause
con cl;red
con hi z 5 sh; int 14
msg The electrolyte free surface height
pause
con cl;red
vec z 1 sh
set prop off;set por off
con prps z 1 col 14;va;1;197
msg Metal Vectors above the cathode, showing outline of freeze
pause
vec cl;con cl;red
vec z 2 sh
con prps z 2 col 14;va;1;197
msg Metal vectors below the interface, showing outline of freeze
pause
vec cl;con cl;red
vec z 3 sh
msg Electrolyte vectors above the interface
pause
vec cl;red
vec z 4 sh
msg Electrolyte vectors below the anodes
pause
vec cl;red
vec z 5 sh
msg Electrolyte vectors between the anodes
pause
vec cl;red
con jz z 2 sh
int 14
msg Contours of Z-directed current below the interface
pause
con cl;red
con vist z 2 sh;int 14
msg Contours of turbulent viscosity below the interface
pause
con cl;red
con vist z 3 sh;int 14
msg Contours of turbulent viscosity above the interface
ENDUSE
DISPLAY
This case concerns a simplified 4-anode Hall Cell, with no
under-anode gas generation.
The K-Epsilon turbulence model is used.
The constants in the Bx, By and Bz expressions are
-0.012, 0.002 and -0.004 respectively.
ENDDIS
ESTRBEGIN
Debug switches for GROUND
=========================
LG(9) - Interface & free surface adjustment
LG(10) - Current calculation
LG(11) - Interface velocity for transient
LG(14) - Zeroing of pressure & potential inside freeze
IG(17),IG(18) - IZF,IZL First and last slab
IG(19),IG(20) - ISWPF,ISWPL First and last sweep
*************************************************************
The next lines define the cathode current distribution, in
the orientation
^
| IY
|
------> IX
The current is given in A/m**2
*JZC
0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. -5000. -5000. -5000. -5000. -5000. -5000. 0.
0. -5000. -5000. 0. 0. -5000. -5000. -5000. -5000. 0.
0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
*************************************************************
The next lines define the anode locations, in the orientation
^
| IY
|
------> IX
0 means open cell, 1 means cell with anode.
*IANO
0 0 0 0 0 0 0 0 0 0
0 1 1 1 0 0 1 1 1 0
0 1 1 1 0 0 1 1 1 0
0 1 1 1 0 0 1 1 1 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 1 1 1 0 0 1 1 1 0
0 1 1 1 0 0 1 1 1 0
0 1 1 1 0 0 1 1 1 0
0 0 0 0 0 0 0 0 0 0
*************************************************************
The next lines define the freeze locations and heights, in
the orientation
^
| IY
|
------> IX
0 means no freeze, >0 defines the height of freeze in Z-slabs.
*IFRZ
1 1 1 1 1 1 1 1 1 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 0 1
1 0 0 1 1 0 0 0 0 1
2 1 1 1 1 1 1 1 1 1
*************************************************************
The next lines define the heights of the anode underside
above the cathode, in the orientation
^
| IY
|
------> IX
The anode height is in metres.
*HANO
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26 0.26
Alternative anode height definition, making anodes slope
inwards towards the inter-anode gap. To activate, move
*HANO above to right by several spaces, and that below to
start in column 3.
*HANO
.260 .260 .260 .260 .260 .260 .260 .260 .260 .260
.260 .260 .260 .260 .260 .260 .260 .260 .260 .260
.267 .267 .267 .267 .267 .267 .267 .267 .267 .267
.274 .274 .274 .274 .274 .274 .274 .274 .274 .274
.274 .274 .274 .274 .274 .274 .274 .274 .274 .274
.274 .274 .274 .274 .274 .274 .274 .274 .274 .274
.274 .274 .274 .274 .274 .274 .274 .274 .274 .274
.267 .267 .267 .267 .267 .267 .267 .267 .267 .267
.260 .260 .260 .260 .260 .260 .260 .260 .260 .260
.260 .260 .260 .260 .260 .260 .260 .260 .260 .260
ESTREND
************************************************************
IRUNN = 1 ;LIBREF = 498
************************************************************
Group 1. Run Title
TEXT(ESTER demonstration, one-phase, KE-EP )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,10,7)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(14)=EPOT ;NAME(15)=HI
NAME(132)=SURF ;NAME(133)=JIZ
NAME(134)=JIY ;NAME(135)=JIX
NAME(136)=LEN1 ;NAME(137)=VIST
NAME(138)=PRPS ;NAME(139)=NPOR
NAME(140)=EPOR ;NAME(141)=CPOR
NAME(142)=BZ ;NAME(143)=BY
NAME(144)=BX ;NAME(145)=JZ
NAME(146)=JY ;NAME(147)=JX
NAME(148)=FZ ;NAME(149)=FY
NAME(150) =FX
* Solved variables list
SOLVE(P1,U1,V1,W1,EPOT)
* Stored variables list
STORE(FX,FY,FZ,JX,JY,JZ,BX,BY)
STORE(BZ,CPOR,EPOR,NPOR,PRPS,VIST,LEN1,JIX)
STORE(JIY,JIZ,SURF,HI)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(EPOT,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
TERMS (EPOT,N,N,Y,N,Y,N)
UDIFF = T
************************************************************
Group 9. Properties
* Domain material index is 51 signifying:
* etal
SETPRPS(1, 51,ESTR)
DRH1DP =1.0E-10
PRNDTL(EPOT)=-5.0E+06
PRT(EP)=1.314 ;PRT(EPOT)=1.0E+10
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(KE)=1.0E-05 ;FIINIT(EP)=2.0E-08
FIINIT(EPOT)=3. ;FIINIT(HI)=1.001E-10
FIINIT(JIZ)=1.001E-10 ;FIINIT(JIY)=1.001E-10
FIINIT(JIX)=1.001E-10 ;FIINIT(LEN1)=1.001E-10
FIINIT(VIST)=1.001E-10 ;FIINIT(PRPS)=-1.
FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.
FIINIT(CPOR)=1. ;FIINIT(BZ)=1.001E-10
FIINIT(BY)=1.001E-10 ;FIINIT(BX)=1.001E-10
FIINIT(JZ)=1.001E-10 ;FIINIT(JY)=0.
FIINIT(JX)=0. ;FIINIT(FZ)=1.001E-10
FIINIT(FY)=1.001E-10 ;FIINIT(FX)=1.001E-10
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
This patch is attached to object PREF1
PATCH(PREF1, CELL, -1, 1, 0, 0, 0, 0, 1, 1)
COVAL(PREF1, P1, 1., 0.)
COVAL(PREF1, U1, 0., SAME)
COVAL(PREF1, V1, 0., SAME)
COVAL(PREF1, W1, 0., SAME)
This patch is attached to object PREF2
PATCH(PREF2, CELL, -1, 1, 0, 0, 0, 0, 1, 1)
COVAL(PREF2, P1, 1., 0.)
COVAL(PREF2, U1, 0., SAME)
COVAL(PREF2, V1, 0., SAME)
COVAL(PREF2, W1, 0., SAME)
This patch is attached to object ANOPTNTL
PATCH(ANOPTNTL, HIGH, -1, 0, 0, 0, 0, 0, 1, 1)
COVAL(ANOPTNTL, EPOT, GRND, GRND)
This patch is attached to object CATHCURR
PATCH(CATHCURR, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
COVAL(CATHCURR, EPOT, FIXFLU, GRND)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
SARAH =1000.
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
ENDIT(EPOT)=GRND1
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.8 )
RELAX(KE ,LINRLX,0.3 )
RELAX(EP ,LINRLX,0.3 )
RELAX(EPOT,LINRLX,0.8 )
KELIN = 0
************************************************************
Group 18. Limits
VARMAX(U1)=1. ;VARMIN(U1)=-1.
VARMAX(V1)=1. ;VARMIN(V1)=-1.
VARMAX(W1)=1. ;VARMIN(W1)=-1.
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = F
NAMGRD =ESTR
GENK = T
PARSOL = F
ISG62 = 1
CSG8 ='ESTR'
CSG10 ='ESTR'
SPEDAT(SET,MATERIAL,51,L,T)
SPEDAT(SET,MATERIAL,52,L,T)
SPEDAT(SET,MATERIAL,100,L,T)
SPEDAT(SET,MATERIAL,198,L,T)
SPEDAT(SET,MATERIAL,199,L,T)
SPEDAT(SET,MATERIAL,151,L,T)
LG( 1) = T
LG( 2) = T
LG( 3) = T
IG( 1) = 2
IG( 2) = 4
IG( 3) = 1
IG( 4) = 1
IG( 6) = 5
IG( 11) = 33
RG( 1) =2250.
RG( 2) =2050.
RG( 3) =1.
RG( 4) =5.0E+06
RG( 5) =200.
RG( 6) =1000.
RG( 7) =-9.81
RG( 9) =35.
RG( 11) =0.5
RG( 12) =100.
RG( 13) =0.9999
RG( 14) =1.
RG( 15) =0.5
RG( 18) =4.444445E-07
RG( 19) =4.878049E-07
RG( 20) =2.
RG( 21) =1.0E-05
RG( 22) =3.
RG( 23) =1.189
RG( 24) =1.0E-03
RG( 26) =-0.012
RG( 27) =2.0E-03
RG( 28) =-4.0E-03
RG( 29) =1.
RG( 31) =9.09E+04
RG( 32) =999.869995
RG( 33) =7800.
RG( 34) =9.08E+05
RG( 35) =2700.
RG( 36) =1.98E+07
RG( 37) =2700.
RG( 38) =2.86E+07
RG( 71) =7800.
RG( 72) =2.17E+06
RG( 73) =8950.
RG( 74) =1.34E+07
RG( 75) =7800.
RG( 76) =6.45E+05
RG( 77) =7800.
RG( 78) =8.7E+05
CG( 2) ='q1'
CG( 5) ='q1'
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(HI ,Y,N,Y,N,Y,Y)
OUTPUT(PRPS,N,N,Y,N,N,N)
OUTPUT(NPOR,N,N,Y,N,N,N)
OUTPUT(EPOR,N,N,Y,N,N,N)
OUTPUT(CPOR,N,N,Y,N,N,N)
OUTPUT(BZ ,N,N,Y,N,N,N)
OUTPUT(BY ,N,N,Y,N,N,N)
OUTPUT(BX ,N,N,Y,N,N,N)
OUTPUT(JZ ,Y,N,Y,N,Y,N)
OUTPUT(FZ ,Y,N,Y,N,Y,N)
OUTPUT(FY ,Y,N,Y,N,Y,N)
OUTPUT(FX ,Y,N,Y,N,Y,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 4 ;IYMON = 3 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = 12345
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.859158,-0.440119,0.261041)
GVIEW(UP,-0.204934,0.171492,0.963635)
GVIEW(DEPTH,3674.429932)
GVIEW(NEARPLANE,0.02)
GVIEW(VDIS,5.176495)
GVIEW(CENTRE,5.,2.,0.275)
> DOM, SIZE, 1.000000E+01, 4.000000E+00, 5.500000E-01
> DOM, MONIT, 3.500000E+00, 1.000000E+00, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_X_2, 3, 1.000000E+00
> GRID, RSET_X_3, 2, 1.000000E+00
> GRID, RSET_X_4, 3, 1.000000E+00
> GRID, RSET_X_5, 1, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Y_2, 3, 1.000000E+00
> GRID, RSET_Y_3, 2, 1.000000E+00
> GRID, RSET_Y_4, 3, 1.000000E+00
> GRID, RSET_Y_5, 1, 1.000000E+00
> GRID, RSET_Z_1, 2, 1.000000E+00
> GRID, RSET_Z_2, 2, 1.000000E+00
> GRID, RSET_Z_3, 1, 1.000000E+00
> GRID, RSET_Z_4, 2, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> DOM, CARBON-BLOCKS,NO
> DOM, ANODES, YES
> DOM, SIZE-X, 3.000000E+00, 2.000000E+00, 1.000000E+00
> DOM, SIZE-Y, 1.200000E+00, 8.000000E-01, 4.000000E-01
> DOM, SIZE-Z-1, 2.000000E-01, 2.900000E-01
> DOM, SIZE-Z-2, 2.500000E-01, 6.000002E-02
> DOM, NUM-X, 3, 2, 1
> DOM, NUM-Y, 3, 2, 1
> DOM, NUM-Z, 2, 2, 1, 2
> DOM, NUM-ANODES, 2, 2
> DOM, TAPPING-GAP, NO
> DOM, FREEZE, YES
> DOM, FREEZE-TOE, 0.000000E+00
> DOM, LORENTZ, ON
> DOM, AUTOMESH, ON
> DOM, INTERFACE, SLOPE_IN_Y
> OBJ, NAME, INIAIR
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 4.500000E-01
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 1.000000E-01
> OBJ, GEOMETRY, cube14
> OBJ, VISIBLE, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, OPAQUE, 50
> OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls
> OBJ, NAME, INIELEC
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E-01
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 2.500000E-01
> OBJ, GEOMETRY, cubet
> OBJ, VISIBLE, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, MATERIAL, 52,Bath
> OBJ, NAME, WALL1
> OBJ, POSITION, 1.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 5.500000E-01
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, WALL2
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 4.000000E+00, 5.500000E-01
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, WALL3
> OBJ, POSITION, 0.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 0.000000E+00, 5.500000E-01
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, WALL4
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 0.000000E+00, 5.500000E-01
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, PREF1
> OBJ, POSITION, 4.000000E+00, 1.600000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 4.000001E-01, 9.999999E-02
> OBJ, GEOMETRY, default
> OBJ, TYPE, CELLTYPE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 49
> OBJ, OPAQUE, 50
> OBJ, PATCHES, PREF1
> OBJ, NAME, PREF2
> OBJ, POSITION, 4.000000E+00, 1.600000E+00, 2.000000E-01
> OBJ, SIZE, 1.000000E+00, 4.000001E-01, 3.000002E-02
> OBJ, GEOMETRY, default
> OBJ, TYPE, CELLTYPE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 49
> OBJ, OPAQUE, 50
> OBJ, PATCHES, PREF2
> OBJ, NAME, ANOPTNTL
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 5.500000E-01
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, default
> OBJ, VISIBLE, NO
> OBJ, TYPE, USER_DEFINED
> OBJ, PATCHES, ANOPTNTL
> OBJ, NAME, CATHCURR
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, default
> OBJ, VISIBLE, NO
> OBJ, GRID, N,N,Y
> OBJ, TYPE, USER_DEFINED
> OBJ, PATCHES, CATHCURR
> OBJ, NAME, ANOD1
> OBJ, POSITION, 1.000000E+00, 4.000000E-01, 2.600000E-01
> OBJ, SIZE, 3.000000E+00, 1.200000E+00, 2.900000E-01
> OBJ, GEOMETRY, default
> OBJ, TYPE, ANODE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 238
> OBJ, OPAQUE, 50
> OBJ, NAME, ANOD2
> OBJ, POSITION, 1.000000E+00, 2.400000E+00, 2.600000E-01
> OBJ, SIZE, 3.000000E+00, 1.200000E+00, 2.900000E-01
> OBJ, GEOMETRY, default
> OBJ, TYPE, ANODE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 238
> OBJ, OPAQUE, 50
> OBJ, NAME, ANOD3
> OBJ, POSITION, 6.000000E+00, 4.000000E-01, 2.600000E-01
> OBJ, SIZE, 3.000000E+00, 1.200000E+00, 2.900000E-01
> OBJ, GEOMETRY, default
> OBJ, TYPE, ANODE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 238
> OBJ, OPAQUE, 50
> OBJ, NAME, ANOD4
> OBJ, POSITION, 6.000000E+00, 2.400000E+00, 2.600000E-01
> OBJ, SIZE, 3.000000E+00, 1.200000E+00, 2.900000E-01
> OBJ, GEOMETRY, default
> OBJ, TYPE, ANODE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 238
> OBJ, OPAQUE, 50
> OBJ, NAME, INTERFAC
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E-01
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, cube
> OBJ, VISIBLE, NO
> OBJ, TYPE, NULL
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, OPAQUE, 50
> OBJ, NAME, FREESURF
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 4.500000E-01
> OBJ, SIZE, 1.000000E+01, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, cube
> OBJ, VISIBLE, NO
> OBJ, GRID, N,N,N
> OBJ, TYPE, NULL
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, OPAQUE, 50
> OBJ, NAME, FREEZE1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 4.000000E+00, 2.000000E-01
> OBJ, GEOMETRY, wedge
> OBJ, ROTATION24, 4
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 192
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, FREEZE2
> OBJ, POSITION, 9.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 4.000000E+00, 2.000000E-01
> OBJ, GEOMETRY, wedge
> OBJ, ROTATION24, 2
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 192
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, FREEZE3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 4.000000E-01, 2.000000E-01
> OBJ, GEOMETRY, wedge
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 192
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, FREEZE4
> OBJ, POSITION, 0.000000E+00, 3.600000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 4.000000E-01, 2.000000E-01
> OBJ, GEOMETRY, wedge
> OBJ, ROTATION24, 3
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, USER
> OBJ, COLOR-VAL, 192
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
************************************************************
MENSAV(S,MENU,ESTR,ESTR)
STOP