【WRF-Hydro】geo_em.d01.nc和Fulldom.nc的创建

一、什么是Fulldom.nc文件

1.更高精度：用经过水文处理的DEM（CON）进行插值形成的，包含了更高精度的地形高程、径流方向、土地利用等数据（注意：单纯的DEM文件必须要经过处理，转化为有用的信息，比如水文处理）

2.水文过程：相比LSM-grid，Routing-grid的分辨率更高（一般来说），作为一个分布式水文模型WRF-Hydro，它使得流域可以划分地更细，从而在每个方块上进行产流过程的计算，最后再汇流。

二、如何生成nc文件

需要的文件有：WPS生成的geo_em.d01.nc文件（粗糙域，包含了分辨率较低的数据）；高精度的、经过处理的DEM文件

重点：WPS坐标系、研究区的划定以及LSM-grid大小的设置；目标区域的DEM文件

1.WPS的geo_em_d01.nc

①位置的确定：确定研究区的范围(33.5N~41N，124.5E~130E)

                         确定研究区中心点的经纬度（ref_lat:37.25N°，ref_lon:127.25E°）

                         确定东西，南北方形的距离（纬向7.5°，832.6km；经向5.5°，526km）

②区域大小的划定：确定单位网格大小(dx=5000,dy=5000,分辨率为5km)

                                确定网格数量（e_sn=832.6/5=167，e_we=526/5=106）

③全球参照系的设定：确定参照系(proj=Lambert,分为低中高纬度选择)

                                    确定纬度准确点（truelat_1=30;truelat_2=60）

                                     确定经度参照（stan_lon=180-127.25=52.75）东经+，西经-

④对应数据的选择：确定网格、次网格大小（LSM-grid=5km）

                                确定分辨率的转化3s=100m；15s=500m；1m=15s*4=2km（即分辨率不可超过2.5m即5km）

                                选择数据（geo_data_res=modis_15s,modis_lai_30s）

注意，这里表红的m表示的是分，1度°=60分min=3600秒s，WPS数据大部分的m代表分，具体要看index文件，若proj=regular_all=分；proj=Lambert=米）

注意：geo_data_res的参数选择可以查看WPS中的geogrid文件夹中的GEOGRID.TBL

⑤运行./geogrid.exe

⑥用ncview geo_em_d01.nc查看文件

2.FULLDOM.nc

①下载GIS-Pre-Processing处理工具

ArcGIS Tools for Preparing WRF-Hydro Routing Grids | Research Applications Laboratory (ucar.edu)

②配置其环境（一般要python3.6）

③使用Build_Routing_Stack.py

python Build_Routing_Stack.py -h 查看说明

其中-R根据你径流网格大小决定（LSM-grid=5km，Routing-grid=500m，-R=10）

④寻找DEM，这里推荐hydroSHED的CON文件，分辨率为3s

 推荐数据下载网站：Download HydroSHEDS Core Products

 我选择直接下载了整个亚洲的3s精度，其实可以下载10°*10°类型的

东经，北纬 n10=10~20        e100=100~110

西经，北纬w100=-100~-90

⑤运行：python Build_Routing_Stack.py -i geo_em.d01.nc -d as_con_3s.tif

⑥查看结果：WRF_Hydro_routing_grids.log和压缩包（三个文件）

                      ncview FULLDOM.nc中的flow direction

  三、常见错误

stan_lon设置成了负的×

√

四、对于flow direction和flow accumulation的理解

DEM高程——从最高到相邻最低——→

flow direction指向——流经的个数

正确的拿来做对比试验很重要！！验证假设！！控制变量法

多看说明书，更重要的是有了理解之后，立刻调整相应的设置，然后看结果，并且记录下来，分析！