大氣科學(xué)WRF模式(Weather Research and Forecasting Model)
WRF案例
準(zhǔn)備一、創(chuàng)建DATA文件夾
準(zhǔn)備二、修改namelist.input內(nèi)容
第一步、運(yùn)行靜態(tài)地理數(shù)據(jù)
一、進(jìn)入WPS文件夾
二、修改namelist.wps
三、生成地理數(shù)據(jù)
第二步、處理氣象數(shù)據(jù)
一、將數(shù)據(jù)放入DATA文件夾中
二、鏈接數(shù)據(jù)
二、啟動(dòng)ungrib.exe
第三步、融合氣象和地理數(shù)據(jù)
第四步、修改WRF中的namelist.input
第五步、鏈接WPS到WRF
第六步、運(yùn)行WRF(時(shí)間較久,耐心等待)
WRF的預(yù)處理系統(tǒng)(WPS)
wps三個(gè)程序組成模塊的作用
(1)geogrid:確定模式區(qū)域并把靜態(tài)地理數(shù)據(jù)插值到格點(diǎn)。
(2)ungrib:從GRIB格式的數(shù)據(jù)中提取氣象要素場(chǎng)。
(3)metgrid:把提取出的氣象要素場(chǎng)水平插值到由geogrid確定的網(wǎng)格點(diǎn)上。
準(zhǔn)備一、創(chuàng)建DATA文件夾
在Build_WRF下創(chuàng)建一個(gè)DATA文件夾用于存放要用的數(shù)據(jù)。當(dāng)然,你也可以自己搞個(gè)名字,自己找個(gè)路徑,只要后面合并的時(shí)候能搞清楚Linux的路徑命令即可。
準(zhǔn)備二、修改namelist.input內(nèi)容
該文件路徑:em_real文件下,文件內(nèi)容:
&share
?wrf_core = 'ARW',
?max_dom = 1,
?start_date = '2016-01-01_00:00:00','2013-06-18_00:00:00','2013-07-24_00:00:00',
?end_date?? = '2016-01-01_18:00:00','2013-06-20_00:00:00','2013-07-24_18:00:00',
?interval_seconds = 21600
?io_form_geogrid = 2,
?/
?
?&geogrid
?parent_id? =? 1, 1, 2,
?parent_grid_ratio =?? 1, 3, 3,
?i_parent_start??? =?? 1, 66, 55,
?j_parent_start??? =?? 1, 59, 43,?
?e_we????????????? =? 155, 235, 460,
?e_sn????????????? =? 135, 182, 400,
?geog_data_res???? = 'modis_30s+10m','2m','3s',
?dx = 25000,
?dy = 25000,
?map_proj = 'lambert',
?ref_lat?? =? 26,
?ref_lon?? =? 80,
?truelat1? =? 30.0,
?truelat2? =? 60.0,
?stand_lon =? 88,
?geog_data_path = '/BIGDATA1/sysu_1/hzh/wrfchem3.9/geog'
?/
?
?&ungrib
?out_format = 'WPS',
?prefix = 'FILE',
?/
?
?&metgrid
?fg_name = 'FILE'
?io_form_metgrid = 2,
?/
?
——————————開始運(yùn)行wrf
第一步、運(yùn)行靜態(tài)地理數(shù)據(jù)
一、進(jìn)入WPS文件夾
cd WPS
二、修改namelist.wps
這里建議直接通過記事本去修改。
該文件位置:WPS\namelist.wps
注意:改動(dòng)的前提是大家按照之前帖子中的namelist.wps已做了覆蓋
改動(dòng)如下圖
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?
簡(jiǎn)單介紹一下:
其中max_dom = 1,代表嵌套層數(shù)為1。
后面的start_date和end_date代表了起始和結(jié)束時(shí)間,代表了2層,如果3層,后面繼續(xù)加。
照理說現(xiàn)在設(shè)置的是1,只需要改第一列的時(shí)間,但是本人是本著改了2個(gè)更保險(xiǎn)。
那為什么
start_date = ‘2020-07-27_00:00:00’
end_date = ‘2020-07-27_18:00:00’
一會(huì)兒跑的數(shù)據(jù)是從 '2020-07-27_00:00:00’到’2020-07-27_18:00:00’的。
修改完成后保存退出。
三、生成地理數(shù)據(jù)
在WPS文件下。
輸入命令:./geogrid.exe
當(dāng)出現(xiàn)下面文字代表成功
———生成geo…….nc文件
如果出現(xiàn)下面的內(nèi)容,需下載對(duì)應(yīng)的地理信息包。
成功后,輸入命令,檢查地理信息文件
ls -lah geo_em.d01.nc
第二步、處理氣象數(shù)據(jù)
一、將數(shù)據(jù)(fnl)數(shù)據(jù)放入DATA文件夾中
根據(jù)namelist將對(duì)應(yīng)的fnl數(shù)據(jù)放入DATA文件
?
二、鏈接數(shù)據(jù)
在WPS文件下。
輸入命令:
ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
?
./link_grib.csh ~/(fnl的路徑)/fnl
(./link_grib.csh ~/data//fnl_20160601_* ./)
——生成AAA/AAB/AAC/AAD文件
二、啟動(dòng)ungrib.exe
在WPS文件下。
輸入命令:
./ungrib.exe
成功則出現(xiàn):
——生成SST文件
檢查文件信息 輸入:
ls -lah FILE*
得到處理好的氣象數(shù)據(jù)
這些文件在WPS根目錄下。
第三步、融合氣象和地理數(shù)據(jù)
在WPS文件下,輸入命令
./metgrid.exe
——生成met文件
此時(shí),WRF的預(yù)處理部分(WPS部分)完成了。
第四步、修改WRF中的namelist.input
跳轉(zhuǎn)至em_real文件夾(后面的操作都是WRF\test\em_real下進(jìn)行的)
cd ~/Build_WRF/WRF/test/em_real/
下面修改namelist.input,文件位置:
\Build_WRF\WRF\test\em_real\namelist.input
修改的是起止時(shí)間以及需要運(yùn)行模擬幾個(gè)小時(shí)的數(shù)據(jù)。一層所以只需要改第一列。
第五步、鏈接WPS到WRF
在WRF\test\em_real下進(jìn)行的
ln -sf ~/Build_WRF/WPS/met_em* .
(ln -sf ../wps/met_em.d01)
./real.exe
——生成wrfinput_d01+rsl文件
查看文件:ls -alh wrfbdy_d01 wrfinput_d01
——————————————
若是運(yùn)行WRF-chem
此時(shí)打開ANTHRO/src,運(yùn)行化學(xué)模式,然后將結(jié)果拷貝至em_real文件
1.???? radm2_MADE.inp——生成排放清單的文件
2.???? 運(yùn)行 ./anthro_emis < radm2_MADE.inp
3.拷貝cp wrfchemi_* ~/zy/em_real/
第六步、運(yùn)行WRF(時(shí)間較久,耐心等待)
后面的操作都是em_real文件下運(yùn)行
./wrf.exe
成功后會(huì)出現(xiàn)
這樣,WRF算是跑完了,結(jié)果也成功了。
天河常用命令
yhi——查看可用節(jié)點(diǎn)
yhq——查看任務(wù)進(jìn)程
后臺(tái)提交(以wrf.sh文件運(yùn)行wrf.exe)
vi wrf.sh——修改此文件 ?:wq保存并退出
文件內(nèi)容:
yhrun -N 2 -n 48 ./wrf.exe
提交
yhbatch -N 2 -p bigdata ./wrf.exe
查看運(yùn)行狀況
tail -f rsl.error.0000
yhalloc——搶占節(jié)點(diǎn)
ssh cn3856(節(jié)點(diǎn)名)
yhbatch -N 1 -p bigdate ./wps.sh? (后臺(tái)運(yùn)行)