PHOTON USE
  p
 
 
 
  gr out x 1 col 10
  text
  1
  1
  The particle is bouncing in a gravity field.
   0.34606E+03 0.26380E+04 CR
  text
  1
  1
  The loss of total energy is due to
   0.34606E+03 0.25077E+04 CR
  text
  1
  1
  the drag force.
   0.34606E+03 0.23775E+04 CR
  geometry read
  T00001
  upause 2
  text
  3
 
  Press  to end
   0.68570E+03 0.14265E+04 CR
  pause
  enduse
 
TEXT(Particle bouncing in grav field: G721
TITLE
    GROUP 3. X-direction grid specification
XULAST=0.001
NY=5
GRDPWR(Y,NY,1.0,1.0)
    GROUP 5. Z-direction grid specification
NZ=10
GRDPWR(Z,NZ,1.0,1.0)
    GROUP 7. Variables stored, solved & named
STORE(P1,V1,W1)
    GROUP 11. Initialization of variable or porosity fields
FIINIT(V1)=0.0;FIINIT(W1)=0.0;FIINIT(P1)=0.0
    GROUP 13. Boundary conditions and special sources
     * Inlet Boundary Condition
PATCH(GXINLET,LOW,1,NX,1,NY,1,1,1,LSTEP)
     * Outlet Boundary Condition
PATCH(GXOUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,LSTEP)
    GROUP 15. Termination of sweeps
     * Number of Sweeps
LSWEEP=1
    GROUP 23. Field print-out and plot control
TRACE=F
#GENTRA
  *-----------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 30 - ISOTHERMAL_PARTICLES
GPTYPE=30
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0.0;GGRAY=-0.01;GGRAZ=0;GBUOYA=F
  * Stochastic model of turbulence
GSTOCH=F
  * Data for isothermal particles
GDRAG=GRND1
  *-------------------------------------------------------
  * GENTRA GROUP 2: Boundary conditions for particles
  *-------------------------------------------------------
  * Inlet-data file-name
GINFIL=Q1
  * Use this space for inlet data id GINFIL=Q1
  
  *__YP__ZP__VP__WP__DIAM_DENSTY_FRATE
   0.9 0.0 -0.001 0.01 0.01 1000.   1.0
  
  * Wall treatment, and rest coefficient if appropriate
GWALLC=3;GWREST=1
  * Porosity threshold
GPOROS=0
  *-------------------------------------------------------
  * GENTRA GROUP 3: Numerical controls
  *-------------------------------------------------------
  * 1st GENTRA sweep; frequency of calls
GSWEP1=1;GSWEPF=1
  * Maximum Lagrangian time-step
GDTMAX=0.1
  * Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=-1000;GTIMMX=100000
  *-------------------------------------------------------
  * GENTRA GROUP 4: Output controls
  *-------------------------------------------------------
  * Trajectory-file identifier, history-file identifier
    and frequency for output
GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1
GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: Provisions for the GENTRA-Earth run
  *-------------------------------------------------------
  * Storage for interphase sources of momentum
STORE(MOMY,MOMZ)
  * Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,5,FIXFLU,GRND)
COVAL(GENPAT,7,FIXFLU,GRND)
  * Relaxation of interphase sources
RELAX(MOMY,LINRLX,7.0000E-01)
RELAX(MOMZ,LINRLX,7.0000E-01)
  * Activation of GENTRA GROUND-station
L(GENSET
  * ----------------------------------------------------------
  * GENTRA GROUP 7: GENTRA test
  * ----------------------------------------------------------
   * GENTRA auto test  RG(51)--RG(58) are reserved
RG(51)= 1.56E+02;RG(52)= 1.57E+03;RG(53)= 9.90E-02
RG(54)= 5.00E-04;RG(55)= 4.53E-01;RG(56)= 1.00E+00
RG(57)= 1.00E+16;RG(58)= 1.00E+01