- 积分
- 900
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-14
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-12-3 15:42:36
|
显示全部楼层
俺是这么改的,lats = (/ 34.0, 39 /)
lons = (/ 114.0, 123.0 /)
loc = wrf_user_ll_to_ij(a, lons, lats, True)
x_start = loc(0,0) - 1
x_end = loc(0,1) - 1
y_start = loc(1,0) - 1
y_end = loc(1,1) - 1
mpres1@ZoomIn = True ; set up map info for zoomed area
mpres1@Xstart = x_start
mpres1@Ystart = y_start
mpres1@Xend = x_end
mpres1@Yend = y_end
u1_zoom = u1(y_start:y_end,x_start:x_end)
v1_zoom = v1(y_start:y_end,x_start:x_end)
tc2_zoom = tc2(y_start:y_end,x_start:x_end) |
|