WRF学习之 ch3 WPS（六）Namelist 描述 GEOGRID、index、METGRID选项，geogrid和Metgrid可用的插值选项，静态数据中土地利用类别，WPS输出场
参考自 https://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.1/users_guide_chap3.html
- Parallelism in the WPS ............................................................. 3-28
- Checking WPS Output .............................................................. 3-29
- WPS Utility Programs ............................................................... 3-30
- Writing Meteorological Data to the Intermediate Format ........... 3-34
- Required Meteorological Fields for Running WRF.....................3-36
- Using MPAS Output for WRF Initial and Boundary Conditions..3-37
- Creating and Editing Vtables .................................................... 3-39
- Writing Static Data to the Geogrid Binary Format ..................... 3-41
- Creating an Urban Fraction Field from NLCD Data .................. 3-44
- Description of Namelist Variables ............................................. 3-46
- Description of GEOGRID.TBL Options ..................................... 3-52
- Description of index Options ..................................................... 3-55
- Description of METGRID.TBL Options ..................................... 3-58
- Available Interpolation Options in Geogrid and Metgrid ............ 3-61
- Land Use and Soil Categories in the Static Data ...................... 3-64
- WPS Output Fields ................................................................... 3-66
Namelist变量的描述
A. SHARE 部分
本节描述了多个WPS程序使用的变量。例如，wrf_core变量指定WPS使用ARW还是NMM内核生成数据，geogrid和metgrid程序都需要该信息。
-
WRF_CORE : 设置为'ARW'或'NMM'的字符串，该字符串告诉WPS正在为输入数据准备哪个动力内核。默认值为'ARW'。
-
MAX_DOM : 一个整数，指定模拟区域/嵌套网格的总数，包括母网格。默认值为1。
-
START_YEAR : 模拟的起始年份。
-
START_MONTH : 模拟的起始月份。
-
START_DAY : 模拟的起始日。
-
START_HOUR : 模拟的起始小时。
-
END_YEAR : 模拟的终止年份。
-
END_MONTH : 模拟的终止月份。
-
END_DAY : 模拟的终止日期。
-
END_HOUR : 模拟的终止时间。
-
START_DATE : 格式为'YYYY-MM-DD_HH:mm:ss'的MAX_DOM字符串列表，指定每个嵌套网格的模拟开始UTC日期。start_date变量是替代使用start_year，start_month，start_day和start_hour的方法。并且，如果两种方法都用于指定开始时间，则start_date变量将优先。无默认值。
-
END_DATE : 格式为'YYYY-MM-DD_HH:mm:ss'的MAX_DOM字符串列表，指定每个嵌套的模拟结束UTC日期。end_date变量是替代使用end_year，end_month，end_day和end_hour的方法。并且，如果两种方法都用于指定结束时间，则end_date变量将优先。无默认值。
-
INTERVAL_SECONDS : 随时间变化的气象输入文件之间的整数秒数。无默认值。
-
ACTIVE_GRID : MAX_DOM逻辑值的列表，为每个网格指定是否应通过geogrid和metgrid处理该网格。默认值为.TRUE.。
-
IO_FORM_GEOGRID : 将写入由geogrid程序创建的区域文件的WRF I/O API格式。可能的选项是：1表示二进制； 2表示NetCDF；3表示GRIB1。给定选项1时，区域文件的后缀为.int；给定选项2时，区域文件的后缀为.nc；给定选项3时，区域文件的后缀为.gr1。默认值为2（NetCDF）。
-
OPT_OUTPUT_FROM_GEOGRID_PATH : 提供绝对路径或者相对路径的字符串，该位置应指向或读取来自geogrid的输出文件。默认值为'./'。
-
DEBUG_LEVEL : 一个整数值，指示应将不同类型的消息发送到标准输出的程度。当debug_level设置为0时，仅一般有用的消息和警告消息将被写入标准输出。当debug_level大于100时，提供更多运行时详细信息的参考消息也将写入标准输出。调试消息和专门用于日志文件的消息永远不会写入标准输出，而总是写入日志文件。预设值为0。
B. GEOGRID 部分
该节指定了用于geogrid程序的变量。geogrid的变量主要定义所有模拟区域域的大小和位置，以及静态地理数据的位置。
-
PARENT_ID : 一个由MAX_DOM个整数组成的列表，对于每个嵌套，网格号为嵌套网格的母网格的编号；对于最外层网格，该数目设置为1。默认值是1。
-
PARENT_GRID_RATIO : 一个由MAX_DOM个整数组成的列表，对于每个嵌套，嵌套比为嵌套网格相对于母网格的比例。无默认值
-
I_PARENT_START : 一个由MAX_DOM个整数组成的列表，对于每个嵌套网格，是在其的母网格的左下角的x坐标的值。对于最粗的网格，这个值指定为1。无默认值。
-
J_PARENT_START : 一个由MAX_DOM个整数组成的列表, 对于每个嵌套网格，是在其的母网格的左下角的y坐标的值。对于最粗的网格，这个值指定为1。无默认值。
-
S_WE : 一个由MAX_DOM个整数组成的列表，其所有值应当都是1，默认的值是1。
-
E_WE : 一个由MAX_DOM个整数组成的列表，对于每个嵌套网格，该值是在这个嵌套网格在东-西方向的大小。对于嵌套网格，e_we必须为该层嵌套的parent_grid_ratio倍加1（即，对于正整数n, e_we = n*parent_grid_ratio+1）。无默认值
-
S_SN : 一个由MAX_DOM个整数组成的列表，其所有值应当都是1，默认的值是1。
-
E_SN : 一个由MAX_DOM个整数组成的列表，对于每个嵌套网格，该值是在这个嵌套网格在东-西方向的大小。对于嵌套网格，e_sn必须为该层嵌套的parent_grid_ratio倍加1（即，对于正整数n, e_sn = n*parent_grid_ratio+1）。无默认值
-
GEOG_DATA_RES : 一个由MAX_DOM个字符串组成的列表，对于每一个嵌套，一个或一组相应的分辨率，用“+”号隔开了用于内插到嵌套的网格中的静力地形数据。对于每个嵌套，该字符串应该包含一个分辨率。这个分辨率与GEOGRID.TBL文件中的对于每个变量的rel_path或abs_path的说明（参见GEOGRID.TBL选项的说明）中的冒号前面的字符串匹配.
如果字符串中的分辨率不匹配GEOGRID.TBL文件中的rel_path或abs_path指定的中的任何此类字符串，那么对于该变量场将使用默认的分辨率的数据，否则就用匹配的。如果有多个分辨率的对应，那么就会使用在GEOGRID.TBL文件中的rel_path或abs_path指定的第一个匹配分辨率。默认值是'default'。 -
DX : x方向上的格距（当地图放大因子是1时候的）。对于ARW，当使用'polar','lambert'，和'mercator'投影时，格距单位是米，当使用'lat-lon'投影是单位是1经度（degrees longitude）。对于NMM，格距的单位是1经度（degrees longitude）。嵌套网格的格距是由给定的parent_grid_ratio和parent_id递归地确定的。无默认值。
-
DY : y方向上的格距（当地图放大因子是1时候的）。对于ARW，当使用'polar','lambert'，和'mercator'投影时，格距单位是米，当使用'lat-lon'投影是单位是1纬度（degrees latitude）。对于NMM，格距的单位是1纬度（degrees latitude）。嵌套网格的格距是由给定的parent_grid_ratio和parent_id递归地确定的。无默认值。
-
MAP_PROJ : 模拟的投影。对于ARW, 可使用的投影有 'lambert', 'polar', 'mercator',和'lat-lon'；对于NMM，投影必须指定为'rotated_ll'。默认的值为'lambert'。
-
REF_LAT : 一个实数值，指定（纬度，经度）位置的纬度部分，该位置在模拟网格中的（i，j）位置已知。对于ARW, 默认情况下（也就是说，当ref_x和ref_y没有指定时)ref_lat给出粗网格中心点的纬度。对于NMM, ref_lat总是给出原点旋转的纬度。没有默认值。
-
REF_LON : 一个实数值，指定（纬度，经度）位置的经度部分，该位置在模拟网格中的（i，j）位置已知。对于ARW, 默认情况下（也就是说，当ref_x和ref_y没有指定时)ref_lat给出粗网格中心点的经度。对于NMM, ref_lon总是给出原点旋转的纬度。对于ARW和NMM，西侧的经度是负的，ref_lon的值应当在[-180,180]范围内。没有默认值。
-
REF_X : 一个实数值，指定(i,j)位置的i，其中模拟区域中的(纬度，经度)位置已知。(i,j)位置通常是针对质量交错网格给出的，其维数通常比非交错网格的维数小1。默认值是(((E_WE-1.)+1.)/2.) = (E_WE/2.)。
-
REF_Y : 一个实数值，指定(i,j)位置的j，其中模拟区域中的(纬度，经度)位置已知。(i,j)位置通常是针对质量交错网格给出的，其维数通常比非交错网格的维数小1。默认值是(((E_SN-1.)+1.)/2.) = (E_SN/2.)。
-
TRUELAT1 : 给定的实数值，对于ARW，是Lambert正形投影中的第一个真纬度，或者是Mercator和极射赤面投影中的唯一的真纬度。对于NMM，不考虑truelat1。无默认值。
-
TRUELAT2 : 给定的实数值，对于ARW，是Lambert兰伯特等角圆锥投影中的第二个真纬度。对于其他所有的投影，都不考虑truelat2。无默认值。
-
STAND_LON : 给定的实数值，对于ARW，该经度是平行于Lambert正形投影和极射赤面投影的y轴。对于常规经纬度投影，该值给出了关于地球地理极点的旋转。对于NMM，不考虑stand_lon。无默认值。
-
POLE_LAT : 对于ARW中的经纬度投影，是在计算经纬度网格的北极点的纬度，其中（指的对于计算网格而言，译者注），-90.0°纬度是全球区域的底部，90.0°纬度是顶部，180.0°经度是在中心。默认值是90.0。
-
POLE_LON : 对于ARW中的经纬度投影，是在计算经纬度网格的北极点的经度，其中，-90.0°纬度是全球区域的底部，90.0°纬度是顶部，180.0°经度是在中心。默认值是0.0。
-
GEOG_DATA_PATH : 一个字符串，给定了地形数据的相对或绝对路径。该路径与指定了GEOGRID.TBL文件中关于的rel_path的位置。没有默认值。
-
OPT_GEOGRID_TBL_PATH : 一个字符串，给出GEOGRID.TBL文件的相对或绝对路径。路径不应包含实际文件名，因为假定GEOGRID.TBL，而应仅提供此文件所在的路径。默认值为 './geogrid/'
C. UNGRIB 部分
目前，这个部分只包含两个变量，它们决定ungrib的输出格式和输出文件的名称。
-
OUT_FORMAT : 一个字符串，其可用的值有'WPS', 'SI', 或者 'MM5'。如果设置为“MM5”，ungrib将以MM5 pregrid程序的格式输出；如果设置为“SI”，ungrib将以grib_prep.exe的格式输出；如果设置为“WPS”，ungrib将以WPS中间格式输出数据。默认值为“WPS”。
-
PREFIX : 将用作ungrib创建的中间格式文件的前缀的字符串；此处，prefix是指中间文件的文件名 PREFIX:YYYY-MM-DD_HH 中的字符串PREFIX。前缀可能包含相对或绝对路径信息，在这种情况下，中间文件将写入指定的目录中。如果要在多个GRIB数据源上运行 ungrib，此选项可能有助于避免重命名中间文件。默认值为'FILE'。
-
ADD_LVLS : 一个逻辑值，决定ungrib是否尝试垂直插值到一组额外的垂直层次（该层次用namelis选项NEW_PLVL和INTERP_TYPE指定）。默认值为.FALSE.。
-
INTERP_TYPE : 一个整数值，指定 ungrib 在垂直插值到新层次时将使用的方法。值为 0 会导致 ungrib 在压力中线性插值，值为 1 会导致 ungrib 在对数压力中线性插值。默认值为 0。
-
NEW_PLVL : 一组实数，指定额外的垂直层次，以 Pa 为单位，当 ADD_LVLS 为真时，ungrib程序将尝试对其进行插值。可以显式地指定新层次的集合，或者，如果层次在气压层上均匀分布，则可以指定三个值：起始气压、终止气压和气压增量。当指定了起始气压、终止气压和气压增量时，气压增量必须是一个负数，用以告诉ungrib程序，说明该值不是目标的压力层，而是要在第一个和第二个之间使用的增量值。没有默认值。
-
PMIN : 一个实数值，指定要从GRIB数据处理的最小压力层，以 Pa 为单位。此选项仅适用于气压层数据集。默认值为 100。
D. METGRID 部分
本节定义了metgrid程序使用的变量。大多数情况下，只是修改fg_name变量，而本节的其他变量较少改动。
-
FG_NAME : 字符串列表，用于指定未解压数据文件的路径和前缀。路径可以是相对的或绝对的，前缀应该包含文件名的所有字符，直到但不包括日期前的冒号。当指定了多个 fg_name 并且在两个或多个输入源中发现相同字段时，最后遇到的源中的数据将优先于该字段的所有先前源。默认值是一个空列表（即没有气象场）。
-
CONSTANTS_NAME : 一个字符串列表给定了时间不变场的未解压数据文件的路径和全部名称。路径可以是相对或者觉得的，文件名应当是完整的文件名； 由于数据是时间不变的，不需要在给定的文件名中加上日期。默认的值是空列表（即，无常数变量场）。
-
IO_FORM_METGRID : 将写入由metgrid 程序创建的输出的WRF I/O API 格式。可能的选项有： 1 表示二进制； 2 为 NetCDF； 3 用于 GRIB1。当给出选项 1 时，输出文件的后缀为 .int；当给出选项 2 时，输出文件的后缀为 .nc；当给出选项 3 时，输出文件的后缀为 .gr1。默认值为 2 (NetCDF)。
-
OPT_OUTPUT_FROM_METGRID_PATH : 一个字符串，给出了应该写入来自metgrid 的输出文件的位置的相对或绝对路径。默认值为当前工作目录（即默认值为'./'）。
-
OPT_METGRID_TBL_PATH : 一个字符串，给出 METGRID.TBL 文件的相对或绝对路径；该路径不应包含实际文件名，因为假定为 METGRID.TBL，而应仅提供此文件所在的路径。默认值为 './metgrid/'。
-
PROCESS_ONLY_BDY: 一个整数，指定初始时间之后的时间段内由metgrid程序处理的边界行和列的数目；对于初始时间，metgrid将始终插入到每个网格点。在 WRF namelist.input 中将此选项设置为 spec_bdy_width 的预期值将加快在metgrid 中的处理速度，但如果网格内部需要插值数据，则不应设置它。如果此选项设置为零，则metgrid 会将气象数据水平插入到模型域中的每个网格点。此选项仅适用于 ARW。默认值为 0。
GEOGRID.TBL选项说明
GEOGRID.TBL文件是一个文本文件，定义了要由geogrid插值的每个数据集的参数。每个数据集都在一个单独的部分中定义，各部分之间用一行等号符号（例如，“ ==============”）分开。在每个部分中都有清单，每个清单都具有 keyword=value 的形式。每个数据集部分中某些关键字是必需的，而其他关键字是可选的。一些关键字与其他关键字互斥。以下，描述了可能的关键字及其期望的值范围。
-
NAME : 指定输出时将分配给插值字段的名称的字符串。无默认值。
-
PRIORITY : 一个整数，指定在table中指定的数据源相对于同一变量场的其他数据源的优先级。如果一个变量场具有n个数据源，则该字段必须有n个单独的表条目，并且必须为每个表条目赋予唯一的priority值，范围为[1，n]。无默认值。
-
DEST_TYPE : 一个字符串，可以是categorical或continuous的，它指示将表部分中给出的数据源中的内插场视为连续场还是分类场。无默认值。
-
INTERP_OPTION : 一个或多个字符串的序列，它们指定了在对变量进行水平内插时要使用的内插方法的名称。可用的插值方法是：average_4pt, average_16pt, wt_average_4pt, wt_average_16pt, nearest_neighbor, four_pt, sixteen_pt, search(r), 和 average_gcell(r)。对于搜索方法（search），可选参数r以源数据网格中的网格点为单位指定最大搜索半径；默认搜索半径为1200点。对于网格单元平均方法（average_gcell），可选参数r指定了应用该方法需要满足的，源数据分辨率与模拟网格分辨率的最小比率；除非指定，否则r = 0.0，这表示选项可用于任何比率。当给出两个或多个方法的序列时，方法应用+号分隔。无默认值。
-
SMOOTH_OPTION : 字符串，给出要在插值后应用于变量场的平滑方法的名称。可用的平滑选项为：1-2-1，smth-desmth和smth-desmth_special（仅适用于ARW）。默认值为null（即不应用任何平滑处理）
-
SMOOTH_PASSES : 如果要对插值变量场执行平滑处理，则smooth_passes指定了应用在变量场的平滑的passes。预设值为1。
-
REL_PATH : 一个字符串，用于指定相对于namelist变量geog_data_path中给定路径的。规范的格式为RES_STRING:REL_PATH，其中RES_STRING是一个字符串，用于以某种独特的方式标识数据的来源或分辨率，并且可以在namelist变量geog_data_res中指定。而REL_PATH是相对于geog_data_path的路径，其中索引并找到数据源的数据图块。如果有多个源或数据源的分辨率，则在table部分中可以给出一个以上的rel_path规定，就像可以为geog_data_res指定多个分辨率（以+符号分隔的顺序）一样。另请参见abs_path。无默认值。
-
ABS_PATH : 一个字符串，它指定数据源的索引和数据切片的绝对路径。规范的格式为RES_STRING:ABS_PATH，其中RES_STRING是一个以某种独特的方式标识数据源或分辨率的字符串，可以在名称列表变量geog_data_res中指定，而ABS_PATH是数据源文件的绝对路径。。如果有多个源或多个分辨率的数据，则在table部分中可以给出一个以上的abs_path规范，就像可以为geog_data_res指定多个分辨率（以+符号分隔的顺序）一样。另请参见rel_path。无默认值。
-
OUTPUT_STAGGER : 一个字符串，它指定要插入字段的交错网格。对于ARW域，可能的值为U，V和M。对于NMM域，可能的值为HH和VV。ARW的默认值为M；NMM的默认值为HH。
-
LANDMASK_WATER : 一个或多个逗号分隔的整数值，给出代表水面类别的场的索引。当在dest_type=categorical field的表部分指定了landmask_water时，L将使用指定类别作为水类别进行计算LANDMASK field。关键词landmask_water和landmask_land是互斥的。默认值为空（即，不会从变量场计算landmask）。
-
LANDMASK_LAND : 一个或多个逗号分隔的整数值，给出代表土地的字段中类别的索引。当在dest_type = categorical field的表部分中指定了landmask_water时，将使用指定类别作为土地类别从该字段计算LANDMASK field。关键词landmask_water和landmask_land是互斥的。默认值为空（即，不会从变量场计算landmask）。
-
MASKED : land 或water，分别表示在土地或水域该变量场无效。如果将masked关键字用于变量场，则将为那些masked类型（land 或water）的网格点分配由fill_missing指定的值。默认值为null（即未屏蔽该变量场）。
-
FILL_MISSING : 用于填充内插变量场中所有丢失或mask的网格点的实数值。默认值为1.E20。
-
HALT_ON_MISSING : yes或no，指示在插值变量场中遇到缺失值时，geogrid是否以错误消息停止。默认值为 no。
-
DOMINANT_CATEGORY : 当指定为字符串时，其作用是使geogrid从分散的categorical场中计算出 dominant category，并输出由dominant_category值指定的dominant catrgory 场。此选项只能用于dest_type = categorical的字段。默认值为null（即，不会从fractional categorical field中计算出主要类别）。
-
DOMINANT_ONLY : 当指定为字符串时，其效果类似于dominant_category关键字的效果：geogrid将根据fractional categorical field计算dominant category，并输出具有由dominant_only值指定的名称的优势类别字段。但是，与dominant_category不同，当使用dominant_only时，分数分类字段将不会出现在土工格栅输出中。此选项只能用于dest_type = categorical的字段。默认值为null（即，不会从分数分类字段中计算出主要类别）。
-
DF_DX : 当给df_dx设置一个字符串值时，geogrid程序会计算气象场的沿x方向导数。在网格的内部，使用中央差分方法，在网格的边界，使用单侧差分；导数场将根据关键字df_dx的字符串值来命名。默认值是null（即，不计算导数场）。
-
DF_DY : 当给df_dy设置一个字符串值时，geogrid程序会计算气象场的沿y方向导数。在网格的内部，使用中央差分方法，在网格的边界，使用单侧差分；导数场将根据关键字df_dy的字符串值来命名。默认值是null。
-
Z_DIM_NAME : 对于输出的3维场，使用该字符串给出垂直纬度（又称z维度）的名称。一个continuous field可以具有多个层次，因此是3维场；而一个categorical field，如果将每个category写成fractional fields，那么可以认为是3维场。无默认值。
-
FLAG_IN_OUTPUT : 提供全局属性名称的字符串，该全局属性将被赋值为1并写入geogrid输出。默认值为空（即，不会为字段写入任何标志）。
-
OPTIONAL : yes或no，说明由geog_data_res选项中指定的分辨率的数据是否是可选的。如果GEOGRID.TBL文件中的条目是可选的，且指定分辨率的数据无法被读取，geogrid将会打印一个消息，说明数据集没有被插值，随后继续运行；否则如果这个条目不是可选的，且指定分辨率的数据无法被读取，geogrid会报错并停止。可以对于相同的变量指定不同的优先级条目的optional关键字，例如，对于一般变量场priority=2可能是optionnal的，而priority=1可能是非-optional的（即，optional=no)。默认值是no
index 选项
与GEOGRID.TBL相关的是index files(索引文件)，他们是与每个静态数据集关联的。一个index文件定义了针对该数据集的参数，而GEOGRID.TBL文件描述了geogrid应如何处理每个数据集。就像在GEOGRID.TBL文件中，索引文件中的规范的格式为keyword=value。下面是可能的关键字及其可能的值。
-
PROJECTION : 指定数据投影方式的字符串，可以是lambert、polar、mercator、regular_ll、albers_nad83或polar_wgs84。没有默认值。
-
TYPE : 一个字符串，其值要么是分类的要么是连续的，用于确定数据文件中的数据应解释为连续场还是离散场。对于由每个可能类别的分数字段表示的类别数据，type值应设置为continuous。没有默认值。
-
SIGNED : yes或no，指示数据文件中的值（始终表示为整数）是否以2的补码形式签名。默认值为"no"。
-
UNITS : 用引号(")括起来的字符串，说明插值变量场的单位；字符串将作为与时间无关的可变属性写入geogrid输出文件。没有默认值。
-
DESCRIPTION : 一个字符串，用引号(")括起来，给出插值变量场的简短描述；该字符串将作为与时间无关的变量写入geogrid输出文件。没有默认值。
-
DX : 一个实数值，给出数据集在x方向上的网格间距。如果投影是lambert, polar, mercator, albers_nad83, 或者polar_wg84其中之一，那么dx就是网格的格距，单位是米；如果投影是regular_ll，那么dx是网格格距，单位是度。没有默认值。
-
DY : 一个实数值，给出数据集在y方向上的网格间距。如果投影是lambert, polar, mercator, albers_nad83, 或者polar_wg84其中之一，那么dy就是网格的格距，单位是米；如果投影是regular_ll，那么dy是网格格距，单位是度。没有默认值。
-
KNOWN_X : 一个实数值，指定与投影中已知的(latitude, longitude)位置相对应的(i,j)位置的i坐标。默认值为1。
-
KNOWN_Y : 一个实数值，指定与投影中已知的(latitude, longitude)位置相对应的(i,j)位置的j坐标。默认值为1。
-
KNOWN_LAT : 一个实数值，指定投影中已知的(latitude, longitude)位置的纬度。没有默认值。
-
KNOWN_LON : 一个实数值，指定投影中已知的(latitude, longitude)位置的经度。没有默认值。
-
STDLON : 一个实数值，指定了在圆锥(conic)和方位(azimuthal)投影中指定与y轴平行的经度。没有默认值。
-
TRUELAT1 : 一个实数值，指定圆锥投影的第一个真实纬度或方位投影的唯一真实纬度。没有默认值。
-
TRUELAT2 : 一个实数值，指定圆锥投影的第二个真实纬度。没有默认值。
-
WORDSIZE : 一个整数，表示数据文件中每个网格点的值的所占的字节数。没有默认值。
-
TILE_X : 一个整数，用于指定源数据的单个碎片在x方向上的格点数（不包括任何外围点(halo points)）。没有默认值。
-
TILE_Y : 一个整数，用于指定源数据的单个碎片在y方向上的格点数（不包括任何外围点(halo points)）。没有默认值。
-
TILE_Z : 一个整数，用于指定源数据的单个碎片在z方向上的格点数；该关键字可作为关键字tile_z_start和tile_z_end的替代品，使用该关键字时，起始z索引假定为1。没有默认值。
-
TILE_Z_START : 一个整数，指定数据文件中数组在z方向的起始索引。如果使用此关键字，则还必须指定tile_z_end。没有默认值。
-
TILE_Z_END : 一个整数，指定数据文件中数组在z方向的结束索引。如果使用此关键字，则还必须指定tile_z_start。无默认值
-
CATEGORY_MIN : 一个整数，对于类别数据（type=categorical），指定在数据集中找到的最小的分类索引。如果使用此关键字，还必须指定category_max 。没有默认值。
-
CATEGORY_MAX : 一个整数，对于类别数据（type=categorical），指定在数据集中找到的最大分类索引。如果使用此关键字，还必须指定类别category_min。没有默认值。
-
TILE_BDR : 一个整数，指定每个数据碎片的在网格点中的外围区宽度(halo width)。默认值为0。
-
MISSING_VALUE : 一个实数，设置了数据集中缺测值。没有默认值。
-
SCALE_FACTOR : 一个实数值，在以整数形式读入数据集碎片中的数据后，数据（通过乘法）缩放的比例。默认值为1。
-
ROW_ORDER : 一个字符串，要么是bottom_top，要么是top_bottom，指定数据集数组的行是从最低索引行到最高索引行（bottom_top）还是从最高索引行到最低索引行（top_bottom）。当使用一些USGS数据集时，这个关键字可能很有用，这些数据集是按top_bottom顺序提供的。默认值为bottom_top。
-
ENDIAN : 字符串，要么是big，要么是little，指定静态数据集数组中的值是按大端字节顺序还是按小端字节顺序。默认值是big。
-
ISWATER : 整数，指定土地利用类别。默认值为16。
-
ISLAKE : 整数，指定了内陆水体的土地利用类别。默认值是-1（即，没有单独的内陆水体类别）。
-
ISICE : 整数，指定了冰的土地利用类别。默认值为24。
-
ISURBAN : 整数，指定了城市用地类别。默认值为1。
-
ISOILWATER : 整数，指定了水的土壤类别。默认值为14。
-
MMINLU : 字符串，用引号(")括起来，告诉WRF的LANDUSE.TBL以及VEGPARM.TBL，使用哪种土地使用分类。默认值为“USGS”。
-
FILENAME_DIGITS : 整数，说明数据碎片名称中使用的数字的位数。可能的值为5或6。默认值为5。
METGRID.TBL选项
METGRID.TBL文件是一个文本文件，它定义了要由metgrid插值的每个气象场的参数。每气象场的参数在一个单独的部分中定义，各部分直接由一行等号分隔（例如，'============='）。在每个部分中，都有定义(specification)，每个规范的形式都是keyword=value。一些关键字是必需的，而另一些关键字是可选的；某些关键字与其他关键字互斥。下面，将描述可能的关键字及其期望的值范围。
-
NAME : 一个字符串，给出对应是部分所属的气象场的名称。名称应该与中间文件中给出的气象场的名称完全匹配（因此，是在生成中间文件时使用的Vtable中给出的名称）。此字段必填。没有默认值。
-
OUTPUT : yes或no，说明是否写到metgrid输出文件中。默认值是yes。
-
MANDATORY : yes或no，说明是否需要该气象场才能成功完成metgrid。默认值是no。
-
OUTPUT_NAME : 一个字符串，给出插值变量输出的名称。当设置output_name值时，将使用table中与具有指定名称的变量相关的插值选项。因此，指定output_name的效果有两个方面：在写入之前为插值变量分配指定的名称，并且使用与变量相关的插值方法，该变量的名称与指定给output_name关键字的值匹配。没有默认值。
-
FROM_INPUT : 用于与fg_name namelist变量中的值进行比较的字符串；如果指定了from_input，则仅当输入的时间变化变量文件的文件名包含from_input的值作为子字符串时，才会使用containing table节。因此，根据正在处理的场的哪个源，可以使用from_input为同一场指定不同的插值选项。没有默认值。
-
OUTPUT_STAGGER : 需要被插值的变量的网格交错。对于ARW，这必须是U、V和M中的一个；对于NMM，这必须是HH和VV中的一个。ARW的默认值为M；NMM的默认值为HH。
-
IS_U_FIELD : yes或no，表示该场是否被用作风U分量场。对于ARW，风U分量场必须内插到U staggering（output_stagger=U）；对于NMM，风U分量场必须内插到V staggering（output_stagger=VV）。默认值为no。
-
IS_V_FIELD : yes 或 no, 指示是否将该场用作风V分量场。对于ARW，风V分量必须内插到V staggering（output_stager=V）；对于NMM，风V分量必须内插到Vstaggering（output_stagger=VV）。默认值为no。
-
INTERP_OPTION : 一种由一个或多个字符串组成的序列，这些字符串说明了需要对变量进行水平内插时使用的方法名称。可用的插值方法有：average_4pt、average_16pt、wt_average_4pt、wt_average_16pt、nearest_neighbor、four_pt、thiven_pt、search(r)和average_gcell(r)。对于搜素方法（search），可选参数r指定了最大搜索半径（以源数据网格中的网格点为单位）；默认搜索半径为1200格点。对于网格细胞平均方法（average_gcell），可选参数r指定将应用该方法的源数据分辨率与模拟网格分辨率的最小比率；除非指定，否则r=0.0，并且该选项用于任何比率。当给定两个或多个方法的序列时，这些方法应以+号分隔。默认值为nearest_neighbor。
-
INTERP_MASK : 要使用插值mask的场的名称，以及该场中表示屏蔽点的值和可选的关系符号，< 或者 > 。规范采用表单域（？）？maskval），其中field是字段的名称？是可选的关系符号（
），maskval是实值。如果字段中的对应点分别等于、大于或小于无关系符号、>符号或<symbol的maskval值，则不会在插值中使用源数据点。默认值为“无遮罩”。
The name of the field to be used as an interpolation mask, along with the value within that field which signals masked points and an optional relational symbol, < or >. A specification takes the form field(?maskval), where field is the name of the field, ? is an optional relational symbol (< or >), and maskval is a real value. Source data points will not be used in interpolation if the corresponding point in the field field is equal, greater than, or less than, the value of maskval for no relational symbol, a > symbol, or a < symbol, respectively. Default value is no mask. -
INTERP_LAND_MASK : 插值到水点时用作插值掩码的字段的名称（由静态LANDMASK字段确定），以及该字段中表示陆点的值和可选的关系符号
。规范采用表单域（？）？maskval），其中field是字段的名称？是可选的关系符号（ ），maskval是实值。默认值为“无遮罩”。he name of the field to be used as an interpolation mask when interpolating to water points (determined by the static LANDMASK field), along with the value within that field which signals land points and an optional relational symbol, < or >. A specification takes the form field(?maskval), where field is the name of the field, ? is an optional relational symbol (< or >), and maskval is a real value. Default value is no mask. -
INTERP_WATER_MASK : 插值到岸点（由静态LANDMASK字段确定）时用作插值掩码的字段的名称，以及该字段中表示水点的值和可选关系符号
。规范采用表单域（？）？maskval），其中field是字段的名称？是可选的关系符号（ ），maskval是实值。默认值为“无遮罩”。The name of the field to be used as an interpolation mask when interpolating to land points (determined by the static LANDMASK field), along with the value within that field which signals water points and an optional relational symbol, < or >. A specification takes the form field(?maskval), where field is the name of the field, ? is an optional relational symbol (< or >), and maskval is a real value. Default value is no mask. -
FILL_MISSING : 一个实数，指定分配给模型网格点的值，例如，由于缺少或不完整的气象数据而没有得到插值。默认值为1.E20。A real number specifying the value to be assigned to model grid points that received no interpolated value, for example, because of missing or incomplete meteorological data. Default value is 1.E20.
-
Z_DIM_NAME : 对于三维气象场，一种字符串，给出输出字段要使用的垂直维度的名称。默认值为num\u metgrid\u levels。For 3-dimensional meteorological fields, a character string giving the name of the vertical dimension to be used for the field on output. Default value is num_metgrid_levels.
-
DERIVED : 是或否，指示是否从其他插值字段而不是从输入字段中插值字段导出字段。默认值为“否”。Either yes or no, indicating whether the field is to be derived from other interpolated fields, rather than interpolated from an input field. Default value is no.
-
FILL_LEV : fill\u lev关键字可以在一个表节中多次指定，它指定如果字段的某个级别不存在，则该级别应该如何填充。关键字的泛型值采用以下形式DLEVEL:字段（SLEVEL），其中DLEVEL指定要填充的字段中的级别，field指定要从中复制级别的源字段，SLEVEL指定要使用的源字段中的级别。DLEVEL可以是整数或字符串all。字段可以是另一个字段的名称、字符串常量或字符串索引。如果将字段指定为const，则SLEVEL是用于填充的常量值；如果将字段指定为垂直索引，则不能指定（SLEVEL），并使用源字段的垂直索引值；如果DLEVEL为'all'，然后使用level\u template关键字指定的字段中的所有级别来填充字段中相应的级别，一次填充一个级别。没有默认值。The fill_lev keyword, which may be specified multiple times within a table section, specifies how a level of the field should be filled if that level does not already exist. A generic value for the keyword takes the form DLEVEL:FIELD(SLEVEL), where DLEVEL specifies the level in the field to be filled, FIELD specifies the source field from which to copy levels, and SLEVEL specifies the level within the source field to use. DLEVEL may either be an integer or the string all. FIELD may either be the name of another field, the string const, or the string vertical_index. If FIELD is specified as const, then SLEVEL is a constant value that will be used to fill with; if FIELD is specified as vertical_index, then (SLEVEL) must not be specified, and the value of the vertical index of the source field is used; if DLEVEL is 'all', then all levels from the field specified by the level_template keyword are used to fill the corresponding levels in the field, one at a time. No default value.
-
LEVEL_TEMPLATE :一种字符串，给出一个字段的名称，从该字段中可以获得一系列垂直级别并用作模板。此关键字与fillèlev规范一起使用，该规范在其规范的DLEVEL部分使用all。没有默认值。 A character string giving the name of a field from which a list of vertical levels should be obtained and used as a template. This keyword is used in conjunction with a fill_lev specification that uses all in the DLEVEL part of its specification. No default value.
-
MASKED : Either land, water, or both. 要么是陆地，要么是水，要么两者兼而有之。将“MASKED”设置为“land”或“water”表示字段不应分别插值到WRF land或water points；但是，将“MASKED”设置为“both”表示字段应仅使用源数据中的land points插值到WRF land points，并仅使用源数据中的water points插值到WRF water points。当某个字段被屏蔽或无效时，静态LANDMASK字段将用于确定该字段应插值到哪些模型网格点；无效点将被分配FILL\u MISSING关键字给定的值。源数据点是陆地还是水由使用INTERP\u land\u MASK和INTERP\u water\u MASK选项指定的掩码确定。默认值为空（即该字段对陆点和水点都有效）。Setting MASKED to land or water indicates that the field should not be interpolated to WRF land or water points, respectively; however, setting MASKED to both indicates that the field should be interpolated to WRF land points using only land points in the source data and to WRF water points using only water points in the source data. When a field is masked, or invalid, the static LANDMASK field will be used to determine which model grid points the field should be interpolated to; invalid points will be assigned the value given by the FILL_MISSING keyword. Whether a source data point is land or water is determined by the masks specified using the INTERP_LAND_MASK and INTERP_WATER_MASK options. Default value is null (i.e., the field is valid for both land and water points).
-
MISSING_VALUE : 实数，缺测值。
-
VERTICAL_INTERP_OPTION : 一个字符串，指定垂直插值到缺失点时应使用的垂直插值方法。目前，此选项尚未实现。没有默认值。A character string specifying the vertical interpolation method that should be used when vertically interpolating to missing points. Currently, this option is not implemented. No default value.
-
FLAG_IN_OUTPUT : 一个字符串，给出一个全局属性的名称，如果要输出插值字段（output=yes），它将被赋值为1并写入metgrid输出。默认值为空（即，不会为字段写入任何标志）。A character string giving the name of a global attribute which will be assigned a value of 1 and written to the metgrid output if the interpolated field is to be output (output=yes). Default value is null (i.e., no flag will be written for the field).
在Geogrid and Metgrid中可用的插值选项
通过GEOGRID.TBL以及METGRID.TBL文件，用户可以控制对源数据进行插值的方法，无论是geogrid中的静态场还是metgrid中的气象场。事实上，可以给出插值方法的列表，在这种情况下，如果不可能使用列表中的第i种方法，那么将会使用第（i+1）中方法，直到找到可以使用的方法，或者在列表中没有可以使用的方法。例如，要对变量场使用四点双线性插值方案，可以指定interp_option=four_pt。但是，如果变量中有缺测值区域存在，这可能会阻止使用four_pt选项，那么我们可以通过指定interp_option=four_pt+average_4pt来请求尝试使用简单的“四点平均”方法。下面，从概念上描述了WPS中的每个可用插值选项；有关每个方法的详细信息，请用户参考文件WPS/geogrid/src/interp_options.F中的源代码。
- four_pt : 四点-双线性插值
如上图所示，四点双线性插值法需要四个有效的点aij，围绕着geogrid或metgrid必须插值的点(x,y)。直观地说，该方法通过线性插值到a11和a12之间以及a21和a22之间的点(x,y)的x坐标，然后使用这两个插值值线性插值到y坐标。
- sixteen_pt : 十六点重叠抛物线插值(Sixteen-point overlapping parabolic interpolation)
The sixteen_pt overlapping parabolic interpolation method requires sixteen valid source points surrounding the point (x,y), as illustrated in the figure above. The method works by fitting one parabola to the points ai1, ai2, and ai3, and another parabola to the points ai2, ai3, and ai4, for row i, ; then, an intermediate interpolated value pi within row i at the x-coordinate of the point is computed by taking an average of the values of the two parabolas evaluated at x, with the average being weighted linearly by the distance of x from ai2 and ai3. Finally, the interpolated value at (x,y) is found by performing the same operations as for a row of points, but for the column of interpolated values pi to the y-coordinate of (x,y).
-
average_4pt : Simple four-point average interpolation
The four-point average interpolation method requires at least one valid source data point from the four source points surrounding the point (x,y). The interpolated value is simply the average value of all valid values among these four points. -
wt_average_4pt : Weighted four-point average interpolation
The weighted four-point average interpolation method can handle missing or masked source data points, and the interpolated value is given as the weighted average of all valid values, with the weight wij for the source point aij, , given by
Here, xi is the x-coordinate of aij and yj is the y-coordinate of aij.
-
average_16pt : Simple sixteen-point average interpolation
The sixteen-point average interpolation method works in an identical way to the four-point average, but considers the sixteen points surrounding the point (x,y). -
wt_average_16pt : Weighted sixteen-point average interpolation
The weighted sixteen-point average interpolation method works like the weighted four-point average, but considers the sixteen points surrounding (x,y); the weights in this method are given by
where xi and yj are as defined for the weighted four-point method, and .
-
nearest_neighbor : Nearest neighbor interpolation
When used for continuous datasets (i.e., datasets that have type=continuous in their index files), the nearest neighbor interpolation method simply sets the interpolated value at (x,y) to the value of the nearest source data point, regardless of whether this nearest source point is valid, missing, or masked. For categorical datasets (i.e., datasets that have type=categorical in their index files), this option actually causes the geogrid program to consider all source pixels that lie within each WRF grid cell, and to find the fraction of the WRF grid cell that is comprised of each category in the source data. -
search : Breadth-first search interpolation
The breadth-first search option works by treating the source data array as a 2-d grid graph, where each source data point, whether valid or not, is represented by a vertex. Then, the value assigned to the point (x,y) is found by beginning a breadth-first search at the vertex corresponding to the nearest neighbor of (x,y), and stopping once a vertex representing a valid (i.e., not masked or missing) source data point is found. In effect, this method can be thought of as "nearest valid neighbor". -
average_gcell : Model grid-cell average
The grid-cell average interpolator may be used when the resolution of the source data is higher than the resolution of the model grid. For a model grid cell Γ, the method takes a simple average of the values of all source data points that are nearer to the center of Γ than to the center of any other grid cell. The operation of the grid-cell average method is illustrated in the figure above, where the interpolated value for the model grid cell – represented as the large rectangle – is given by the simple average of the values of all of the shaded source grid cells.
静态数据中的土地利用和土壤类别
默认土地使用和土壤类别数据集是作为WPS静态文件的一部分默认提供，它们包含了与WRF run 目录下的VEGPARM.TBL和SOILPARM.TBL文件中描述的USGS类别相匹配的类别。下表提供了24种土地利用类型和16种土壤类型的描述。
Table 1: USGS 24-category Land Use Categories
| Land Use Category | Land Use Description |
|---|---|
| 1 | Urban and Built-up Land |
| 2 | Dryland Cropland and Pasture |
| 3 | Irrigated Cropland and Pasture |
| 4 | Mixed Dryland/Irrigated Cropland and Pasture |
| 5 | Cropland/Grassland Mosaic |
| 6 | Cropland/Woodland Mosaic |
| 7 | Grassland |
| 8 | Shrubland |
| 9 | Mixed Shrubland/Grassland |
| 10 | Savanna |
| 11 | Deciduous Broadleaf Forest |
| 12 | Deciduous Needleleaf Forest |
| 13 | Evergreen Broadleaf |
| 14 | Evergreen Needleleaf |
| 15 | Mixed Forest |
| 16 | Water Bodies |
| 17 | Herbaceous Wetland |
| 18 | Wooden Wetland |
| 19 | Barren or Sparsely Vegetated |
| 20 | Herbaceous Tundra |
| 21 | Wooded Tundra |
| 22 | Mixed Tundra |
| 23 | Bare Ground Tundra |
| 24 | Snow or Ice |
Table 2: IGBP-Modified MODIS 20-category Land Use Categories
| Land Use Category | Land Use Description |
|---|---|
| 1 | Evergreen Needleleaf Forest |
| 2 | Evergreen Broadleaf Forest |
| 3 | Deciduous Needleleaf Forest |
| 4 | Deciduous Broadleaf Forest |
| 5 | Mixed Forests |
| 6 | Closed Shrublands |
| 7 | Open Shrublands |
| 8 | Woody Savannas |
| 9 | Savannas |
| 10 | Grasslands |
| 11 | Permanent Wetlands |
| 12 | Croplands |
| 13 | Urban and Built-Up |
| 14 | Cropland/Natural Vegetation Mosaic |
| 15 | Snow and Ice |
| 16 | Barren or Sparsely Vegetated |
| 17 | Water |
| 18 | Wooded Tundra |
| 19 | Mixed Tundra |
| 20 | Barren Tundra |
Table 3: 16-category Soil Categories
| Soil Category | Soil Description |
|---|---|
| 1 | Sand |
| 2 | Loamy Sand |
| 3 | Sandy Loam |
| 4 | Silt Loam |
| 5 | Silt |
| 6 | Loam |
| 7 | Sandy Clay Loam |
| 8 | Silty Clay Loam |
| 9 | Clay Loam |
| 10 | Sandy Clay |
| 11 | Silty Clay |
| 12 | Clay |
| 13 | Organic Material |
| 14 | Water |
| 15 | Bedrock |
| 16 | Other (land-ice) |
WPS 输出场
下面给出了写入geogrid程序输出文件的全局属性和变量列表。这个列表是运行ncdump程序对一个典型的geo_em.d01文件进行的输出的一个简化版本。
ncdump输出
netcdf geo_em.d01 {
dimensions:
Time = UNLIMITED ; // (1 currently)
DateStrLen = 19 ;
west_east = 73 ;
south_north = 60 ;
south_north_stag = 61 ;
west_east_stag = 74 ;
land_cat = 21 ;
soil_cat = 16 ;
month = 12 ;
num_urb_params = 132 ;
variables:
char Times(Time, DateStrLen) ;
float XLAT_M(Time, south_north, west_east) ;
XLAT_M:units = "degrees latitude" ;
XLAT_M:description = "Latitude on mass grid" ;
float XLONG_M(Time, south_north, west_east) ;
XLONG_M:units = "degrees longitude" ;
XLONG_M:description = "Longitude on mass grid" ;
float XLAT_V(Time, south_north_stag, west_east) ;
XLAT_V:units = "degrees latitude" ;
XLAT_V:description = "Latitude on V grid" ;
float XLONG_V(Time, south_north_stag, west_east) ;
XLONG_V:units = "degrees longitude" ;
XLONG_V:description = "Longitude on V grid" ;
float XLAT_U(Time, south_north, west_east_stag) ;
XLAT_U:units = "degrees latitude" ;
XLAT_U:description = "Latitude on U grid" ;
float XLONG_U(Time, south_north, west_east_stag) ;
XLONG_U:units = "degrees longitude" ;
XLONG_U:description = "Longitude on U grid" ;
float CLAT(Time, south_north, west_east) ;
CLAT:units = "degrees latitude" ;
CLAT:description = "Computational latitude on mass grid" ;
float CLONG(Time, south_north, west_east) ;
CLONG:units = "degrees longitude" ;
CLONG:description = "Computational longitude on mass grid" ;
float MAPFAC_M(Time, south_north, west_east) ;
MAPFAC_M:units = "none" ;
MAPFAC_M:description = "Mapfactor on mass grid" ;
float MAPFAC_V(Time, south_north_stag, west_east) ;
MAPFAC_V:units = "none" ;
MAPFAC_V:description = "Mapfactor on V grid" ;
float MAPFAC_U(Time, south_north, west_east_stag) ;
MAPFAC_U:units = "none" ;
MAPFAC_U:description = "Mapfactor on U grid" ;
float MAPFAC_MX(Time, south_north, west_east) ;
MAPFAC_MX:units = "none" ;
MAPFAC_MX:description = "Mapfactor (x-dir) on mass grid" ;
float MAPFAC_VX(Time, south_north_stag, west_east) ;
MAPFAC_VX:units = "none" ;
MAPFAC_VX:description = "Mapfactor (x-dir) on V grid" ;
float MAPFAC_UX(Time, south_north, west_east_stag) ;
MAPFAC_UX:units = "none" ;
MAPFAC_UX:description = "Mapfactor (x-dir) on U grid" ;
float MAPFAC_MY(Time, south_north, west_east) ;
MAPFAC_MY:units = "none" ;
MAPFAC_MY:description = "Mapfactor (y-dir) on mass grid" ;
float MAPFAC_VY(Time, south_north_stag, west_east) ;
MAPFAC_VY:units = "none" ;
MAPFAC_VY:description = "Mapfactor (y-dir) on V grid" ;
float MAPFAC_UY(Time, south_north, west_east_stag) ;
MAPFAC_UY:units = "none" ;
MAPFAC_UY:description = "Mapfactor (y-dir) on U grid" ;
float E(Time, south_north, west_east) ;
E:units = "-" ;
E:description = "Coriolis E parameter" ;
float F(Time, south_north, west_east) ;
F:units = "-" ;
F:description = "Coriolis F parameter" ;
float SINALPHA(Time, south_north, west_east) ;
SINALPHA:units = "none" ;
SINALPHA:description = "Sine of rotation angle" ;
float COSALPHA(Time, south_north, west_east) ;
COSALPHA:units = "none" ;
COSALPHA:description = "Cosine of rotation angle" ;
float LANDMASK(Time, south_north, west_east) ;
LANDMASK:units = "none" ;
LANDMASK:description = "Landmask : 1=land, 0=water" ;
float XLAT_C(Time, south_north_stag, west_east_stag) ;
XLAT_C:units = "degrees latitude" ;
XLAT_C:description = "Latitude at grid cell corners" ;
float XLONG_C(Time, south_north_stag, west_east_stag) ;
XLONG_C:units = "degrees longitude" ;
XLONG_C:description = "Longitude at grid cell corners" ;
float LANDUSEF(Time, land_cat, south_north, west_east) ;
LANDUSEF:units = "category" ;
LANDUSEF:description = "Noah-modified 21-category IGBP-MODIS landuse" ;
float LU_INDEX(Time, south_north, west_east) ;
LU_INDEX:units = "category" ;
LU_INDEX:description = "Dominant category" ;
float HGT_M(Time, south_north, west_east) ;
HGT_M:units = "meters MSL" ;
HGT_M:description = "GMTED2010 30-arc-second topography height" ;
float SOILTEMP(Time, south_north, west_east) ;
SOILTEMP:units = "Kelvin" ;
SOILTEMP:description = "Annual mean deep soil temperature" ;
float SOILCTOP(Time, soil_cat, south_north, west_east) ;
SOILCTOP:units = "category" ;
SOILCTOP:description = "16-category top-layer soil type" ;
float SCT_DOM(Time, south_north, west_east) ;
SCT_DOM:units = "category" ;
SCT_DOM:description = "Dominant category" ;
float SOILCBOT(Time, soil_cat, south_north, west_east) ;
SOILCBOT:units = "category" ;
SOILCBOT:description = "16-category top-layer soil type" ;
float SCB_DOM(Time, south_north, west_east) ;
SCB_DOM:units = "category" ;
SCB_DOM:description = "Dominant category" ;
float ALBEDO12M(Time, month, south_north, west_east) ;
ALBEDO12M:units = "percent" ;
ALBEDO12M:description = "Monthly surface albedo" ;
float GREENFRAC(Time, month, south_north, west_east) ;
GREENFRAC:units = "fraction" ;
GREENFRAC:description = "MODIS FPAR" ;
float LAI12M(Time, month, south_north, west_east) ;
LAI12M:units = "m^2/m^2" ;
LAI12M:description = "MODIS LAI" ;
float SNOALB(Time, south_north, west_east) ;
SNOALB:units = "percent" ;
SNOALB:description = "Maximum snow albedo" ;
float SLOPECAT(Time, south_north, west_east) ;
SLOPECAT:units = "category" ;
SLOPECAT:description = "Dominant category" ;
float CON(Time, south_north, west_east) ;
CON:units = "" ;
CON:description = "Subgrid-scale orographic convexity" ;
float VAR(Time, south_north, west_east) ;
VAR:units = "" ;
VAR:description = "Subgrid-scale orographic variance" ;
float OA1(Time, south_north, west_east) ;
OA1:units = "" ;
OA1:description = "Subgrid-scale orographic asymmetry" ;
float OA2(Time, south_north, west_east) ;
OA2:units = "" ;
OA2:description = "Subgrid-scale orographic asymmetry" ;
float OA3(Time, south_north, west_east) ;
OA3:units = "" ;
OA3:description = "Subgrid-scale orographic asymmetry" ;
float OA4(Time, south_north, west_east) ;
OA4:units = "" ;
OA4:description = "Subgrid-scale orographic asymmetry" ;
float OL1(Time, south_north, west_east) ;
OL1:units = "" ;
OL1:description = "Subgrid-scale effective orographic length scale" ;
float OL2(Time, south_north, west_east) ;
OL2:units = "" ;
OL2:description = "Subgrid-scale effective orographic length scale" ;
float OL3(Time, south_north, west_east) ;
OL3:units = "" ;
OL3:description = "Subgrid-scale effective orographic length scale" ;
float OL4(Time, south_north, west_east) ;
OL4:units = "" ;
OL4:description = "Subgrid-scale effective orographic length scale" ;
float VAR_SSO(Time, south_north, west_east) ;
VAR_SSO:units = "meters2 MSL" ;
VAR_SSO:description = "Variance of Subgrid Scale Orography" ;
float LAKE_DEPTH(Time, south_north, west_east) ;
LAKE_DEPTH:units = "meters MSL" ;
LAKE_DEPTH:description = "Topography height" ;
float URB_PARAM(Time, num_urb_params, south_north, west_east) ;
URB_PARAM:units = "dimensionless" ;
URB_PARAM:description = "Urban_Parameters" ;
// global attributes:
:TITLE = "OUTPUT FROM GEOGRID V4.0" ;
:SIMULATION_START_DATE = "0000-00-00_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 74 ;
:SOUTH-NORTH_GRID_DIMENSION = 61 ;
:BOTTOM-TOP_GRID_DIMENSION = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 73 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 74 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 60 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 61 ;
:GRIDTYPE = "C" ;
:DX = 30000.f ;
:DY = 30000.f ;
:DYN_OPT = 2 ;
:CEN_LAT = 34.83001f ;
:CEN_LON = -81.03f ;
:TRUELAT1 = 30.f ;
:TRUELAT2 = 60.f ;
:MOAD_CEN_LAT = 34.83001f ;
:STAND_LON = -98.f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:corner_lats = 28.17127f, 44.36657f, 39.63231f, 24.61906f, 28.17842f, 44.37617f, 39.57812f, 24.57806f, 28.03771f, 44.50592f, 39.76032f, 24.49431f, 28.04485f, 44.51553f, 39.70599f, 24.45341f ;
:corner_lons = -93.64893f, -92.39661f, -66.00165f, -72.64047f, -93.80048f, -92.59155f, -65.83557f, -72.5033f, -93.65717f, -92.3829f, -65.9313f, -72.68539f, -93.80841f, -92.57831f, -65.76495f, -72.54843f ;
:MAP_PROJ = 1 ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:grid_id = 1 ;
:parent_id = 1 ;
:i_parent_start = 1 ;
:j_parent_start = 1 ;
:i_parent_end = 74 ;
:j_parent_end = 61 ;
:parent_grid_ratio = 1 ;
:FLAG_MF_XY = 1 ;
:FLAG_LAI12M = 1 ;
:FLAG_LAKE_DEPTH = 1 ;
}
The global attributes corner_lats and corner_lons contain the lat-lon location of the corners of the domain with respect to different grid staggerings (mass, u, v, and unstaggered). The locations referred to by each element of the corner_lats and corner_lons arrays are summarized in the table and figure below.
全局属性corner_lats和corner_lons包含相对于不同网格交错（mass、u、v和unstaggered）的域角的lat lon位置。下表和下图总结了角点阵列和角点阵列的每个元素所指的位置。
| Array index | Staggering | Corner |
|---|---|---|
| 1 | Mass | Lower-left |
| 2 | Upper-left | |
| 3 | Upper-right | |
| 4 | Lower-right | |
| 5 | U | Lower-left |
| 6 | Upper-left | |
| 7 | Upper-right | |
| 8 | Lower-right | |
| 9 | V | Lower-left |
| 10 | Upper-left | |
| 11 | Upper-right | |
| 12 | Lower-right | |
| 13 | Unstaggered | Lower-left |
| 14 | Upper-left | |
| 15 | Upper-right | |
| 16 | Lower-right |
In addition to the fields in a geogrid output file (e.g., geo_em.d01.nc), the following fields and global attributes will also be present in a typical output file from the metgrid program, run with the default METGRID.TBL file and meteorological data from NCEP's GFS model.
除了geogrid输出文件中的字段（例如geo_em.d01.nc）之外，以下字段和全局属性也将出现在metgrid程序的典型输出文件中，以默认值运行金属网格.TBLNCEP GFS模型的文件和气象数据。
geo_em.d01.nc
netcdf met_em.d01.2016-04-07_00\:00\:00 {
dimensions:
Time = UNLIMITED ; // (1 currently)
DateStrLen = 19 ;
west_east = 73 ;
south_north = 60 ;
num_metgrid_levels = 27 ;
num_st_layers = 4 ;
num_sm_layers = 4 ;
south_north_stag = 61 ;
west_east_stag = 74 ;
z-dimension0132 = 132 ;
z-dimension0012 = 12 ;
z-dimension0016 = 16 ;
z-dimension0021 = 21 ;
variables:
char Times(Time, DateStrLen) ;
float PRES(Time, num_metgrid_levels, south_north, west_east) ;
PRES:units = "" ;
PRES:description = "" ;
float SOIL_LAYERS(Time, num_st_layers, south_north, west_east) ;
SOIL_LAYERS:units = "" ;
SOIL_LAYERS:description = "" ;
float SM(Time, num_sm_layers, south_north, west_east) ;
SM:units = "" ;
SM:description = "" ;
float ST(Time, num_st_layers, south_north, west_east) ;
ST:units = "" ;
ST:description = "" ;
float GHT(Time, num_metgrid_levels, south_north, west_east) ;
GHT:units = "m" ;
GHT:description = "Height" ;
float HGTTROP(Time, south_north, west_east) ;
HGTTROP:units = "m" ;
HGTTROP:description = "Height of tropopause" ;
float TTROP(Time, south_north, west_east) ;
TTROP:units = "K" ;
TTROP:description = "Temperature at tropopause" ;
float PTROPNN(Time, south_north, west_east) ;
PTROPNN:units = "Pa" ;
PTROPNN:description = "PTROP, used for nearest neighbor interp" ;
float PTROP(Time, south_north, west_east) ;
PTROP:units = "Pa" ;
PTROP:description = "Pressure of tropopause" ;
float VTROP(Time, south_north_stag, west_east) ;
VTROP:units = "m s-1" ;
VTROP:description = "V at tropopause" ;
float UTROP(Time, south_north, west_east_stag) ;
UTROP:units = "m s-1" ;
UTROP:description = "U at tropopause" ;
float HGTMAXW(Time, south_north, west_east) ;
HGTMAXW:units = "m" ;
HGTMAXW:description = "Height of max wind level" ;
float TMAXW(Time, south_north, west_east) ;
TMAXW:units = "K" ;
TMAXW:description = "Temperature at max wind level" ;
float PMAXWNN(Time, south_north, west_east) ;
PMAXWNN:units = "Pa" ;
PMAXWNN:description = "PMAXW, used for nearest neighbor interp" ;
float PMAXW(Time, south_north, west_east) ;
PMAXW:units = "Pa" ;
PMAXW:description = "Pressure of max wind level" ;
float VMAXW(Time, south_north_stag, west_east) ;
VMAXW:units = "m s-1" ;
VMAXW:description = "V at max wind" ;
float UMAXW(Time, south_north, west_east_stag) ;
UMAXW:units = "m s-1" ;
UMAXW:description = "U at max wind" ;
float SNOWH(Time, south_north, west_east) ;
SNOWH:units = "m" ;
SNOWH:description = "Physical Snow Depth" ;
float SNOW(Time, south_north, west_east) ;
SNOW:units = "kg m-2" ;
SNOW:description = "Water equivalent snow depth" ;
float SKINTEMP(Time, south_north, west_east) ;
SKINTEMP:units = "K" ;
SKINTEMP:description = "Skin temperature" ;
float SOILHGT(Time, south_north, west_east) ;
SOILHGT:units = "m" ;
SOILHGT:description = "Terrain field of source analysis" ;
float LANDSEA(Time, south_north, west_east) ;
LANDSEA:units = "proprtn" ;
LANDSEA:description = "Land/Sea flag (1=land, 0 or 2=sea)" ;
float SEAICE(Time, south_north, west_east) ;
SEAICE:units = "proprtn" ;
SEAICE:description = "Ice flag" ;
float ST100200(Time, south_north, west_east) ;
ST100200:units = "K" ;
ST100200:description = "T 100-200 cm below ground layer (Bottom)" ;
float ST040100(Time, south_north, west_east) ;
ST040100:units = "K" ;
ST040100:description = "T 40-100 cm below ground layer (Upper)" ;
float ST010040(Time, south_north, west_east) ;
ST010040:units = "K" ;
ST010040:description = "T 10-40 cm below ground layer (Upper)" ;
float ST000010(Time, south_north, west_east) ;
ST000010:units = "K" ;
ST000010:description = "T 0-10 cm below ground layer (Upper)" ;
float SM100200(Time, south_north, west_east) ;
SM100200:units = "fraction" ;
SM100200:description = "Soil Moist 100-200 cm below gr layer" ;
float SM040100(Time, south_north, west_east) ;
SM040100:units = "fraction" ;
SM040100:description = "Soil Moist 40-100 cm below grn layer" ;
float SM010040(Time, south_north, west_east) ;
SM010040:units = "fraction" ;
SM010040:description = "Soil Moist 10-40 cm below grn layer" ;
float SM000010(Time, south_north, west_east) ;
SM000010:units = "fraction" ;
SM000010:description = "Soil Moist 0-10 cm below grn layer (Up)" ;
float PSFC(Time, south_north, west_east) ;
PSFC:units = "Pa" ;
PSFC:description = "Surface Pressure" ;
float RH(Time, num_metgrid_levels, south_north, west_east) ;
RH:units = "%" ;
RH:description = "Relative Humidity" ;
float VV(Time, num_metgrid_levels, south_north_stag, west_east) ;
VV:units = "m s-1" ;
VV:description = "V" ;
float UU(Time, num_metgrid_levels, south_north, west_east_stag) ;
UU:units = "m s-1" ;
UU:description = "U" ;
float TT(Time, num_metgrid_levels, south_north, west_east) ;
TT:units = "K" ;
TT:description = "Temperature" ;
float PMSL(Time, south_north, west_east) ;
PMSL:units = "Pa" ;
PMSL:description = "Sea-level Pressure" ;
// global attributes:
:TITLE = "OUTPUT FROM METGRID V4.0" ;
:SIMULATION_START_DATE = "2016-04-07_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 74 ;
:SOUTH-NORTH_GRID_DIMENSION = 61 ;
:BOTTOM-TOP_GRID_DIMENSION = 27 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 73 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 74 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 60 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 61 ;
:GRIDTYPE = "C" ;
:DX = 30000.f ;
:DY = 30000.f ;
:DYN_OPT = 2 ;
:CEN_LAT = 34.83001f ;
:CEN_LON = -81.03f ;
:TRUELAT1 = 30.f ;
:TRUELAT2 = 60.f ;
:MOAD_CEN_LAT = 34.83001f ;
:STAND_LON = -98.f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:corner_lats = 28.17127f, 44.36657f, 39.63231f, 24.61906f, 28.17842f, 44.37617f, 39.57812f, 24.57806f, 28.03771f, 44.50592f, 39.76032f, 24.49431f, 28.04485f, 44.51553f, 39.70599f, 24.45341f ;
:corner_lons = -93.64893f, -92.39661f, -66.00165f, -72.64047f, -93.80048f, -92.59155f, -65.83557f, -72.5033f, -93.65717f, -92.3829f, -65.9313f, -72.68539f, -93.80841f, -92.57831f, -65.76495f, -72.54843f ;
:MAP_PROJ = 1 ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:grid_id = 1 ;
:parent_id = 1 ;
:i_parent_start = 1 ;
:j_parent_start = 1 ;
:i_parent_end = 74 ;
:j_parent_end = 61 ;
:parent_grid_ratio = 1 ;
:NUM_METGRID_SOIL_LEVELS = 4 ;
:FLAG_METGRID = 1 ;
:FLAG_EXCLUDED_MIDDLE = 0 ;
:FLAG_SOIL_LAYERS = 1 ;
:FLAG_SNOW = 1 ;
:FLAG_PSFC = 1 ;
:FLAG_SM000010 = 1 ;
:FLAG_SM010040 = 1 ;
:FLAG_SM040100 = 1 ;
:FLAG_SM100200 = 1 ;
:FLAG_ST000010 = 1 ;
:FLAG_ST010040 = 1 ;
:FLAG_ST040100 = 1 ;
:FLAG_ST100200 = 1 ;
:FLAG_SLP = 1 ;
:FLAG_SNOWH = 1 ;
:FLAG_SOILHGT = 1 ;
:FLAG_UTROP = 1 ;
:FLAG_VTROP = 1 ;
:FLAG_TTROP = 1 ;
:FLAG_PTROP = 1 ;
:FLAG_PTROPNN = 1 ;
:FLAG_HGTTROP = 1 ;
:FLAG_UMAXW = 1 ;
:FLAG_VMAXW = 1 ;
:FLAG_TMAXW = 1 ;
:FLAG_PMAXW = 1 ;
:FLAG_PMAXWNN = 1 ;
:FLAG_HGTMAXW = 1 ;
:FLAG_MF_XY = 1 ;
:FLAG_LAI12M = 1 ;
:FLAG_LAKE_DEPTH = 1 ;
}
