TALK=T;RUN(1,1) 2ph580

L($W574)
TEXT(2 PH BACKWARD STEP,WITH GRAVITY    :W580
TITLE
  DISPLAY
  This case is like W574; but the density ratio is increased to
  1000:1, and gravity causes additional separation of the two phases
  ENDDIS
RHO2=1000.0;VALUE(INLET,P2,RHO2*R2IN*UIN)
COVAL(OUTLET,P2,1.0*RHO2/RHO1,0.0)
PATCH(GRAVITY,PHASEM,1,NX,1,NY,1,1,1,1)
COVAL(GRAVITY,V2,FIXFLU,-9.81*(1.0-RHO1/RHO2))
  photon use
  msg press return to use q1ear
  pause
  use q1ear
  enduse