NC数据打开 提取 合并的问题

lizhiboys

大家好 在grads连续打开两个NC文件  下面设置经纬度和层次后 D X.2  X代表变量  为什么不能画图？而D X.1 是可以的  

这个问题其实是个衍生出来的问题，因为最近要把365个文件（每个代表一天 NC格式）放在一个文件下面（让其时次为1-365） 用循环提取数据的生成grd的过程中发现 从第二个开始后面的数据都有问题，所以发现了上面的问题！
希望那个大家能帮帮忙，通过何种方式能把365个NC放在一起？  fortran可以直接读取NC格式的数据吗？我在论坛里也搜过了 发现帮助不大 如果这个问题解决了 可以极大的方便以后论坛里的初学者！再次感谢大家！

喵~

是不是循环哪里不对？？？

lqouc

gs发上来吧。

夏朗的芒果

LZ把gs传上来吧，大伙一块瞧瞧。

lizhiboys

'reinit'
i=1
while(i<=7)
'sdfopen I:\hou\'i'.nc'
'set gxout fwrite'
'set fwrite I:\hou\'i'.grd'
'set x 1 1440'
'set y 1 720'
'set t 1'
'set z 1'
'd sst.'i''
'disable fwrite'
i=i+1
endwhile

;

lizhiboys

喵~ 发表于 2013-5-15 16:29

                               
登录/注册后可看大图

是不是循环哪里不对？？？

已上传GS 我制作了7个文件来 做检查的

lizhiboys

lqouc 发表于 2013-5-15 16:41

                               
登录/注册后可看大图

gs发上来吧。

已上传GS 我制作了7个文件来 做检查的

lizhiboys

这样生成的是一个一个的GRD   我还做了另一种 把上面那个GS里的输出路径改成
'set fwrite I:\hou\SST.grd'
结果发现生成的grd配上CTL后 只有第一个时次可以画图 后面都是 undefined grid value

lizhiboys

夏朗的芒果 发表于 2013-5-15 16:43

                               
登录/注册后可看大图

LZ把gs传上来吧，大伙一块瞧瞧。

已发   大家帮我看看吧

lqouc

lizhiboys 发表于 2013-5-15 16:59

                               
登录/注册后可看大图

已上传GS 我制作了7个文件来 做检查的

'disable fwrte'后面加一个reinit，就不需要文件号了