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.
REAL(XSIZ,YSIZ,ZSIZ); XSIZ=100.; YSIZ=10.; ZSIZ=20.
   ** Start position of submarine, m.
REAL(XPOS,YPOS,ZPOS); 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)
  ************************************************************
  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