天河中运行WFR模式进行植被覆盖变化的气候效应模拟-----WPS部分
首先建立自己的文件夹cuihn,将wps文件夹复制到自己的文件夹下,开始修改namelist.wps:
- geogrid.exe ,建立地面数据集;修改namelist.wps中的 &geogrid 变量);具体修改如下:
&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '1990-03-01_00:00:00','1990-03-01_00:00:00',
end_date = '1990-11-01_00:00:00','1990-11-01_00:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 50,
j_parent_start = 1, 30,
e_we = 190, 271,
e_sn = 120, 181,
geog_data_res = 'usgs_lakes+default', 'usgs_lakes+default',
dx = 30000,
dy = 30000,
map_proj = 'lambert',
ref_lat = 28.25,
ref_lon = 110,
truelat1 = 24.0,
truelat2 = 32.0,
stand_lon = 110.0,
geog_data_path = '/BIGDATA1/igsnrr_xzzhang_1/wrfdata/geog_data/WPS_GEOG/'
/
命令为:
打开:vi namelist.wps
编辑:光标移至文档最后输入 i
保存:先按Esc退出编辑模式,然后光标移至文档结尾输入: :wq
(保存并退出);:q
(退出);:q!
(强制退出) 详细学习LINUX
编辑退出后,执行:
yhrun -n 1 -N 1 -p bigdata ./geogrid.exe >& log.geogrib
得到nc文件,如下所示(本次选取2层嵌套,因此生成geo_em.d01.nc和geo_em.d02.nc):
在matlab中打开研究区和domain1,domain2检查:
2.ungrib.exe(解码大尺度环流数据,建立大气强迫场;修改namelist.wps中的 &share 变量)
首先与大气数据建立链接,然后运行ungrib.exe;由于第一个实验使用NCEP2数据,有两种文件,需要先链接flx文件,运行ungrib.exe,生成前缀名为FLX的文件,再链接pgb文件,运行ungrib.exe,生成前缀名为PGB的文件
在此之前,需要设定所使用 Vtable 文件,对于不同的气象数据,需要不同的Vtable来描述(NCEP2用Vtable.NCEP2;NCEP_FNL用Vtable.GFS)。WRF所支持的气象数据集对应的 Vtable 保存在WPS代码目录的 ungrib/Variable_Tables/ 中,链接相应的Vtable至wps中的ungrib.exe运行目录下:
ln -sf ungrib/Variable_Tables/Vtable.NCEP2 Vtable
因此第一次的namelist.wps设置为:
&ungrib
out_format = 'WPS',
prefix = 'FLX',
Linux(1链接&2运行)命令如下:
./link_grib.csh /BIGDATA1/igsnrr_xzzhang_1/wrfdata/ncep2/flx.ft06.1990*.tar
yhrun -n 1 -N 1 -p bigdata ./ungrib.exe >& log.ungrib
ps:链接完成wps运行目录中出现如下图文件才是链接成功:
第二次的namelist.wps设置为:
&ungrib
out_format = 'WPS',
prefix = 'PGB',
Linux命令如下:
./link_grib.csh /BIGDATA1/igsnrr_xzzhang_1/wrfdata/ncep2/pgb.anl.1990*.tar
yhrun -n 1 -N 1 -p bigdata ./ungrib.exe >& log.ungrib
另外,在最开始链接到气象数据时一直出现-bash: ./link_grib.csh: 错误,仔细研究后是权限不足,在XFTP6中手动修改了link_grib.csh文件的权限(右键->修改权限,发现原本只有读取权限,勾选三个执行->确认,可以成功链接!)
该部分结束生成两种6小时一个的气象文件,文件贯穿namelist.wps中设置的模拟时间,文件如下:
- metgrid.exe(融合地面数据与大气强迫场,建立适用于模拟区域的边界数据集,生成met文件);运行前要修改namelist 【牢记重新定位 geo_em 文件,不同年份,用不同的】
namelist.wps中的相关设置为:
&metgrid
fg_name = 'PGB','FLX'
io_form_metgrid = 2,
执行的命令:
yhrun -n 1 -N 1 -p bigdata ./metgrid.exe >& log.met
会给各个嵌套区都生成形如 met_em.d0x.yyyy-mm-dd_hh:mm:ss.nc 的数据文件,这些文件将会用来驱动WRF运行,如: