8RK2DHU III I II mil 1 1 II 1 1 II HI 1 1 1 II nil II niiiiinit II I II I III n II 1 1 II 1 1 II iiin 



RUNT VERSION FEB 74 B 17:1? Qst^milb 



PROGRAM bRK?D (INPUT. 



C***LATbST REVISinN ***** 2 



C***PRaGRAM BRK?D CO^IPUTEi 



C CYLINDER ON OR "^EAR T 



C INFINITE DEPTH 



3 COMMON RI12(?5. 25).RK 



1FI(25,6), RI(25.25), 



2RMU(3>3.10) . RLAM(3.3 



3 DELW(2i.6.10) . XQK 



3 COMMON/aNE/X(25).Y(25 



l.FEINt2£). ^TT'Jdt/). 

 3 CQMM0N/0NE2/CC3(25).S 



3 COMMON /TWO/ N.NNW, N 



3 COMMON/THREE/ WAVELd 



3 C0MMCN/SIX/i<N{5),CN('j 



3 COMMON ISt^Etl AREA. 



3 COMMON / EIGHT/L8LMU( 



1LBL(10»3). DEGO.IO) 

 3 C0MM0r4 /NINE/ L8L-<AR( 



3 ■ CGMMQN/TeN/r)ELT(2.3). 

 3 REAL K 



3 DATA XN/.2635603197ia 



1 7.0^581C005B'5<5.1 

 3 DATA CN/. 521755610563 



1 .0036U75>Jb79?2. 

 3 PlE'ATAN2(').,-2. ) 



7 TP«2.*PIE 



iO GAMMA'O. 57721566 



C***BFGIN READING INPUT OAT 

 12 3000 FORMAT (6410) 



C*****kEAD LABLES FOR PRIiO 

 12 READ 3000. (((LSLMUd 



36 READ 3000. (((LSLA.KI 



63 READ 3000» ( (LBLFB 



103 READ 3000. ((LBLHWa 



123 READ 30UJ. (( LBL 



1<.3 READ 3C00. ((L3LSAR 



163 READ 3000. ((LBLH^R 



203 READ 3000, (LBLRd.L) 



220 READ 20. TITLE 



226 20 FORMAT (8A1CJ) 

 226 PRINT 30. TITLE 



23't 30 FORMAT (HI. 3Aln///) 

 234 READ 50. N. NW, NwAVE 



256 50 FORMAT (5110. 615) 



C N « NUMBER OF <:i 



C THE HJLL 



C NW • NUMBER 3F P 



C TO Bt COMPUTE 



C WAVE HEIGHT 4 



C WHICH IS Pe>?F 



C NWAVtL=N UMBER 



C BE PERFOR'iEO 



C ISYM • 1 FOR SYM 



C « ANYTHING 



OUTPUT. PUNCH. r4Pe5«TNP')T,T4PF6 -OUTPUT) 

 7 AUGUST 19TJ 



THE FIRST-ORDER ?i=SPTNSF OF AN OSCILLATING 

 HE FREE SMRFirF OF AN ITEAL FLUn OF 



56(25,25). PnT(25.'5). 40W( 25, 25 ). PF (??.*:) , 



RJ(25.25). t!K(25,4), »L(?5.4), 

 ,10). F3(3.10), OFLP'CIO) . HWR(2'^,6,10). 

 25,10) 



),)fB(25).YB(?5).ANG(?5).DEL<25).>/V(?5) 

 RNllRM(25.3), jr(5) 

 S3(23) 

 WA7EL. ISYM. TSi<T0. sir. Pie.GAMMA.M.TK.To 



0). wN(io). mL(n).TL 



) 



R. D, ROE, GEF, BTTTLE(3). TITLt(R) 

 3,3.3). L8LAM{^,3.^), LBLFB ( 3, 3 ) . LBLHW-^ ( 7. 3) , 



3.3). L8LHWR(=.3). l'»LR(5.3) 



F3R(^,lO).PH4S(2.lO).F^RNO(2.lo).PH^'>n(2.1o) 



.i..'il3<.a3t'5"?1';7.?.'=9',42'^7710'.), 



2.640B0084<>2''ft/ 



. . 39666681 10B3.. 075 qA244P6F' 17. 



.00002336997?39/ 



A AND PRINTING FCHT CHECK 



OU 



.J 

 • J 

 (J 

 (J 

 (J 

 (J 

 (J 





3) . 





S 





,^), I = 1.3) 









= 





.^). T ■= 1.3) 









B 





>3) 









= 





7) 









■ 





. '0) 









B 





.3) 









a 





.5) 



1.3) 



I. ISYM. ISKTP. LC. J*" 



PAIGHT LINE SFGM^NTS TO 3f USED TO FTT 

 NOTE. THERE "UST =5 F N + 1 OFF-^tT POINTS 

 OINTS JN FREF SUPPACE WHERE WAVE HEIGHT t<; 

 D. THIS IS IN ADDITION TO THE COMPUTATIPM HE 

 .0 WAVELENGTHS ON pTTHPP SIDE OF THE BODY 

 ORMEO AUTOMATICALLY 

 OP WAVELENGTHS 4^ i^HTCH COMPUTATIONS «pF TO 



METRIC SECTION 



FLS6 FOR NON-<;ymmptpic SECTION 



Table D-2. Listing of program BRK2D 



III 



