Use data from 'GLOBAL SURFACE SUMMARY OF DAY DATA '
ftp://ftp.ncdc.noaa.gov/pub/data/globalsod/

1. Get data
	* the path should be dos path to use uncomp.exe
	copy all-1994.tar xxx.tar and all-2005.tar
	Unzip them
	
	use Generate_name.f get run1.bat
	run2.bat
	run run1
	run run2.bat
	
	get ascii format data
	jan1994.txt
	feb1994.txt
	....
	nov2000.txt 
	
	copy stnlist-sorted.txt
	rename stnlist-sorted.txt stnlist.txt
	
	use data-selection.for abstract data
	get output station_wmo.txt
	
2. Covert to metric (00: optional)
	Use conver-data.f to covert data to metric
	output data is put in to ./metric/
	
3. Select for SiB-DHM use (01: selected station)
Background: 
	There is data source from Chinese Meteorological Bureau (CMB).
	There is data source from WMO.
	The stations from CMB is save in station_info.
	The stations from wmo is save in station_wmo.txt
	Now, select all the stations which are NOT included in station_info but included in station_wmo.txt

How to do:
	get station_info file
	Use select-stn.for
	output data is put into ./selmet/
	xxxxxx.txt data file in metric (selected stations)
	selstn.txt selected station information
	allstn.txt all stations information
	
Check results:
	Arc: generate selstn		
	Generate: point
	ID,X,Y,{Angle},{Scale}: &r selstn.txt
	ID,X,Y,{Angle},{Scale}: END	
	
	Arc: generate allstn		
	Generate: point
	ID,X,Y,{Angle},{Scale}: &r allstn.txt
	ID,X,Y,{Angle},{Scale}: END	
	
	Use ArcCatalog to set projection Geographic, Sphere of Radius 6370997 meters
	
4. Put the selected stations into SIB_DHM 
	Copy 	./selmet/xxxxxx.txt to SIB_DHM (e.g. \ATM_INPUT_DATA\gauge_15item\)
	Copy allstn.txt to SIB_DHM (e.g. \GEO_INPUT_DATA\)
	allocate allstn.txt as station information file
	allocate interpolate directory for new stations
	e.g. set stn_info_file='../GEO_INPUT_DATA/allstn.txt'
	set ETPATH ='../ATM_INPUT_DATA/ETPATH2/'
	set ATM_ITP='../ATM_INPUT_DATA/Interpolate2/'
	
	f77 interpolation.f -assume byrecl -o interpolation.exe
	run interpolation.exe to generate interpolation results