用grads合并多个NC文件转化为grd文件出错

bobobiu

小白求教！我想将一个月的500hpa高度场数据的nc文件合并成一个grd文件，一个月31天，一天4个时次，总共124个文件，但是和合成的grd有问题，画出来一堆乱线条，但是又不知道gs文件哪里出了错.求大神告知
'reinit'
'set gxout fwrite'
'set fwrite C:\Users\lenovo\Desktop\hgt1.grd'
i=1
while(i<=124)
'sdfopen C:\Users\lenovo\Desktop\datafiles\d'i'.nc'
'set t 1'
'set x 1 360'
'set y 1 181'
'd hgt_p0_l100_gll'
'close 1'
i=i+1
endwhile
'disable fwrite'

river

一堆乱线？可否上图？每个步骤出错都可能导致问题，你怎么就肯定是合并的脚本有问题？你编的描述文件呢？描述文件编的有问题也可能会画图不对。你的缺测值设置的是什么？

bobobiu

dset C:\Users\lenovo\Desktop\hgt1.grd
title height
options byteswapped
undef 1E+20
xdef 360 linear  0     1
ydef 181 linear -90    1
zdef 1   linear  0     1
tdef 124 linear  00:00Z01JAN2016 6hr
vars 2
HGT_P0_L100_GLL0=>hgt_p0_l100_gll  0  t,y,x  Geopotential height
initial_time0_encoded=>initial_time0_e  0  t  initial time encoded as double
endvars

ctl我是用q命令查询之后写的
grd文件我在meteoinfo里面看不到时次

bobobiu

river 发表于 2017-9-17 12:52
一堆乱线？可否上图？每个步骤出错都可能导致问题，你怎么就肯定是合并的脚本有问题？你编的描述文件呢？描 ...

我已经贴上来了 麻烦指教

风之精灵

试试将'sdfopen C:\Users\lenovo\Desktop\datafiles\d'i'.nc' 修改成'sdfopen C:\Users\lenovo\Desktop\datafiles\d'%i%'.nc'呢

junyang517

利用fwrite处理数据，曾经也好几个问题困扰过我。据我的经验，分享给你：
1.undef设置不正确，GrADS的fwrite输出的默认undef值好像是-9.99E8。
2.可能是网格点不匹配。你查看一下二进制文件属性，看看其大小是否等4*124*2*360*181。如果不等于，你在看看nc格式文件的ctl。用sdfopen打开nc格式的文件，qctlinfo查看，审核他的x、y定义，需要一致。

river

bobobiu 发表于 2017-9-17 13:24
我已经贴上来了 麻烦指教

options byteswapped
这句什么意思？不是随便乱加的