利用探空资料和micaps第五类数据格式绘制历史探空图

流觞曲水

program&data.rar (82.37 KB, 下载次数: 375, 售价: 1 贡献)

2013-7-28 00:49 上传

点击文件名下载附件        探空图可以用grads画，偶然想到也可以利用micaps里自带的tlogp图的功能来画，只要数据格式对，应该也可以实现，micaps画的探空图看着更漂亮舒服。这里利用怀俄明大学共享的探空数据，如图1，网址：http://weather.uwyo.edu/upperair/seasia.html，点开可以查到1973年至今全球任意一站的探空数据(前提是该站有数据)，利用这些数据可以画和micaps一样漂亮的温度对数压力图。

图1：

  

       打开网页后，在Type of plot里选择“Text.list”，其他按自己需求选择。研究了一下怀俄明大学的这个探空数据，如图2，共11列，从左到右依次为：气压(hPa)、位势高度(位势米，gpm)、温度(℃)、露点温度(℃)，相对湿度(%)、比湿(g/kg)、风向(0-360度，顺时针，0度为北风)、风速(0.5m/s)、位温(K氏温标)、假相当位温(K)和静力总温度(K)。行数上，每一层为一行，从地面直到气球爆炸时的高度。

图2：

       再看micaps第五类数据格式的数据，如图3或图4的右栏，文件头第一行分别为diamond 5(第5类数据格式)和数据说明；第二行分别为年、月、日、时次(北京时)、该文件里总站点数；第三行分别为：区站号、经度、纬度、测站海拔高度、该站的数据个数(行数×列数)。如图3或图4的左栏，micaps探空数据只有6列，从左到右依次为：气压、位势高度(10gpm)、温度、露点温度、风向和风速，并且micaps数据除了地面气压(第一行数据)、1000hPa、925hPa、850hPa、700hPa、500hPa、400hPa、300hPa、250hPa、200hPa、150hPa和100hPa外(若测站海拔较高，就没有1000hPa、925hPa···)，其它层的位势高度、风向和风速都是缺测的(缺测值用9999表示)。行数上，只有从地面直到100hPa的高度。micaps第五类数据格式用到的只是怀俄明大学探空数据的第1、2、3、4、7、8列的数据，而且位势高度、温度和露点温度的精度要低一级。显然，怀俄明大学的这个探空数据要比micaps的探空数据的信息量要全，包括层次、要素和数据精度。
图3北京(左：micaps数据，右：程序结果)：
  
图4上海(左：micaps数据，右：程序结果)：
   

       具体步骤：如图2所示，将地面到100hPa的数据复制到txt里，命名为tlogp.txt，别忘记在最后一行末尾按下回车。现编写一个fortran小程序，将怀俄明大学的探空数据转换成和micaps的要素存放顺序、类型一样、micaps可读取的数据格式(即micaps第五类数据格式)，缺测值用9999表示。程序及说明如图5、6。
图5、6：


      
        这里利用2008年8月8日08时北京(54511)和1998年8月8日08时上海(58362)的探空数据试验，程序运行结果即图3、图4的右栏。对比图3或图4的左右栏，micaps数据里缺测层的位势高度、风向和风速，程序结果里都有。搞到这两个时次micaps的探空数据作图再来验证，一样！如图7(北京程序结果出图)和图8(北京实况结果出图)，图9(上海程序结果出图)和图10(上海实况结果出图)。方便的也可以用当天的high文件里的探空数据检验该程序。
图7(北京程序结果出图)：
  
图8(北京实况结果出图)：


图9(上海程序结果出图)：

图10(上海实况结果出图)：


       层次资料越多，tlogp图越光滑好看，根据需要也可将100hPa层以上的加进去，但探空图一般画到100hPa就够了。fortran程序及数据有需要的可下载，程序上面也已贴出。

chenhongfengye

请教版主一个问题，怀俄明大学共享的探空数据和这个网址http://www.ncdc.noaa.gov/oa/climate/igra/是不是一样的啊？

river

谢谢分享，再出一个将怀俄明大学探空资料转换成grads可读资料来画探空图就更完美了

mofangbao

这个真不错，值得好好学习一下

kongfeng0824

感谢分享！！！！！！！！！

又是那隻貓

太棒了 收藏

mbt125

师兄，牛b啊

jqh99999

流觞曲水原来是师兄啊，今天才知道，好迟钝啊。。。

走走

图画得很漂亮

岳焕

这个很好，正在找这个数据呢