226 



54050 

 54100 

 54150 

 54200 

 54250 

 54300 

 54350 

 54400 

 54450 

 54500 

 54550 

 54600 

 54650 

 54700 

 54750 

 54800 

 54850 



54 900 

 54950 

 55000 

 55050 



55 100 

 55150 

 55200 

 55250 

 55300 

 55350 

 55400 

 55450 

 55500 

 55550 

 55600 

 55650 

 55700 

 55750 

 55800 

 55850 

 55900 

 55950 

 56000 

 56050 

 56100 



56 150 

 56200 

 56250 

 56300 

 56350 

 56400 

 56450 

 56500 

 56550 

 56600 

 56650 

 56700 

 56750 

 56800 

 56850 

 56900 

 56950 

 57000 



SUBROUTINE BLOCK( n, X , maxb 1 k . bl k ) 

 DIMENSION X(4100) 

 INTEGER BLK.STOP 



STOP = (N-BLK) 

 NM = 1 



DO 100 IJ=2,stop,blk 



SUM=0.0 



MJ=IJ+BLK-1 



KJ=Id 



DO 130 KKd=kj ,mj 

 SUM=SUM + X(kkj ) 



X(nm) = SUM/BLK 



NM-NM+1 

 Cont 1 nue 



SUM=0.0 



DO 140 kkj = 1 J .n 



SUM = SUM + X(kkj ) 

 X(nm) = SUM/(n-mJ ) 



MAXBLK = nni 



Return 

 End 



SUBROUTINE SURF2(ab.d1 . cl2 , d3 , d4 , angi , ang2 , ang3 . ang4 , Icut, 

 wvnr.blk.or lent ) 



DIMENSION ab(2050) .d1(2050) ,d2(2050) ,d3(2050) ,d4(2050) . 



1 ang 1 ( 400 ) . ang2 ( 400 ) , ang3 ( 400 ) . ang4 ( 400 ) , 



2 wvnr(2050) 

 INTEGER blk 



C Estimates direction using only absolute gage and one 

 C differential channel at a time. 



pi = 3 . 14159 



DO 300 J =2. icut 



z = 1./wvnr(j) * sqrt(d1(j)/ab(j)) 

 IF(Z. LE . 1 . )ANG1( j )=0RIENT-acos(Z)*180./pi 

 Z = 1./wvnr(J) * sqrt(d1( j )/ab( j )) 

 IF(Z.LE. 1 . )ANG2( j ) =ORI ENT-as i n( Z ) * 180 . /p1 

 Z = 1./wvnr(j) * sqrt(d3( j )/ab( j ) ) 

 IF(Z.LE. 1 . )ANG3( j ) =0RI ENT-acos( Z ) * 180 . /p i 

 Z = 1./wvnr(j) * sqrt(d4(j)/ab(j)) 

 IF(Z.LE. 1 . )ANG4( j ) =0RIENT-as in( Z ) * 180. /pi 



IF(ang1( j ) . 1 t O. )ang1( j )=ang1( j ) +360.0 

 IF(ang2( J ) . 1 t O. )ang2( j )=ang2( J ) + 360.0 

 IF(ang3(J ) . 1 t .0. )ang3( j )=ang3( j ) +360.0 

 IF(ang4(J ) . 1 t .0. )ang4( j )=ang4( j )+360.0 



