WPS运行---利用fnl数据做气象场

肉嘟嘟的科研笔记

已于 2024-06-03 22:20:11 修改

阅读量422

点赞数 5

分类专栏： WRFChem学习文章标签：笔记

于 2024-04-15 17:39:13 首次发布

本文链接：https://blog.csdn.net/weixin_64799803/article/details/134106436

版权

WRFChem学习专栏收录该内容

6 篇文章 0 订阅

订阅专栏

WPS内容参考

https://www.cnblogs.com/jiangleads/articles/14211982.html

1. fnl数据下载

https://rda.ucar.edu/datasets/ds083.2/

GRIB2 6 HOURLY FILES 2007.12.06 to current

Complete File List

选择自己需要的年份数据点进去（本文以2022年6月数据为例）

选择python脚本下载

将该脚本放在选定的文件夹下

python rda-download.py 运行该脚本即可下载需要的数据

2. 新建case文件夹运行WPS

mkdir case

cd case

cp -rf WPS/* . ## 根据你的文件夹位置设定

3. 根据选择的模拟时间设置namelist.wps share部分

&share
 wrf_core = 'ARW', 
 max_dom = 1, !单层嵌套
 start_date = '2022-06-01_00:00:00','2006-08-16_12:00:00',!模拟的开始时间
 end_date   = '2022-06-30_18:00:00','2006-08-16_12:00:00',!模拟的中止时间，与你下载的fnl时间一致
 interval_seconds = 21600 ! 时间间隔6h=6*3600
 io_form_geogrid = 2,! 输出格式为nc
/

4. geogrid 设置

&geogrid
 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  31,
 j_parent_start    =   1,  17,
 e_we              =  270,
 e_sn              =  200,  97,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the HGT_M, GREENFRAC, 
 ! and LU_INDEX/LANDUSEF fields have changed in WPS v3.8. The HGT_M field
 ! is now interpolated from 30-arc-second USGS GMTED2010, the GREENFRAC 
 ! field is interpolated from MODIS FPAR, and the LU_INDEX/LANDUSEF fields 
 ! are interpolated from 21-class MODIS.
 !
 ! To match the output given by the default namelist.wps in WPS v3.7.1, 
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'gtopo_10m+usgs_10m+nesdis_greenfrac+10m','gtopo_2m+usgs_2m+nesdis_greenfrac+2m',
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = 'default','default',
 dx = 27000,
 dy = 27000,
 map_proj = 'lambert',
 ref_lat   =  36,
 ref_lon   = 105,
 truelat1  =  36,
 truelat2  =  36,
 stand_lon = 105,
 geog_data_path = '/WPS_GEOG/'
/

运行WPS

新建一个case文件，复制WPS和WRF目录里的内容进行运行

mkdir case
cd case
mkdir WPS_case !新建文件夹复制WPS内容并运行
cd WPS_case
cp -rf ../../WPS/* . !复制WPS内容并运行
vi namelist.wps      !编辑wps
./geogrid.exe >& geogrid.log !运行geogrid并将运行记录输出到geogrid.log
./link_grib.csh fnl_data/fnl_2022*
ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
./ungrib.exe >& ungrib.log !运行ungrib并将运行记录输出到ungrib.log
./metgrid.exe  >& metgrid.log 行metgrid并将运行记录输出到metgrid.log

将WRF目录下的文件复制到新的文件夹下，并将上述WPS结果met文件链接到相应case文件夹中并修改对应的namelist.input部分运行real.exe


cp -rf ../WRFV3/run/* . !新建文件夹复制WRF内容
ln -sf ../WRFV3/main/*.exe . !链接WRF文件夹下的*.exe内容
ln -sf WPS_case/met_em.d01.2022-06-* .!链接met数据所在路径

ncdump -h met_em.d01.2022-06-01_00\:00\:00.nc | grep -i "num_metgrid_levels" !查看met数据的层数，方便后边修改wrf模拟部分的namelist.input
vi namelist.input
./real.exe !准备运行先运行real输出wrfinput和wrfbdy