IF (SENSE SWITCH 3) lOO.lOl 



100 TYPE 103»MAX 



103 FORMAT ( '!fHMAX= » I 4 ) 



101 IF (MAX-500) 399,15.15 

 399 CALL MOVE (X»Y,A.FK) 



IF (CXY) 15»15,396 



396 GO TO (397,395»15) . MIT 

 395 WRITE 200, MAX 



200 FORMAT ( 32HCURVATURE APPROXIMATED FOR MAX =,I4) 



397 ANGLE=A»57. 29577951 

 XXX = X 



YYY = Y 



PCTDIF = ABSF( (C(3)-E(l)-E(2)-E(3) )/C(3! )*100. 



CALL ROUND (XXXtlOO. ) 



CALL ROUND ( YYY.IOO. ) 



CALL ROUND {ANGLE»100.) 



CALL ROUND (PCTDIF»10.) 



PUNCH 1 2 »M AX, XXX » YYY 9 ANGLE. PCTDIF 



12 FORMAT ( I7.2F9.2,F11.2»F10,1) 

 GO TO (3»15) ,NGO 



15 PUNCH 13 



13 FORMAT (12HRAY STOPPED.) 

 RETURN 



END 



8 6 



SUBROUTINE SURFCE (X.Y.A.FK) 



DIMENSION S(3»3).EM(4.3).E(3) ,YVW(3) .CMATt 50.51) .C(4) ,XLABL(12) 



COMMON S.EM,E.YVW,CMAT,C,XLABL.D.TT»CXY.IT.NGO»AMM»ANN,MAX,MIT 



I=X + 1. 



J=Y+1. 



FI = I 



FJ = J 



XL=X+1.-FI 



YL=Y+1.-FJ 



IF (MAX-1) 1.1,4 

 4 IF (ZI-FI) 1.2.1 

 2 IF (ZJ-FJ) 1,3.1 

 1 ZI = FI 



ZJ = FJ 



C(1)=CMAT( I.J) 



C(2)=CMAT( I+l.J) 



C(3)=CMAT( I+1,J+1) 



C(4)=CMAT( I,J+1) 



DO 318 11=1,3 



YVW( II ) = 0. 



DO 318 L=1.4 

 318 YVW(II) = YVWdl )+C(L)*EM(L,II ) ~i 



DO 319 II=1»3 



E(II) = 0. 



RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 

 RAYN 



07 

 08 

 09 

 10 

 11 

 12 

 13 

 14 

 15 

 16 

 17 

 18 

 19 

 20 

 21 

 22 

 23 

 24 

 25 

 26 

 27 

 28 

 29 

 30 



SURFCEOl 

 SURFCE02 

 SURFCE03 

 SURFCE04 

 SURFCE05 

 SURFCE06 

 SURFCE07 

 SURFCE08 

 SURFCE09 

 SURFCElO 

 SURFCEll 

 SURFCE12 

 SURFCE13 

 SURFCE14 

 SURFCE15 

 SURFCE16 

 SURFCE17 

 SURFCE18 

 SURFCE19 

 SURFCE20 

 SURFCE21 

 SURFCE22 

 SURFCE23 

 SURFCE24 

 SURFCE25 



48 



