摘要：MeteoInfo特别是MeteoInfoLab的推广需要写大量详细的帮助文档和示例程序，MeteoInfo原先的网站使用最原始的编写.html文件的方式来更新，效率实在太低，最近学习了一下Sphinx (http://sphinx-doc.org/) - Python文档创建工具，利用了reStr 阅读全文

摘要：读取CloudSAT HDF Swath数据，绘图分上下两部分，上面是时间和高度维的Radar Reflectivity Factor二维图，下面是卫星轨迹图。示例程序：# Add filef = addfile('D:/Temp/hdf/2010128055614_21420_CS_2B-GEOP... 阅读全文

摘要：在MeteoInfoLab界面中，图形的大小会随着它所在的窗口的大小改变而改变，在需要精确控制图中一些要素的位置的时候会比较困难，这时可以用figure函数的一些参数来控制图形版面大小。figure函数是创建一个新的图形窗体，其中的参数figsize=None，即缺省情况下图形窗体中的图形不指定大小... 阅读全文

摘要：打开包含站点填图的站点数据文件（比如micaps 1）之后，用文件对象的smodeldata函数获取StationModel数据对象，然后用stationmodel函数绘制站点填图图层。脚本程序：#Set data foldersbasedir = 'D:/MyProgram/Distributio... 阅读全文

摘要：天气现象符号分布图实际就是散点图，可以用scatterm函数绘制，但之前需要创建天气符号图例，用weatherspec函数。如果只需要绘制某些天气现象（比如雾、霾），可以给出相应的天气符号序号列表（可以在网上查找）。脚本程序：#Set data foldersbasedir = 'D:/MyProg... 阅读全文

摘要：绘制Stramline流线图的函数是streamline，需要两个变量（U/V分量或者风向/风速）。脚本程序：f = addfile('D:/Temp/GrADS/model.ctl')u = f['U'][0,[500],[10,60],[60,140]]v = f['V'][0,[500],[1... 阅读全文

摘要：尝试编写MeteoInfoLab脚本计算垂直螺旋度，结果未经验证。脚本程序：print 'Open data files...'f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc')f_vwnd = addfile('D:/Temp/nc/vwnd.2011.nc')... 阅读全文

摘要：尝试了用MeteoInfoLab编写计算水平螺旋度的脚本，结果未经验证。脚本程序：print 'Open data files...'f_uwnd = addfile('D:/Temp/nc/uwnd.2011.nc')f_vwnd = addfile('D:/Temp/nc/vwnd.2011.n... 阅读全文

摘要：SeaWiFS HDF Grid数据读取，特别是涉及到了文件的众多属性数据的读取，数据取对数后绘图。脚本程序：#Add data filef = addfile('D:/Temp/hdf/S1999001.L3m_DAY_CHL_chlor_a_9km.hdf')#Get data variable... 阅读全文

摘要：对于全球数据来说，经度要么是-180 - 180，要么是0 - 360，都会存在边界数据不连续的问题。比如0 - 360的数据，怎么得到 -20 - 30度的连续格点数据就是个问题（跨越了数据的经度边界），在MeteoInfoLab中可以用DimArray或MIArray的join函数来将两个数组合... 阅读全文

摘要：在MeteoInfoLab中增加了创建netCDF文件并写入数据的功能，这里利用合并多个netCDF文件为一个新的netCDF文件为例。1、创建一个可写入的netCDF文件对象（下面用ncfile表示），用addfile函数，第一个参数是文件名，第二次参数'c'表示创建新的netCDF文件。ncfi... 阅读全文

摘要：读取HYSPLIT输出的轨迹数据文件和相应时间的气象数据文件，生成轨迹图层，循环每条轨迹的节点，读出该节点的经度、纬度、气压、时间，通过对气象数据插值获得该节点的气象数据。脚本程序：#-----------------------------------------------------# Aut... 阅读全文