ERA5再分析资料驱动WRFV4.6，温度湿度明显偏低

zhoulm177 · 发表于 6 小时前

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

使用ERA5再分析资料驱动WRFV4.6，模拟海南岛的天气，海南岛上2米温度偏低5度以上，湿度偏低40%以上。请帮忙分析。
namelist.wps

&share
wrf_core = 'ARW',
max_dom = 3,
start_date = 'STARTYEAR-STARTMONTH-STARTDAY_STARTHOUR:00:00','STARTYEAR-STARTMONTH-STARTDAY_STARTHOUR:00:00','STARTYEAR-STARTMONTH-STARTDAY_STARTHOUR:00:00',
end_date = 'ENDYEAR-ENDMONTH-ENDDAY_ENDHOUR:00:00','ENDYEAR-ENDMONTH-ENDDAY_ENDHOUR:00:00','ENDYEAR-ENDMONTH-ENDDAY_ENDHOUR:00:00',
interval_seconds = 3600
io_form_geogrid = 2,
opt_output_from_geogrid_path = './',
/
&geogrid
parent_id = 0, 1, 2,
parent_grid_ratio = 1, 5, 5,
i_parent_start = 1, 39, 29,
j_parent_start = 1, 39, 29,
e_we = 100, 111, 251,
e_sn = 100, 111, 231,
geog_data_res = 'default', '30s', '30s',
! geog_data_res = 'default','default','default',
dx = 25000, ! 25km
dy = 25000, ! 25km
map_proj = 'mercator',
ref_lat = 19.50,
ref_lon = 110.00,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = 110.0,
geog_data_path = '/media/E/geog_high_res_mandatory/WPS_GEOG'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE','SFILE'
io_form_metgrid = 2,
!constants_name = 'ERA5_INVARIANT'
/

复制代码

namelist.input

&time_control
run_days = 0,
run_hours = RUNHOURS,
run_minutes = 0,
run_seconds = 0,
start_year = STARTYEAR, STARTYEAR, STARTYEAR,
start_month = STARTMONTH, STARTMONTH, STARTMONTH,
start_day = STARTDAY, STARTDAY, STARTDAY,
start_hour = STARTHOUR, STARTHOUR, STARTHOUR,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = ENDYEAR, ENDYEAR, ENDYEAR,
end_month = ENDMONTH, ENDMONTH, ENDMONTH,
end_day = ENDDAY, ENDDAY, ENDDAY,
end_hour = ENDHOUR, ENDHOUR, ENDHOUR,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 3600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 30, 15,
frames_per_outfile = 12, 18, 12,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/
&domains
time_step = 150,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 100, 111, 251,
e_sn = 100, 111, 231,
e_vert = 40, 40, 40,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 25000, 5000, 1000,
dy = 25000, 5000, 1000,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 39, 29,
j_parent_start = 1, 39, 29,
parent_grid_ratio = 1, 5, 5,
parent_time_step_ratio = 1, 5, 5,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 4, 4, 4,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 15, 9, 3,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
!sf_ocean_physics = 0, 0, 0,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
usemonalb = .false.,
! num_land_cat = 21,
! isurban = 13,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 2,
km_opt = 4,
diff_6th_opt = 0,
diff_6th_factor = 0.10,
base_temp = 290.0,
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.02, 0.02, 0.02,
khdif = 0, 0, 0,
kvdif = 0,
epssm = 0.3,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 2,
scalar_adv_opt = 1, 2, 2,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&dfi_control
dfi_opt = 0
dfi_nfilter = 7
dfi_write_filtered_input = .false.
dfi_write_dfi_history = .false.
dfi_cutoff_seconds = 3600
dfi_time_dim = 1000
dfi_bckstop_year = 2000
dfi_bckstop_month = 01
dfi_bckstop_day = 24
dfi_bckstop_hour = 11
dfi_bckstop_minute = 00
dfi_bckstop_second = 00
dfi_fwdstop_year = 2000
dfi_fwdstop_month = 01
dfi_fwdstop_day = 24
dfi_fwdstop_hour = 12
dfi_fwdstop_minute = 30
dfi_fwdstop_second = 00
/

复制代码

运行脚本：

