K1(KST)=JST 



K2(KST)=JJ 



K3(KST)=IPNUM 



JST=JJ 



SUM=0.0 



IPNUM=X(JST) 



IF(KST.EQ.l) GO TO 116 



MST=KST-1 



IF(IAVE.NE.K3(MST)) GO TO 116 



K2(MST)=K2(KST) 



GO TO 35 



116 KST=KST+1 

 IF(KST.GT.LTPS) GO TO 117 



35 CONTINUE 

 GO TO 118 



117 WRITE(06,119) 



119 FORMATC YOU HAVE TOO MANY PROVINCES-INCREASE DIMENSION OR POEL') 

 GO TO 500 

 C NOW OUTPUT FINAL PROVINCE BOUNDARIES OF AT LEAST MIN. SIZE 



118 K3(KST)=99 

 KST=KST-1 



C SMOOTHE SMALL PROVINCE SEGMENTS 



CALL PRVFIX(K1,K2,K3,MIN,KST) 



J2=l 



DO 36 I=1,KST 



11=1+1 



IF(K3(II).EQ.K3(I)) GO TO 36 



JST=K1(J2) 



JJ=K2(I) 



IPNUM=K3(J2) 



ISDLAT=AY(JST) 



A=ISDLAT 



SELAT=ABS(AY(JST)-A)*60. 



ISDLNG=AX(JST) 



A=ISDLNG 



SELNG=ABS(AX(JST)-A)*60. 



IDLAT=AY(JJ) 



A=IDLAT 



EELAT=ABS(AY(JJ)-A)*60. 



IDLNG=AX(JJ) 



A=IDLNG 



EMLNG=ABS(AX(JJ)-A)*60. 



IDIFF=K2(I)-K1(I) 

 C COMPUTE RMS LEVEL OF HIGH PASSED DATA FOR THIS PROVINCE 

 C ALSO AVERAGE SPECTRAL SLOPE AND Y-INTERCEPT 



C3=ABS(JJ-JST+1) 



RMS=0.0 



SLOPEX=0.0 



XINTER=0.0 



DO 37 IE=JST,JJ 



SL0PEX=SL0PEX+(TZ(IE,11)/C3) 



XINTER=XINTER+(TZ(IE,12)/C3) 



156 



