1. Orinal Map
	Get FAO Irrigation area map
	Copy yr_dem10 as the basemap

2. Get FAO irrigation zone
	Use ArcMap
	Add yr_dem10
	Add gmia
	Set Data Frame: (Up: -240500 , Down: -1500500, Left: -457500, Right: 1752500)
	Export gmia, Extend: Data Frame, Spatial reference: Data Frame
	Cell size: 200
	Get gmia1, Add gmia1
	Raster Calculator: yrgmia = con ( yr_dem10 , gmia1 )
	Export (Save ) yrgmia to yrgmia1
	irriratio.asc = gridascii ( yrgmia1 )
	
3. Get SiB2 Land use
	Use ArcMap
	Add yr_dem10
	Add sib22_01a
	Raster Calculator: sib2yr = con ( yr_dem10 , sib22_01a )
	sib2y.asc = gridascii ( sib2yr )
	
4. Get FAO Soil map
	Put data at \FAOsoil\
	import auto fescntll sane
	Use ArcCatalog 	define the projection Geographic GCS_Sphere_ARC_INFO
	Use ArcMap add it and project it to LAMBERT_AZIMUTHAL
	Export it and use the same coordinate system as the data frame, get sanel.shp
	Use ArcMAP, covert covergae sane (or sanel.shp) to raster
	Input features: sane, Field: SUNM, Output cell size 0.01, Output raster: sane01
	Use raster calculator: 	soilyr1 = con ( [yr_dem10] , [sane01 - sane01] )
	soily.asc = gridascii ( soilyr1 )
	
	At \FAOsoil\SOILX\
	Compile Soil_texture.f
	f77 Soil_texture.f
	Soil_texture.exe
	Get Texture.DAT