运行和验证
操作步骤
- 下载WRFDA验证需要的观察结果文件，并且使用SFTP工具将其上传到“/path/to/WRFTEST/test2/”目录下（这里以2019年12月26日数据为例）。
- 执行以下命令，进入test2目录创建“namelist.input”文件，并且新增以下内容。
cd /path/to/WRFTEST/test2 vi namelist.input
&wrfvar1 var4d=false, print_detail_grad=false, / &wrfvar2 / &wrfvar3 ob_format=1, / &wrfvar4 use_gpsrefobs=false, / &wrfvar5 / &wrfvar6 max_ext_its=1, ntmax=50, orthonorm_gradient=true, / &wrfvar7 cv_options=3, / &wrfvar8 / &wrfvar9 / &wrfvar10 test_transforms=false, test_gradient=false, / &wrfvar11 / &wrfvar12 / &wrfvar13 / &wrfvar14 / &wrfvar15 / &wrfvar16 / &wrfvar17 / &wrfvar18 analysis_date="2019-12-26_00:00:00.0000", / &wrfvar19 / &wrfvar20 / &wrfvar21 time_window_min="2019-12-25_23:00:00.0000", / &wrfvar22 time_window_max="2019-12-26_01:00:00.0000", / &time_control run_days = 0, run_hours = 3, run_minutes = 0, run_seconds = 0, start_year = 2019, 2019, 2019, start_month = 12, 12, 12, start_day = 26, 26, 26, start_hour = 00, 00, 00, start_minute = 00, 00, 00, start_second = 00, 00, 00, end_year = 2019, 2019, 2019, end_month = 12, 12, 12, end_day = 26, 26, 26, end_hour = 12, 12, 12, end_minute = 00, 00, 00, end_second = 00, 00, 00, interval_seconds = 10800, input_from_file = .true.,.true.,.true., history_interval = 60, 60, 60, frames_per_outfile = 1, 1, 1, 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 = 50, time_step_fract_num = 0, time_step_fract_den = 1, max_dom = 1, e_we = 701, 691, 94, e_sn = 601, 703, 91, e_vert = 30, 30, 30, p_top_requested = 5000, num_metgrid_levels = 32, num_metgrid_soil_levels = 4, dx = 9000, 3000, 3333.33, dy = 9000, 3000, 3333.33, grid_id = 1, 2, 3, parent_id = 0, 1, 2, i_parent_start = 1, 309, 30, j_parent_start = 1, 238, 30, parent_grid_ratio = 1, 3, 3, parent_time_step_ratio = 1, 3, 3, feedback = 1, smooth_option = 0 / &physics mp_physics = 3, 3, 3, ra_lw_physics = 1, 1, 1, ra_sw_physics = 1, 1, 1, radt = 30, 30, 30, sf_sfclay_physics = 1, 1, 1, sf_surface_physics = 2, 2, 2, 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 = 3, num_soil_layers = 4, num_land_cat = 21, sf_urban_physics = 0, 0, 0, / &fdda / &dynamics w_damping = 0, diff_opt = 1, 1, 1, km_opt = 4, 4, 4, diff_6th_opt = 0, 0, 0, diff_6th_factor = 0.12, 0.12, 0.12, base_temp = 290. damp_opt = 0, zdamp = 5000., 5000., 5000., dampcoef = 0.2, 0.2, 0.2 khdif = 0, 0, 0, kvdif = 0, 0, 0, non_hydrostatic = .true., .true., .true., moist_adv_opt = 1, 1, 1, scalar_adv_opt = 1, 1, 1, / &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, /
- 执行以下命令，开始da_wrfvar.exe处理。
mpirun --allow-run-as-root -N 128 -hostfile $HOSTFILE -x OMP_NUM_THREADS=1 -x OPENBLAS_NUM_THREADS=1 -x PATH -x LD_LIBRARY_PATH -mca pml ucx -x UCX_NET_DEVICES=mlx5_0:1 -mca btl ^vader,tcp,openib,uct -x UCX_TLS=self,sm,rc --bind-to core ./da_wrfvar.exe
- 作业完成后执行以下命令，查看da_wrfvar.exe处理过程生成的日志“rsl.error.0000”来判断是否正常结束，如果出现“WRF-Var completed successfully”字段则表示正常结束。
less rsl.error.0000
taskid: 0 hostname: XA320V2-7 module_io_quilt_old.F 2931 T Namelist logging not found in namelist.input. Using registry defaults for variables in logging. Ntasks in X 8 , ntasks in Y 12 ************************************* Parent domain ids,ide,jds,jde 1 701 1 601 ims,ime,jms,jme -4 95 -4 57 ips,ipe,jps,jpe 1 88 1 50 ************************************* DYNAMICS OPTION: Eulerian Mass Coordinate alloc_space_field: domain 1 , 298744960 bytes allocated WRF NUMBER OF TILES FROM OMP_GET_MAX_THREADS = 1 Tile Strategy is not specified. Assuming 1D-Y WRF TILE 1 IS 1 IE 88 JS 1 JE 50 WRF NUMBER OF TILES = 1 *** WRF-Var completed successfully ***