TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 10/02/03
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
An In-Form (MOVOB...) statement prescribes the motion of a
submarine in a density-stratified ocean, the latter also being
defined by an In-Form statement.
The submarine is represented by a rectangular so as to allow use
of a coarse grid.
ENDDIS
PHOTON USE
p
m1
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 1 time step
upause 2
p
m2
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 2 time step
upause 2
p
m3
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 3 time step
upause 2
p
m4
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 4 time step
upause 2
p
m5
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 5 time step
upause 2
p
m6
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 6 time step
upause 2
p
m7
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 7 time step
upause 2
p
m8
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 8 time step
upause 2
p
m9
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 9 time step
upause 2
p
m10
up z; vi -1 -1 1
gr ou x 1; gr ou y 1; gr ou z m
con rho1 x 1 fi;.0001
vec y 1 sh; vec z m sh; red
surf obid x .99; surf obid y .99; surf obid z .99
msg Velocity vectors at 10 time step
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 365
************************************************************
Group 1. Run Title
TEXT(MOFOR by In-Form: Underwater navigation )
INFORM1BEGIN
Initial data of problem
-----------------------
** Height of the pycnocline, m.
REAL(HIG); HIG=140.
** Parameter that controls the thickness of the pycnocline, m.
REAL(ETH); ETH=20.
** Density near the bottom, kg/m^3
REAL(RHOL); RHOL=1027.
** Density near the surface, kg/m^3
REAL(RHOH); RHOH=1023.
** The size of the simulation area, m.
REAL(XLAST,YLAST,ZLAST); XLAST=1000.; YLAST=50.; ZLAST=200.
Auxiliary data
--------------
** Sizes of submarine, m.
XSIZ=100.; YSIZ=10.; ZSIZ=20.
** Start position of submarine, m.
XPOS=0.; YPOS=0.; ZPOS=HIG-0.5*ZSIZ
INFORM1END
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,3.000000E+02,10)
* Modify regions
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,50,10,25,1.000000E+02)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(148) =RHO1 ; NAME(149) =OBID
NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 )
* Stored variables list
STORE(PRPS,OBID,RHO1)
store(#fxu,#fxv,#fxw)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(u1 ,Y,Y,Y,N,N,Y)
SOLUTN(v1 ,Y,Y,Y,N,N,Y)
SOLUTN(w1 ,Y,Y,Y,N,N,Y)
************************************************************
Echo InForm settings for Group 7
INFORM7BEGIN
** Definition of the VR moving objects by In-Form
** Velocity of submarine movement, m/c.
char(usubm,times)
usubm=3.
times=tim
(MOVOB of BLOCK is POS(:times:*:USUBM:&0&0&0&0&0))
store(obol)
(stored var obol is old(obid))
INFORM7END
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
SETPRPS(1, 0)
ENUL = 1.792000E-06
ENUT = 1.000000E-02
DRH1DP = 1.000000E-05
DVO1DT = 3.410000E-03
EMISS = 0.000000E+00 ;SCATT = 0.000000E+00
************************************************************
Echo InForm settings for Group 9
INFORM9BEGIN
** For the stratification of a water.
CHAR(FORM)
FORM=:RHOL:+0.5*(TANH((ZG-:HIG:)/:ETH:)+1.)*(:RHOH:-:RHOL:)
(PROPERTY of RHO1 is FORM)
INFORM9END
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = 0.000000E+00
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
No PATCHes used for this Group
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 60
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
CSG1 ='M'
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
SPEDAT(SET,MOFOR,VELCON,C,1.E6)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
ISG62 = 0
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
store(prps,vpor,tmp1)
EX (P1 ) = 6.683000E+00 ;EX (U1 ) = 1.1570000E-01
EX (V1 ) = 5.751000E-03 ;EX (W1 ) = 9.265000E-03
EX (RHO1) = 1.026000E+03 ;EX (OBID) = 1.280000E-02
EX (PRPS) = 0.000000E+00 ;EX (OBOL) = 1.536000E-02
EX (#FXW) = 1.175000E-11 ;EX (#FXV) = 1.226000E-11
EX (#FXU) = 2.304000E-02
************************************************************
Group 21. Print-out of Variables
OUTPUT(V1 ,Y,Y,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='M'
************************************************************
Echo InForm settings for Group 25
GVIEW(P,0.000000E+00,-9.995500E-01,2.999551E-02)
GVIEW(UP,0.000000E+00,2.999551E-02,9.995500E-01)
> DOM, SIZE, :XLAST:, :YLAST:, :ZLAST:
> DOM, MONIT, 4.650000E+02, 2.400000E+01, 8.653362E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_Y_1, 10, 1.500000E+00
> GRID, RSET_Z_1, -25, 7.000000E-01
> OBJ, NAME, WESTFREE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E+01, 2.000000E+02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, GROUND
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+03, 5.000000E+01, 0.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, EASTFREE
> OBJ, POSITION, 1.000000E+03, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E+01, 2.000000E+02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, NORTFREE
> OBJ, POSITION, 0.000000E+00, 5.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.000000E+03, 0.000000E+00, 2.000000E+02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, SOUTFREE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+03, 0.000000E+00, 2.000000E+02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, HIGHFREE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+02
> OBJ, SIZE, 1.000000E+03, 5.000000E+01, 0.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, BLOCK
> OBJ, POSITION, :XPOS:, :YPOS:, :ZPOS:
> OBJ, SIZE, :XSIZ:, :YSIZ:, :ZSIZ:
> OBJ, CLIPART, cubet
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, -1
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, INI_PRESS, 0.000000E+00
> OBJ, SCAL_FIXF, 0.000000E+00
STOP