ncl处理grd数据格式 但是图是空白的

尽头的尽头

有可能是缺测值的问题

subtropical

尽头的尽头 发表于 2019-7-11 21:19
有可能是缺测值的问题

但是也设置了缺测了，除了缺测还有很多数值是正常的，但是也没画出来，真的不知道怎么解决了。。。。

尽头的尽头

subtropical 发表于 2019-7-12 08:32
但是也设置了缺测了，除了缺测还有很多数值是正常的，但是也没画出来，真的不知道怎么解决了。。。。

先用grads直接画图看是否正常

subtropical

尽头的尽头 发表于 2019-7-12 10:22
先用grads直接画图看是否正常

用grads画图是正常的样子

兰溪之水

subtropical 发表于 2019-7-11 19:25
谢谢兰溪大神指点 我看了之后 试着写了 也显示生成了nc文件，但是我用ncl画图的时候还是不对
这是写一个 ...

不要自己设缺省值，grads出来已经带有缺省值。

subtropical

兰溪之水 发表于 2019-7-13 12:10
不要自己设缺省值，grads出来已经带有缺省值。

谢谢大神 我把设置缺测那句去掉了 但结果还是一样。。。。

尽头的尽头

subtropical 发表于 2019-7-12 19:22
用grads画图是正常的样子

那应该是缺测值的问题，可以用grads写成nc数据，这样会避免一些麻烦

subtropical

尽头的尽头 发表于 2019-7-15 09:48
那应该是缺测值的问题，可以用grads写成nc数据，这样会避免一些麻烦

谢谢 我试着写了
这是写一个时次一个变量nc的gs：
'reinit'
'open D:\work2019\wrf_dat_read\wrf02d2018072612.ctl'
'set gxout fwrite'
'set x 1 747'
'set y 1 606'
'set lev 850'
'set t 1'
'define dbz1 = dbz'
'set sdfwrite d:\work2019\wrf_dat_read\dbz1.nc'
'sdfwrite dbz1'
'clear sdfwrite'
;
这是ncl读取脚本：
begin
min_lat =  28.5
max_lat =35.0
min_lon = 113.0
max_lon = 121.0
f1=addfile("d:/work2019/wrf_dat_read/dbz1.nc","r")
print(f1)
   lat=f1->latitude
lon=f1->longitude
;time   = f1->time
dbz = f1->dbz1(:,:)
print(dbz({30:35},{114:120}))

结果还是有问题，，，，T_T

尽头的尽头

subtropical 发表于 2019-7-15 10:12
谢谢 我试着写了
这是写一个时次一个变量nc的gs：
'reinit'

那就直接grads画图吧

subtropical

尽头的尽头 发表于 2019-7-15 21:02
那就直接grads画图吧

好吧。。。还是想弄清楚为什么不能写出来读，之前也遇到了这个问题，一直没解决。。。。