中国自动站与CMORPH融合的逐时降水量0.1°网格数据用grads无法画图

一顾留白

下载了融水数据想要批量累加降水，但是画不出图来，以下是我的ctl文件：
DSET   ^SURF_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-201001%d2%h2.grd   
UNDEF  -999.0
OPTIONS   little_endian  template
TITLE  China Merged Analysis Daily Working File         
XDEF   700 LINEAR  70.05 0.1      
YDEF   440 LINEAR  15.05 0.1      
ZDEF     1 LEVELS 1  
TDEF 999999  LINEAR 00Z21Jan2010 1hr      
VARS 2                  
crain      1 00  CH08   combined analysis (mm/Hour)
gsamp      1 00  CH02   gauge sampling
ENDVARS
file:///C:/Users/linl/AppData/Roaming/Tencent/Users/1245975149/TIM/WinTemp/RichOle/D%7DQU)$E%601H%60QQ~LXNOAOWVQ.png
file:///C:/Users/linl/AppData/Roaming/Tencent/Users/1245975149/TIM/WinTemp/RichOle/D%7DQU)$E%601H%60QQ~LXNOAOWVQ.png


输入命令 d crain
显示Cannot contour grid - all undefined values
print后图像显示

一顾留白

已解决，需要建立循环，读对时次才能读出数据

ldwloveyyz

nsw704406795

GRD数据文件如果不是和ctl文件在同一文件夹内，GRD数据要用绝对路径打开

一顾留白

是的呢~谢谢啦~

Joan0125

求问楼主这个怎么解决的呀 我和同学也碰到了这个问题
这个循环怎么写呀

一顾留白

Joan0125 发表于 2019-4-26 21:15
求问楼主这个怎么解决的呀 我和同学也碰到了这个问题
这个循环怎么写呀

'open D:\HUANAN\2011.7.15-17\CPAS_MERG_HLY_WORKING.lnx.ctl'
'define rain24=0.0'
i=1
while(i<=49)                    （这个根据时次改一下，要和ctl里的对应起来）
'set t 'i''
'rain24=rain24+crain.1(t='i')'
i=i+1
endwhile
  
'set lon 104 122'
'set lat 18 30'

'set grads off'
'set grid off'
'set gxout shaded'


'set rgb 16 255 255 255'
'set rgb 17 173 255 47'
'set rgb 18 154 205 50'
'set rgb 19 0 255 0'
'set rgb 20 34 139 34'
'set rgb 21 255 255 0'
'set rgb 22 255 165 0'
'set rgb 23 255 0 0'
'set rgb 24 238 130 238'


*'set ccols  0 16 17 18 19 20 21 22 23 24'
*'set cterp on'
*'set clevs  0 50 70 90 110 130 150 250'

'set ccols  0 16 17 18 19 20 21 22 23 24'
'set clevs  0 80 120 160 200 250 300 350'
'd rain24'
'cbar.gs'
'set mpdset cn cnriver'
'printim D:/2011.gif white'
'reinit'

这个是我的程序，应该可以直接运行的，你可以试一下

Joan0125

谢谢楼主！这是gs文件么
请问要做24h累计降水怎么弄呀实在头疼
在哪可以添加读取多个文件呀 老师给的是每小时的

一顾留白

Joan0125 发表于 2019-5-2 18:38
谢谢楼主！这是gs文件么
请问要做24h累计降水怎么弄呀实在头疼
在哪可以添加读取多个 ...

这个就是求取累计降水的gs文件，具体时次可以修改i的值，循环累加，另外ctl文件中也要注意修改t的值

小汪小汪

您好！请问您那里还有中国自动站与CMORPH融合的逐时降水量0.1°网格数据吗?