grads處理臺站數(shù)據(jù)精講

1、Nanjinguniversityofinfmationscience&technology阿木阿木1廢話不說，開始講解廢話不說，開始講解第一步：處理站臺第一步：處理站臺id文件：文件：文件處理成三列數(shù)據(jù)，第一列為站臺號，第二列為緯度，第三列為經(jīng)度注意，如：站號緯度經(jīng)度5013652.97122.525024652.35124.725034951.70124.335035351.72126.65上述數(shù)據(jù)中50136連同后面的空格要保證一

2、共有8位（這點不容質(zhì)疑）第二步：處理數(shù)據(jù)第二步：處理數(shù)據(jù)一般要將所要的要素，如日照時間，單個要素放到一個txt文件中，數(shù)據(jù)的順序是：同一個時間點371個站數(shù)據(jù)，然后放下一個時間點的371個要素值，如此直到最后時刻。第三步：處理成對應(yīng)格式數(shù)據(jù)第三步：處理成對應(yīng)格式數(shù)據(jù)grd文件文件案例程序（371個站18263天的數(shù)據(jù)）：說明：這里的rct_50years_371_5vars50year_371id.txt為第一步生成的臺站id文件rct

3、_50years_371_5vars371_sd_1956_2005.txt為第二步生成的數(shù)據(jù)文件注意：1、數(shù)據(jù)必須設(shè)成實型的，站號必須為8位字符型（與第一步照應(yīng)）；2、數(shù)據(jù)存放依然是先站，后時次，即最外層為時間循環(huán)3、每一行數(shù)據(jù)：站號—緯度—經(jīng)度—相對時次—記錄標志—nflag—要素數(shù)據(jù)4、程序中tim為相對時次，指的是相對于一天的世界時，如果數(shù)據(jù)本身是日數(shù)據(jù)，那么time就取0.0，一般都是這樣，除非是有08時、20時等多個時刻，才

4、考慮time不為0.0，此處設(shè)為0.0恒定，元數(shù)據(jù)為日數(shù)據(jù)5、程序中nlev標志一個時刻記錄開始，結(jié)束時賦值為0，對應(yīng)程序；6、nflag為有沒有地面報，有為1，否則為0，此處為地面觀測相對濕度，為17、循環(huán)中每個時次的371個站數(shù)據(jù)輸完后，要加一行標志記錄（紅色）8、生成的文件為二進制grd文件，所以為write（30），沒有格式說明；直接加文件號9、第一二三步數(shù)據(jù)的順序一定要與站號順序?qū)?yīng)。programex00implicitno

5、neintegernlevnflagijreallat(371)lon(371)timsd(371)acter8stid(371)tim=0.0nflag=1open(10file=D:ProgramFilesFTRAN6.0txt_gridrct_50years_371_5vars50year_371id.txt)open(20file=D:ProgramFilesFTRAN6.0txt_gridrct_50years_371_5va

6、rs371_sd_1956_2005.txt)Nanjinguniversityofinfmationscience&technology阿木阿木3enddodoi=1nx1lon(i1)=lon(i)1.0enddodok=118263doi=1nxdoj=1nys(ij)=1enddoenddoenddowrite(1)send第六步：寫出對應(yīng)插值第六步：寫出對應(yīng)插值grd的ctldsetc:371yue371_sd_1956_20

7、05.grddtypestationstnmapc:371yuesd50.mapundef32766titlesunhourtdef18263linear01jan19561dyvars1sd099sunhour(0.1h)endvars第七步：生成數(shù)據(jù)地圖文件第七步：生成數(shù)據(jù)地圖文件臺站的數(shù)據(jù)文件，需要地圖文件才能讀取方法：在grads啟動顯示圖形窗口后，輸入!stnmap回車輸入：c:371yue371_sd_1956_2005.c