爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: MeteoInfo

MeteoInfoLab脚本汇总贴

  [复制链接]

新浪微博达人勋

发表于 2016-12-1 10:05:08 | 显示全部楼层
这个是我WRFOUT的ctl
dset F:\xuexi\WRFOUT\monthlydata\2009010131.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.7.1 MODEL
pdef   58  43 lcc  42.500  120.500   29.500   22.000  60.00000  30.00000  120.50000  30000.000  30000.000
xdef  195 linear  107.26420   0.13513513
ydef  110 linear   34.83622   0.13513513
zdef   19 levels  
tdef  241 linear 00Z01JAN2009     180MN      
VARS  100
变量太多,省略了,而后土层也省略了!
这个是我提取数据后的ctl
dset F:\xuexi\WRFOUT\comparativeanalysis\monthly\monthlyvary.grd
undef -9.99e+08
title 1L
xdef  281 linear  114.15895   0.04504504
ydef  150 linear   39.03433   0.04504504
zdef   1 levels 1
tdef  1 linear 00Z01JAN2009     1mon      
VARS  7
LH             0  99  LATENT HEAT FLUX AT THE SURFACE (W m-2)
T2             1  0  TEMP at 2 M (K)
VEGFRA         1  0  VEGETATION FRACTION (-)
GRDFLX         1  0  GROUND HEAT FLUX (W m-2)
LAI            1  0  LEAF AREA INDEX (m-2/m-2)
ALBEDO         1  0  ALBEDO (-)
HFX            1  0  UPWARD HEAT FLUX AT THE SURFACE (W m-2)
ENDVARS
然后两个ctl我给上传一下 2009010131.ctl (11.48 KB, 下载次数: 0)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-1 10:06:10 | 显示全部楼层
王老师这个提取的数据直接用meteoinfolab可以打开,数据我验证过提取的应该没问题,但是这个坐标轴的问题我没有解决,麻烦王老师了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-12-1 10:26:26 | 显示全部楼层
ruanjinshuai 发表于 2016-12-1 10:06
王老师这个提取的数据直接用meteoinfolab可以打开,数据我验证过提取的应该没问题,但是这个坐标轴的问题我 ...

注意到2009010131.ctl中有pdef行(也就是投影设置),说明此时的数据还是Lambert投影的。monthvary.ctl中已经没有pdef行了,数据已经投影为等经纬度,其实数据还是矩形,只是模式区域范围外的格点都是缺测值,在绘制等值线的时候会空缺。

在绘制monthvary.ctl文件时可以尝试将 axesm() 中投影的参数设置得和原始Lambert投影一样,看看效果。

或者用MeteoInfoLab对原始数据进行处理得到月变化之类的结果数据(不会改变数据的投影)。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-1 19:59:26 | 显示全部楼层
好的,谢谢王老师
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-5 15:29:57 | 显示全部楼层
必须收藏学习啊!赞!赞!赞!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-6 22:45:24 | 显示全部楼层
王老师,您好,请问meteoinfo能批处理像FNL、GLDAS这样的grib数据吗?我这是一年的数据,但是分了2920个小grb文档,如果可以的话,麻烦王老师给我发一个帖子就好,我没找到,麻烦王老师了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-6 23:05:30 | 显示全部楼层
王老师,就是刚才那个问题。我找到你之前回答这个问题http://bbs.06climate.com/forum.php?mod=viewthread&tid=17632,问题是那个微博已经删除了,没有了!网上有把这种grib格式用grads处理生成的ctl文件,但是这种ctl无法被meteoinfo所识别,我用grads检验发现此ctl是有数据的,所以我认为meteoinfo无法识别这种方式生成的ctl,但我不知道如何处理这个grib格式的数据,请教一下王老师,麻烦您了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-12-6 23:15:24 | 显示全部楼层
ruanjinshuai 发表于 2016-12-6 22:45
王老师,您好,请问meteoinfo能批处理像FNL、GLDAS这样的grib数据吗?我这是一年的数据,但是分了2920个小g ...

你希望能怎样处理?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-6 23:50:24 | 显示全部楼层
因为我想做一个时间序列的分析,在时间上把它们连在一起
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-6 23:51:52 | 显示全部楼层
然后,再用meteoinfo画图,比如我想画一个月均的空间分布图,这如何实现?麻烦您了,王老师!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表