NSEC=SEC/GRIDKM 

 IF(NSEC.GT,NP) NSEC=NP 

 C RECREATE DATA ARRAY 



DO 3645 ICHG=1,NP 

 LOCCHG=ICHG+NFIR 

 DATA(ICH6)=SDATA(L0CCHG) 

 IF(LOCCHG.EQ.NSEC) GO TO 3646 



3645 CONTINUE 



C INTERPOLATE NEW LAT & LON 



3646 SDIFF=FLOAT(NSEC)/FLOAT(NP) 

 FDIFF=FLOAT(NFIR)/FLOAT(NP) 

 FLAT=DECDEG(ISDLAT,SELAT) 

 FLON=DECDEG(ISDLNG,SELNG) 

 SLAT=DECDEG(IDLAT,EELAT) 

 SLON=DECDEG( IDLNG.EMLNG) 

 FNLAT=( (SLAT -FLAT )*FDIFF )+FLAT 

 SNLAT=( (SLAT -FLAT )*SDIFF)+FLAT 

 FNLON=( (SLON-FLON )*FDIFF )+FLON 

 SNLON=( (SLON-FLON )*SDIFF )+FLON 



C RECONSTRUCT INTO DEGREES AND MINUTES 



CALL DEGMIN(FNLAT,ISDLAT,SELAT) 



CALL DEGMIN(FNL0N,ISDLN6,SELNG) 



CALL DEGMIN(SNLAT,IDLAT,EELAT) 



CALL OEGMIN(SNLON,IDLNG,EMLNG) 



NP=NSEC-NFIR 



CALL ERASE 



GO TO 2451 

 C 



C INTERACTIVE ROUTINE Fu MODIFY REGRESSION FIT TO SPECTRUM 



C 



3650 WRITE (6, 3651) 



3651 F0RMAT(T70,' ARE YOU SATISFIED WITH THE REGRESSION FIT?') 

 READ(5,3652)RESPON 



3652 FORMAT (Al) 

 IF(RESPON.EQ.'Y') GO TO 3656 

 WRITE (6, 3653) 



3653 F0RMAT(T70,' ENTER LOG(LOWER FREQUENCY ),LOG(UPP£R FREQUENCY)') 

 READ(5,*) CT0FF1,CT0FF2 



CALL ERASE 

 LPASS=1 

 GO TO 4911 

 3656 CT0FF1=CT1SV 

 CTOFF2=CT2SV 

 CALL ERASE 

 LPASS=0 



WRITE (14, 399 )NP,IPNUM,ISDLAT,SELAT,ISDLNG,SELN6,IDLAT,EELAT, 

 1IDLNG,EMLNG,B1SAVE,B0SAVE 

 399 F0RMAT(I5,I8,4(1X,I4,F6.2),F6.3,E9.3) 

 C PLOT PHASE SPECTRUM 



3658 IF(KPHA.NE.l) GO TO 581 

 CALL PL0T(15.0,0.0,-3) 

 X=ALOG10(D£LK/XORG)*3.0125 

 Y=ABS(PHASE(2)/30.0) 



194 



