(grads)关于底图画上地形的问题

樱花下的一场雨

   现在数据都有，cn，cnworld等地形文件都有，但是不知道怎么编写ctl和gs，使地图成为西南地区的地图。而且本身就是格点！

aiai1992

我觉得你调用地图文件后，draw map 一下就好了，你想显示哪个区域就设一下经纬度，如果你想画某个省的底图，你看一下论坛的maskout ，应该能够解决你的问题

樱花下的一场雨

aiai1992 发表于 2014-4-24 18:11
我觉得你调用地图文件后，draw map 一下就好了，你想显示哪个区域就设一下经纬度，如果你想画某个省的底图 ...

地图是set mpdset cnworld cnriver 之后draw map 再设置我要的x 21 34，y 97 110 。就可以了么？但是我现在发现一个问题，就是为什么只有中间一部分有数据，周围都是白的，你看看图

樱花下的一场雨

aiai1992 发表于 2014-4-24 18:11
我觉得你调用地图文件后，draw map 一下就好了，你想显示哪个区域就设一下经纬度，如果你想画某个省的底图 ...

你看看为啥

aiai1992

樱花下的一场雨 发表于 2014-4-24 22:21
你看看为啥

周围是空白就是因为你先draw map 了，你先画的底图，之后画的数据，你把顺序换一下试试

樱花下的一场雨

aiai1992 发表于 2014-4-25 08:19
周围是空白就是因为你先draw map 了，你先画的底图，之后画的数据，你把顺序换一下试试

'reinit'
'open d:/1/700/quqv700.ctl'
'set grid off'
'set grads off'
'set gxout stream'
'set cterp on'
'set csmooth on'
'set cthick 6'
'd hdivg(qu,qv)'
'set mpdset cnworld'
'set map 1 3 1'
'draw map'
'set mpdset cnriver'
'set map 4 1 10'
'draw map'
'set mpdset chongqing sichuan guizhou yunnan xizang'
'draw map'
'draw title water fiux divergence at 700hpa'
'printim d:/1/700/qdiv700.png white'
;

这是变完之后的，还是一样

aiai1992

数据CTL你检查了吗？

樱花下的一场雨

aiai1992 发表于 2014-4-25 21:16
数据CTL你检查了吗？

我知道为什么了额，hdivg这个函数本来就很不稳定

飞鹰

我想请问高人一下你的底图是在吗画的？

river

樱花下的一场雨 发表于 2014-4-26 09:41
我知道为什么了额，hdivg这个函数本来就很不稳定

怎么叫这个函数不稳定呢，有没有去看过这个函数求散度时的原理？如果用fortran来求，使用的什么方法？他们都是一样的原理，用了差分，无论是钱差还是后差，最终在边界上都是没办法差分出值的，所以才会缺一块。你在求散度的时候求一个大范围的，然后画图的时候再限制到你要的范围就不会出现现在这种情况了