GrADS读取若干nc文件转换并合并为二进制文件

速冻。八爪鱼

求助一下论坛各位大神：
现在我有若干NCEP的nc文件，比如mean daily air temperature和geopotential height。分别为air.1948.nc和hgt.1948.nc。
之前已经通过grads分别转换成二进制文件，命名为air.grd和hgt.grd。
现需要将两个文件合并，想到是否可以用转换成grd文件的方法打开两个nc文件分别display。测试gs文件如下：

"reinit"

"sdfopen /home/air.1948.nc"

"sdfopen /home/hgt.1948.nc"

"set gxout fwrite"

"set fwrite /home/combine.1948.grd"

"set x 1 144"

"set y 1 73"

tt=1

while (tt <= 366)

"set t "tt""

zz=1

while (zz <= 17)

"set z "zz""

"d air.1"

"d hgt.2"

zz=zz+1

endwhile

tt=tt+1

endwhile

"disable fwrite"

;

写了相对应的ctl，结果air变量是正确的，而hgt变量不正确。一开始认为是变量储存方式不对，后来发现新输出的combine.grd和air.grd文件大小一致，不知道怎么回事。

顺便求教grads里面二进制数据的存储方法。谢谢。

edwardli

我的建议用wgrib(2)

速冻。八爪鱼

edwardli 发表于 2014-8-30 07:33
我的建议用wgrib(2)

是wgrib这个？还是wgrib2？

沙颖凯

两个文件分开用不也行么。
GrADS的存放顺序是 时次 -> 变量 -> 层次。http://bbs.06climate.com/forum.php?mod=viewthread&tid=10541

river

能说一下合并成一个的原因吗？一共就两个文件，没必要费这个劲儿吧

hahahaha

zz=1
while (zz <= 17)
"set z "zz""
"d air.1"
zz=zz+1
endwhile
zz=1
while (zz <= 17)
"set z "zz""
"d hgt.2"
zz=zz+1
endwhile
改成这样，应该先把一个变量的层次都写完，再写下一个变量，存储顺序是x-y-z-v-t

速冻。八爪鱼

river 发表于 2014-8-30 08:17
能说一下合并成一个的原因吗？一共就两个文件，没必要费这个劲儿吧

总共有八个文件要合并，每个文件从1948-2014，所以把一年的合并了以后再用ctl批量处理。

速冻。八爪鱼

沙颖凯 发表于 2014-8-30 08:14
两个文件分开用不也行么。
GrADS的存放顺序是 时次 -> 变量 -> 层次。http://bbs.06climate.com/forum.php ...

好，谢谢指点grads储存，主要老板要合并我也没办法。

river

速冻。八爪鱼 发表于 2014-8-30 09:46
总共有八个文件要合并，每个文件从1948-2014，所以把一年的合并了以后再用ctl批量处理。

老板要求那就另说了。看一下6楼的回复，需要先写完一个变量，再写另一个

速冻。八爪鱼

hahahaha 发表于 2014-8-30 09:19
zz=1
while (zz

了解了谢谢