GrADS 读 FY2G PRE AWX 出图时出现问题

听海落雪

我想用 FY2G有一小时降水估计资料画图，参考了论坛里的贴子，编写出了ctl和gs。在GrADS打开ctl，直接 'd pre'，能够出正常的图，但是运行我编写的gs文件出图明显有问题。希望大神们不吝赐教。

ctl 如下：

dset F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001_OTG_2017%m2%d2_%h200.AWX
title 2017 FY2G_PRE
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1001 linear 55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 48 linear 00z23JUN2017 1hr
vars 1
pre 0 -1,40,1 FY2G_PRE(m)
endvars



gs如下：

'reinit'

'open F:\Data\GrADS\2017\20170624\PRE\FY2G_PRE_001.ctl'

'set grads off'
'set grid off'

'set lon 96 106'
'set lat 25 33'

****** 1小时降水 ******
'set rgb 54 255 255 255'  
'set rgb 55 220 250 211'  
'set rgb 56 166 242 143'  
'set rgb 57 61 186 61'   
'set rgb 58 97 184 255'   
'set rgb 59 0 0 225'      
'set rgb 60 250 0 250'   
'set rgb 61 128 0 64'

********* 一页四图 *********

p=1
while(p<=4)

'run page q'p''
'set t 'p''

'set parea 0.5 9.0    1.0 8.0'
'set xlopts 1 5 0.26'
'set ylopts 1 5 0.26'
'set xlab on'
'set ylab on'
'set xlabs 96`3.|98`3.|100`3.|102`3.|104`3.|106`3.`1E'
'set ylabs 25`3.|27`3.|29`3.|31`3.|33`3.`1N'

'set gxout shaded'
'set csmooth on'
'set clevs 0 0.1 1.5 7 15 40 50'
'set ccols 54 55 56 57 58 59 60 61'
'd smth9(pre)'

'set mpdset country1'
'set map 15 1 2'
'draw map'
'set mpdset xnqy_sj'
'set map 15 1 2'
'draw map'
'set mpdset pxdq_dj'
'set map 8 1 5'
'draw map'

'query w2xy 102.2147369385 27.9398202896'
x1=subwrd(result,3)
y1=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.16'
'draw mark 1  'x1' 'y1'  0.25'

'query w2xy 103.4 32.7'
x11=subwrd(result,3)
y11=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.28'
'draw string 'x11' 'y11'  `2  'p+13'BJZ23JUN'

p=p+1
endwhile

************************************************

'printim F:\Data\GrADS\2017\20170624\pic\PRE_2314-2317.png white'

river

你把这句'set clevs 0 0.1 1.5 7 15 40 50'先注释掉看看

听海落雪

river 发表于 2017-7-7 21:13
你把这句'set clevs 0 0.1 1.5 7 15 40 50'先注释掉看看

还是不行，颜色变了。
之前看到一篇网页，是这样写的“……再说CTT，PRE等一些因存储格式导致grads不能顺利读出的产品，原因有老师在之前帖子中已分析。由于强迫症同类产品想出一致的图，我的做法比较简单粗暴，该写AWX成grads可读即可。给出matlab处理的部分脚本……”


http://bbs.06climate.com/forum.php?mod=viewthread&tid=46113

rabin_xu

AWX格式前面有说明信息，数据只能重新读写一下，过滤掉说明信息，否则grads读取不了

听海落雪

rabin_xu 发表于 2017-8-19 15:07
AWX格式前面有说明信息，数据只能重新读写一下，过滤掉说明信息，否则grads读取不了

您说是 CTL 文件？

rabin_xu

你可以看看3.1，那里比较清楚

rabin_xu

还要扣积分啊，那自己去这个连接下载吧
http://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=10

arctichare

简单粗暴的那个帖子是我写的。。。脚本不知道放哪儿了，现在FY2G PRE AWX读取也出岔子了{:eb303:}{:eb303:}{:eb303:}

arctichare

应该是整出来了，matlab脚本留在这里备份一下

list=dir(['./','FY2*.AWX']);   %必须与CTA对应
k=length(list);
for i=1:k
     str= strcat ('./', list(i).name);
end

file=list(i).name;
fid = fopen(file,'rb');
test = fread(fid,[1001 503],'ubit16');% 读取TBB数据
data=test(:,3:503);
for m = 1:1001
    for n = 1:501
        if(data(m,n)==64546)
        data(m,n)=0;
        end        
    end
end

听海落雪

rabin_xu 发表于 2017-8-21 08:30
你可以看看3.1，那里比较清楚

好的，谢谢