MeteoInfo

摘要： LaTeX是排版常用的语法，科学计算软件中也常用它来写数学公式（比如MatLab, Matplotlib等），MeteoInfo通过调用JMathLaTeX库也可以实现这样的功能。LaTeX的语法介绍可以参考此网页：http://matplotlib.org/users/mathtext.html下... 阅读全文

摘要： subplot命令可以自动等间距分配多个坐标系（Axes），命令中有三个参数，前两个定义了行数和列数，第三个指定了当前的坐标系（Axes），绘图命令是作用在当前坐标系中的。脚本程序：def f(t): return exp(-t) * cos(2*pi*t)t1 = arange(0., 5.... 阅读全文

摘要： 绘图的时候首先要有坐标系（Axes），可以用axes命令创建，如果没有创建在绘图时会自动创建一个。参数里的position是用来置顶坐标系的图形（figure）中的位置的，通过位置置顶，可以将多个坐标系组合在一个figure中。位置是以figure为单位坐标（0 - 1之间），前两个参数是左下角的x... 阅读全文

摘要： 气象数据基本为多维数据（通常是4维，空间3维加时间维），只让数据中一维可变，其它维均固定即可提取一维数据。比如此例中固定了时间维、高度维、纬度维，只保留经度维可变：hgt = f['hgt'][0,[500],[4],[180,360]]。固定某一维可以用序号（从0开始），比如此例中的时间维为0，既... 阅读全文

摘要： 绘制格点数据的散点图，用scaterm函数。脚本程序：f = addfile('D:/Temp/GrADS/model.ctl')ps = f['PS'][0,(10,60),(60,140)]axesm()mlayer = shaperead('D:/Temp/map/country1.shp')... 阅读全文

摘要： 读取风场U/V变量数据，可以从U/V计算出风速：speed = sqrt(u*u+v*v)。quiverm函数用来绘制风场矢量图，参数中包括U/V变量，如果要绘制彩色风场还需要第三个变量，这里是风速speed，10是颜色的分级数，size=8指定了风场箭头的基准大小。脚本程序：f = addfile... 阅读全文

摘要： Maskout通常有两种类型：Maskout图形和Maskout数据。这里是Maskout图形的示例。需要用shaperead读取地图数据形成图层作为Maskout图层（这里是中国的行政区域china.shp），然后用masklayer函数对数据图层进行屏蔽，第一个参数是用来Maskout的图层，第... 阅读全文

摘要： FLEXPART是一个类似HYSPLIT的扩散模式，它输出的netcdf文件参照了WRF，可惜全局属性没有写全，比如只有一个投影名称（例如Lambert），没有相关的投影参数：中央经度，标准纬度等等。必须查阅WRF的头文件才能重建投影（为什么不照猫画虎把属性写全呢？）。数据的经纬度坐标是有的，但在L... 阅读全文

摘要： 此例中的降水文本文件下载自“中国气象科学数据共享服务网”（http://cdc.nmic.cn/sksj.do?method=ssrjscp），其实是ESRI的文本格点数据格式。对于文本数据，MeteoInfoLab中有asciiread函数来读取数据，第一个参数是文件名，然后可以有一些可选参数，h... 阅读全文

摘要： MODIS的气溶胶光学厚度（AOD）产品应用很广，数据可以在Giovanni上下载：http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html。有HDF, netCDF和文本格式，这里示例打开netCDF格式数据，提取AOD数据并绘图。脚本程... 阅读全文

摘要： 用ncep数据计算水汽通量散度的脚本。需要air, uwnd, vwnd和rhum变量。数据是4维数据，需要固定时间维和高度维，数据中纬度维的数据是反向的，因此读取时需要特殊的设置（::-1）。脚本中用到几个内置的函数：cdiff, hdivg和magnitude，和GrADS中同名函数的作用一样。... 阅读全文