PHOTON USE
  p
 
 
 
 
  gr ou x 1
  MSG Velocity vectors
  vec x 1 sh
  msg
  msg Press return to plot pressure contours
  pause
  cont p1 x 1 fil;.01
  msg
  msg Type e to End
  ENDUSE
 
#$B515
TEXT(POT FLOW OVER ELLIPSE;ALGEBR. GRID: B516
TITLE
  CHAR(ANSW)
mesg(Press return to continue
readvdu(answ,char,y)
  DISPLAY
    This case is similar to 514 except that the grid file now
  represents a half elliptic cylinder of aspect ratio 2:1 with its
  long axis in the z-direction. The grid is generated by
  interpolation with constant z-coordinates and is non-orthogonal
  ENDDIS
 
IZ1=5;IZ2=17
REAL(BB,CC);AA=0.4;BB=0.2
   ** Reset corner points for ellipse.
GSET(P,C,0,1,6.0000E-01)
GSET(P,D,0,1,1.4000E+00)
GSET(P,G,0,0,1.4000E+00)
GSET(P,H,0,0,6.0000E-01)
   ** Reset Y grid lines.
GSET(L,L1,A,B,NY,1.5)
GSET(L,L5,F,E,NY,1.5)
   ** Reset Z grid lines over ellipse surface
GSET(L,L3,C,D,12,S1.77)
   -- Defining a curve for L7.  Again, using S1.35... to approximate
GSET(V,ELLI,S,G,SPLINE)
DO II=1,11
CC=3.1416*:II:/12
GSET(V,0.0,BB*SIN(CC),1.0+AA*COS(CC))
ENDDO
GSET(V,ELLI,E,H)
 
GSET(L,L7,G,H,12,S1.35CRVELLI)
   ** Set Z grid lines of remaining frame
GSET(L,L2,B,C,4,1.0)
GSET(L,L4,D,E,4,1.0)
GSET(L,L6,F,G,4,1.0)
GSET(L,L8,H,A,4,1.0)
   ** Define frame.
GSET(F,F1,A,-,B,C.D,E,-,F,G.H)
   ** Match frame on to I1
GSET(M,F1,+J+K,1,1,1,TRANS)
   ** Copy I1 to I2
GSET(C,I2,F,I1,+,1,0,0)
NONORT=T