1、第第8 8章章 站點資料的使用站點資料的使用站點資料站點資料是指資料在四維空間中的分布是不規(guī)則的，即離散型離散型數(shù)據(jù)資料。GrADS提供了部分支持站點資料的功能。但實際應用中，常常將站點資料插值為格點資料以后再使用繪圖。8.1 使用站點資料繪圖使用站點資料繪圖 一般站點資料多以文本形式存放（即ASCII數(shù)據(jù)文件），需要按照GrADS的要求轉換為二進制站點資料文件。臺站號，緯度，經(jīng)度，tim，nlev，nflag，地面物理量數(shù)值，高空物理量數(shù)值.臺站號，緯度，經(jīng)度，tim，0，nflag160臺站第一時次結束記錄行第一時次第二時次臺站號，緯度，經(jīng)度，tim，nlev，nflag，地面物理量數(shù)值，

2、高空物理量數(shù)值.臺站號，緯度，經(jīng)度，tim，0，nflag160臺站第二時次結束記錄行.站點資料的排放順序 臺站號必須是8個字符 tim為一個確定的時刻值，是本報時次的相對值，一般取0.0； nlev表示物理量數(shù)據(jù)的組數(shù)，包括一個地面報組加上各高度層數(shù)，取0表示一個時次記錄的結束； nflag取1表示有地面觀測變量，取0表示記錄中沒有地面觀測變量。臺站號、緯度、經(jīng)度、tim 、nlev、nflag、物理量數(shù)值頭記錄數(shù)據(jù)內(nèi)容每個地面變量組具有相同的數(shù)據(jù)量大小，缺測點用缺測值代替，地面變量以浮點的形式寫出，變量排列順序在各個時刻的報告中必須形同，記錄順序在描述文件中給出。每一層中的高空變量組須將所

3、有的變量寫出，缺測的用缺測值代替，所以每層高空變量組在各層和整個文件的各個時刻報告中都有相同的數(shù)據(jù)量。高空變量組寫出順序如下：level 浮點數(shù)值寫出地球坐標系中z維方向的值variables 本層中左右變量報告例子：站點溫度 54398 116.617 40.133 54399 116.283 39.983 54412 116.633 40.733 54419 116.633 40.300 54421 117.117 40.650 54431 116.633 39.917 54499 116.217 40.220 54501 115.683 39.967 54505 116.120 39.9

4、20 54513 116.200 39.950 54514 116.250 39.870 54596 116.000 39.700 54597 115.733 39.733 54410 116.133 40.600 54406 115.967 40.450 54424 117.100 40.150 54433 116.483 39.950 54594 116.330 39.750 54416 116.867 40.383 54511 116.467 39.800 st.txt文件 27.90 28.80 24.50 28.20 26.30 27.90 27.60 26.90 28.90 29.

5、30 29.60 27.10 25.10 21.90 24.90 28.10 29.10 28.90 27.80 29.00t1508.txt文件1）將十進制數(shù)據(jù)文件按GrADS要求轉換為二進制站點數(shù)據(jù)文件 program station parameter (mt=20) ! mt means nums of station character stid(mt)*8 integer nflag,nlev real lat(mt),lon(mt),tim,t2(mt) open(30,file=e:beijingobservest.txt) open(50,file=e:beijingobse

6、rvet1508.txt)open(40,file=e:beijingobservest1508.dat,form=binary) 編寫Fortran程序： do i=1,mt read(30,*)stid(i),lon(i),lat(i) write(*,*)stid(i),lat(i),lon(i) end do close(30) read(50,*) (t2(i),i=1,mt) close(50) tim=0.0 nflag=1 nlev=1 do i=1,mt write(40)stid(i),lat(i),lon(i),tim,nlev,nflag,t2(i) end do nl

7、ev=0 write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,nflag close(40) end 2）經(jīng)Fortran程序轉換后得到二進制站點資料文件st1508.dat，編寫其對應的數(shù)據(jù)描述文件st1508.ctldset e:beijingobservest1508.dat dtype station stnmap e:beijingobservest1508.mapundef 999.9 title temp duration tdef 1 linear 14Jul2002 1hr vars 1 t2 0 99 temp duration

8、 endvars 該文件是即將要創(chuàng)建的，與數(shù)據(jù)文件對應的站點映射文件3）編輯完數(shù)據(jù)描述文件后，利用“stnmap”命令生成映射文件“st1508.map”ga_! stnmap。出現(xiàn)提示信息Enter stnctl filename: 輸入剛才編輯完成的數(shù)據(jù)描述文件的路徑及名稱e:beijingobservest1508.ctl 回車即可在e:beijingobserve目錄下生成st1508.map4）編輯gs文件繪圖reinitopen e:beijingobservest1508.ctlset lon 115.4 117.4set lat 39.4 41.0set grads offse

9、t grid offset mpdset beij*set gxout stnmark（）set rbcols 3 3 3set cmark 3set digsiz 0.1set stid on （站點標號資料值一起顯示）d t2;采用set gxout stnmark8.2 站點資料處理為格點資料繪圖站點資料處理為格點資料繪圖 以降水資料為例，利用站點資料繪制降水量分布場（等值線）資料說明：多臺站記錄的2007年7月8日24h累積降水量資料，保存為文本資料。54452 119.81 40.80 367 5.754454 120.35 40.35 16 13.354455 120.70 40.

10、58 10 24.354237 121.65 42.03 145 .156018 95.30 32.90 4068 6.750950 125.25 45.70 150 .0150963 128.73 45.97 110 3.554096 131.15 44.38 498 15.856645 99.41 26.41 2345 23.456664 101.26 26.63 1242 .756684 103.28 26.42 2110 1.156739 98.48 25.12 1649 3.056768 101.53 25.02 1773 .0156786 103.83 25.58 1900 .15

11、6886 103.77 24.53 1708 1.056969 101.58 21.50 633 .156973 102.83 23.61 1310 .150953 126.77 45.75 143 .0150867 127.35 46.08 9999 .1。rain24.txt1）將十進制數(shù)據(jù)文件按GrADS要求轉換為二進制站點數(shù)據(jù)文件real vec(1092)real lat(1092),lon(1092)integer h(1092)character*5 stid(1092)open(1,file=e:rain0rain24.txt,status=old)do 20 k=1,1092

12、20 read(1,*) stid(k),lat(k),lon(k),h(k),vec(k) close(1) OPEN (3,FILE=e:rain0rain24.grd,FORM=BINARY) TIM=0.0 NLEV=1 NFLAG=1 DO 40 I=1,1092 WRITE(3) STID(I),LAT(I),LON(I) # ,TIM,NLEV,NFLAG,vec(i) 40 continueNLEV = 0WRITE(3) STID(I-1),LAT(I-1),LON(I-1),TIM,NLEV,NFLAG close(3)2）經(jīng)Fortran程序轉換后得到二進制站點資料文件r

13、ain24.grd，編寫其對應的數(shù)據(jù)描述文件station.ctldset e:rain0rain.grddtype stationstnmap e:rain0rain.mapundef -999.0title autumn rain tdef 1 linear 08:00Z07JUL2007 1hrvars 1p 0 99 rainfall dataENDVARS3）編輯完數(shù)據(jù)描述文件后，利用“stnmap”命令生成映射文件“rain.map”ga_! stnmap。出現(xiàn)提示信息Enter stnctl filename: 輸入剛才編輯完成的數(shù)據(jù)描述文件的路徑及名稱e:rain0statio

14、n.ctl 回車即可在e:rain0目錄下生成rain.map4)準備格點數(shù)據(jù)資料及對應的數(shù)據(jù)描述文件 該格點資料以備站點資料插值使用 格點資料設置的分辨率決定了站點資料插值后的分辨率 格點資料的精度和范圍應比站點資料大 格點資料的描述時間應與站點資料一致DSET E:rain0PRC5.DATTITLE PRCpUNDEF -9.99E33XDEF 141 LINEAR 70 0.5YDEF 81 LINEAR 15 0.5ZDEF 1 LINEAR 1 1tdef 1 linear 08:00Z07JUL2007 1hrVARS 1g 0 99 afdENDVARS格點數(shù)據(jù)資料“prc5.dat”的數(shù)據(jù)描述文件5）將二進制站點資料“rain24.grd”插值到格點文件“prc5.dat”的網(wǎng)格點上，并繪圖reinitopen e:rai