求助，grads做的平均输出成一个grd出错

心如止水的饭团

因为我wrf输出的是grd格式，每小时的数据，我现在想给做成每月的。也就是每小时的数据平均成每月的。主要是自己的fortran编程不行，所以想到了用grads的ave函数，然后把计算好的平均输出成一个grd文件，但是现在发现输出的文件数据有问题。
首先是每小时文件的ctl
dset ^201001/wrfd01_%y4-%m2-%d2_%h2.dat
options  byteswapped template
undef 1.e30
title  OUTPUT FROM WRF V3.4.1 MODEL
pdef  182 172 lcc  24.000  106.000   91.500   86.500  44.00000   4.00000  115.00000  45000.000  45000.000
xdef  578 linear   45.65011   0.20270270
ydef  387 linear  -17.81141   0.20270270
zdef   20 levels
   0.00250
   0.00800
   0.01500
   0.02350
   0.03400
   0.04650
   0.06100
   0.07850
   0.09950
   0.12500
   0.15550
   0.19200
   0.23550
   0.28800
   0.35150
   0.42700
   0.51750
   0.62650
   0.75750
   0.91450
tdef    774 linear 00Z01JAN2010      60MN
VARS   1
tk            20  0  Temperature (K)
end vars


这是我写的gs
'reinit'
'open wrf201001.ctl'
'mtk = ave(tk,t=1,t=774)'
'set gxout fwrite'
'set fwrite monthly/201001.grd'
'set x 1 182'
'set y 1 172'
'd mtk'
'disable fwrite'
;

如果是直接在grads里面做平均，然后画图是可以画的，但是把数据输出来，然后画图就有问题了
这是删除投影后的输出文件的ctl：
dset 201001.grd
options little_endian
undef -9.9900000E+08
title  OUTPUT FROM WRF V3.4.1 MODEL
xdef  182 linear   45.65011   0.20270270
ydef  172 linear  -17.81141   0.20270270
tdef    1 linear 00Z01JAN2010      1mon
zdef    1 levels 1000
VARS   1
tk            0  0  Temperature (K)
ENDVARS

现在我能猜到的原因应该是pdef的影响，请问这个有什么方法可以解决这个影响的吗?

river

输出的数据已经是等经纬投影的了，描述文件里不需要投影参数了

心如止水的饭团

river 发表于 2017-1-13 21:09
输出的数据已经是等经纬投影的了，描述文件里不需要投影参数了

我把投影参数删除后
dset 201001.grd
options little_endian
undef -9.9900000E+08
title  OUTPUT FROM WRF V3.4.1 MODEL
xdef  182 linear   45.65011   0.20270270
ydef  172 linear  -17.81141   0.20270270
tdef    1 linear 00Z01JAN2010      1mon
zdef    1 levels 1000
VARS   1
tk            0  0  Temperature (K)
ENDVARS
得到的图片也不对，得到的图片在一楼

river

心如止水的饭团 发表于 2017-1-13 21:30
我把投影参数删除后
dset 201001.grd
options little_endian

你自己看看原始资料的经纬度范围，和你计算平均的经纬度范围一致吗？你自己设置的'set x 1 182' 'set y 1 172' 你能对应上相应的经纬度吗？