- 积分
- 15624
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-8
- 最后登录
- 1970-1-1
|
2金钱
本帖最后由 李文帅 于 2021-11-11 11:18 编辑
根据“GrADS站点文件作图详细解决方案.pdf”,利用Grads站点作图,脚本如下:
'reinit'
'sdfopen D:\Research_Room\Data\ThirdGroup\nc\multiple\ERA5_2018.12_thpw.nc'
'open D:\Research_Room\program\station_plot\PM10\campaign201812\china_sites_2018120100.ctl'
'open D:\Research_Room\program\station_plot\PM10\campaign201812\grid_1.0.ctl'
'set grid off'
'set grads off'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set map 1 1 2'
'set mpdset hires'
'set lon 108 128'
'set lat 30 45'
'set ylabs 30N|35N|40N|45N'
'set xlabs 108E|113E|118E|123E|128E'
'set gxout shaded'
'set t 1'
'd oacres(g.3,pm10.2)'
'set t 8'
'set lev 1000'
'set gxout contour'
'set clab masked'
'set clopts 2 6 0.1'
'set ccolor 1'
'set cthick 5'
'd u.1'
;
这是能出图的,
但只要更改一下气象数据,即将第2行的“ERA5_2018.11_thpw.nc”改为“ERA5_2018.12_thpw.nc”,
就画不出来了站点数据了,如下图:
这两个气象数据只有月份不一样,其他参数(层数,经纬度,变量)完全一样,这是为什么?
|
最佳答案
查看完整内容
后面站点数据和网格文件的ctl中时间改为11月应该就行了。
当然最正确的做法就是在画完了站点插值之后,close这两个站点文件再打开nc文件画图,这样两类数据的维度互不影响。
|