c	revised by tang for SiB2_DHM
c	tangqh@iis.u-tokyo.ac.jp
c	2005/9/18
        program sod

        integer i, j, m, q, ct, nv, stnmsc 
        integer stnch(1), stn(1), hdrct, obsct, blkch(100)
        integer days, dayf, dayrun, ranglo, ranghi
        integer stnf, yrf, mof, daf
        integer pckyrs, pckmos, pckdas, pckyrf, pckmof, pckdaf, done

        character*1 ch, invch, morech, limch
        character*2 cntry, ctrych(50), statch(50)
        character*3 moabb(12), mopick
        character*4 call, yrabb(15), yrpick
        character*5 lat, elev 
        character*6 lon
        character*19 city
        character*80 filler, invfile, datfile
        character*116 dataf
        character*132 header
		
	integer latmin,latmax,lonmin,lonmax
	integer lati,loni,elevi

	latmin=-18000 
	latmax=18000 
	lonmin=-18000
	lonmax=18000

	nstn = 0
	        
        open (unit=10, file='stnlist.txt', form='formatted',
     +   access='sequential', status='old')
        open (unit=100, file='stn_arc1.txt', form='formatted',
     +   access='sequential', status='unknown')

        do 15 j=1,14
           read (unit=10, fmt=5010) filler
 5010   format(A80)
 15     continue

 110   read (unit=10, fmt=5041, end=999) stnmsc, call, city,   
     +   cntry, lati, loni, elevi   
     
 5041   format(I6.6,1X,A4,3X,A19,1X,A2,1X,I5,1X,I6,1X,I5)
 5042   format(I,3f15.5)

		if (city(1:5).ne.'BOGUS') then

		if (lati.le.latmax.and.lati.ge.latmin) then
		if (loni.le.lonmax.and.loni.ge.lonmin) then

			if (lati.ne.0.and.loni.ne.0) then
				nstn = nstn +1
				write(100,*) nstn, loni/100., lati/100.
			endif

		endif !(search box)			   
		endif !(search box)

		endif !(city.ne.'BOGUS CHINESE      ')	 

 25     continue        
        goto 110
 999    continue              

        close (unit=10)
	close (unit=100)
	
	end