```
PHOTON USE
p

gr ou x 1
MSG Velocity vectors
vec x 1 sh
msg
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)
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
```