使用说明¶
整个系统的工作目录如图5所示,其中:App 为格点化、格式转换等相关程序和脚本目录;BaseData 为雷达基数据目录;BufrData 为转换成BUFR格式的数据目录;GridData 为格点化数据目录;ImageData 为反射率组网拼图图片目录;MosaicExe 为雷达拼图程序目录;MREF 为反射率拼图数据目录;QCData 为质控后的雷达数据目录;Temp 为临时文件目录;Install 为相关安装程序目录;Script 则为相关程序脚本目录。
整个系统的运行逻辑和步骤如下:
主程序( /App/main.exe )¶
该程序启动后会保持在后台运行,根据相关设置读取雷达基数据,进行质量控制(质控后的雷达数据存储至../QCData目录下),然后对质控后的雷达数进行据格点化处理(格点雷达数据存储至 ../GridData 目录下)。
程序相关的配置文件( /App/config.ini )说明如下:
[PATH]
BaseDataPath=c:\Test\BaseData #雷达基数据路径
GridDataPath=c:\Test\GridData #雷达格点数据路径
QCDataPath=c:\Test\QCData #雷达质控数据路径
MREFPath=c:\Test\MREF #雷达拼图数据路径
BUFRPath=c:\Test\BufrData #雷达BUFR格式数据路径
[RADAR]
RadarTotal=9 #雷达站数量
雷达组网拼图程序( /MosaicExe/RadarMosaicRT.exe )¶
该程序启动后会在后台运行,根据相关设置将已经格点化的雷达数据进行组网拼图,并将生成的拼图数据导出至/MREF目录下存储。
程序相关的配置文件( /MosaicExe/SysConfigParams.ini )说明如下:
#拼图省份
Fujian
#参与拼图的雷达站数量
9
#站名 站号 经度 纬度 海拔 型号 探测范围 格点数据路径
福州 9591 119.5319 25.9819 642.6 SA 230 C:\Test\GridData\z9591\
厦门 9592 118.0789 24.4819 185.1 SA 230 C:\Test\GridData\z9592\
宁德 9593 119.6781 26.5642 923.0 SA 230 C:\Test\GridData\z9593\
泉州 9595 118.5869 24.9075 531.0 SA 230 C:\Test\GridData\z9595\
漳州 9596 117.4 23.5994 250.0 SA 230 C:\Test\GridData\z9596\
龙岩 9597 117.1819 25.0578 1506.7 SA 230 C:\Test\GridData\z9597\
三明 9598 117.6381 26.2633 302.6 SA 230 C:\Test\GridData\z9598\
建阳 9599 118.1008 27.35 291.3 SA 230 C:\Test\GridData\z9599\
南昌 9791 115.8989 28.5908 94.9 SA 230 C:\Test\GridData\z9791\
#临时文件路径
C:\Test\Temp\
#系统运行主路径
C:\Test\
#自动站数据路径
C:\Test\
#图片导出路径
C:\Test\ImageData\
#经向分辨率
0.01
#纬向分辨率
0.01
#拼图产品垂直层数
9
#拼图产品各层高度(单位:米)
1000
1500
2000
2500
3000
3500
4000
4500
5000
#拼图产品时间宽度(多少分钟以内的数据都认为是该时刻参与拼图的数据)
4
#拼图产品时间分辨率
6
#延迟时间(当某个时次的数据未到齐时的最长等待时间)
2
雷达数据BUFR格式转换程序( /App/BufrConverter.sh )¶
由于BUFR格式的转换涉及Linux平台,因此该功能通过调用shell脚本(.sh)来实现。BufrConverter.sh 分别完成雷达反射率拼图(process_mosaic)和雷达径向速度数据(process_vel)的BUFR格式转换功能,生成的BUFR格式数据( Ref_xxx.bufr、Vel_xxx.bufr )导出至 /BUFRData 目录下存储。其中,该脚本将在main.exe程序执行时自动被调用(监测到有新的雷达反射率组网拼图数据生成时触发),但需将 bash.exe(C:\Windows\System32\bash.exe) 程序拷贝到 App 目录下。
程序相关的配置文件( /App/namelist.input )说明如下:
&vel !雷达径向速度BUFR格式转换
InputPath='/mnt/c/Test/QCData', !输入数据路径(质控后雷达数据)
OutputPath='/mnt/c/Test/BufrData', !输出数据路径(转换后BUFR数据)
DataTime='201807110012', !数据时间(世界时)
TimeGap=4, !时间间隔(分钟)
RadarCount=9, !雷达站数量
RadarType='SA','SA','SA','SA','SA','SA','SA','SA','SA', !雷达型号
RadarNum='Z9591','Z9592','Z9593','Z9595','Z9596','Z9597','Z9598','Z9599','Z9791', !雷达站号
/
&mosaic !雷达反射率拼图BUFR格式转换
InputPath='/mnt/c/Test/MREF', !输入数据路径(反射率拼图数据)
OutputPath='/mnt/c/Test/BufrData', !输出数据路径(转换后BUFR数据)
FileName='Z_OTHE_RADAMOSAIC_20180711001200.bin,!输入拼图数据名称
/
Note
注意:由于脚本是在Linux环境下执行,因此涉及路径的变量一定要加上 /mnt/ 。时间间隔则表示BUFR数据生成时刻前后多少分钟内的雷达数据被采用。如果要重新配置待处理的雷达站信息,除了修改相应的配置文件以外,还需要修改 ModuleRadarPosition.f90(/Script/VelBufr) ,具体修改方式可参考文件内的代码注释。
雷达数据格点化程序( /App/3DGrid.exe )¶
该程序在 main.exe 运行过程中自动被调用,主要完成雷达数据格点化功能。
程序相关的配置文件( /App/SglMosaicParams.ini )说明如下:
#参与格点化的雷达站数量
9
#站名 站号 经度(度、分、秒) 纬度(度、分、秒) 海拔 型号 探测范围(径向速度、反射率) 距离库长 方位角分辨率 地形遮挡标志位 方位屏蔽标志位
福州 9591 119 31 55 25 58 55 642.6 SA 230 460 1000 0.95 1 1
厦门 9592 118 04 44 24 28 55 185.1 SA 230 460 1000 0.95 1 1
宁德 9593 119 40 41 26 33 51 923 SA 230 460 1000 0.95 1 1
泉州 9595 118 35 13 24 54 27 531 SA 230 460 1000 0.95 1 1
漳州 9596 117 24 00 23 35 58 250 SA 230 460 1000 0.95 1 1
龙岩 9597 117 10 55 25 03 28 1506.7 SA 230 460 1000 0.95 1 1
三明 9598 117 38 17 26 15 48 302.6 SA 230 460 1000 0.95 1 1
建阳 9599 118 06 03 27 21 00 291.3 SA 230 460 1000 0.95 1 1
南昌 9791 115 53 56 28 35 27 94.9 SA 230 460 1000 0.95 1 1
#经向分辨率
0.01
#纬向分辨率
0.01
#格点化产品垂直层数
9
#格点化产品各层高度(单位:米)
1000
1500
2000
2500
3000
3500
4000
4500
5000
Note
注意:雷达原始数据、格点数据和质控数据都将在各自的根目录下根据站号分别存储。
程序运行日志( /App/Log )¶
程序运行产生的各种日志文件均放置在Log文件夹下( /App/Log ),并按日期建立子目录存放。主程序的日志名为 Main.log ,而雷达反射率拼图和径向速度的BUFR格式转换日志分别对应于 RefBufr.***.log 和 VelBufr.***.log 。
系统关联说明¶
本系统的核心内容是为RAP/HRRR模式系统的雷达资料同化提供符合系统格式要求(BUFR格式)并经过质控的雷达数据。因此,本系统属于独立于RAP/HRRR模式系统及其他子系统的独立子系统,主要负责雷达数据的实时收集、读取、质控和格式转换工作。系统可单独运行在服务器端,由ECFLOW在每个时次(cycle)设置数据监控任务(task)实现运行管理的目的,即通过时间触发的数据监控任务实时监测本系统处理生成的雷达BUFR格式数据的到达情况。其中,数据到达监测需在资料同化(DA)步骤之前执行,若该时次雷达资料同化所需数据已存在则直接运行资料同化步骤,若不存在则暂停,等待数据达到一段时间;若等待时间超过设置的阈值则提示数据缺失错误,此时可查询相关日志了解错误原因。
