mm5说明书(中文版)
中尺度模式（Mesoscale Model 5 v3）用户手册
一、概述
1.mm5模式系统的结构
第五代中尺度模式mm5是近年来由美国大气研究中心（NCAR）和美国滨州大学（PSU）在mm4基础上联合研制发展起来的中尺度数值预报模式，已被广泛应用于各种中尺度现象的研究。Mm5在以往的模式基础上作了许多变化，主要有以下几点：1）复合区域嵌套功能，2）菲静力部分扩展3）四位数据同化功能以及较多的物理过程参数化，能够方便、广泛地应用于各种计算平台。这些变化使得许多工作在这一模式系统下建立起来。
图是整个mm5模式系统的结构框图，它表现了模式的模块次序、数据流程以及各模块主要功能的简短说明。TERRAIN和REGRID模块用来处理在麦卡托或兰博托或极射赤面投影下，地形数据和等压面气象数据从规则经纬网格点到高分辨可变中尺度区域的水平插值。由于插值不能提供全面的中尺度信息，因此插值数据必须加大，RAWINS/little_r就是用连续扫描Cressman客观分析方法和复合二次曲面技术来处理水平网格观测资料和无线电探空资料。INTERP模块处理MM5系统中气压坐标到sigma坐标的垂直插值，接近地面的sigma平面与地形相似，高水平sigma面与等压面近似。MM5模块是系统的核心部分，包含气象过程的主控程序，主要求解大气运动基本方程组。 INTERB模块与INTERP模块作用相反，主要是把MM5模块计算结果从sigma坐标插值到气压坐标中。
2．Mm5模式的水平和垂直格点
介绍模式的格点构造是非常有用的，模式系统通常是从等压面上获得、分析数据的，但是这些资料在进入模式之前不得不被插值到模式的垂直坐标中。垂直坐标是地形伴随的，也就是在底层水平网格伴随地形，而上层表面是平坦的。中间层是随着气压的减小趋向顶层气压逐渐变得平坦（如图）。用来定义模式水平层：
p是气压，p
t 是顶层气压，p
s
是表面气压。
从上图可以看出：在顶层等于0，在底层等于1，模式的每一水平层由值来定义，模式的垂直分辨率由0到1之间的数目决定，通常边界层的分辨率高于顶层分辨率，水平层数尽管原则上没有限制，但通常在10到40层之间变化。
水平格点采用Arakawa-Lamb B型跳点网格，如图，从图中可以看出T,q等量被定义在格点平面的中央，而向东风速u,向北风速v被配置到两边。模式中中央点被当作叉X，两边点被当作点·，因而水平风速被定义在点上。
所以上层变量被定义在模式垂直层的中央，作为半层由一条虚线表示出来（图），垂直速度在格水平层上表示为实线，在定义sigma坐标时所有层都要被列出，包括0层和1层，因此模式层的数量比sigma层少一层。
3.区域嵌套
mm5具有区域嵌套的功能，最多可以同时运转9个区域的嵌套和相互作用。图时其中的一种配置，对于双向交互作用的区域，它们的嵌套比率为3：1，所谓的双向就是嵌套区域边界条件的输入来自于粗糙区域，而嵌套区域对粗糙区域又有反馈。
从图可以看出嵌套区域被给定一个嵌套级别（如区域2，3），而且嵌套区域之间也可以相互覆盖。区域4是第三级别，说明他的网格尺寸和时步相当于区域1的1/9。每一个次区域都有一个母区域，而且被完全嵌入，区域2，3的母区域是1，区域4的是3。在模拟试验中嵌套区域可能随时被转动或停止，也就是无论何时母区域被停止，所有他的派生嵌套区域就被停止。在模拟试验中移动一个区域也是有可能的，条件是没有母区域作用嵌套网格和没有粗糙网格。
单向嵌套在模式中也是可能的，这是模式首先被运行产生一个输出，一旦单
向嵌套区域位置被定下来，边界条件文档也就生成了。典型的边界条件文档是一小时产生一次（依赖于粗糙区域输出的频率），这些数据被时间插值供给嵌套区域，因此单向嵌套不同于双向嵌套：1、没有反馈2、在边界处分辨率是粗糙临时的。单向嵌套区域也可能被初始化在增强数据分辨率和地形分辨率的条件下。
4.侧边界条件
在运行任何区域数值天气预报模式中，都要求有侧边界条件，mm5模式中所有的四个边界条件被指定为水平风速、温度、压强和湿度场，如果需要可能指定一些微物理场（比如云）作为边界条件。因此在运行模式前，边介值必须被设定。
边界值来源于对未来时刻的分析或是前粗糙区域的模拟值（单向嵌套），或是来源于其他母区域的实时预报值。对于实时预报侧边界条件最终依赖于全球模式预报。当然通过观察分析侧边界条件可能被被增强，同样的方法用于初始条件。上层空气分析值，边界值只能12小时一次，而对于模式产生的边界值可能更高频率比如6小时一次甚至1小时一次。
通过线性插值模式可以利用离散时间分析值，这些分析值能够全面地说明模式网格输出的行与列的状态。双向嵌套边界条件与单向是相似的，但每粗糙区域
一个时步他就会更新，而且没有松弛区域。
5. 四维数据同化
拓展时间上的的数据可以输入到模式中，四维数据同化只是模式中的一个选项。本质上四维数据同化允许模式在强迫条件下执行，也就是张弛朝着观测场和分析场。这样做得好处是一段张弛逼近后，模式将某种程度上适合那一时间间隔内的所有数据，而保留值接近动力平衡。这样做是优于一个简单天气分析初始化过程，因为一段时间内的附加数据将有效增加数据的分辨率。在一个位置上的观测数据通过模式将被带入下一阶段，从而可以填充下一时段数据的空缺。
四维数据同化的两个主要应用是：1、动力初始化2、四维数据设置。动力初始化就是在实时预报中利用前面的预报值是初始条件最优化。结果表明相对于静态初始化附加数据对预报是有益的。第二个应用，四维数据设置是一种产生动态平衡分析的方法，在一个拓展是段张弛逼近数据同化过程中这种方法包含流程中的实际连续性和地转风与热量风的平衡。
数据同化的两种方法依赖于数据是否是格点的还是单个观察的。考虑模式网格分析场的形成，在一个给定时间间隔格点数据将用来张弛逼近模式点对点数据，对于一个大尺度信息这是非常有用的。对于小尺度信息，非天气图定时资料或特殊平台象机床、飞行器等不能被利用，但个观测资料能够被用来张弛逼近模式。这些观测资料被给定一个时间窗和影响半径。每个格点上的观测资料的权重依赖于他们在时间上和空间上到观测站的距离。
6. 地表资料的分类
关于地表分类模式有三个选项设置，这些设置在TERRAIN模块的存档资料中，它们是13、16、24种，（植被、沙漠、城市、水、冰）。模式每个格点单元被分派一种地表资料，从而决定了表面特性如反射率、粗糙长度、长波反射率、热容量和湿度等。此外如果积雪覆盖着一项被设置，表面特性相应的被改动。列表中
的值相应季节变化。
7. 地图投影及地图比例因子
模式系统有几种投影方式可供选择，兰博托投影适合中纬度，极射赤面投影适合高纬度，麦卡托投影适合低纬度，模式的x,y方向不对应东西或南北，麦卡托投影除外，因此观测到的风俗不得不被旋转到模式中。模式的u,v在于观测支队比前需要旋转。这些变换在模式的前处理中就被计算出。
地图的比例因子m被定义为：
m = (distance on grid) / (actual distance on earth)
它的值通常接近于1，虽纬度变化，模式投影保持小区域的形状，以便每处dx=dy，但是网格长度是变化的，对应的区域在平面上表示为球面。无论那儿水平梯度被用到时，地图比例因子都会在模式的方程中被计算出。
8. 运行模式系统所需要的数据
因为mm5模式系统首先是为实时资料研究设计的，因而它需要以下资料来运转：
地形资料和地表分类资料
网格大气资料至少有以下几种：海平面气压、风、温度、相对湿度、位
势高度，气压分层：表面气压、1000、850、700、500、400、300、250、
200、150、100mb
观测资料
二、模式中提供了一些基本的地形、地表分类和全球覆盖但分辨率可变的植
被数据，关于怎样从NCAR上获得资料，请查阅网站: 模式分模块介绍
1、TERRAIN模块
1.1用途
mm5模式系统中的第一个模块是TERRAIN，这个模块水平插值规则经纬度上的地形高度和植被（地表分类资料）到指定的中尺度区域上，如图。如果地表模式LSM在MM5中应用时，附加资料如土壤类型、植被分布率、每年深层土壤温度将被产生。
1.1.1TERRAIN模块的任务
TERRAIN模块主要完成两个任务：
建立中尺度区域：粗糙区域和划分网格区域（移动嵌套区域除外）
为所有的中尺度区域产生地形数据场，这将首先被REGRID模块用到，
以后将被MM5和NESTDOWN模块使用。这一模块仍然计算一些模式需要
的常数：经纬度、地图比例因子、柯里奥利参数。
1.1.2TERRAIN模块流程
TERRAIN模块包括四部分：如图
原始数据输入
从规则经纬度的原始数据到中尺度网格的插值
嵌套区域的调整和反馈
地形高度、地表分类资料和其他地形数据以MM5格式的输出
输入数据
原始数据
作为输入到TERRAIN模块的可利用数据包括地形高度、地表分类资料/植被，海陆分界线、土壤类型、植被分布率和深层土壤温度。所有这些数据由六种分辨率：1度、30，10，5，2分、30秒。下面是一个通常可用数据的列表：原始高度资料有六种分辨率：1度、30，10，5，2分、30秒，所有低分辨率（1度到2分）数据产生于30秒的USGS资料
原始地表分类资料/植被有三种类型：
i）13种类，全球覆盖资料 1度、30，10分（3个文件夹）
ii）17种类北美覆盖资料 1度、30，10，5，2分、30秒（6个文件夹）iii）25种类全球覆盖资料 1度、30，10，5，2分、30秒（6个文件夹）两种类型的海陆分界线文件定义水体
i)17种类北美覆盖资料 1度、30，10，5，2分、30秒（6个文件夹）
ii)25种类全球覆盖资料 1度、30，10，5，2分、30秒（6个文件夹）对于mm5中的LSM，土壤类型、植被分布率和深层土壤温度是必需的，其原始数据为：
i)17种类 6个分辨率的全球土壤数据（6个文件夹）
ii)12个月，10分，全球植被分布（1个文件夹）
iii)1度全球每年深层土壤温度（1个文件夹）
数据格式
既然原始数据有不同的来源那么他们就有不同的格式，这些数据要被转化成一种标准格式以适应TERRAIN模块的使用。数据排列和格式如下：
纬度与纬度之间从北到南一度分割，数据点从西到东排列，通常开始经度为0度（或180度）
我们用两个字符存储地形高度和深层土壤温度，（最大值<215或32768）图一个字符存储其它数据（值<100）图，这将使得读数据更加方便。
所有原始数据文件都是直接路径
所有数据都被假定在格点中央上是有效的，因此有360*180个1度数据点，360*2*180*2个30分数据点，360*120*180*120个30秒数据点。
输入数据的来源和文件大小
高度
地形高度数据
分辨率数据来源覆盖范围大小（bytes）1 deg. km)USGS Global129,600
30 min. km)USGS Global518,400
10 min. km)USGS Global4,665,600
5 min. km)USGS Global18,662,400
2 min. km)USGS Global116,640,000
Tiled 30 sec. km)**GTOPO30
by
. Geological
Survey's
Global
(33 tiles:
40o lon. x 50o lat.
or
57,600,000
or
51,840,000
for each of
EROS Data Center in late 199660o lon. x 30o lat.)
tiles
30 sec. km)USGS Global***1,866,240,000
植被和地表分类资料
1)来自于PSU/NCAR全球13种类数据
PSU/NCAR Land-use Data
分辨率数据来源覆盖范围大小（bytes
1 deg. km)PSU/NCAR Global842,400
30 min. km)PSU/NCAR Global3,369,600
10 min. km)PSU/NCAR Global30,326,400
2）来自于USGS北美17种类数据
17-category Vegetation Data
数据来源覆盖面大小(bytes) 1 deg. km)Simple Biosphere model0o-90oN, 60o-180oW183,600
30 min. km)Simple Biosphere model0o-90oN, 60o-180oW734,400
10 min. km)Simple Biosphere model0o-90oN, 60o-180oW6,609,600 5 min. km)Simple Biosphere model0o-90oN, 60o-180oW26,438,400 2 min. km)Simple Biosphere model0o-90oN, 60o-180oW165,240,000 30 sec. km)Simple Biosphere model0o-90oN, 60o-180oW155
3）来自于USGS全球25种类数据
25-category USGS Vegetation Data
分辨率数据来源覆盖面大小(bytes)
1 deg. km)USGS Global1,620,000
30 min. km)USGS Global6,480,000
10 min. km)USGS Global58,320,000
5 min. km)USGS Global233,280,000
2 min. km)USGS Global1,458,000,000
30 sec. km)USGS Global933,120,000
海陆分界线
1）来自于SiB Vegetation 北美海陆分界线数据
Land-Water Mask Data
分辨率数据来源覆盖面大小(bytes) deg. km)SiB Vegetation0o-90oN, 60o-180oW10,800
0 min. km)SiB Vegetation0o-90oN, 60o-180oW43,200
0 min. km)SiB Vegetation0o-90oN, 60o-180oW388,800
min. km)SiB Vegetation0o-90oN, 60o-180oW1,555,200 min. km)SiB Vegetation0o-90oN, 60o-180oW9,720,000
0 sec. km)SiB Vegetation0o-90oN, 60o-180oW155,520,000
2）来自于USGS全球海陆分界线数据
USGS Land-Water Mask Data
分辨率数据来源覆盖面大小(bytes)
1 deg. km)USGS Vegetation Global64,800
30 min. km)USGS Vegetation Global259,200
10 min. km)USGS Vegetation Global2,332,800
5 min. km)USGS Vegetation Global9,331,200
2 min. km)USGS Vegetation Global58,320,000 30 sec. km)USGS Vegetation Global933,120,000
土壤
Global 17-category Soil Data
分辨率数据来源覆盖面大小(bytes) 1 deg. km)FAO+STATSGO Global1,101,600 30 min. km)FAO+STATSGO Global4,406,400 10 min. km)FAO+STATSGO Global39,657,600 5 min. km)FAO+STATSGO Global158,630,400 2 min. km)FAO+STATSGO Global991,440,000 30 sec. km)FAO+STATSGO Global933,120,000
地表分类指数
Global Monthly Vegetation Fraction Data
分辨率数据来源覆盖面大小(bytes)
10 min. km)AVHRR Global27,993,600
土壤温度
Global Annual Deep Soil Temperature Data
分辨率数据来源覆盖面大小(bytes)
1 deg. km)ECMWF analysis Global129,600
数据信息
如果一个用户哟不同的数据来源，那么这些数据必须被转化成标准格式和直接路
径的文件，此外下面的信息必须提供给TERRAIN模块，通过一个数据声明在
或 in , 和.
数据的分类
水分类的ID号
数据分辨率
初始经纬度
文件的总计数量
在一个纬度内的数据点数量
连接到FORTRAN单元的文件名地表分类和土壤种类
13种(PSU/NCAR)地表分类和物理参数
Landuse
Integer Identification Landuse
Description
Albedo(%)
Moisture
Avail. (%)
Emissivity
(% at 9 m m)
Roughness
Length (cm)
Thermal Inertia
(cal cm-2 k-1 s-1/2) Sum Win Sum Win Sum Win Sum Win Sum Win
1Urban land181851088885050 2Agriculture172330609292155 3Range-grassland1923153092921210 4Deciduous forest1617306093935050 5Coniferous forest1212306095955050 6
Mixed forest and
wet land
1414357095954040 7Water881001009898.01.01 8Marsh or wet land1414507595952020 9Desert25252585851010
10Tundra1570509092921010
11Permanent ice55709595959555
Tropical or sub
12
1212505095955050 tropical forest
13Savannah2020151592921515
17种SIB植被分类和物理参数Vegetation
Integer Identification Vegetation
Description
Albedo(%)
Moisture
Avail. (%)
Emissivity
(% at 9 m m)
Roughness
Length (cm)
Thermal Inertia
(cal cm-2 k-1 s-1/2) Sum Win Sum Win Sum Win Sum Win Sum Win
1Evergrn. Broadlf.1212505095955050 2Broadlf, Decids.1617306093935050 3Decids. Evergrn.1414357095954040 4Evergrn. Needlf.1212306095955050 5Decids. Needlf.1617306093935050 6Grnd. Tree Shrb.2020151592921515 7Ground only192315309292.12.10 8Broadlf.
DT. Shrb1570509092921010 11Bare Soil25252585851010 12Agricltr. or C3 Grs172330609292155 13Perst. Wetland1414507595952020 14Dry Coast Cmplx1923153092921210 15Water881001009898.01.01 16Ice cap & Glacier55709595959555 17No data
25种USGS植被分类和物理参数
Vegetation
Integer Identification Vegetation
Description
Albedo(%)
Moisture
Avail. (%)
Emissivity
(% at 9 m m)
Roughness
Length (cm)
Thermal Inertia
(cal cm-2 k-1 s-1/2) Sum Win Sum Win Sum Win Sum Win Sum Win
1Urban1818101088885050 2Drylnd Crop. Past.172330609292155 3Irrg. Crop. Past.182350509292155 4Mix. Dry/
Mosaic182325409292145 6Crop./Wood Mosc1620356093932020 7Grassland192315309292.12.10 8Shrubland2225102088881010 9Mix Shrb./Grs.2024152590901110 10Savanna2020151592921515 11Decids. Broadlf.1617306093935050 12Decids. Needlf.1415306094935050 13Evergrn. Braodlf.1212505095955050