【气象水文】WRF示例运行与后处理软件

【目的】用于备忘WRF的安装过程，记录所遇到问题与解决方式；

【安排】该内容分为三个部分：WRF模型安装、示例运行、可视化环境配置。

【操作系统与软件】Windows系统与其子系统Ubuntu系统

【参考】

示例运行：https://blog.csdn.net/qq_40505953/article/details/124080415

ncview安装：服务器安装ncview-专业气象研究-气象家园_气象人自己的家园

【安装步骤】

1、WRF模型环境配置安装参考该栏目上一篇博文；

2、WRF案例运行-Single Domain Case

（1） 静态地理背景数据下载

静态地理背景数据下载地址

# 服务器Build_WRF文件夹下新建文件夹WPS_GEOG，将下载的数据压缩包放入其中，并解压；

 # 解压过程中，可能会遇到许可限制不能创建文件夹，解决方式为获取root权限后再解压，如下：

首先，临时的获取到root权限，使用root权限更改root用户的密码。然后利用su命令键入我们刚刚更改过的root密码。

键入当前的用户的密码之后，在输入两遍root用户的新密码，这样就把root用户的密码给更改了。

具体参考：Linux 下获取Root权限的几种方法-蒲公英云

# 由于存在多个压缩包，文件夹内循环解压是最好的解决方式。笔者采用如下：

for tar in *.tar.bz2; do tar xvf $tar; done

 # 其他格式的压缩包解压方式参考：linux 批量解压gz文件夹,linux 批量解压gz bz2文件_Z198的博客-CSDN博客

（2）更改 namelist.wps文件

# cd到安装WPS的文件夹，利用命令行输入：

nano namelist.wps

# 更改后的 namelist.wps文件

&share
 wrf_core = 'ARW',
 max_dom = 1,
 start_date = '2000-01-26_12:00:00',
 end_date   = '2000-01-26_00:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,
 parent_grid_ratio =   1,
 i_parent_start    =   1,
 j_parent_start    =   1,
 e_we              =  100,
 e_sn              =  100,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the MAXSNOALB and ALBEDO12M
 ! fields have changed in WPS v4.0. These fields are now interpolated
 ! from MODIS-based datasets.
 !
 ! To match the output given by the default namelist.wps in WPS v3.9.1,
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = 'default',
 dx = 27000,
 dy = 27000,
 map_proj = 'mercator',
 ref_lat   =  41,
 ref_lon   =  86,
 truelat1  =  30.0,
 truelat2  =  60.0,
 stand_lon =  86,
 geog_data_path = '/home/jiang/WRF/BuildWRF/WPS_GEOG/'
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE'
 io_form_metgrid = 2,
/

 注意geog_data_path = '/home/jianghu/WRF/BuildWRF/WPS_GEOG/'一定要换成自己的路径，然后保存+退出（ctrl+s,ctrl+x)。

参考：WRF案例运行-Single Domain Case_GISer凌的博客-CSDN博客

(3)运行geogrid.exe

(base) jiang@LAPTOP-TIPVR7MR:~/WRF/BuildWRF/WPS$ ./geogrid.exe

结果输出如下：

(base) jiang@LAPTOP-TIPVR7MR:~/WRF/BuildWRF/WPS$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1
  Processing XLAT and XLONG
  Processing MAPFAC
  Processing F and E
  Processing ROTANG
  Processing LANDUSEF
  Calculating landmask from LANDUSEF ( WATER =  17  21 )
  Processing HGT_M
  Processing SOILTEMP
  Processing SOILCTOP
  Processing SCT_DOM
  Processing SOILCBOT
  Processing SCB_DOM
  Processing ALBEDO12M
  Processing GREENFRAC
  Processing LAI12M
  Processing SNOALB
  Processing CON
  Processing VAR
  Processing OA1
  Processing OA2
  Processing OA3
  Processing OA4
  Processing OL1
  Processing OL2
  Processing OL3
  Processing OL4
  Processing VAR_SSO

  Optional fields not processed by geogrid:
    LAKE_DEPTH (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/lake_depth/')
    URB_PARAM (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/NUDAPT44_1km/')
    FRC_URB2D (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/urbfrac_nlcd2011/')
    IMPERV (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/nlcd2011_imp_ll_9s/')
    CANFRA (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/nlcd2011_can_ll_9s/')
    EROD (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/erod/')
    CLAYFRAC (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/clayfrac_5m/')
    SANDFRAC (priority=1, resolution='default', path='/home/jiang/WRF/BuildWRF/WPS_GEOG/sandfrac_5m/')

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of geogrid.        !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

3、数据后处理：ncl安装和ncview安装，以及linux下的图形展示配置

（1）miniconda 安装 

# minconda是安装后面两个图形显示工具的前提，不然其他可实现的安装过程极其繁琐，笔者花费大量时间已经为大家踩过坑了，希望后者能走捷径！

wget https://mirrors.bfsu.edu.cn/anaconda/miniconda/Miniconda3-py37_4.8.3-Linux-x86_64.sh
sh Miniconda3-py37_4.8.3-Linux-x86_64.sh

# 安装后，添加环境变量：

export PATH="~/miniconda3/bin:"$PATH
source ~/.bashr

# 安装mininconda后，换源：

conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/pkgs/free/ 
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/bioconda/ 
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/msys2/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/menpo/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/peterjc123/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/pkgs/main/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/conda-forge/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/cloud/pytorch/
conda config --add channels https://mirrors.bfsu.edu.cn/anaconda/pkgs/r/
conda config --set show_channel_urls yes

参考：Linux安装miniconda_新客草草的博客-CSDN博客_linux安装miniconda

（2）ncl 安装

# conda 环境下创建ncl环境并且激活：

conda create -n ncl_stable -c conda-forge ncl
source activate ncl_stable

# 添加环境变量

#适合本机的环境变量
export NCARG_ROOT=/home/jiang/miniconda3/envs/ncl_stable
export PATH=$NCARG_ROOT/bin:$PATH

（3）ncview 安装 （仅支持Ubuntu系统，默认安装版本为ncview2.1.8）

sudo apt-get install ncview

# 参考：ncview简易安装与运行_钢筋火龙果的博客-CSDN博客_安装ncview

（4）图形配置：

# 安装x11

sudo apt-get install x11-apps

# 添加环境变量DISPLAY，指定windows主机实现图形显示：根据IPv4本地机IP修改：

vim ~/.bashrc
export DISPLAY=192.168.***.***:0.0
source ~/.bashrc

# Windows10本机安装Xming

下载地址：https://sourceforge.net/projects/xming/?source=typ_redirect

配置：打开Xming的安装目录，找到文件“X0.hosts”（刚安装的Xming改文件名一般为X0，也可能X1），以文本形式打开这个文件，在第二行补充服务器IP地址：

localhost
127.0.0.***

# 更改后，重启Xming

参考：Centos服务器下ncview 安装_栩栩913的博客-CSDN博客_centos安装ncview

# 尝试 ncview geo_em.d01.nc