站点资料使用

1、第第8 8章章 站点资料的使用站点资料的使用站点资料站点资料是指资料在四维空间中的分布是不规则的，即离散型离散型数据资料。GrADS提供了部分支持站点资料的功能。但实际应用中，常常将站点资料插值为格点资料以后再使用绘图。8.1 使用站点资料绘图使用站点资料绘图 一般站点资料多以文本形式存放（即ASCII数据文件），需要按照GrADS的要求转换为二进制站点资料文件。台站号，纬度，经度，tim，nlev，nflag，地面物理量数值，高空物理量数值.台站号，纬度，经度，tim，0，nflag160台站第一时次结束记录行第一时次第二时次台站号，纬度，经度，tim，nlev，nflag，地面物理量数值，

2、高空物理量数值.台站号，纬度，经度，tim，0，nflag160台站第二时次结束记录行.站点资料的排放顺序 台站号必须是8个字符 tim为一个确定的时刻值，是本报时次的相对值，一般取0.0； nlev表示物理量数据的组数，包括一个地面报组加上各高度层数，取0表示一个时次记录的结束； nflag取1表示有地面观测变量，取0表示记录中没有地面观测变量。台站号、纬度、经度、tim 、nlev、nflag、物理量数值头记录数据内容每个地面变量组具有相同的数据量大小，缺测点用缺测值代替，地面变量以浮点的形式写出，变量排列顺序在各个时刻的报告中必须形同，记录顺序在描述文件中给出。每一层中的高空变量组须将所

3、有的变量写出，缺测的用缺测值代替，所以每层高空变量组在各层和整个文件的各个时刻报告中都有相同的数据量。高空变量组写出顺序如下：level 浮点数值写出地球坐标系中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）将十进制数据文件按GrADS要求转换为二进制站点数据文件 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）经Fortran程序转换后得到二进制站点资料文件st1508.dat，编写其对应的数据描述文件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 该文件是即将要创建的，与数据文件对应的站点映射文件3）编辑完数据描述文件后，利用“stnmap”命令生成映射文件“st1508.map”ga_! stnmap。出现提示信息Enter stnctl filename: 输入刚才编辑完成的数据描述文件的路径及名称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）将十进制数据文件按GrADS要求转换为二进制站点数据文件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）经Fortran程序转换后得到二进制站点资料文件r

13、ain24.grd，编写其对应的数据描述文件station.ctldset e:rain0rain.grddtype stationstnmap e:rain0rain.mapundef -999.0title autumn rain tdef 1 linear 08:00Z07JUL2007 1hrvars 1p 0 99 rainfall dataENDVARS3）编辑完数据描述文件后，利用“stnmap”命令生成映射文件“rain.map”ga_! stnmap。出现提示信息Enter stnctl filename: 输入刚才编辑完成的数据描述文件的路径及名称e:rain0statio

14、n.ctl 回车即可在e:rain0目录下生成rain.map4)准备格点数据资料及对应的数据描述文件 该格点资料以备站点资料插值使用 格点资料设置的分辨率决定了站点资料插值后的分辨率 格点资料的精度和范围应比站点资料大 格点资料的描述时间应与站点资料一致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格点数据资料“prc5.dat”的数据描述文件5）将二进制站点资料“rain24.grd”插值到格点文件“prc5.dat”的网格点上，并绘图reinitopen e:rai