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

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

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

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

4、考慮time不為0.0，此處設(shè)為0.0恒定，元數(shù)據(jù)為日數(shù)據(jù)5、程序中nlev標(biāo)志一個(gè)時(shí)刻記錄開始，結(jié)束時(shí)賦值為0，對(duì)應(yīng)程序；6、nflag為有沒有地面報(bào)，有為1，否則為0，此處為地面觀測(cè)相對(duì)濕度，為17、循環(huán)中每個(gè)時(shí)次的371個(gè)站數(shù)據(jù)輸完后，要加一行標(biāo)志記錄（紅色）8、生成的文件為二進(jìn)制grd文件，所以為write（30），沒有格式說明；直接加文件號(hào)9、第一二三步數(shù)據(jù)的順序一定要與站號(hào)順序?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第六步：寫出對(duì)應(yīng)插值第六步：寫出對(duì)應(yīng)插值grd的ctldsetc:371yue371_sd_1956_20

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