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
 TEXT(MOFOR by In-Form: 2D motion of 2 objects
libref=360
TITLE 
  DISPLAY
  This input file shows how MOFOR, the moving-frame-of-reference 
  feature, can be used without a .MOF file, when suitable In-Form 
  statements are used to define its motion, as described
  here.
  
  Two cylinders move through a fluid. The first follows a parabolic 
  trajectory from left to right, the second a horizontal trajectory 
  from right to left.

  '(MOVOB' statements dictate these motions of by use of the POS()
  function, of which the arguments set: 
  X, Y and Z coordinates and rotation angles about X, Y and Z axes,
  as functions of the time variable, TIM.

  The variable OBID (= object identifier) marks the positions, for
  ease of display by VIEWER, PHOTON or WinPHOTON.
  
  ENDDIS
  VRV USE
    * Domain scaling factors
    SCALE  1.000000E+00 1.000000E+00 1.000000E+00
    * Settings for current slice
    PROBE  1.070000E+01  3.500000E+00  5.000000E-02; PROBE ON
    SLICE Z
    * View and up directions
    VIEW  0.000000E+00  0.000000E+00  1.000000E+00
    UP  0.000000E+00  1.000000E+00  0.000000E+00
    * View centre
    VIEW CENTRE  1.100001E+01  3.500000E+00  5.000000E-02
    * View size
    VIEW SIZE  1.038916E+01
    * View perspective 
    VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
    DOMAIN ON
    VARIABLE  Velocity
    VARIABLE RANGE  0.000000E+00  1.400000E+01
    VECTOR ON; CONTOUR OFF; SURFACE OFF
    VECTOR REFERENCE  5.000000E+00
    ANIMATE START 1 END 20 INTERVAL 1                                               
  ENDUSE
  photon use
  p
  m1;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m2;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m3;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m4;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m5;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m6;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m7;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m8;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m9;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m10;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m11;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m12;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m13;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m14;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m15;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m16;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m17;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m18;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m19;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  upause 2
  p
  m20;
 
  gr ou z 1
  ve z 1 sh
  surf obid z 1
  enduse
 ************************************************************
 IRUNN   =       1 ;LIBREF =     360
 ************************************************************
  Group 1. Run Title
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.000000E+00,2.000000E+00,20)
    * Modify regions
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,110,35,1,.01)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(149) =OBID ; NAME(150) =PRPS
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(PRPS,OBID)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN
   ** Definition of the VR moving objects by In-Form
char(vel,gravt,times)
vel=10.; gravt=9.81; times=tim
(MOVOB of CYL1 is POS(:times:*:vel:,:times:*:vel:-0.5*:gravt:*:t$
imes:^2,0,0,0,0))                                                   
(MOVOB of CYL2 is POS(-:times:*:vel:,0,0,0,0,0))
store(obol)
(stored var obol is old(obid))
  INFORM7END
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  0)
 ENUT    = 0.000000E+00
 DRH1DP  = 1.000000E-05
 DVO1DT  = 3.410000E-03
 EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
 ************************************************************
  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  =      20
 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
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 CSG1    ='M'
 SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 DISTIL  =    T  ;NULLPR =    F
 NDST    =       0
 DSTTOL  = 1.000000E-02
 EX    (P1  ) =  2.451000E+01 ;EX    (U1  ) =  3.334000E+00
 EX    (V1  ) =  1.756000E+00 ;EX    (OBID) =  3.896000E-02
 EX    (PRPS) =  0.000000E+00 ;EX    (OBOL) =  3.948000E-02
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(V1  ,Y,Y,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      54 ;IYMON  =      18 ;IZMON  =       1
 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
  INFORM25BEGIN
  DEBUG=f
  MOFOR=t
  formula=f
  INFORM25END
 
 GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        2.200000E+01, 7.000000E+00, 1.000000E-01
> DOM,    MONIT,       1.070000E+01, 3.500000E+00, 5.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
 
> OBJ,    NAME,        RELIEF
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 2.000000E-01, 1.000000E-01
> OBJ,    CLIPART,     cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PRESSURE_RELIEF
> OBJ,    PRES_RELIEF,   1.000000E+03, 0.000000E+00
 
> OBJ,    NAME,        CYL1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cylinder
> 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
 
> OBJ,    NAME,        CYL2
> OBJ,    POSITION,    2.100000E+01, 3.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cylinder
> 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
 
> OBJ,    NAME,        WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.200000E+01, 0.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
libref=360
STOP