实用标准文档
值得注意的是，由于场库资料中整数和浮点数都是以64位格式存储的，也就是8个字节。因此，在程序的声明部分注意场库调用参数的存储长度；或者在编译程序的时候，选用64位编译方式编译。当然，也可以在编译的时候，用编译选项 -qintsize=8 和 -qrealsize=8两个选项来进行编译。
index_dic是模式字典文件全名，在调用前必须指定。在IBM SP上，场库在SP:/gpfs/fs2/nwp/gdbt213目录下用index_dic(1:28)='/u/nwp/gdbt213/dic/index_dic'；在场库在盘阵库/asm2/nwpgdb/t213gdb目录下用index_dic(1:32)='/u/nwp/gdbt213 /dic/index_nfs_dic’。
3.2.2 读取T213场库资料的程序设计
读取T213场库资料的基本思路是首先设置好相应的气象要素场的具体调用参数；然后利用T213场库提供的通用接口程序读取场库资料，保存在内存变量中；再将内存变量中的相应气象要素场数据按不同的要求输出为不同的数据文件格式。根据这一思想，可以得到读取T213场库资料的程序流程图（见图10）。 实现读取T213场库资料的程序名称为readt213.F。程序开始调用子程序parse_table来变量列表Vtable文件中的变量信息。Vtable的格式如下：
GRIB Code| Level Code| Level 1 | Level 2 | REGRID Name | REGRID Units | REGRID Description |
-----+------+------+------+----------+----------+------------------------------------------+
80 | 1 | 0 | | SST | K | Sea Surface Temperature | 1 | 1 | 0 | | PSFC | Pa | Surface Pressure | 2 | 102 | 0 | | PMSL | Pa | Sea-level Pressure |
7 | 1 | 0 | | HGT | m | Height | Terrain Height
211 | 1 | 0 | | SKINTEMP | K | Skin temperature (can use for SST also) | 229 | 105 | 2 | | T | K | Temperature | At 2 m 230 | 105 | 2 | | RH | % | Relative Humidity | At 2 m 202 | 105 | 10 | | U | m s{-1}| U | At 10 m 203 | 105 | 10 | | V | m s{-1} | V | At 10 m 7 | 100 | | | HGT | m | Height | 11 | 100 | * | | T | K | Temperature | 33 | 100 | * | | U | m s{-1} | U | 34 | 100 | * | | V | m s{-1} | V | 52 | 100 | * | | RH | % | Relative Humidity | -----+------+------+------+----------+----------+---------------------------------------------+
变量列表Vtable文件中，第一行是说明信息，调用parse_table子程序时，会忽略这些信息。而parse_table子程序识别这些信息与有用的变量信息的重要标志就是第一列包含中画线“-”的第二行。也就是说，程序将忽略第一次出现中画线的那一行以前的所有信息，然后读取第一次出现中画线到第二次出现中画线之间的所有变量信息。如果在第二次出现中画线的那一行后面还有信息，那么将忽略第二次出现中画线的那一行之后的所有信息，直到第三次出现中画线的那一行，同时必须有第四次出现中画线的行。第三和第四次之间有没有变量信息不
文案大全
实用标准文档
重要，重要的是第三次和第四次必须同时配套出现。
文案大全
实用标准文档
T213场库 调用子程序parse_table， 读取变量列表文件Vtable 读取要从T213场库中读取的变量列表 读取参数列表文件namelist.input 读取时间参数及输出格式参数 初始化计算时间 时间循环：当计算时间在起止时间之内 文件打开标志赋值flopen＝.true. 变量循环：当1≤loopj≤变量总数 读取SST 场数据，并按要求格式输出 要读取的变量为地形？ 否 是 要读取的变量为SST？ 否 call proc_sst 读取场库气象场 是 打开地形文件:t213_terrain.dat 读取地形数据、 设置场的相关信息 关闭地形文件 气象场是否为高斯格点场? 是 call GAUSS2LATLON 否 1 文案大全
实用标准文档 1 flopen 为真(.true.)? 是 =1 输出类型＝? =2 创建文本格式 输出文件 创建AVN第4类格式输出文件 =3 否 创建AVN第3类格式输出文件 call proc_landsea 读取landsea文件，截取计算区域范围中的水陆标志信息，并按要求格式输出 否 设置输出区域范围 变量为地形？ 是 call get_sec1 设置地形场相关输出信息 call get_sec2 call get_gridinfo =1 输出类型＝? =2 =3 call out_pregrid call out_hinterp call out_mm5 结束：变量循环 结束：时间循环 结束 文案大全
图1. 读取T213场库资料程序流程图
实用标准文档
在变量列表Vtable文件中，这些变量信息设计程表格形式，其中第一列是变量的GRIB码，这些GRIB码是T213场库中规定的编码（详见附录2）。当要使用WRF SI的prep_grib.exe程序时，还必须将T213场库的GRIB编码转换成AVN的编码形式。第二列是气象要是变量所处的垂直层的GRIB编码。其中100指气压层，105指特定的垂直层，102是指海平面层，1是指包括海面与陆面的地球表面层。第三、四列是相应于具体的次类型的垂直层的值。比
如垂直层类型为100，则第三列就是去各个等压面的值，如1000、950、850等等，则第三列用通配符“*”表示，而第四列为空；如果垂直层类型为105，则第三列则为具体指定的垂直层面的值，如10米高度的风速则用10米等高面的值，第三列为10，而第四列为空；如果垂直层类型为102或1，则第三列取为0，第四列为空。第五列为变量名称，按照GRIB编码的标准进行变量命名。第六列为变量单位，采用GRIB编码的单位标准。第七列为变量文字描述信息，不超过42个字符。在表格中，每一列都是用竖线“|”分隔。第一列的起始处不要加分隔符“|”，但最后一列也就是第七列必须以分隔符结束。后面可以跟一些说明文字信息，但这些信息会被忽略掉。 程序读完变量列表信息后，紧接着便打开参数设置文件“namelist.input”，读取用户调用场库的起止时间、输出区域的范围以及输出文件的格式。例如： &filetimespec
START_YEAR = 2003 START_MONTH = 04 START_DAY = 07 START_HOUR = 00 START_MINUTE = 00 START_SECOND = 00 END_YEAR = 2003 END_MONTH = 04 END_DAY = 10 END_HOUR = 00 END_MINUTE = 00 END_SECOND = 00 INTERVAL = 21600 START_X = 1 START_Y = 1 DELTA_X = 320 DELTA_Y = 160 /
! OUT_TYPE IS A OUTPUT TYPE, THE VALID VALUE IS: ! 1 : PREGRID FORMAT ! 2 : HINTERP FORMAT ! 3 : MM5 FORMAT &typespec OUT_TYPE=2 /
文案大全
WRF模式简易操作 - 中文指南设计 - 图文
