micaps资料处理成nc资料后的一些问题

766168584

最近由于需要用micaps数据并利用grads作图，果断用了Meteoinfo这个强大的好软件，将我的T639micaps资料（海平面气压场）和实况气压场转换成了NC格式数据资料，但是却没法将这两个nc资料画在同一张图上。就像问题截图那样能放在一起就好了



。
经过认真研究发现，这两个数据不但经纬度的范围不同，而且格点的分辨率也不一样，看了很多帖子，似乎方法不对，也或许自己在这个里面大海捞针一样，帖子太多，我搜索了很多关于grads  NC之类的关键词查找。最后还是没法解决，希望那个高手能够助我一臂之力，给我指点个方向，或者给我说说我应该仔细看谁的帖子都行，由于本人Fortran之类的都不太熟悉，但是一直在努力看。。希望能悉心指导，在此表示感谢。下面是我fwrite的程序


'reinit'
'sdfopen I:\wjh\wjh\639dm13061620.nc'
'set fwrite I:\wjh\wjh\639dm13061620.grd'
'set gxout fwrite'

'set lon 50 120'
'set lat 20 70'
'set z 1'
'set t 1'

  'd dust'
'disable fwrite'

但是这个程序运行后，grd文件0kb。



我也用了其他方法，用grads2ASCII这个方法转换了nc成grd格式资料，但是却没法画图。这几个问题一直困扰我。现在贴出这个转换后我写的ctl文件
dset 639dm13061720.grd
title NCEP data
undef -9.9E33
xdef 72 linear 0.0 2.5
ydef 36 linear 0 2.5
zdef 1 levels 1000
tdef 1 linear 12Z17JUN2013 1mo
vars 1
dust 0 99 Surface Pressure
endvars


希望能得到帮助。

mofangbao

楼主现在的实际情况是怎样？能分别出图吗？你为啥要用frwite输出数据？如果只是为了在同一个图片中输出两张图，而不是把两种数据叠加一起画，那是不要求分辨率什么之类的保持一致的，只要合理设置parea和vpage即可，当成是两次画图，只是画在屏幕的不同地方。画完第一张，关掉数据，然后再画第二张，你可以参考下这个帖子：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=5505

766168584

mofangbao 发表于 2014-11-5 09:12
楼主现在的实际情况是怎样？能分别出图吗？你为啥要用frwite输出数据？如果只是为了在同一个图片中输出两张 ...

谢谢清风，终于解决了！

一个不小心

看其他文献的时候发现有的图只有等值线没有底图，原来是用的MicAPS资料啊，长见识了

mofangbao

一个不小心 发表于 2014-11-5 20:57
看其他文献的时候发现有的图只有等值线没有底图，原来是用的MicAPS资料啊，长见识了

想多了，数据和底图是两码事，micaps数据也可以加底图的，别的数据也可以不加底图的

一个不小心

mofangbao 发表于 2014-11-6 09:15
想多了，数据和底图是两码事，micaps数据也可以加底图的，别的数据也可以不加底图的

我一直以为其他的软件实现不了呢，从没试过不要地图的画法