运行和验证

操作步骤

  1. 下载WRFDA验证需要的观察结果文件,并且使用SFTP工具将其上传到“/path/to/WRFTEST/test2/”目录下(这里以2019年12月26日数据为例)。

  2. 执行以下命令,进入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,
     /

  3. 执行以下命令,开始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

  4. 作业完成后执行以下命令,查看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 ***