1、Chapter 1: Overview翻译:by 澳洲的牛牛 laiwf 由于译者水平有限,其中一定存在翻译不妥的地方,希望大家能帮忙指正IntroductionThe Advanced Research WRF (ARW) modeling system has been in development for the past few years. The current release is Version 3, available since April 2008. The ARW is designed to be a flexible, state-of-the-art atmosp
2、heric simulation system that is portable and efficient on available parallel computing platforms. The ARW is suitable for use in a broad range of applications across scales ranging from meters to thousands of kilometers, including: Idealized simulations (e.g. LES, convection, baroclinic waves) Param
3、eterization research Data assimilation research Forecast research Real-time NWP Coupled-model applications Teaching简介Advanced Research WRF (ARW)模式系统在过去的数年中得到了发展。最近公布了第三版，从 2008 年 4 月开始可供使用。ARW 是灵活的，最先进的大气模拟系统，它易移植，并且有效的应用于各种操作系统。ARW 适用于从米到成千上万公里尺度的各种天气系统的模拟，它的功能包括： 理想化模拟（如，LES，对流，斜压波） 参数化研究 数据同化研究 预
4、报研究 实时数值天气预报 耦合模式应用 教学The Mesoscale and Microscale Meteorology Division of NCAR is currently maintaining and supporting a subset of the overall WRF code (Version 3) that includes: WRF Software Framework (WSF) Advanced Research WRF (ARW) dynamic solver, including one-way, two-way nesting and moving n
5、est. The WRF Preprocessing System (WPS) WRF Variational Data Assimilation (WRF-Var) system which currently supports 3DVAR capability Numerous physics packages contributed by WRF partners and the research community Several graphics programs and conversion programs for other graphics toolsAnd these ar
6、e the subjects of this document.The WRF modeling system software is in the public domain and is freely available for community use.NCAR 的中尺度以及微尺度气象部门最近维护以及支持整个 WRF（第三版）的代码的子集，其中包括： WRF 软件框架 Advanced Research WRF (ARW)动力求解方法，包括单向，双向嵌套以及移动嵌套 预处理系统 WRF 多种数据同化系统（WRF-Var），该系统还支持 3 维同化能力 为 WRF 合作伙伴以及研究论坛提
7、供的数值物理包. 一些画图程序和转换适合其他画图工具的程序这也是本文档的主题。WRF 模式系统软件现在是公开以及免费使用的。 As shown in the diagram, the WRF Modeling System consists of these major programs:The WRF Preprocessing System (WPS) WRF-Var ARW solver Post-processing 2) interpolating terrestrial data (such as terrain, landuse, and soil types) to the s
8、imulation domain; and 3) degribbing and interpolating meteorological data from another model to this simulation domain. Its main features include:GRIB 1/2 meteorological data from various centers around the worldMap projections for 1) polar stereographic, 2) Lambert-Conformal, 3) Mercator and 4) lat
9、itude-longitudeNestingUser-interfaces to input other static data as well as met data这个程序的主要用于实时数值模拟。其中包括：1）定义模拟区域；2）插值地形数据（如地势，土地类型，以及土壤类型）到模拟区域；3）从其他模式结果中细致网格以及插值气象数据到此模拟区域。它的主要特点包括： GRIB 1/2 格式全球格点气象数据 地图投影有 1）极地投影，2）兰伯特-保角投影，3）麦卡托投影以及 4）经纬度投影嵌套 由用户接口输入其他数据及 met 数据WRF-VarThis program is optional,
10、but can be used to ingest observations into the interpolated analyses created by WPS. It can also be used to update WRF models initial condition when WRF model is run in cycling mode. Its main features are as follows. It is based on incremental variational data assimilation technique Conjugate gradi
11、ent method is utilized to minimized the cost function in analysis control variable space Analysis is performed on un-staggered Arakawa A-grid Analysis increments are interpolated to staggered Arakawa C-grid and it gets added to the background (first guess) to get final analysis at WRF-model grid Con
12、ventional observation data input may be supplied both in ASCII or “PREPBUFR” format via “obsproc” utility Multiple radar data (reflectivity w* C/ E. N y“ UWPS 是由三个单独的程序geogrid，ungrib 和 metgird 组成。当然，也包括了很多其它的应用程序，这些程序放在 util 目录下。下面是对这三个主要程序的一个简单描述，更详细的内容将在后边的章节进一步介绍。程序geogrid 的目的是确定模拟区域，及把各种地形数据集插值到
13、模式格点上。模拟区域的确定是通过设置 namelist.wps 文件中的与 “geogrid”有关的参数来实现的。除了计算经纬度和地图每个格点的比例因子外，geogrid 还会根据默认值来插值土壤类型、地表利用类型、地形高 度、年平均深层土壤温度、月季植被覆盖、月季反照率、最大的积雪反照率及斜坡的类别。可以通过 WRF 的官方网站来下载这些场的全球数据集，而且这些资料可 以被看成是不随时间改变的，因此只需下一次就可以了。很多数据集只有在某一特定精度上才能使用，但是其它的则可在 30“,2, 5, 和 10中的任何一个精度使用；其中“的意思是弧秒而则是代表弧分（具体见下面的表格）。尽管一个与模拟
14、区域分辨率接近的地形数据被应用后，这个插值后 的地形场可能更具代表性，但是用户仍是无需下载所有可用的分辨率的数据（SJ 认为还是全下载的好）。尽管如此，那些希望应用拥有可以覆盖大范围的格点空间 的区域的用户，还是愿意下载所有分辨率的静态地形数据。除了插值默认的静态数据，geogrid 程序还可以插值进更多的连续的且不同种类的地形到模拟区域。可以通过应用表格文件GEOGRID.TBL 来实现插 值新的或额外的数据集到模拟区域。GEOGRID.TBL 文件定义了所有可以被 geogrid 生成的地形场；它描述了插值一个地形场所需的方法，及所需数 据所放的具体位置。由 geogrid 生成的文件的格
15、式是 WRF I/O API，因此可以通过选择 NetCDF I/O格式使 geogrid 生成 NetCDF 格式的输出文件，以便更方便的用一些外部软件ncview，NCL 和最新版本的 RIP4 来实现可视化（画出 地形图）。ungrib 程序ungrib 程序读 GRIB 文件，“抽出”数据，然后把它们用一个简单的格式写出来，这种格式就是“过渡”格式（格式是 SJ 的个人翻译，具体参考英文文献，及 writing data to the intermediate format 一节以获得关于这个格式的更多细节）。 GRIB 文件包含随时间变化的气象要素场，而且它是从其它区域或全球模式如
16、NCEP 的 NAM 或 GFS 模式而来的。ungrib 程序可以读 GRIB1 的数据，当然，如果编译时选择了 GRIB2 的选项，则就可以读 GRIB2 的数据。GRIB 文件包含的要素场多于启动 WRF 所需的数量。两种不同格式的 GRIB 用了不同的编码来确定变量和在 GRIB 文件中的层次。ungrib 用这些编码表格 Vtable（variable tables）来确定那些场需要从 GRIB 文件里提取出来并写成过渡格式。关于编码的细节可以在 WMO GRIB 文档中找到。对于相同 GRIB 模式输出文件的各种 Vtables 可以在 WPS 主目录下的/ungrib/Varia
17、ble_Tables/里。 WPS为 NAM104 和 212 格点、NAM AWIP 格式、GFS、NCEP/NCAR 再分析、RUC（气压坐标数据和混合坐标数据），AFWA 的 AGRMET 地表模式输出，ECMWF 及其它数据 集。用户可以以其它 Vtable 做模板来为其它的模式输出创造自己的 Vtable；要进一步了解关于 Vtable 中要素场的细节，可以参考 creating and editing Vtables。Ungrib 可以用三个用户可选格式中任何一个来写过渡数据，这三个格式是：WPS一个新的格式，它包含了对接下来的程序有用的额外信息；SIWRF 之前使用过的 过渡格式
18、；MM5用来向 MM5 模式输入 GRIB2 数据。尽管 WPS 格式是被推荐使用的，但是这三个中的任何一个格式都可以被用来启动 WRF。程序气metgrid 程序的作用是把 ungrib 程序提取出的气象要素场水平插值到 geogrid 确定的模拟区域上。这个插值后的数据可以被 WRF 的real 程 序所识别并吸收。metgird 插值的那些数据的时间段可以通过设置 namelist.wps 中 share 记录部分来调整，而且每个模拟区域（最外围区和 嵌套区）的时间都要单独设置。与 ungrib 程序一样，metgird 所处理的数据也是随时间改变的，因此每次做新的模拟时，都要运行 me
19、tgird 程 序。METGRID.TBL 文件是用来控制如何把气象要素场进行插值的。METGRID.TBL 文件为每个要素场都提供了一个区间，在这个区间里，可能会确定诸如要素场的插值方式、 作为标记插值以及要素场所要插值的网格（如 ARW 的 U，V；NMM 的 H，V）。由 metgrid 生成的文件的格式是 WRF I/O API，因此可以通过选择 NetCDF I/O 格式使 metgrid 生成 NetCDF 格式的输出文件，以便更方便的用一些外部软件ncview，NCL 和最新版本的 RIP4 来实现可视化（画出 地形图）。运行 WPSMe运行 WRF Preprocessing
20、System（WPS）有如下三个步骤：1，利用 geogrid 模块确定一个模式的粗糙区域（最外围的范围），及其它嵌套区域2，利用 ungrib 把模拟期间所需的气象要素场从 GRIB 资料集中提取出来3，利用 metgird 把上述的气象要素场（第二步所做的工作）水平插值到模式区域（第一步所做的工作）中当多个模拟在同一区域重复进行时，只需要做一次第一步的工作即可（也就是说 geogrid.exe 所做出的地形资料geo_em.d0*.nc 可以重复 利用）；因此，只有随时间改变的数据才需要在每次模拟时用第二、三步来处理。类似的，如果在多次模拟中，气象数据是类似的，但是地形区域却不断改变的话，
21、 那第二步是可以省略的。下面是各个步骤的详细说明：Step 1: 确定模式区域如 果 WPS 安装成功，会在 WPS 根目录下出现三个可执行程序geogrid.exe, ungrib.exe 和 metgird.exe（原程序在各自的同名子目录下）的链接。除了这三个可执行程序的链接，还有一个 namelist.wps 文 件。下面是一个 WPS 主目录下各组件的列表：drwxr-xr-x 2 4096 arch-rwxr-xr-x 1 1672 clean 气-rwxr-xr-x 1 3510 compile-rw-r-r- 1 85973 compile.output-rwxr-xr-x 1
22、 4257 configure-rw-r-r- 1 2486 configure.wpsdrwxr-xr-x 4 4096 geogridlrwxrwxrwx 1 23 geogrid.exe - geogrid/src/geogrid.exedrwxr-xr-x 3 4096 metgridlrwxrwxrwx 1 23 metgrid.exe - metgrid/src/metgrid.exe-rw-r-r- 1 1101 namelist.wps-rw-r-r- 1 1987 namelist.wps.all_options-rw-r-r- 1 652 namelist.wps.nmm-
23、rw-r-r- 1 4786 READMElrwxrwxrwx 1 21 ungrib.exe - ungrib/src/ungrib.exedrwxr-xr-x 3 4096 utilMete模式的最外层区域和其它嵌套区域都是在 namelist.list 中 geogrid 记录里设置的，另外，在 share 记录里也有需要设置的参数。下面是关 于这两个部分的一个模板，并且如果想获得更多关于每个参数的可能设置及其目的的信息，可以参考 description of namelist variablesr$ g% N i* ainterval_seconds = 21600,io_form_g
24、eogrid = 2/或者 start_year, start_month, start_day, start_hour, end_year, end_month, end_day 和 end_hour）以及气象数据文件的时间间隔（interval_seconds）。在“ungrib”部分，参数 out_format 是用来设置过渡数据文件的格式的；metgird 程序可以读任何一个被ungrib 支持的格式，因此虽然WPS是被推荐的格式, 但是WPS，SI, 和MM5中的任何一个都可以作为选项被选择；另外，用户可以通过设置 prefix 参数来确定过渡数据文件的保存路径和前缀。例如，如果 p
25、refix 被 设成ARGRMET，那过渡数据的文件名就是AGRMET:YYYY-MM-DD_HH，其中 YYYY-MM-DD_HH 是文件中数据的真实的日期。在合理的修改了 namelist.wps 文件后，Vtable 文件也是需要事先被提供的，并且 GRIB 文件也要被链接（拷贝）到 ungrib.exe 所期望的目录下（一 般是 WPS 的主目录，如果 WPS 被安装成功的话）。如果 WPS 要处理气象数据，就要在 Vtable 的辅助下进行才行，即使Vtabel 仅仅被象征性地以 Vtable 来命名并链接到 ungrib.exe 所期望的目录下（一般是 WPS 的主目录）。例如，如
26、果是来自 GFS 模式的 GRIB 数据，可以如此处理：$ ln -s ungrib/Variable_Tables/Vtable.GFS Vtable程序将尝试以 GRIBFILE.AAA,GRIBFILE.AAB, , GRIBFILE.ZZZ 这样的文件名来读取 GRIB 数据。为了能简化把 GRIB 数据链接到到合适的目录，并以这些名字来命名，程序提供了一个 shell 的脚本 link_grib.csh。这个脚本通过读取一个以 GRIB 数据文件名为列表的命令行参数来完成上述工作。例如，如果 GRIB 数据被下载到一个 /data/gfs 的目录下，那这些文件可以用下面的方法来进行链
27、接：$ ls /data/gfs/ j. ! * G% Z-rw-r-r- 1 42728372 gfs_080324_12_00-rw-r-r- 1 48218303 gfs_080324_12_06$ ./link_grib.csh /data/gfs/gfs*上述操作完成后，WPS 目录下的列表将是如下所示：drwxr-xr-x 2 4096 archM-rwxr-xr-x 1 1672 clean-rwxr-xr-x 1 3510 compile-rwxr-xr-x 1 4257 configure-rw-r-r- 1 2486 configure.wps-rw-r-r- 1 1957
28、004 geo_em.d01.ncdrwxr-xr-x 4 4096 geogridlrwxrwxrwx 1 23 geogrid.exe - geogrid/src/geogrid.exe-rw-r-r- 1 11169 geogrid.loglrwxrwxrwx 1 38 GRIBFILE.AAB - /data/gfs/gfs_080324_12_06-rwxr-xr-x 1 1328 link_grib.csh“ s D9 H3 v: d4 Kdrwxr-xr-x 3 4096 metgrid-rw-r-r- 1 1094 namelist.wps-rw-r-r- 1 1987 nam
29、elist.wps.all_options-rw-r-r- 1 1075 namelist.wps.global-rw-r-r- 1 652 namelist.wps.nmmdrwxr-xr-x 4 4096 ungriblrwxrwxrwx 1 21 ungrib.exe - ungrib/src/ungrib.exedrwxr-xr-x 3 4096 utillrwxrwxrwx 1 33 Vtable - ungrib/Variable_Tables/Vtable.GFS修改完 namelist.wps 并链接了合适的 Vtable 和正确的 GRIB 文件后，ungrib.exe 就可
30、以被执行以生成过渡数据格式的气象数据文件了。Ungrib 可按如下命令进行操作：气象,数值$ ./ungrib.exe o8 Cconstants_name = LANDSEA, SOILHGT/第二个功能从多个源数据里获得资源进行插值在有两个或两个以上补充数据（随时间变化）且需要被合并到一起后进行 real.exe处理的情况下是非常有 用的。可以通过 v0 x! 4 = Gaussian (global only!)! 5 = Polar stereographic x(real : nlats ! Number of latitudes north of equator! (for Gau
31、ssian grids)real : xfcst ! Forecast hour of datareal : xlvl ! Vertical level of data in 2-d arrayreal : startlat, startlon ! Lat/lon of point in array indicated by! startloc stringreal : dx, dy ! Grid spacing, kmreal : xlonc ! Standard longitude of projection6real : truelat1, truelat2 ! True latitud
32、es of projectionreal : earth_radius ! Earth radius, kmreal, dimension(nx,ny) : slab ! The 2-d array holding the datalogical : is_wind_grid_rel ! Flag indicating whether winds are! relative to source grid (TRUE) or, U8 V$ o! relative to earth (FALSE)character (len=8) : startloc ! Which point in array
33、 is given by# s: I, F/ / p) n9 ! startlat/startlon; set either q. o7 o; l( c! Gaussianwrite(unit=ounit) hdate, xfcst, map_source, field, &units, desc, xlvl, nx, ny, iprojwrite(unit=ounit) startloc, startlat, startlon, &! Polar stereographicelse if (iproj = 5) then! q 7 b, K/ I9 Iunits, desc, xlvl, n
34、x, ny, iprojxlonc, truelat1, earth_radius! 3) WRITE WIND ROTATION FLAGwrite(unit=ounit) is_wind_grid_rel! 4) WRITE 2-D ARRAY OF DATAwrite(unit=ounit) slab生成并编辑 Vtables尽管已经提供了很多 Vtables，它们也对应了很多普遍使用的数据集，但是仍然不可能把所有与 ungrib 可使用的 GRIB 格式的数据集的相应 Vtables 全部列举。当一个新的源数据需要被 ungrib.exe 处理时，用户可能需要自己生成一个新的 Vtab
35、le，或者以某个已经有的 Vtable 为模板来修改。无论怎样，对 Vtable 中各个要素场的意义和用法的基本了解都是有帮助的。每个 Vtable 可能会包含 7 个或 11个要素场， 这取决于所使用的源数据是 GRIB1 格式还是 GRIB2 格式。每个 Vtable 中的要素场无非是以下三种中的一种：一，描述在GRIB 文件中数据如何被确 认的要素场；二，描述数据如何被 ungrib.exe 和 metgrid.exe 确认的要素场；三，确定 GRIB2 格式的要素场。每个被 ungrib.exe 提取出的变量在 Vtable 中都会对应一或更多行，同时多行数据又会根据层次的不同分割如地
36、面高度和高空层。要素场在 Vtable 文件中是对应一行或全部行，取决于场和层次的具体情况。第一个要素组描述在 GRIB 文件中数据如何被定义的量的模板如下：气Param| Type |Level1|Level2|- - - - GRIB1 参数规定了气象要素场的 GRIB 编码，是在数据集中一系列独一无二的数字。尽管如此，不同的数据集可能会用不同的 GRIB 编码标示相同的要素场例 如，GFS 数据集中的高层气温的 GRIB 码是 11，但是在 ECMWF 数据中则是 130。为了找到某个要素场的 GRIB 码，可以使用 g1print.exe 和 g2print.exe 应用程序。在 GR
37、IB 码中，“Level Type“, “From Level1“,和“From Level2“是用来确定一个要素场可以在那些层次上找到。与“GRIB1 Param“一样，g1print.exe 和 g2print.exe 可以用来找出多层次要素场的值。多层要素场的意义取决于“Level Type“，总结如下：当分层的要素场（Level Type 112）被确定后，层次的起止点是有单位的，这个单位是决定于要素场本身；可以用 g1print.exe 和g2print.exe 应用程序来看这些值。气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,
38、HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛 0 R/ 4 4 M6 t( E% yVtable 的第二个要素场组是来描述在 metgrid 和 real 程序中数据如何被确认的，下面是个模板：3 ?( Q& f. x0 K, M( E) “ T) b| metgrid | metgrid | metgrid |气| Name | Units | Description |这些场中最重要的要数 是“metgrid Name“，它决定了变量的名字，这些名字在 ungrib 写过渡格式文件时被分配到各个气象要素场。这个名字同样要与 METGRID.TBL 文件中的相匹 配，以便于 m
39、etgrid 程序觉得要素场如何被水平插值。“metgrid Units“和“metgrid Description“分别确定了要素场的单位和一个简单的描述；这里需要特别注意的是，当一个要素场没有被描述时，那 ungrib 是不会把这个要素场 写进过渡文件的。最后一组提供了 GRIB2 的具体描述，具体如下：|GRIB2|GRIB2|GRIB2|GRIB2|Discp|Catgy|Param|Level|- 气尽管在一个 Vtable 文件中有 GRIB2 要素场并不能妨碍这个 Vtable 被 GRIB1 数据使用，而且 GRIB2 要素场仅仅在需要使用 GRIB2 数据集 时才会被用到。例
40、如，Vtable.GFS 文件包含了 GRIB2 的要素场，但是却可以被 1（GRIB1）的和 0.5（GRIB2）的 GFS 数据集使 用。由于 WPS 提供了很多与 GRIB2 格式相对于的 Vtables，所以在此就不一一列举了。用 Geogrid 二进制格式写静态数据geogrid.exe 用来插值的地形静态数据以标准的 2-D 和 3-D 数值形式写成一个简单的二进制格式。当用户有了新的静态数据资源可以通过把数据写成二进制形式让 WPS 来使 用。geogrid 格式支持单层或多层连续的场，场的分类（作为区域的分类）及为每种分类给出的一小部分的场。根据二进制格式的含义，下面给出了一个
41、 30-second 的 USGS 土地利用数据集及的模板：对 于一个主导分类的分类要素场，数据首先要存储在一个规则的 2-D 整型数组中，每个整数都给出了与格点相对应的主导分类的说明。考虑到这个数组的作用，数据 是由低（空）到高（空）或由南向北一行行的写入文件中。例如，在上边的图中，nm 的数组的各个元素是按照 x11, x12, ., x1m, x21, . , x2m, . , xn1, . , xnm 的顺序写入的。尽管可以通过调整 “index”文件中的参数endian=little 来设置数据集中排列顺序为 little-endian，在写 入文件时，每个元素是按 big-endi
42、an 的顺序排列的（例如，对于 4-字节整数 ABCD，字节 A 被存储在最低的地址中，而字节 D 被存储在最高的字节 中）。而且虽然在一个数组中用所需的最小数字来代表一个完整的数值范围是最好的，但在这里一个文件中必须用相同的为位数来存储各个元素。当写入二进制数据 时，不用头文件，记录标志，或者多余的位数。例如，一个两位的 10001000 数组文件的大小是 2,000,000。因为 Fortran 无格式的写记录 标志，因此它是不可能直接写出一个 geogrid 的二进制文件的，但是可以通过用 C 或 Fortran 在写或读时调用 read_geogrid.c 和 write_geogri
43、d.c（在 geogrid/src 目录下）来实现。与主导分类类似的是连续或真实的要素场。像主导分类一样，单层连续场作为 2- D 数组首先被组织，然后一行行写入二进制文件。尽管如此，因为一个连续长可以包含非整型或负值，所以每个文件中的存储表示有点复杂。数组中所有的元素都要 先被转成整型。首先通过扫描所有的元素，找出有需要的精度的，然后把剩下的部分通过四舍五入来移除。例如，要求精度达到小数点后三位，则-2.71828 需要被舍入成-2.718。然后所有的的数组元素都被转化成整数，如果有负数，则会有第二次转化：如果用一个字节储存一个元素，则每个负数都会添加 28； 如果是两个字节则是 216；三
44、个字节是 224；四个字节是232。这对无需转化成正数的过程很重要。最后的正的整数数组按照主导分类的格式被写入文件。多 层的联系场类似单层的。对于一个n m r 数组，先按照上边的方法转化成正的整数场，然后从低到高开始写每个层次。分类场可以被看成是一个多层连续场的一部分，它的层次是 k（1 k r）。当把一个场按照 geogrid 二进制格式写入一个文件时，用户应该遵守 geogrid 程序的命名规则，即希望数据文件的名字格式是 xstart- xend.ystart-yend，其中 xstart, xend, ystart,和 yend 分别是 5 位正整数，数值的起止 x 指针、起止 y
45、指针都包含在文件中；指针起始为 1，而不是 0。所以，一个 8001200 的数组 （如 800 行，1200 列）可以被命名为 00001-01200.00001-00800。当一个数据被分成多个部分，每个部分都是一个规则的矩形数组， 每个数组都是一个单独的文件。在这种情况下，数组相对的位置是由每个数组的文件名中的 x 和 y 指针决定的。尽管如此，值得注意的是，数据集中的每个部分都要 有相同的 x-和 y-维数，并且每一块数据不能相互重叠；并且每块的起止都要在索引的范围内。例如，全球的 30-second 的 USGS 地形数据集被分成多 个12001200 的数组，每个数组都代表了一个
46、1010的区域；这个文件的最西南角的位置是（90S, 180W），命名为 00001-01200.00001-01200，而最东北角的位置是（90N, 180E），命名为 42001-43200.20401-21600。如果数据集被分成多个小的部分，而且格点的数字不能在 x-方向上被平分，那最后剩下 的几列数据就要用标志数（index file 文件中参数 missing_value 确定）来填补。例如，一个 x-方向上有 2456 个格点的数据集，如果分成 3 个部分的话，在规定范围内每块应 该是 1820,8211640,和 16412460，而2457-2460 之间就要用标志数（类似缺
47、测值）来填充。因为起 止索引必须是 5 位数字，所以一个场是不能在 x-或 y-方向上超过 99999 个数据格点的，如果超过了，用户要把这个数据集分成多个 geogrid 能识别的 小的部分。例如，一个大的全球数据集可以把东半球和西半球分成多个小的部分。除了二进制数据文件，geogrid 要求每一个数据集都有一个额外的说明文 件，即“index”文件，因此，两个数据集是不能放在同一个目录下的。其实，当要处理数据集时，index 文件是最先被 geogrid 检索的，因为这个 文件里包含了所有可能与用到的数据相关的有用信息，下面是一个index 的模板：type = continuoussigned = yesprojection = regular_lldx = 0.00833333known_x = 1.0known_y = 1.0known_lat = -89.99583known_lon = -179.99583Meteotile_x = 1200tile_z = 1tile_bdr=3description=“Topography height“用户可以通过参考 index file options 来了解所有可能出现在 index 文件中的参数及相关的意义。
