MeteoInfoLab脚本汇总贴

ruanjinshuai · 发表于 2016-11-30 21:53:37

恩恩，解决了，谢谢王老师

ruanjinshuai · 发表于 2016-11-30 22:24:55

王老师您好，请教您一下，我用WRFPUT转出来的数据画出T2的月均分布图为F:\xuexi\WRFOUT\comparativeanalysis\T2\T2-d02.png；然后，我先把这套数据用grads提取出来做了一个月平均，然后再用meteoinfo画图，结果其坐标轴不对，如图所示F:\xuexi\WRFOUT\comparativeanalysis\T2\T2-d02monthly.png。

ruanjinshuai · 发表于 2016-11-30 22:27:00

图没上去，这个是图

ruanjinshuai · 发表于 2016-11-30 22:29:45

王老师，就是第二张图的这种扇形图形始终改不掉！我添加了proj=f.proj,还是没有改变！麻烦您了！

MeteoInfo · 发表于 2016-11-30 22:41:17

ruanjinshuai 发表于 2016-11-30 22:29
王老师，就是第二张图的这种扇形图形始终改不掉！我添加了proj=f.proj,还是没有改变！麻烦您了！

或许数据本身就是这样

ruanjinshuai · 发表于 2016-12-1 09:29:34

王老师，这没懂，首先我用grads画图不会出现这个所谓扇形的图形，然后我用grads提取数据再用meteoinfo画图会出现这个情况，但本身这个数据直接用meteoinfo画图也不会出现这个扇形，所以我不知道这个情况是为啥？难道是我在提取数据出现的问题吗？

MeteoInfo · 发表于 2016-12-1 09:46:08

ruanjinshuai 发表于 2016-12-1 09:29
王老师，这没懂，首先我用grads画图不会出现这个所谓扇形的图形，然后我用grads提取数据再用meteoinfo画图 ...

原始数据应该是Lambert投影（矩形），GrADS在提取数据的时候会将数据投影到等经纬度（扇形）。

也就是说原始数据在Lambert投影的地图中绘制是矩形，在等经纬度的投影中绘制是扇形。

ruanjinshuai · 发表于 2016-12-1 09:49:17

那老师你的意思是说我在meteoinfo中把投影方式改为Lambert投影吗？有个疑问为啥用grads用提取出来的数据画图也是矩形，

MeteoInfo · 发表于 2016-12-1 09:50:59

ruanjinshuai 发表于 2016-12-1 09:49
那老师你的意思是说我在meteoinfo中把投影方式改为Lambert投影吗？有个疑问为啥用grads用提取出来的数据画 ...

要不你把数据文件贴出来看看

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)

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