对于全球数据来说，经度要么是-180 - 180，要么是0 - 360，都会存在边界数据不连续的问题。比如0 - 360的数据，怎么得到 -20 - 30度的连续格点数据就是个问题（跨越了数据的经度边界），在MeteoInfoLab中可以用DimArray或MIArray的join函数来将两个数组合并为一个，参数分别是另一个数组和合并的维的序号，比如下面例子中的二维数组（y, x），经度维是第二维（序号为1），也就是说两个数组沿着经度维合并。还需要给合并后的数组的经度维设置新的经度值。

脚本程序：

f = addfile('D:/Temp/GrADS/model.ctl')

psv = f['PS']

ps = psv[0,[10,50],[340,360]]

ps1 = psv[0,[10,50],[0,30]]

ps = ps.join(ps1, 1)

n = ps.dimlen(1)

ps.setdimvalue(1, linspace(-20,30,n))

axesm()

mlayer = shaperead('D:/Temp/map/country1.shp')

geoshow(mlayer, edgecolor=(0,0,255))

layer = contourfm(ps, 20)

title('Pressure')

colorbar(layer, aspect=12)

#savefig('D:/Temp/test/map_test.png', 600, 400)

最新文章

1. Getting Started With Hazelcast 读书笔记(第一章)
2. [问题2014A06] 复旦高等代数 I（14级）每周一题（第八教学周）
3. JS常用的腳本庫--包括在線編輯器
4. 利用JDK的中Proxy动态代理实现Spring的AOP技术
5. Android的那些轮子
6. linux下利用openssl来实现证书的颁发（详细步骤）
7. Java中的ExceptionInInitializerError异常及解决方法
8. Gray Code 解答
9. Kafka学习(一)配置及简单命令使用
10. 20190409-层叠の层叠上下文、层叠水平、层叠顺序、z-index、伪元素层叠
11. python pip 安装模块步骤
12. 一些面试题（关于string的）
13. rocketmq自启动配置
14. P3147 [USACO16OPEN]262144
15. RxJava + Retrofit
16. Java不定参数
17. POJ 1260 Pearls （斜率DP）题解
18. Django-工程搭建
19. 最短路问题（dijkstral 算法）（优化待续）
20. leetcode680

热门文章

1. 用Nginx将web请求引导到本机两个tomcat
2. 亚马逊DRKG使用体验
3. Node.js向MongoDB中插入并查询数据
4. turtle空间坐标系
5. 关于JAVA的一些零碎小知识
6. springboot中关于Long类型返回前端精度丢失问题处理
7. 配置静态 IP、网卡命名规范
8. Java基础一篇过（二）泛型
9. 刷题[FBCTF2019]Event
10. CCNP：重发布及实验