日志
【转】中国地区降水资料的下载及fortran读取
热度 1 ||
1.数据集信息
数据集中文名称：中国地面降水月值0.5°×0.5°格点数据集
数据集代码：SURF_CLI_CHN_PRE_MON_GRID_0.5
数据集版本：V2.0
数据集建立时间：20120801
2.数据来源：该数据集的数据来源包括2个部分：由国家气象信息中心基础资料专项收集、整理的1961年至最新的全国国家级台站（基本、基准和一般站）的降水月值资料；由GTOPO30数据（分辨率为0.05°×0.05°）经过重采样生产的中国陆地0.5°×0.5°的数字高程模型DEM。
3.数据集实体
3.1.数据集实体内容说明
3.1.1.数据集实体文件名称：
中国地面降水月值0.5°×0.5°格点数据集文件命名由数据集代码（SURF_CLI_CHN_PRE_MON_GRID_0.5）、年份、月份标识（YYYYMM）组成。
具体形式：SURF_CLI_CHN_PRE_MON_GRID_0.5 -YYYYMM.TXT
3.1.2.数据集实体文件的内容描述：
数据集实体包括1961年1月至最新的逐月数据文件，每个文件中包括的前6行为头文件信息，其中：
第一行"ncols 128"：表示实体数据有128列；
第二行"nrows 72"表示实体数据有72行；
第三行"xllcorner 72"表示实体数据左下方对应最小经度是72°E；
第四行"yllcorner 18"表示实体数据左下方对应最小纬度是18°N；
第五行"cellsize 0.5"表示网格是0.5°×0.5°的；
第六行"NODATA_value -9999.0"表示中国区域以外的值用-9999.0表示。
从第七行开始是对应网格的月降水量，保留了1位小数。
经度单位：度，纬度单位：度，格点降水单位：mm
3.1.3.特征值说明：中国区域以外的值用-9999.0表示。
3.2.数据存储信息
3.2.1.存储格式和读取：数据集存储格式为ARCGIS标准格式，文本文件，固定长记录，按行读取。
3.2.2.数据集在介质中的放置
存储介质及数量：存储介质：光盘，1张
存储目录结构：
datasets：存放数据集实体文件。共包括1961年1月到2012年5月间的617个月值网格点降水数据文件。
metadata：元数据文档(SURF_CLI_CHN_PRE_MON_GRID_0.5_META_C.doc)。
description：说明文档(SURF_CLI_CHN_PRE_MON_GRID_0.5_DOCU_C.DOC)。
documents：存放了数据集附属文件。
数据总量：45.7MB
3.3.时间属性
时间范围：1961年1月－最新
时间分辨率：逐月
3.4.空间属性
3.4.1.地理范围
地理范围描述：中国
最西经度：72°E
最东经度：136°E
最北纬度：54°N
最南纬度：18°N
3.4.2.台站信息：无
3.4.3.空间分辨率：0.5°×0.5°网格点
3.4.4.垂直范围：无
3.4.5.投影方式：无
3.5.观测仪器：祥见《地面观测规范》的相关内容
3.6.数据处理方法：
步骤1：利用1961年1月－2012年5月降水月值数据，经格式转换，生成符合ANUSPLIN软件要求的数据；
步骤2：利用GTOPO30数据，经过分幅解压、格式转换重采样等，生成符合要求的中国范围0.5°×0.5°的数字高程模型数据dem；
步骤3：运行ANUSPLIN软件中的selnot.exe模块，按照广义交叉验证值最优原则选择初始结点；
步骤4：选择经度、纬度、海拔高度3个独立变量，运行splinb.exe模块，生成表面系数文件等；
步骤5：运行lapgrd.exe模块，生成每月的0.5°×0.5°的地面降水网格数据；
步骤6：对降水格点数据集进行交叉验证和误差分析。
3.7.数据质量状况
3.7.1质量控制方法：交叉验证质量检验法
3.7.2质量状况：经交叉验证，质量状况良好。
3.8.数据完整性：包括1961年1月－2012年5月间共617个月数据文件，数据完整。
4.引用文献：无
5.数据集制作及技术支持
5.1.数据集制作人
姓名：赵煜飞 朱江 许艳
单位：中国气象局国家气象信息中心资料服务室
电话：68407255
5.2.数据集文档编撰者
姓名：赵煜飞
单位：中国气象局国家气象信息中心资料服务室
电话：68407255
5.3.技术支持
单位：中国气象局国家气象信息中心资料服务室
电话：68407499
传真：68407499
e-mail：cdc@cma.gov.cn
邮政编码：100081
单位地址：北京海淀区中关村南大街46号
在人人网上有好友分享了一篇较为详细的对本数据的下载，这里一并贴出并表示感谢。
不会处理站点资料的童鞋请进！这里有福利！新增资料处理方法！！
众所周知，由于站点资料本身的特点，其处理难度很大。在这里我为大家提供一种替代资料，可以避开站点资料的处理工作。同时下面提供的网站里面还含有许多其他中国地区的气象资料，需要的同学请自取。
一、资料下载
1.注册：
进入“中国气象科学数据共享服务网”网址：http://cdc.nmic.cn/home.do（建议用IE浏览器打开，否则可能无法下载（[url=http://bbs.06climate.com/forum.php?mod=viewthread&tid=1105]自动回复：请不要使用迅雷等下载工具，点我查看下载帮助[/url]）资料），新用户需要注册，点击注册后进入下面的页面。
点击“一般会员注册”，填写个人信息后即可登录下载资料。
2.登录并下载
登录后进入首页
下一步点击数据集检索，进入下面的页面
左边一排是数据类型，选择你需要的数据即可进入下载。有一点需要说明，由于我们是一般会员，所以只能下载共享级别为0级和1级的资料。下面以下载中国降水格点资料为例来具体说明下载步骤。
点击最左面一栏里面的“气象服务产品”出现如下界面
这样可下载的数据就在右边显示出来了。
找到“中国地面降水月值0.5?x0.5?格点数据集（v2.0）”
点击“检索”则可以进入下载界面（这里需要注意一点，我们一次只能下载不大于50M的资料，如果你选择的时间段过长会无法下载（[url=http://bbs.06climate.com/forum.php?mod=viewthread&tid=1105]自动回复：请不要使用迅雷等下载工具，点我查看下载帮助[/url]），这时建议你分时间段多次下载）：
选择你需要的年月后点击“查询”，进入如下界面（最下面有个全选）：
选择你需要的文件就可以直接下载了。数据说明文件有关于数据的详细解释。
二、存在的问题：
资料的质量没有验证，对于有能力处理站点资料的同学建议依然使用站点资料。
三、资料读取及画图方法
1.FORTRAN读取资料的代码
program prec
parameter(xx=128,yy=72)
real prec(xx,yy)
open(10,file='SURF_CLI_CHN_PRE_MON_GRID_0.5-198301.txt',status='old',form='formatted')
do j=1,6
read(10,*)
enddo
read (10,100) prec
close(10)
100 format(128f8.1)
open(11,file='rain.grd',form='binary')
do j=72,1,-1
write(11) (prec(i,j),i=1,128)
enddo
close(11)
end
生成“rain.grd”文件之后就可以用GrADS画图了
2.ctl描述文件
dset rain.grd
undef -9999.0
title CHINA PRECIPITATION GRID DATA
xdef 128 linear 72.000 0.500
ydef 72 linear 18.000 0.500
zdef 1 levels 1000
tdef 1 linear JAN1983 1mo
vars 1
rain 0 99
endvars
