L(116)
TEXT(Comp Flow Over Double-Wedge Airfoil 
TITLE
     Here the compressibility is introduced by the setting
  of POTCMP to T. Necessary information is carried through to
  GREX3 by way of PORIA and PORIB. More sweeps are needed because
  the compressibility effect is introduced by updating the
  porosities at each sweep.
     Only the high porosity is modified in GREX3, by reference
  to the W velocity alone. The validity is thus restricted
  to "slender-body" flows in which the W velocity greatly
  exceeds the other components.
POTCMP=T
REAL(GAMMA,SONVEL);GAMMA=1.4;SONVEL=4.0
PORIA=(GAMMA-1.0)/((GAMMA+1.0)*SONVEL**2)
PORIB=1.0/(GAMMA-1.0)
LITER(POT)=-100;RESREF(POT)=1.E-20;LSWEEP=200;NPLT=1;RESFAC=0.01
RELAX(HPOR,LINRLX,0.5)
IYMON=NY/2;IZMON=NZ/2;TSTSWP=-1
IZPRF=18;IZPRL=22;IYPRL=10;NZPRIN=1;NYPRIN=1
OUTPUT(POT,Y,Y,Y,Y,Y,N);OUTPUT(HPOR,Y,N,N,N,Y,N)
OUTPUT(W1,Y,N,N,N,Y,N);OUTPUT(V1,Y,N,N,N,Y,N)