MAIN          OS/360  -ORTRAN H EXTENDED


ISN 0065


ISN
ISN ISN
ISN
ISN ISN
ISN
ISN ISN
ISN ISN
ISN ISN
ISN
ISN
ISN
ISN ISN
ISN ISN
ISN


0066 0067
0068 0069 0071 0072 0073
0074 0075
0076 0077 0078
0079
0080 0081 0082 0083
0084
0085 0086 0087


   ISN 0088 ISN 0089 ISN 0090 ISN 0091 ISN 0092 ISN 0093 ISN 0094 ISN 0095 ISN 0096 ISN 0097 ISN 0098 ISN 0099 ISN 0100 ISN 0101 ISN 0102 ISN 0103 ISN 0104 ISN 0105

FIGURE 4.12 -


          CALL LINE (DV,-DXDZ)
    C THIS PORTION WILL CONSTRUCT ARROWHEADS ON THE VECTORS WHICH ARE IN
    C  THE PLANE OF VECTOR-Z AXIS, 30% OF VECTOR LENGTH, AND 30 DEGREE
    C  INCLUDED ANGLE.  TANGENTOF 1/2 ANGLE{ 15) IS** 0.26795
          A=DK*D X+D Z*D Z
          B=SQRT (A*DX*DX+A*Al
          C=O
          IF(3.EQ.O.O0) C=1.
          VUI=A/( B+C)
          VU2= DX*DY/(B+Cj
          VU3=-DY*DZ/( B+C)
          RI=-.3*DY
          R 2=- * 3* (-D X)
          R3=-.3*DZ
          D=0.3*( SQR.T(A+DX*DX) ) *0,26795
          P1X=RI+VUI *D P1Y=R2+VU2  D PIZ=R3+VU3*D
          P2X=RI+VU1 *(-D,)
          P2Y=R2 +VU2 *(-D) P2Z=R3+VU3*(-D)
          CALL LINE (PIX.PIY,PIZ)
          CALL MOVE (-PlX,-PIYl,-PIZ)
          CALL LINE (P2X9P2YvP2Z)
      5   CONTINUE
    C THIS S--TION WILL PLOT 5 CUTS AND A VIEW OF ALL PLANE VELOCITY VECTORS
    C FOR BOT+ THE R-THETA PLANE AND THE X-R PLANE
    C DRAW PLOT OF R-THETA CUTS
    6     CALL FACT (1.0)
          CALL CLEAR. (4,5.9)
          CALL XY
          CALL ORIG (18.6915.2)
          CALL TXSIZ (0.0875)
          CALL TXPLT (-1.6-1.844,0.+1lO0&2@*)
          CALL TXPLT (-1 .6,-l,266,0., 1,
          CALL TXPLT (-1.6,-0.588,0.,+I,'0.4@*)
          CALL TXPLT (-1.6,1.156,0,9+I1,C*2@') CALL TXPLT (-.6,1.73490.*, Y/R@") CALL TXPLT (-1.6,2.312.0.i+1,0.4°) CALL TXPLT (0.*02.70,O.,Oi'X=2-10@) CALL TXPLT (3.0,2.70,0.,09'X=10-18@') CALL TXPLT (6*0,2.7090#s0*'X=18--26')
          CALL. TXPLT (0.0,-030,0.,0,'X=26-34@*)"
          CALL TXPLT (3,0,-0*30,09.090X=34-42@0).
          CALL TXPLT (-6.0,-030q,0,sX=ALL@').
          DO 10 J=192
continued


LEVEL 2",1 ( JAN 75 )


DATE 7,6,251/15