fortran插值micaps小时降水量值增大

zoe_lover

参考清风大哥的《fortran+grads站点作图详解》时，用第一类小时降水量资料插值，得到的降水量级有几百毫米每小时，明显与实际降水不符，以下是资料和gs，以及过程。麻烦请问哪里出来问题？

站点数据读入

zoe_lover

'reinit'
'open e:\zhandian\test\grid.ctl'
'open e:\zhandian\test\sta.ctl'
'set  lon 100 110'
'set lat 30 35'
'set grads off'
'set grid off'
'set mpdset  hires'
'define rgrid=oacres(g,rain.2)'
'set gxout shaded '
'set cmin 10'

*'set rgb 21 249 242 242'
*'set rgb 22 223 217 217'
*'set rgb 23 197 192 192'
*'set rgb 24 171 167 167'
*'set rgb 25 145 142 142'
*'set rgb 26 119 117 117'
*'set rgb 27 93 92 92'
*'set rgb 28 67 67 67'
*'set rgb 29 41 42 42'
*'set rbcols 21 22 23 24 25 26 27 28 29'
*'set clevs  10 20 30 40 50 60 70'

'd rgrid'
'set gxout contour'
'set cmin 10'
'cnbasemap rgrid'
'cbarn'
'printim e:\zhandian\test\fop.png white'
；

river

你在你的fortran程序里加上print 看看 你读进去的数据是不是正常的

zoe_lover

river 发表于 2018-3-10 23:15
你在你的fortran程序里加上print 看看 你读进去的数据是不是正常的

C:\Users\Administrator\Desktop\数据.png