MeteoInfoLab脚本汇总贴

ruanjinshuai · 发表于 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)

ruanjinshuai · 发表于 2016-12-1 10:06:10

王老师这个提取的数据直接用meteoinfolab可以打开，数据我验证过提取的应该没问题，但是这个坐标轴的问题我没有解决，麻烦王老师了！

MeteoInfo · 发表于 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对原始数据进行处理得到月变化之类的结果数据（不会改变数据的投影）。

ruanjinshuai · 发表于 2016-12-1 19:59:26

好的，谢谢王老师

oydl1985 · 发表于 2016-12-5 15:29:57

必须收藏学习啊！赞！赞！赞！

ruanjinshuai · 发表于 2016-12-6 22:45:24

王老师，您好，请问meteoinfo能批处理像FNL、GLDAS这样的grib数据吗？我这是一年的数据，但是分了2920个小grb文档，如果可以的话，麻烦王老师给我发一个帖子就好，我没找到，麻烦王老师了！

ruanjinshuai · 发表于 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格式的数据，请教一下王老师，麻烦您了！

MeteoInfo · 发表于 2016-12-6 23:15:24

ruanjinshuai 发表于 2016-12-6 22:45
王老师，您好，请问meteoinfo能批处理像FNL、GLDAS这样的grib数据吗？我这是一年的数据，但是分了2920个小g ...

你希望能怎样处理？

ruanjinshuai · 发表于 2016-12-6 23:50:24

因为我想做一个时间序列的分析，在时间上把它们连在一起

ruanjinshuai · 发表于 2016-12-6 23:51:52

然后，再用meteoinfo画图，比如我想画一个月均的空间分布图，这如何实现？麻烦您了，王老师！

		自动登录	找回密码
密码			立即注册