大气再分析资料的预处理1 JRA资料的下载与预处理
第一章 再分析资料的预处理
1 JRA资料的预处理
1.1 JRA资料简介
The JapanMeteorological Agency (JMA) conducted JRA-55, the second Japanese globalatmospheric reanalysis project. It covers 55 years, extending back to 1958,coinciding with the establishment of the global radiosonde observing system.Compared to its predecessor, JRA-25, JRA-55 is based on a new data assimilationand prediction system (DA) that improves many deficiencies found in the firstJapanese reanalysis. These improvements have come about by implementing higherspatial resolution (TL319L60), a new radiation scheme, four-dimensionalvariational data assimilation (4D-Var) with Variational Bias Correction (VarBC)for satellite radiances, and introduction of greenhouse gases with time varyingconcentrations. The entire JRA-55 production was completed in 2013, andthereafter will be continued on a real time basis.
Specificearly results of quality assessment of JRA-55 indicate that a large temperaturebias in the lower stratosphere has been significantly reduced compared toJRA-25 through a combination of the new radiation scheme and application ofVarBC (which also reduces unrealistic temperature variations). In addition, adry land surface anomaly in the Amazon basin has been mitigated, and overallforecast scores are much improved over JRA-25.
Most ofthe observational data employed in JRA-55 are those used in JRA-25.Additionally, newly reprocessed METEOSAT and GMS data were supplied by EUMETSATand MSC/JMA respectively. Snow depth data over the United States, Russia andMongolia were supplied by UCAR, RIHMI and IMH respectively.
The DataSupport Section (DSS) at NCAR has processed the 1.25 degree version of JRA-55with the RDA (Research Data Archive) archiving and metadata system. The modelresolution data has also been acquired, archived and processed as well,including transformation of the TL319L60 grid to a regular latitude-longitudeGaussian grid (320 latitudes by 640 longitudes, nominally 0.5625 degree). AllRDA JRA-55 data is available for internet download, including completesubsetting and data format conversion services.
Both daily and monthly data areencoded in Gridded binary (GRIB) Edition 1 format (WMO 2011). (Daily datainclude 3-, 6- and 24-hourly data, which vary by category).
1.2 JRA55资料下载
数据下载说明页面:http://jra.kishou.go.jp/JRA-55/index_en.html
JRA55的数据列表参考文件JRA-55_handbook_LL125_en.pdf
手册下载地址为:http://jra.kishou.go.jp/JRA-55/document/JRA-55_handbook_LL125_en.pdf
方法一:手动下载,适用于数据不多的情况
在浏览器中输入ftp://ds.data.jma.go.jp,输入用户名和密码即可登录下载所需数据
方法二:在LINUX下或者CMD命令行中下载,适用于数据稍多的情况
$ ftp ds.data.jma.go.jp # log-on to the JMA DataDissemination System (JDDS)
???????? # user ID
******** # password
binary # changes transfer mode to binary
cd JRA-55/Hist/Monthly/anl_surf125
get anl_surf125.195801
get anl_surf125.195802
...
bye
方法三:在FTP专业下载软件FileZilla下,适用于大量数据的下载
连接主机ftp://ds.data.jma.go.jp
输入账户密码,选择所需要的数据和本地下载路径,即可批量快速下载
注意1: ”JRA55是GRIB1类型的资料,下载数据完成后,将批量描述的CTL和索引文件.idx也一起下载,否则自己按照“E8料认识与应用(GrADS)-_兰溪整合版(grib1_grib2).pdf”中的方法生成CTL文件和对应的索引文件。
注意2:下载资料前需要自己申请JRA官网下载数据的账号和密码,账号申请的地址为:http://jra.kishou.go.jp/cgi-bin/secaccess.cgi?pr=2 ,注意申请时的说明,一般第二天就会邮件回复,给一个账户和登录密码。
1.3 JRA55资料预处理
第一步:确认数据变量的层次数目(可做可不做):在CMD或者LINUX中命令中输入wgrib{文件名},查看数据的层次等信息,以确认层次的数目(注意wgrib.exe需要事先下载好并将其路径添加到环境变量中去)
第二步:修改.ctl文件中的路径,保存
第三步:编写合并多个GRIB文件的.gs文件,并用GRADS运行,生成合并后的.grd文件
注意易错:由于网络不稳定等原因,所有文件下载完成后,首先需要检查文件总个数是否正确(一个时次一个文件),其次检查每个文件大小是否相同,即是否每个文件都完整下载了,以防后续处理出错。
实例:在FTP上下载了JRA55总共696个时次的月平均位势高度场的资料到本地文件夹,
同时还将anl_p125_hgt.monthly.ctl和anl_p125_hgt.monthly.idx文件一同下载到本地。(JRA55资料的处理方法类似于FNL资料的处理方法)
STEP1:根据自己的实际修改anl_p125_hgt.monthly.ctl文件中前两行的默认路径,以及需要的时次,最终得到的CTL文件的内容为:
dset ~\anl_p125_hgt.%y4%m2
index ~\anl_p125_hgt.monthly.idx
undef9.999E+20
titleanl_p125_hgt
* produced by grib2ctl v0.9.12.5p41
dtypegrib 255
options yrev template
ydef145 linear -90.000000 1.25
xdef288 linear 0.000000 1.250000
tdef 696 linear 00Z01Jan1958 1mo
*根据下载的文件数进行调整
* z has 37 levels, for prs
zdef37 levels
1000975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
vars1
HGTprs37 7,100,0 ** (profile) Geopotential height [gpm]
ENDVARS
重要说明:因为官网上下载的索引文件(包括CTL)是针对所有文件时次的,即索引文件开始的时次已经确定,所以下载的时候,最好全部下载!否则自己修改CTL并重新生成索引文件!或者,不必下载到最后一个时次,但必须从最开始时次(1958年1月或者1958年1月1日….)的文件下载,此时对CTL文件时间维的维数稍加修改即可,但是起始时间维不变!此为本人理解,若有错误望大家指正!!!
STEP2:自行编写.gs脚本合并所有的GRIB文件,这里合并所有层次、时次和二维空间场的数据:
'reinit'
'open ~/anl_p125_hgt.monthly.ctl'
*自行修改路径
'setgxout fwrite'
'setfwrite ~/JRA55.hgt.mon.order.grd'
it=253
*这里是从1979年1月开始到2015年12月
while(it<=696)
'setlat -90 90'
'set lon 0 358.75'
*这里易将'set lon 0 358.75'错写为'set lon 0 360'从而导致数据排列发生混乱
*'setx 1 288'
*'sety 1 145'
'sett 'it''
iz=1
while(iz<=37)
'setz 'iz''
'dHGTprs'
iz=iz+1
endwhile
it=it+1
endwhile
'disablefwrite'
;
数据量大的话慢慢等,去喝杯咖啡,但是你不会就此变成洛伦兹。。。。。。
结果:2.56G! 大功告成!为这个巨大的文件配个CTL就可以做后续工作了!
最后别急!一定要对合并的GRD文件做一个正确性验证!指定一个时次,用原始单个文件的数据和合并后的数据分别画图进行验证!!
大气再分析资料的预处理1 JRA资料的下载与预处理相关推荐
- 美国NCEP/NCAR再分析资料简介
NCEP/NCAR再分析数据集是由美国气象环境预报中心(NCEP)和美国国家大气研究中心(NCAR)联合制作的,他们采用了当今最先进的全球资料同化系统和完善的数据库,对各种来源(地面.船舶.无线电探空 ...
- 关于NCEP/NCAR再分析资料学习
在NECP/NCAR官网上下载再分析资料https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html 这里 ...
- Sentinel-2数据下载及预处理
Sentinel-2数据下载及预处理 Sentinel-2简介: 欧空局(ESA)仅发布了哨兵2号(S2)的L1C级多光谱数据(MSI),L1C级数据是经过几何精校正的正射影像,并没有进行辐射定标和大 ...
- c语言预处理命令12个,C语言编译预处理和预处理命令
1.预处理概念:编译程序时,编译器将对程序文件作至少两个阶段的编译预处理,预处理程序预处理阶段和编译程序预处理阶段. 预处理程序预处理阶段:C预处理程序先于C编译程序运行.预处理程序从前向后根据预处理 ...
- c语言提供了三种预处理命令,9、C语言之预处理命令
预处理命令 基本概念:ANSI C标准规定可以在C源程序中加入一些"预处理命令",以改进程序设计环境,提高编程效率. 这些预处理命令是由ANSI C同一规定的,但是它们不是C语言本 ...
- c语言有效的预处理命令,9、C语言之预处理命令(示例代码)
预处理命令 基本概念:ANSI C标准规定可以在C源程序中加入一些"预处理命令",以改进程序设计环境,提高编程效率. 这些预处理命令是由ANSI C同一规定的,但是它们不是C语言本 ...
- c语言输入输出预处理命令,C语言第十一讲,预处理命令.
C语言第十一讲,预处理命令. 一丶预处理简介 什么是预处理,预处理就是预先处理. 我们知道,程序会经过编译,连接形成可执行文件 这些在编译之前对源文件进行简单加工的过程,就称为预处理(即预先处理.提前 ...
- 遥感多光谱数据下载与预处理(一、数据选择 下载)
首先说明本人并非专业大牛,不是教程贴只是记录一下学习过程和大家交流,过程有不严谨不合规范不对的地方欢迎各位大神指正. 本人目前做过接触过最多的是多光谱遥感数据,也是与无人机.雷达.高光 ...
- 珞珈一号数据下载和预处理
文章目录 数据网址 数据下载 数据预处理 数据网址 直接搜索http://59.175.109.173:8888/app/login.html即可进入数据注册页面,只需简单注册登录就可以进入下载页面了 ...
最新文章
- 第一章课后习题(Java)
- android ScrollView 控制行数
- 网络通信之 字节序转换原理与网络字节序、大端和小端模式
- 水星mw310r虚拟服务器,水星MW310R静态网络(IP)设置教程
- You have mail in /var/spool/mail/root
- MySQL数据库 基本操作语句
- SDK 和 API 有什么区别
- OpenID Connect:OAuth 2.0协议之上的简单身份层
- session和cache的区别
- XLNet 和BERT的区别是什么?
- 基于STC89C52单片机的智能灯光毕业设计论文
- 火灾自动报警系统 300W广播功率放大器
- VUE中toast的使用与开发
- IIP3/IMD/IM3/P1dB
- 2022年低压电工复训题库及在线模拟考试
- 关于Java你不知道的那些事之等等与equals的区别
- 使用公网NAT网关SNAT功能访问互联网
- MathJax使用LaTeX语法编写数学公式
- 鉴于现有的形式化方法(简称“桥”系列)有其优点和缺点,特发明新的“秤”系列
- 物联网大数据平台的主要功能和特点