real.exe 报错 参考压力不是单调递减的，这往往是由非常高的地形造成的

yk7sn7

namelist.input
&time_control
run_days                            = 1,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017, 2017,
start_month                         = 02,   02,   02,
start_day                           = 01,   01,   01,
start_hour                          = 12,   12,   12,
end_year                            = 2017, 2017, 2017,
end_month                           = 02,   02,   02,
end_day                             = 02,   02,   02,
end_hour                            = 12,   12,   12,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,  60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 7200,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
/

&domains
time_step                           = 216,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 100,   232,   316,
e_sn                                = 82,    142,   250,
e_vert                              = 33,    33,    33,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 39200, 13066.667,  4355.556,
dy                                  = 39200, 13066.667,  4355.556,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     13,    90,
j_parent_start                      = 1,     8,     15,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
physics_suite                       = 'CONUS'
mp_physics                          = 3,    3,    3,
cu_physics                          = 1,    1,    0,
ra_lw_physics                       = 1,    1,    1,
ra_sw_physics                       = 1,    1,    1,
bl_pbl_physics                      = 1,    1,    1,
sf_sfclay_physics                   = 1,    1,    1,
sf_surface_physics                  = 5,    5,    5,
radt                                = 30,    30,    30,
bldt                                = 0,     0,     0,
cudt                                = 5,     5,     5,
icloud                              = 1,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
hybrid_opt                          = 2,
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                            = 3,
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,     
gwd_opt                             = 1,
/

&bdy_control
spec_bdy_width                      = 5,
specified                           = .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

报错文件的内容：
taskid: 0 hostname: ubuntu
module_io_quilt_old.F        2931 T
Ntasks in X            1 , ntasks in Y            1
*************************************
Configuring physics suite 'conus'


         mp_physics:      3*     3*     3*
         cu_physics:      1*     1*     0*
      ra_lw_physics:      1*     1*     1*
      ra_sw_physics:      1*     1*     1*
     bl_pbl_physics:      1*     1*     1*
  sf_sfclay_physics:      1*     1*     1*
sf_surface_physics:      5*     5*     5*


