好文档 - 专业文书写作范文服务资料分享网站

WRF模式简易操作---中文指南

天下 分享 时间: 加入收藏 我要投稿 点赞

v1.0 可编辑可修改 会忽略这些信息。而parse_table子程序识别这些信息与有用的变量信息的重要标志就是第一列包含中画线“-”的第二行。也就是说,程序将忽略第一次出现中画线的那一行以前的所有信息,然后读取第一次出现中画线到第二次出现中画线之间的所有变量信息。如果在第二次出现中画线的那一行后面还有信息,那么将忽略第二次出现中画线的那一行之后的所有信息,直到第三次出现中画线的那一行,同时必须有第四次出现中画线的行。第三和第四次之间有没有变量信息不重要,重要的是第三次和第四次必须同时配套出现。

29

v1.0 可编辑可修改 T213场库 调用子程序parse_table, 读取变量列表文件Vtable 读取要从T213场库中读读取参数列表文件 读取时间参数及输出格式参数 初始化计算时间 时间循环:当计算时间在起止时间之内 文件打开标志赋值flopen=.true. 变量循环:当1≤loopj≤变量总数 读取SST 场数据,并按要求是 要读取的变量为地形 否 要读取的变量为SST 否 call proc_sst 读取场库气象场 是 打开地形文件: 读取地形数据、 关闭地形文件 气象场是否为高斯格点场 是 call GAUSS2LATLON 否 1 30

v1.0 可编辑可修改 1 flopen 为真是 =1 输出类型= =2 创建文本格式 创建AVN第4类=3 否 创建AVN第3类call 读取landsea文件,截取计算区域范围中的水陆标志否 设置输出区域范围 变量为地形 是 call 设置地形场相关call call =1 输出类型= =2 =3 call call call out_mm5 结束:变量循环 结束:时间循环 结束

图1. 读取T213场库资料程序流程图

31

v1.0 可编辑可修改

在变量列表Vtable文件中,这些变量信息设计程表格形式,其中第一列是变量的GRIB码,这些GRIB码是T213场库中规定的编码(详见附录2)。当要使用WRF SI的程序时,还必须将T213场库的GRIB编码转换成AVN的编码形式。第二列是气象要是变量所处的垂直层的GRIB编码。其中100指气压层,105指特定的垂直层,102是指海平面层,1是指包括海面与陆面的地球表面层。第三、四列是相应于具体的次类型的垂直层的值。比如垂直层

类型为100,则第三列就是去各个等压面的值,如1000、950、850等等,则第三列用通配符“*”表示,而第四列为空;如果垂直层类型为105,则第三列则为具体指定的垂直层面的值,如10米高度的风速则用10米等高面的值,第三列为10,而第四列为空;如果垂直层类型为102或1,则第三列取为0,第四列为空。第五列为变量名称,按照GRIB编码的标准进行变量命名。第六列为变量单位,采用GRIB编码的单位标准。第七列为变量文字描述信息,不超过42个字符。在表格中,每一列都是用竖线“|”分隔。第一列的起始处不要加分隔符“|”,但最后一列也就是第七列必须以分隔符结束。后面可以跟一些说明文字信息,但这些信息会被忽略掉。

程序读完变量列表信息后,紧接着便打开参数设置文件“”,读取用户调用场库的起止时间、输出区域的范围以及输出文件的格式。例如: &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

32

v1.0 可编辑可修改 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 /

获取调用场库的起止时间后,就可以初始化计算时间。初始化时间时,调用子程序

build_hdate。每一次调用后都会返回一个格式为

“yyyy-mm-dd_hh:mm:ss”的时间串。以起始时间为开始,在场库中找符合条件的时次,进行场库资料的读取,直到用户定义的调用终止时间为止。

在每一个要读取场库资料的有效时次内,遍历从变量列表中获取得每一个变量。但由于地形、海温和水陆标志信息特殊,要进行特殊处理。另外,由于T213场库中,有一些要素场保存的量仍为高斯格点场资料,故要求将其内插转换到经纬网格点上。这部分功能是调用子程序GAUSS2LATLON来完成的(这个子程序是在陈启英、姚明明提供的转换程序的基础上修改而成)。

最后模式的输出是根据用于定义的输出格式来输出的。在文件中定义的

33

WRF模式简易操作---中文指南

v1.0可编辑可修改会忽略这些信息。而parse_table子程序识别这些信息与有用的变量信息的重要标志就是第一列包含中画线“-”的第二行。也就是说,程序将忽略第一次出现中画线的那一行以前的所有信息,然后读取第一次出现中画线到第二次出现中画线之间的所有变量信息。如果在第二次出现中画线的那一行后面还有信息,那么将忽略第二次出现中画线的那一行之后的所有信息,直到第三次出现中画线的那一行,同时
推荐度:
点击下载文档文档为doc格式
7scvq1yw0y0088t3x4ji0cqsi0v0jd00p7r
领取福利

微信扫码领取福利

微信扫码分享