WRF学习笔记1：根据官网模拟Matthew飓风

官网地址：ARW Online Tutorial

官网将该案例研究用于多种不同的模型运行，包括在模型运行中添加 SST 数据的示例

一、ungrib

1.下载并了解数据

Download the Matthew case study data and place in the DATA/ directory

GRIB2，1deg，6h，27pressure(1000-10hPa ; excluding surface)

2016-10-06_00 to 2016-10-08_00 

了解数据：WPS/util/g1print.exe and g2print.exe (for GRIB1 and GRIB2 data)

cd Build_WRF/WPS-4.1/util

./util/g2print.exe ../../../../../mnt/e/data/WRF/matthew/fnl_20161006_00_00.grib2 >& g2print.log

会生成一个日志，在文件夹里

2.链接到 GFS Vtable

ln -s ungrib/Variable_Tables/Vtable.GFS Vtable（GFS要根据下载的气象数据修改）

3.链接到下载的气象数据（数据存在E盘）

./link_grib.csh ../../../../../mnt/e/data/WRF/matthew/fnl

4.编辑namelist.wps

max_dom = 1
start_date = '2016-10-06_00:00:00',
end_date = '2016-10-08_00:00:00',
interval_seconds = 21600,
prefix = 'FILE',

5.运行ungrib得到中间文件

./ungrib.exe

文件夹中产生日志：ungrib.log可以排除障碍，理解数据

理解中间文件./util/rd_intermediate.exe FILE:2016-10-06_00

*加入SST数据

原因：1.很长时间（几年，几十年）需要考虑SST           2.沿海城市需要考虑SST

1、下载SST data

2、先删除之前的Vtable

rm -rf GRIBFILE.A* Vtable

之后链接SST的Vtable文件

ln -sf ungrib/Variable_Tables/Vtable.SST Vtable

3、./link_grib.csh ../DATA/matthew_sst/rtg_sst_grb

4、 编辑 namelist.wps

start_date = '2016-10-06_00:00:00',
end_date = '2016-10-08_00:00:00',
interval_seconds =  21600,
prefix = 'SST',

 5、 ./ungrib.exe

这样做出来的数据就会有SST数据

二、geogrid

1、编辑namelist.wps

max_dom = 1 
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 91,
e_sn = 100,
geog_data_res = 'default',
dx = 27000,
dy = 27000,
map_proj = 'mercator',
ref_lat = 28.00,
ref_lon = -75.00,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -75.0,
geog_data_path = 'Your WPS_GEOG data location'

2、确保domain位置正确

ncl util/plotgrids.ncl

版本问题，改成ncl util/plotgrids_new.ncl

这个ncl文件里面输出方式改为PDF，最后文件存在WPS文件夹下面

3、./geogrid.exe

geo_em.d01.nc文件会被覆盖，所以旧的自己存

会有日志 geogrid.log 

Parsed 20 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 = 16 )
Processing HGT_M
Processing SOILTEMP
Processing SOILCTOP
Processing SCT_DOM
Processing SOILCBOT
Processing SCB_DOM
Processing ALBEDO12M
Processing GREENFRAC
Processing SNOALB
Processing SLOPECAT
Processing SLOPECAT
Processing CON
Processing VAR
Processing OA1
Processing OA2
Processing OA3
Processing OA4
Processing OL1
Processing OL2
Processing OL3
Processing OL4
Processing VAR_SSO
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Successful completion of geogrid.         !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

确保出现：

geo_em.d01.nc

4、检查文件内容

用ncdump
    ncdump -h geo_em.d01.nc

三、Metgrid

1、namelist.wps

不需要再更改了，只需要确保正确

2、./metgrid.exe

有个日志 metgrid.log

Processing domain 1 of 1
Processing 2016-10-06_00
FILE
Processing 2016-10-06_06
FILE
Processing 2016-10-06_12
FILE
Processing 2016-10-06_18
FILE
Processing 2016-10-07_00
FILE
Processing 2016-10-07_06
FILE
Processing 2016-10-07_12
FILE
Processing 2016-10-07_18
FILE
Processing 2016-10-08_00
FILE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of metgrid.  !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