#运行的线程数
THREADS=14
echo ============================================= >>/usr/local/wrf/run.log
date>>/usr/local/wrf/run.log
echo 运行开始 >>/usr/local/wrf/run.log
#生成各个namelist文件
#####？？？？如何判断输入的参数是否符合格式？
STARTYEAR=${1:0:4}
STARTMONTH=${1:4:2}
STARTDAY=${1:6:2}
STARTHOUR=${1:8:2}
RUNHOURS=${2:0:2} #运行时长
RUNHOURS=${$RUNHOURS:-72} #如果RUNHOURS为null或者为空时，则RUNHOURS=72
temp1=$STARTYEAR-$STARTMONTH-$STARTDAY
temp2=$STARTHOUR:00:00
STARTTIMESTAMP=$(date +%s -d $temp1\ $temp2) #起始时间戳
echo 起始时间:${1:0:10}
echo 模式时长:$RUNHOURS
ENDTIMESTAMP=$((STARTTIMESTAMP+$(($RUNHOURS*60*60)))) #结束时间戳
ENDTIME=$(date +%Y%m%d%H -d "1970-01-01 UTC $ENDTIMESTAMP seconds"); #seconds不能省略
echo 结束时间:$ENDTIME
ENDYEAR=${ENDTIME:0:4}
ENDMONTH=${ENDTIME:4:2}
ENDDAY=${ENDTIME:6:2}
ENDHOUR=${ENDTIME:8:2}
cp -f /usr/local/wrf/namelistsERA5/20260417_25p5p1_v2/namelist.wps /usr/local/wrf/WPSV4.6.0/namelist.wps
cp -f /usr/local/wrf/namelistsERA5/20260417_25p5p1_v2/namelist.input /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/STARTYEAR/$STARTYEAR/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/STARTMONTH/$STARTMONTH/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/STARTDAY/$STARTDAY/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/STARTHOUR/$STARTHOUR/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/STARTYEAR/$STARTYEAR/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/STARTMONTH/$STARTMONTH/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/STARTDAY/$STARTDAY/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/STARTHOUR/$STARTHOUR/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/RUNHOURS/$RUNHOURS/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/ENDYEAR/$ENDYEAR/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/ENDMONTH/$ENDMONTH/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/ENDDAY/$ENDDAY/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/ENDHOUR/$ENDHOUR/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
sed -i "s/ENDYEAR/$ENDYEAR/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/ENDMONTH/$ENDMONTH/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/ENDDAY/$ENDDAY/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
sed -i "s/ENDHOUR/$ENDHOUR/g" /usr/local/wrf/WRFV4.6.1-GITEE/run/namelist.input
#read -p "按回车键继续..." # 显示提示并等待回车
cd /usr/local/wrf/WPSV4.6.0
#确定GEOGRID.TBL的指向
ln -sf /usr/local/wrf/WPSV4.6.0/geogrid/GEOGRID.TBL.ARW /usr/local/wrf/WPSV4.6.0/geogrid/GEOGRID.TBL
#清除中间文件
rm -f geo_em*
rm -f GRIBFILE.*
rm -f FILE*
rm -f SFILE*
rm -f PFILE*
rm -f met_em.*
rm -f metgrid.log*
#read -p "按回车键继续..." # 显示提示并等待回车
#运行geogrid.exe
./geogrid.exe |&tee rungeogridlog.txt
#链接合适的Vtable文件
ln -sf ./ungrib/Variable_Tables/Vtable.ECMWF Vtable
#链接等压面层GRIB资料，执行ungrib.exe
./link_grib.csh /media/E/data/ERA5/2026020900_48H_era5_forWRF/era5*_pressure.grib
./ungrib.exe |&tee runungriblog.txt # 输出FILE:YYYY-MM-DD_00
#read -p "按回车键继续..." # 显示提示并等待回车
#链接单层GRIB资料，改prefix为SFILE，重新ungrib.exe
./link_grib.csh /media/E/data/ERA5/2026020900_48H_era5_forWRF/era5*_single.grib
sed -i "s/prefix = 'FILE'/prefix = 'SFILE'/g" /usr/local/wrf/WPSV4.6.0/namelist.wps
./ungrib.exe |&tee runungriblog.txt # 输出SFILE:YYYY-MM-DD_00
#read -p "按回车键继续..." # 显示提示并等待回车
#运行metgrid.exe
./metgrid.exe |&tee runmetgridlog.txt
cd /usr/local/wrf/WRFV4.6.1-GITEE/run
rm -f wrfinput_d0? wrfbdy_d01
rm -f wrflowinp_d0? wrfvar_output
rm -f met_em.d0?.*
rm -f rsl.* rsl.error.* rsl.out.*
ln -sf /usr/local/wrf/WPSV4.6.0/met_em.* .
#read -p "按回车键继续..." # 显示提示并等待回车
#运行real.exe
date>>/usr/local/wrf/run.log
echo 开始运行real.exe >>/usr/local/wrf/run.log
#mpirun -np $THREADS ./real.exe |&tee runreallog.txt
./real.exe |&tee runreallog.txt
#read -p "按回车键继续..." # 显示提示并等待回车
#运行wrf.exe
date>>/usr/local/wrf/run.log
echo 开始运行wrf.exe >>/usr/local/wrf/run.log
mpirun -np $THREADS --bind-to core ./wrf.exe |&tee runwrflog.txt

复制代码

VTABLE为 Vtable.ECMWF

		自动登录	找回密码
密码			立即注册