1,S0ATA(4500),BDATA(2050),CDATA(2050), PHASE (2050), SMPH( 2050) 



COMMON/BOUND/I, IPNUM.ISDLAT,SELAT,ISDLNG,SELN6,IDLAT,EELAT, 

 IIDLNG.EMLNG 



CHARACTERS ITH(8) ,ITG(8) 



CHARACTER*! RESPON 

 C READ(5,*)NSETS,JC0DE,IUNIT,JUNIT,IFE0F,GRIDNM 



NSETS=1000 



JC0DE=1 



IUNIT=5 



JUNIT=2 



IFE0F=1 

 C GRID SPACING IS READ FROM DATA FILE 



IF(JCODE.6E.O) CALL PLOTS (STO. 2500, 2) 



CALL FACTOR (.65) 



JSET=1 



IF(JCODE.GE.O) CALL PL0T(1.0,0.0,-3) 

 C REA0(IUNIT,*)NP,KX0RG,KY0RG,AN0RM1,(CUT(I),H(I),N(I),I=1,2),IMEAN, 



C 1 KPHA,IPHA,IORGN,ILIST,IADJ 

 C NUMBER OF POINTS IS READ FROM DATA FILE, MAXIMUM IS HELD TO 2000 



NPTS=2000 



KX0RG=-2 



KY0RG=-7 



AN0RM1=2.0 



CUT(1)=-1.0 



H(l)=.2 



N(l)=3 



CUT(2)=.08 



H(2)=.2 



N(2)=3 



IMEAN=1 



KPHA=0 



IPHA=0 



I0RGN=1 



ILIST=0 



IA0J=1 

 245 CONTINUE 



APHA=IPHA 



X0RG=10.0**KX0RG 



Y0RG=10.0**KY0RG 



I0UT=7 



IF(ABS(AN0RM1).LT.1.0E-20)AN0RM=1.0 



CALL PRVOUT(NP,DATA,ITG,XINTER,SLOPEX,GRIDNM) 



GRIDKM=GRIDNM*1.852 

 C XINTER=10**(AL0610(XINTER)-(ALOG10(l./GRIDKM)+.5)) 



2451 CT0FF2=.5*(1./GRIDKM) 

 CT0FF1=.01*(1./GRIDKM) 

 CT0FF1=AL0G10(CT0FF1) 

 CT0FF2=AL0G10(CT0FF2) 

 CT1SV=CT0FF1 

 CT2SV=CT0FF2 



C SAVE DATA ARRAY IN SDATA ARRAY 



DO 2452 1=1, NP 



2452 SDATA(I)=DATA(I) 



186 