产生文件

met_em.d01.2016-10-06_00:00:00.nc
met_em.d01.2016-10-06_06:00:00.nc
met_em.d01.2016-10-06_12:00:00.nc
met_em.d01.2016-10-06_18:00:00.nc
met_em.d01.2016-10-07_00:00:00.nc
met_em.d01.2016-10-07_06:00:00.nc
met_em.d01.2016-10-07_12:00:00.nc
met_em.d01.2016-10-07_18:00:00.nc
met_em.d01.2016-10-08_00:00:00.nc

检查文件ncdump -h met_em.d01.2016-10-06_00:00:00.nc

四、WRF

1、找到地方

cd ..    

cd WRF-4.1.2/test/em_real或者cd run跑WRF

2、链接

找到刚才生成的WPS下的met_em.d0*链接到文件夹里

ln -s ../../../WPS-4.1/met_em.d0* .（注意这里的点）或者ln -sf ../../../WPS-4.1/met_em.d01.2016-10* .

3、编辑namelist.input

run_days = 0,
run_hours = 48,
run_minutes = 0,
run_seconds = 0,
start_year = 2016,
start_month = 10,
start_day = 06,
start_hour = 00,
end_year = 2016,
end_month = 10,
end_day = 08,
end_hour = 00,
interval_seconds = 21600
input_from_file = .true.,
history_interval = 180,
frames_per_outfile = 1,
restart = .false.,
restart_interval = 1440,
time_step = 150,
max_dom = 1,
e_we = 91,
e_sn = 100,
e_vert = 45,
num_metgrid_levels = 32
dx = 27000,
dy = 27000,

4、./real.exe

wrfinput_d01（旧的自己保存）
wrfbdy_d01

5、./wrf.exe 

wrfout_d01_2016-10-06_00:00:00
wrfout_d01_2016-10-06_03:00:00
wrfout_d01_2016-10-06_06:00:00
wrfout_d01_2016-10-06_09:00:00
wrfout_d01_2016-10-06_12:00:00
wrfout_d01_2016-10-06_15:00:00
wrfout_d01_2016-10-06_18:00:00
wrfout_d01_2016-10-06_21:00:00
wrfout_d01_2016-10-07_00:00:00
wrfout_d01_2016-10-07_03:00:00
wrfout_d01_2016-10-07_06:00:00
wrfout_d01_2016-10-07_09:00:00
wrfout_d01_2016-10-07_12:00:00
wrfout_d01_2016-10-07_15:00:00
wrfout_d01_2016-10-07_18:00:00
wrfout_d01_2016-10-07_21:00:00
wrfout_d01_2016-10-08_00:00:00
wrfrst_d01_2016-10-07_00:00:00
wrfrst_d01_2016-10-08_00:00:00

wrfirst文件将包含07_00和08_00的时间信息，以便可以从这些点重新启动模型。

6、快速检查结果

ncdump -h wrfout_d01_2016-10-06_00:00:00
ncdump -v Times wrfout_d01_2016-10-06_00:00:00 (查看文件中有哪些预测时间)

一个小技巧：

 ncdump -h wrfout_d01_2016-10-06_00:00:00 >xx.log

可以生成一个日志在文本文件查看

可以将数据储存在其他文件夹中，查看时打开文件夹再输入命令

文件存储小技巧：

可以将生成的wrfout文件存储在E盘中，每次运行wrf，test/em_real中之前的文件就可能被替换，可以在test/em_real中运行

cp wrfout_d01_2016* ../../../../../../mnt/e/data/WRF/matthew/wrfout_d01

将文件存在其他地方，之前的文件可以批量删除

rm -rf wrfout_d01_2016*

五、ncl画图

总结：

WPS+WRF

WPS：先编辑namelist.wps，链接下载的气象数据，之后运行ungrib、geogrid、metgrid

WRF：先编辑namelist.input，链接WPS生成的数据，之后运行real.exe、wrf.exe