一、总过程
1.总参考:WPS(使用高精度WPS静态数据)生成的geo_em.nc数据;
2.整合具体方面的数据:
①静态数据:
geo_em.nc+DEM(.tif)+GIS预处理工具+hydro.TBL=Fulldom.nc(径流方面等等)、LDASOUT(地形方面);
geo_em.nc+R语言工具+soil.TBL=wrf_input.nc(初始状态)、soilproperties(土壤方面);
②驱动数据:
geo_em.nc+ESMF(regrided)+气象数据(FORCING,比如GRIB2文件)=DOMAIN1(气象数据FORCING)
3.构建文件夹DOMAIN,FORCING,TBL,exe,RESTART
4.设置hydro.namelist;hrldas.namelist
①输入输出的文件路径、类型(例如IGRID=1)
②物理参数(ZLVL风速、土壤层数、地下水;径流、地面反照率选项)
③冷热启动(RESTART)
④模拟的时间步长(FORCING_TIMESTEP = 3600)
5.运行mpirun -np 4 wrf_hydro_Noah.exe(使用四核,越大越快)
二、常见错误
注1:应当从命令行查看;或者diag.00000查看
注2:一下以namelist为主
注3:warning对于模型运行可忽略;可能对后续数据提取有影响
1.文件名称、路径设置错误:(bad file describe)
(OUTDIR = "./result",没有result文件夹!改成"./"就好)
2.变量设置错误:
①缺少该设置,加上就好(常常是注释掉了)
②设置错误(常常是手贱打多了)
(hydro,hrldas里面相关的soil设置,比如说NSOIL)
(时间步长KHOUR=2 > FORCING的数据量)(diag.00000)
( ./result× ./√)(diag.00000)
3.系统错误
①mpirun -np 4(一直用4核,打开多个程序就会冲突)
三、尚未解决
1.可能还是namelist设置问题?
2.python版本问题?(LDASOUT的gis生成)
3.各种warning(非必要,但可能后面数据缺失)
①缺少了也能跑
②没有径流点
③LDASIN版本问题
④读取径流深度
四、解决问题的总体思路
1.提出大致猜测:可能是namelist、nc、FORCING问题
2.简单增加理解:通过说明书理解各变量,文件的意义
3.总结:模型整体的运作过程,文件的大致意义
4.提出具体猜测:(状态好的时候要先写下总体框架)
总分析:例子成功=模型本身没问题=②配置出问题=①生成的nc有问题
①wps配置问题=geo_em有问题OR=GIS处理文件的问题=LDASOUT的问题
②仅仅是namelist配置的问题
5.控制变量、逐个验证、并记录结果:(知道geo_em,wrf这种主体数据替换不可验证)
①自己的DOMAIN,FORCING,namelist
②正确的DOMAIN,FORCING,namelist
③怀疑是LDASOUT问题?用自己的LDASOUT替换正确的
怀疑是FORCING问题?用自己的FORCING替换正确的
怀疑是namelist问题?用自己的替换正确的
怀疑DOMAIN和namelist要配套?用自己的局部替换正确的