(* = option overrides suite setting)
*************************************
REAL_EM V4.1.2 PREPROCESSOR
*************************************
Parent domain
ids,ide,jds,jde            1         100           1          82
ims,ime,jms,jme           -4         105          -4          87
ips,ipe,jps,jpe            1         100           1          82
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             304825256  bytes allocated
d01 2017-02-01_12:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2017-02-01_12:00:00  Input data is acceptable to use: met_em.d01.2017-02-01_12:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-01_12:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d01 2017-02-01_12:00:00 Timing for input          0 s.
d01 2017-02-01_12:00:00          flag_soil_layers read from met_em file is  1
Max map factor in domain 1 =  1.04. Scale the dt in the model accordingly.
Using sfcprs3 to compute psfc
d01 2017-02-01_12:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d01 2017-02-01_12:00:00 forcing artificial silty clay loam at   47 points, out of   8019
d01 2017-02-01_12:00:00 Total post number of sea ice location changes (water to land) =     14
d01 2017-02-01_12:00:00 Timing for processing          0 s.
d01 2017-02-01_12:00:00 Timing for output          0 s.
d01 2017-02-01_12:00:00 Timing for loop #    1 =          0 s.
d01 2017-02-01_18:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2017-02-01_18:00:00  Input data is acceptable to use: met_em.d01.2017-02-01_18:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-01_18:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d01 2017-02-01_18:00:00 Timing for input          0 s.
d01 2017-02-01_18:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
d01 2017-02-01_18:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d01 2017-02-01_18:00:00 forcing artificial silty clay loam at   15 points, out of   8019
d01 2017-02-01_18:00:00 Timing for processing          0 s.
d01 2017-02-01_18:00:00 Timing for output          0 s.
d01 2017-02-01_18:00:00 Timing for loop #    2 =          0 s.
d01 2017-02-02_00:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2017-02-02_00:00:00  Input data is acceptable to use: met_em.d01.2017-02-02_00:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-02_00:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d01 2017-02-02_00:00:00 Timing for input          0 s.
d01 2017-02-02_00:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
d01 2017-02-02_00:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d01 2017-02-02_00:00:00 forcing artificial silty clay loam at   15 points, out of   8019
d01 2017-02-02_00:00:00 Timing for processing          0 s.
d01 2017-02-02_00:00:00 Timing for output          0 s.
d01 2017-02-02_00:00:00 Timing for loop #    3 =          0 s.
d01 2017-02-02_06:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2017-02-02_06:00:00  Input data is acceptable to use: met_em.d01.2017-02-02_06:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-02_06:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d01 2017-02-02_06:00:00 Timing for input          0 s.
d01 2017-02-02_06:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
d01 2017-02-02_06:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d01 2017-02-02_06:00:00 forcing artificial silty clay loam at   15 points, out of   8019
d01 2017-02-02_06:00:00 Timing for processing          0 s.
d01 2017-02-02_06:00:00 Timing for output          0 s.
d01 2017-02-02_06:00:00 Timing for loop #    4 =          0 s.
d01 2017-02-02_12:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2017-02-02_12:00:00  Input data is acceptable to use: met_em.d01.2017-02-02_12:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-02_12:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d01 2017-02-02_12:00:00 Timing for input          0 s.
d01 2017-02-02_12:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
d01 2017-02-02_12:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d01 2017-02-02_12:00:00 forcing artificial silty clay loam at   15 points, out of   8019
d01 2017-02-02_12:00:00 Timing for processing          0 s.
d01 2017-02-02_12:00:00 Timing for output          0 s.
d01 2017-02-02_12:00:00 Timing for loop #    5 =          0 s.
*************************************
Nesting domain
ids,ide,jds,jde            1         232           1         142
ims,ime,jms,jme           -4         237          -4         147
ips,ipe,jps,jpe            1         232           1         142
INTERMEDIATE domain
ids,ide,jds,jde           11          93           6          58
ims,ime,jms,jme            6          98           1          63
ips,ipe,jps,jpe            9          95           4          60
*************************************
d01 2017-02-02_12:00:00  alloc_space_field: domain            2 ,              54676188  bytes allocated
d01 2017-02-02_12:00:00  alloc_space_field: domain            2 ,            1094399432  bytes allocated
d02 2017-02-01_12:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d02 2017-02-01_12:00:00  Input data is acceptable to use: met_em.d02.2017-02-01_12:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-01_12:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d02 2017-02-01_12:00:00 Timing for input          0 s.
d02 2017-02-01_12:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
d02 2017-02-01_12:00:00 No average surface temperature for use with inland lakes
Assume non-CLM input
d02 2017-02-01_12:00:00 forcing artificial silty clay loam at  189 points, out of  32571
d02 2017-02-01_12:00:00 Total post number of sea ice location changes (water to land) =     27
d02 2017-02-01_12:00:00 Timing for processing          1 s.
d02 2017-02-01_12:00:00 Timing for output          0 s.
d02 2017-02-01_12:00:00 Timing for loop #    1 =          1 s.
*************************************
Nesting domain
ids,ide,jds,jde            1         316           1         250
ims,ime,jms,jme           -4         321          -4         255
ips,ipe,jps,jpe            1         316           1         250
INTERMEDIATE domain
ids,ide,jds,jde           88         198          13         101
ims,ime,jms,jme           83         203           8         106
ips,ipe,jps,jpe           86         200          11         103
*************************************
d02 2017-02-01_18:00:00  alloc_space_field: domain            3 ,             111788028  bytes allocated
d02 2017-02-01_18:00:00  alloc_space_field: domain            3 ,            2498590856  bytes allocated
d03 2017-02-01_12:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d03 2017-02-01_12:00:00  Input data is acceptable to use: met_em.d03.2017-02-01_12:00:00.nc
metgrid input_wrf.F first_date_input = 2017-02-01_12:00:00
metgrid input_wrf.F first_date_nml = 2017-02-01_12:00:00
d03 2017-02-01_12:00:00 Timing for input          0 s.
d03 2017-02-01_12:00:00          flag_soil_layers read from met_em file is  1
Using sfcprs3 to compute psfc
----- ERROR: The reference pressure is not monotonically decreasing
             This tends to be caused by very high topography
(i,j) =           50           1 , topography =    6867.82275      m
k =           12 , reference pressure =    38316.8906      Pa
k =           13 , reference pressure =    38381.2266      Pa
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1236
In the dynamics namelist record, reduce etac from   0.200000003
-------------------------------------------

application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

astiny

垂直层搞密点，或者etac改小点。

yk7sn7

astiny 发表于 2022-2-8 16:18
垂直层搞密点，或者etac改小点。

谢谢，大佬，我换了一下第三层嵌套区域的位置，就成功了。下次再遇到这种问题，会优先考虑您说的方法再试试

猫老师007

astiny 发表于 2022-2-8 16:18
垂直层搞密点，或者etac改小点。

请问垂直层加密只需要修改e_vert这个参数吗