>>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>

   This demonstration shows  how  the  two-phase  option  of
   PHOENICS,  library  case  W249,  is  extended by PLANT to
   permit simulation  of  the  relative  movement  of  water
   (rain, or snow) and air. ( Choose the non-shelter variant
   when prompted).
 
   PLANT is  used  to  introduce  patchwise discontinuity in
   interphase friction distribution:  the particulate  phase
   fragments  in the bottom part of the domain are much less
   influenced by air,  than the ones of the upper part.  The
   inspection of PHOTON pictures show what happened near the
   ground:  the gravitational sedimentation is increased and
   the air becomes purer.
 
    PLANT information :
     * Data input groups used: 10
     * Ground groups planted : 10
     * Headings used  : PRPT??
     * Functions used : None
     * Commands used  : PLACE
 
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
 LOAD(W429)
TEXT(Rain and snow : effect of interphase friction.
IF(.NOT.ROOF) THEN
LSWEEP=400
mesg(Reduce interface friction at the bottom half of domain ? (y/n)
readvdu(ans,char,y)
if(:ans:.eq.y) then
 LG(1)=T
 NAMSAT=MOSG
 CFIPS=GRND

  PLANTBEGIN  
   INTFRC=9.81*(:RHO2:-:RHO1:)/(:FALLVEL:)*MASS1*LIQ
   PLACE(0.0,24.0,5.0,10.0) /LG(1)
   INTFRC=9.81*(:RHO2:-:RHO1:)/(7.*:FALLVEL:)*MASS1*LIQ
   PLACE(0.0,24.0,0.0,5.0) /LG(1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
   The above  two  statements  provide  what is required for
   none-shelter (ROOF=F) situation:  the first,  apllied  to
   the  upper  domain part,  extended from 5.0 m up to 10.m,
   has got the coefficient of 7 times larger than the second
   one  has.  The  lattter is applied for the bottom part of
   the domain extended up to 5.0 m  above  the  ground.  The
   both are activated if logical switch LG(1) equals T.
 
   Please note the use of PLACE command with  the  arguments
   being  physical  distances  rather  than  grid extents in
   terms of cell numbers.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  PLANTEND

endif
ENDIF
tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=2.000E+02; EX(AU1)=2.000E+00; EX(BU2)=1.979E+00
EX(CV1)=5.349E-03; EX(DV2)=2.511E+00; EX(LIQ)=1.079E-05
EX(EPOR)=1.000E+00; EX(NPOR)=1.000E+00
 LIBREF=606
STOP