geo_em.nc文件全局变量解释+Un/Stagger网格解释

// global attributes:全局属性用于存储关于整个netCDF文件的元数据信息,描述了文件的一些基本特征和属性
		:TITLE = "OUTPUT FROM GEOGRID V3.9" ;表示这是GEOGRID V3.9的输出
		:SIMULATION_START_DATE = "0000-00-00_00:00:00" ;模拟起始时间
		:WEST-EAST_GRID_DIMENSION = 100 ;模拟域的网格信息网格尺寸：东西
		:SOUTH-NORTH_GRID_DIMENSION = 105 ;模拟域的网格信息网格尺寸：南北
		:BOTTOM-TOP_GRID_DIMENSION = 0 模拟域的网格信息网格尺寸：垂直
		:WEST-EAST_PATCH_START_UNSTAG = 1 ;模拟域在东西方向未stagier网格的起止索引,这里是从1到99
		:WEST-EAST_PATCH_END_UNSTAG = 99 ;WRF模型中,水平网格分为unstagger和stagger两种,unstagger网格用于标量变量(如温度、压强等),stagger网格用于矢量变量(如风场)。两种网格在空间位置有一个半格距的偏移
		:WEST-EAST_PATCH_START_STAG = 1 ;模拟域在东西方向stagger网格的起止索引,这里是从1到100
		:WEST-EAST_PATCH_END_STAG = 100 ;
		:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
		:SOUTH-NORTH_PATCH_END_UNSTAG = 104 ;
		:SOUTH-NORTH_PATCH_START_STAG = 1 ;
		:SOUTH-NORTH_PATCH_END_STAG = 105 ;
		:GRIDTYPE = "C" ;模拟域的网格信息网格类型，使用Arakawa C类型的水平网格,这是WRF模型常用的网格存储方式
		:DX = 27000.f ;模拟域的网格信息网格间距:经
		:DY = 27000.f ;模拟域的网格信息网格间距:纬
		:DYN_OPT = 2 ;正方体等面积(Square-Square)投影,这是一种保面积的lambertian等面积投影,常用于中尺度模拟
		:CEN_LAT = 32.5f ;模拟域中心位置经纬度
		:CEN_LON = 118.6f ;模拟域中心位置经纬度
		:TRUELAT1 = 30.f ;真实纬度
		:TRUELAT2 = 60.f ;真实纬度
		:MOAD_CEN_LAT = 32.5f ;模式输出在大气数据(MOAD)中使用的中心纬度
		:STAND_LON = 118.f ;标准经线
		:POLE_LAT = 90.f ;地理坐标系北极点的位置,这里设置为经纬度(90N, 0E),即标准的地理坐标系
		:POLE_LON = 0.f ;
		:corner_lats = 19.46492f, 44.11765f, 43.89582f, 19.3251f, 19.44859f, 44.09171f, 43.86773f, 19.30737f, 19.35179f, 44.24067f, 44.01835f, 19.21228f, 19.33549f, 44.21468f, 43.9902f, 19.1946f ;模拟域角点坐标：纬度
		:corner_lons = 106.5773f, 101.4457f, 135.9629f, 130.4113f, 106.4573f, 101.2744f, 136.1329f, 130.5309f, 106.5945f, 101.4097f, 136.0018f, 130.3927f, 106.4747f, 101.2381f, 136.1721f, 130.5121f ;模拟域角点坐标：经度
		:MAP_PROJ = 1 ;投影方式
		:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;地表分类数据集
		:NUM_LAND_CAT = 21 ;下垫面种类总数
		:ISWATER = 17 ;水体类别编号
		:ISLAKE = 21 ;湖面类别编号
		:ISICE = 15 ;冰类别编号
		:ISURBAN = 13 ;城市类别编号
		:ISOILWATER = 14 ;"土壤水分"类别的编号为14
		:grid_id = 1 ;表示这是第1个网格的设置
		:parent_id = 1 ;如果是嵌套网格,则指定这个网格来自哪个母域,这里表明来自母域1
		:i_parent_start = 1 ;
		:j_parent_start = 1 ;
		:i_parent_end = 100 ;
		:j_parent_end = 105 ;
		:parent_grid_ratio = 1 ;
		:sr_x = 1 ;x和y方向的缩放比例,这里都为1,即无缩放
		:sr_y = 1 ;
		:FLAG_MF_XY = 1 ;等于1表示使用数据的地图投影坐标而非地理坐标
		:FLAG_LAI12M = 1 ;等于1表示使用月平均的叶面积指数(LAI)数据
		:FLAG_LAKE_DEPTH = 1 ;湖泊深度数据集
}

在大气模式中,unstagger和stagger两种水平网格是用于存储不同变量的两种网格布局方式。

1. Unstagger网格

unstagger网格也被称为A-grid或矩阵网格,是最简单的网格存储方式。标量变量(如温度、压强、密度等)通常存储在unstagger网格的格点上。unstagger网格的格点分布如下:

```
* * * * *
* * * * *  
* * * * *
* * * * *
* * * * *
```

2. Stagger网格 

stagger网格也被称为C-grid,矢量变量(如风场u、v分量)则存储在这种偏移的网格点上。stagger网格在x方向和y方向上相对unstagger网格各偏移了1/2个网格距离,如下所示:

```
      u
* * * * *
  v   v
* * * * *  
  v   v
* * * * *
  v   v  
* * * * *
      u
```

其中u风分量存储在x方向的stagger点,v风分量存储在y方向的stagger点.

采用stagger网格存储矢量场的优点是:
1) 可以很好地表示差分方程,从而获得较高的数值计算精度。
2) Stagger网格可避免奇异耦合错误,提高数值计算的稳定性。

因此,大气模式像WRF通常会同时使用unstagger网格和stagger网格,前者存储标量场,后者存储矢量场,以获得较高的模拟精度和数值稳定性。GEOGRID文件中定义了两种网格在模拟区域内的网格索引范围。