DIMENSION ITME(25I .IBOyR{25) tIBOyA(25l . IDRGRI2SI .IDRGA125I 

 Dl WENS I ON XRB(25] ,YRB(25I i XRD ( 25 I .YRO( 25 ) ■eUOYAI25l .DROGAI25) 

 DIMENSION XR(25) .VRI25I ,T1ME(3I . IXHR ( 3 I . IMIN ( 3 ) 

 CNUMBER OF READINGS 



READ(2 »1 )N. INDAY 



1 FORMAT ( 12.121 

 C READ INPUT DATA 



DO 10 1=1 iN 



10 READ (2. 2) ITME( I) .IBOYRl I 1 .IBOYA( I 1 . 1 DROR I I I .IDRGA( 1 I 



2 FORMATI I4.1X.2I I5.I3.1XI I 

 C OUTPUT INPUT DATA 



INTHR=ITMEI 11 



IHR«INTHR 



WRITE(3.3) 



3 FORMAT(<.ax.'iHBUOY.3X.'.HBUOY.5X.'.HOROG.3X,4HDROG./.35X.7HTIMEI2l.5X 

 1.5HRANGE.2X.5HANGLE.'>X.5HRANGE.2X.5HANGLE./I 



DO l* 1 = 1. N 

 IF( ITMEI I l-IHRIll.12.12 



11 IHR"ITMEIII 



DTKE= ITMEI 11+1 INDAY+1)»1.E0'> 

 GO TO 13 



12 DTME=ITME( I 1 + INDAY-l.EO'. 



13 WRITEO.'.IDTME.ISOYRI I) .IBOYAl I) . lORGR I I ) . 1 DRGAI I ) 



4 FOPHATI35X.F8.0.4X.I5.3X.13.5X.I5.3X.I3./I 



14 CONTINUE 



C CONVERT INPUT ANGLES TO RADIANS 

 DEGRO=57. 29578 

 DO 20 I=1.N 



BUOYAI I IsIBOYAI I 1/DEGRD 

 20 0ROGA( I )=IDRGAI I 1/DEGRD 

 C CONVERT INPUT VECTORS TO COMPONENT FORM 

 DO 30 1=1. N 



XRBl I 1-IBOYRl I I»SIN1BU0YA( I ) 1 

 YRBl I I'IBOYRI I)»C0S1BU0YA1 I 1 1 

 XROl I 1=IDRGRI II«SIN(0R0GA( 1 1 1 

 YRDI I l = [0RGR( I I'COSIDROGAI I ) I 

 XRI I l=XRD( I l-XRBI I I 

 30 YRI 1 l»YRDI n-YRB( I I 

 WRITE (3.7) 



7 FORMAT I 14X.lHX.6X.lHY.21X.lHX.6X.lHyi 



WRITE 13.8) 1 IT^EI I 1 .XRBl I ) .YRBl I I .XRDII > lYRDl I l.I = l.Nl 



8 FORMAT 1 I6.5x.F9.2.2X.F9.2.19X.F9.2.2X.F9,2./ ) 

 C WRITE OUTPIT HEADINGS 



WRITE13.5I 



5 FORMATI1X.1H1.15X.7HTIME1 I ) . / . 15X .7HT IME 1 2 ) .6X .5HVEL0C .5X. 5HANGLE . 

 1/ ) 



C COMPUTATION OF VELOCITY 



NN=N-1 



DO 126 K-l.NN 



XDELT = XR1K-.1 )-XRUl 



YDELT = YR1K, + 1 )-YR(x) 

 XMAGT=SQRTI IX0ELT**2 1+(YDELT»*2) ) 

 C COMPUTE ANGLE INFORMATION 



IFIYOELTl 119.115.119 



115 IFIXDELTI 117.116.118 



116 ANGLE = 0. 

 GO TO 123 



117 ANGLE = 270. 

 GO TO 123 



118 ANGLE = 90. 

 GO TO 123 



119 ANGLE » ATAN(XDELT/YDELT)»DEGRD*.5 

 IF (YDELTl 120.999.121 



120 ANGLE = ANGLE + 180. 

 GO TO 123 



121 IFIXDELT) 122.999.123 



122 ANGLE = ANGLE * 360. 



123 C0NV»91. 440183/60. 

 IXMRI 1 )=ITME1X ) /lOO. 



1XHR121 = ITMEIK-H)/100. 

 IMIN11)»ITMEU )-( IXHRl 1 )»1.F02 1 

 IMIN12)-ITMElx*l)-lIXHRI2)»l.Fn21 

 199 I0EL=IXHR12 )-IXHRl 1 ) 



IFl IDEL)202.201.200 

 2 00 DELT=I 1 IMIN121*601-IMIN111 1 * I 60* I I DEL-1 I 1 

 GO TO 203 



201 DELT=1MIN12I-IMIN11] 



GO TO 203 



202 IXHR12 )»IXHRt21+24 

 GO TO 199 



203 VEL0C=1XMAGT/DELT1«C0NV 

 COMPUTE DATE-TIME GROUP 



TIMEll) = ITME(K) *■ INDAY*1,E04 

 TIME12I = ITME1(C+1I + INDAY»1.E04 

 WRITE (3.61 TIMEdl.TIMECl .VELOC.AMGLE 

 6 FORMAT ( 15X.FS.0./.15X.F8.0.5X.F8.2.5X.F5,01 



126 CONTINUE 



999 STOP 

 END 



Figure 71. IBM 1130 FORTRAN program used for drogue velocity and true bearing calculation. 



70 



