clear #delimit; set mem 150m; global datapath ~/data/china/china1990/gis1990/historical/; shp2dta using $datapath/CH_CENSUS1990_PARTA.shp, database($datapath/CH_CENSUS1990_PARTA.dta) coordinates($datapath/coord1.dta) genid(id) replace;