77 IPrI+1 



00 79 J=IP»6 

 IFCGGUDCJ) )79»70,78 

 7« SUM = SUM+P(L.I»J)*CnsrCWNCL)*OCI. J)*COSf CTHZRO(K)-PSICI» J))) 

 X -P(L.J»I)*SINF(WNCL)*0(I#J)*COSFcTHZRO(K)-PSi(I»J))) 



79 CONTINUE 

 76 CONTINUE 

 75 TEM(K) = (A\/EP(L)+7.0*SUM)/TE'^MS 



EmaX(L} = TEm( 1 •) 



THMA<(L)=RTU*CTWnPl-THZRO(n) 



no 7a K=2»NZERDS 



TE(EMAXCL)-TE^(K))73»7a.7/i 



73 EMAXCL)=TEM(K) 

 ThMAX(L)=RTD*CTWOPi-THZRO(K)) 



thmax is bearing From magnITic north 



74 continue 

 SUM=0.0 



do 80 I=l»5 

 IFCGOODCI) }80>80»8l 



81 IP=I+1 



00 83 J=IP»6 

 IFCGOQDCJJ )93#83*82 



82 SUM=SUM+PCL/ I.J)*PCL»I»J)+PCL»J»I)*P(L»J»I) 



83 CONTINUE 



B-8 



