nc文件用fwrite提取数据后画图不一致

樱花阿喵

从1960-2014年的nc格式的逐日降水数据中用grads提取grd格式的1995-2014年的6-8月数据，结果用nc文件画出的图是正常的，用提取的文件画出的图只有粗边界，grd文件的ctl文件都是按照nc文件里的数据描述文件写的，现在正在写毕业论文，还请大家帮忙看一下
下面是提取数据的gs文件：
'reinit'
'sdfopen e:\pr\historical\pr_Amon_ACCESS-CM2_historical_1x1regrid_196001-201412.nc'
'set gxout fwrite'
'set fwrite e:\pr\historical\access-cm2_9514summer.grd'
tt=420
while(tt<660)
'set t 'tt+6''
'd pr'
'set t 'tt+7''
'd pr'
'set t 'tt+8''
'd pr'
tt=tt+12
endwhile
'disable fwrite'
;


下面是提取文件的ctl文件：
dset e:\pr\historical\access-cm2_9514summer.grd
title month precipitation
undef 1e+20
xdef 51 linear 60 1
ydef 31 linear 20 1
zdef 1 linear 0 1
tdef 60 linear 12Z16JUN1995 1mo
vars 1
pr=>pr  0  t,y,x  Precipitation
endvars

樱花阿喵

问题已经解决了，感谢大家的关注。
我在网上找到了一篇帖子，https://www.cnblogs.com/jiangleads/p/8284326.html
大致内容是：画图出错原因，GrADS用fwrite输出时，默认的缺测值是-9.9900000E+08，如果输入源的缺测值与之不一致，则会出现两种缺测值的情况。解决的办法：可以在fwrite 声明后面设置默认缺测值设置。

carlye2002

樱花阿喵 发表于 2021-2-26 10:10
问题已经解决了，感谢大家的关注。
我在网上找到了一篇帖子，https://www.cnblogs.com/jiangleads/p/82843 ...

你好，怎么声明缺测值呢？在grads手册中没找到相关内容

樱花阿喵

carlye2002 发表于 2021-3-8 17:27
你好，怎么声明缺测值呢？在grads手册中没找到相关内容

set undef  xx ,我发的那个网址里也有示例的

泡沫学天气

请问一下可以提取为TXT格式的数据嘛？