WRF-Chem模

liwenjia

大家好，想咨询大家一个问题。我现在在用WRF-Chem进行CO2的浓度模拟，边界条件和初始条件用的是CT2022，所有都比较正常，生成的日志文件显示排放也读进去了，但是不知道为什么除了边界在变，中间的区域根本没有变化，我试过将部分地区排放增加与之前对比，发现没有变化，应该就是排放没读进去，想请问一下大家是什么原因啊，是我namelist设置的有问题吗

liwenjia

namelist设置如下
&time_control
run_days                            = 10,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2019, 2017, 2017, 2017,
start_month                         = 09,   07,   07,   07,
start_day                           = 01,   04,   04,   04,
start_hour                          = 00,   00,   00,   00,
start_minute                        = 00,   00,   00,   00,
start_second                        = 00,   00,   00,   00,
end_year                            = 2019, 2017, 2017, 2017,
end_month                           = 09,   07,   07,   07,
end_day                             = 11,   12,   12,   12,
end_hour                            = 00,   18,   18,   18,
end_minute                          = 00,   00,   00,   00,
end_second                          = 00,   00,   00,   00,

interval_seconds                    = 86400,
input_from_file                     = .true.,.true.,.true.,.true.,
history_interval                    = 60,  60,   60,  60,
frames_per_outfile                  = 1, 1, 24,
restart                             = .false.,
write_hist_at_0h_rst = .true.
restart_interval                    = 86400,

io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,


auxinput5_inname                    = 'wrfchemi_d<domain>_d<date>,   
auxinput5_interval_m                = 60, 720, 720,
io_form_auxinput5                   = 2,
frames_per_auxinput5                = 1000,  12, 12,

force_use_old_data                  = T,
debug_level                         = 0,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 80, 451, 571, 391,
e_sn                                = 80, 451, 551, 351,  
e_vert                              = 35, 35, 35, 35,
dx                                  = 12000, 3000, 600, 120,
dy                                  = 12000, 3000, 600, 120,
p_top_requested                     = 5000,
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
grid_id                             = 1,     2,     3,    4,
parent_id                           = 1,     1,     2,    3,
i_parent_start                      = 1,     299,   231,  236,
j_parent_start                      = 1,     141,   197,  205,
parent_grid_ratio                   = 1,     5,     5,    5,
parent_time_step_ratio              = 1,     5,     5,    5,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 10,    -1,
cu_physics                          = 5,    -1,
ra_lw_physics                       = 4,    -1,
ra_sw_physics                       = 4,    -1,
bl_pbl_physics                      = 2,    -1,
sf_sfclay_physics                   = 2,    -1,
sf_surface_physics                  = 2,    -1,
radt                                = 15,    15,
bldt                                = 0,     0,
cudt                                = 0,     0,
icloud                              = 1,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,
fractional_seaice                   = 1,
/

&fdda
/

&dynamics
w_damping                           = 1,
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,     
chem_adv_opt                        = 1,      1,      1,     
/

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

&grib2
/

&chem
kemit                               = 8,
chem_opt                            = 16, 16, 16, 0,
photdt                              = 30, 30, 30,
chemdt                              = 2., 0.6666, 0.3333,
!frames_per_emissfile                = 1, 1, 1,
io_style_emissions                  = 2,
emiss_inpt_opt                      = 16,  16,  16,
emiss_opt                           = 16,  16,  16,
chem_in_opt                         = 0,
phot_opt                            = 1,        1,    1,
gas_drydep_opt                      = 0,        1,    1,
aer_drydep_opt                      = 0,        1,    1,
bio_emiss_opt                       = 16,    16,    3,
vprm_opt                            = "VPRM_table_TROPICS",
dust_opt                            = 0,  
dmsemis_opt                         = 0,  
seas_opt                            = 0,  
gas_bc_opt                          = 16,        16,    16,
gas_ic_opt                          = 16,        16,    16,
aer_bc_opt                          = 1,        1,    1,
aer_ic_opt                          = 1,        1,    1,
gaschem_onoff                       = 1,        1,    1,
aerchem_onoff                       = 1,        1,    1,
wetscav_onoff                       = 0,        0,    0,
cldchem_onoff                       = 0,        0,    0,
vertmix_onoff                       = 1,        1,    1,
chem_conv_tr                        = 0,        0,    0,
biomass_burn_opt                    = 0,        0,    0,
plumerisefire_frq                   = 30,       30,   30,
aer_ra_feedback                     = 0,        0,  0,
have_bcs_chem                       = .true., .false., .false.,
ne_area                             = 100,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

liwenjia

rsl.out.0000部分如下
taskid: 0 hostname: node11
Quilting with   1 groups of   0 I/O tasks.
Ntasks in X            1 , ntasks in Y            1
  Domain # 1: dx = 12000.000 m
WRF V4.3.3 MODEL
No git found or not a git repository, git commit version not available.
*************************************
Parent domain
ids,ide,jds,jde            1          80           1          80
ims,ime,jms,jme           -4          85          -4          85
ips,ipe,jps,jpe            1          80           1          80
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             901944448  bytes allocated
  med_initialdata_input: calling input_input
   Input data is acceptable to use: wrfinput_d01
CURRENT DATE          = 2019-09-01_00:00:00
SIMULATION START DATE = 2019-09-01_00:00:00
Timing for processing wrfinput file (stream 0) for domain        1:    0.04561 elapsed seconds
Max map factor in domain 1 =  1.00. Scale the dt in the model accordingly.
D01: Time step                              =    60.0000000      (s)
D01: Grid Distance                          =    12.0000000      (km)
D01: Grid Distance Ratio dt/dx              =    5.00000000      (s/km)
D01: Ratio Including Maximum Map Factor     =    4.97684622      (s/km)
D01: NML defined reasonable_time_step_ratio =    6.00000000
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          41  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20  CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19  CATEGORIES
*********************************************************************
*             PROGRAM:WRF-Chem V4.3.3 MODEL
No git found or not a git repository, git commit version not available.
*                                                                   *
*    PLEASE REPORT ANY BUGS TO WRF-Chem HELP at                     *
*                                                                   *
*              wrfchemhelp.gsd@noaa.gov                             *
*                                                                   *
*********************************************************************
WARNING: Users interested in the GHG options should check the comments/references in header of module_ghg_fluxes
Warning: the VPRM parameters may need to be optimized depending on the season, year and region!
The parameters provided here should be used for testing purposes only!
  call madronich phot initialization
start_domain_em: numgas =    0
Timing for Writing wrfout_d01_2019-09-01_00:00:00 for domain        1:    0.30561 elapsed seconds
mediation_integrate: med_read_wrf_chem_emissions: Open file wrfchemi_00z_d01
HOURLY EMISSIONS UPDATE TIME        0.0       0.0
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2019-09-01_00:00:00
d01 2019-09-01_00:00:00  Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain        1:    0.01329 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     80 JS      1 JE     80
WRF NUMBER OF TILES =   1
d01 2019-09-01_00:00:00  ----------------------------------------
d01 2019-09-01_00:00:00  W-DAMPING  BEGINS AT W-COURANT NUMBER =    1.00000000
d01 2019-09-01_00:00:00  ----------------------------------------
Timing for main: time 2019-09-01_00:01:00 on domain   1:   10.89764 elapsed seconds
Timing for main: time 2019-09-01_00:02:00 on domain   1:    1.06758 elapsed seconds
Timing for main: time 2019-09-01_00:03:00 on domain   1:    1.11553 elapsed seconds
Timing for main: time 2019-09-01_00:04:00 on domain   1:    1.12049 elapsed seconds
Timing for main: time 2019-09-01_00:05:00 on domain   1:    1.12393 elapsed seconds
Timing for main: time 2019-09-01_00:06:00 on domain   1:    1.12714 elapsed seconds
Timing for main: time 2019-09-01_00:07:00 on domain   1:    1.12679 elapsed seconds
Timing for main: time 2019-09-01_00:08:00 on domain   1:    1.13178 elapsed seconds
Timing for main: time 2019-09-01_00:09:00 on domain   1:    1.13217 elapsed seconds
Timing for main: time 2019-09-01_00:10:00 on domain   1:    1.13837 elapsed seconds
Timing for main: time 2019-09-01_00:11:00 on domain   1:    1.13779 elapsed seconds
Timing for main: time 2019-09-01_00:12:00 on domain   1:    1.14209 elapsed seconds
Timing for main: time 2019-09-01_00:13:00 on domain   1:    1.13917 elapsed seconds
Timing for main: time 2019-09-01_00:14:00 on domain   1:    1.14438 elapsed seconds
Timing for main: time 2019-09-01_00:15:00 on domain   1:    1.14893 elapsed seconds
Timing for main: time 2019-09-01_00:16:00 on domain   1:   10.67142 elapsed seconds
Timing for main: time 2019-09-01_00:17:00 on domain   1:    1.14633 elapsed seconds
Timing for main: time 2019-09-01_00:18:00 on domain   1:    1.15267 elapsed seconds
Timing for main: time 2019-09-01_00:19:00 on domain   1:    1.15964 elapsed seconds
Timing for main: time 2019-09-01_00:20:00 on domain   1:    1.16045 elapsed seconds
Timing for main: time 2019-09-01_00:21:00 on domain   1:    1.16140 elapsed seconds
Timing for main: time 2019-09-01_00:22:00 on domain   1:    1.16382 elapsed seconds
Timing for main: time 2019-09-01_00:23:00 on domain   1:    1.16482 elapsed seconds
Timing for main: time 2019-09-01_00:24:00 on domain   1:    1.17205 elapsed seconds
Timing for main: time 2019-09-01_00:25:00 on domain   1:    1.17565 elapsed seconds
Timing for main: time 2019-09-01_00:26:00 on domain   1:    1.17367 elapsed seconds
Timing for main: time 2019-09-01_00:27:00 on domain   1:    1.17805 elapsed seconds
Timing for main: time 2019-09-01_00:28:00 on domain   1:    1.18237 elapsed seconds
Timing for main: time 2019-09-01_00:29:00 on domain   1:    1.18007 elapsed seconds
Timing for main: time 2019-09-01_00:30:00 on domain   1:    1.18156 elapsed seconds
Timing for main: time 2019-09-01_00:31:00 on domain   1:   10.61524 elapsed seconds
Timing for main: time 2019-09-01_00:32:00 on domain   1:    1.14017 elapsed seconds
Timing for main: time 2019-09-01_00:33:00 on domain   1:    1.14004 elapsed seconds
Timing for main: time 2019-09-01_00:34:00 on domain   1:    1.14214 elapsed seconds
Timing for main: time 2019-09-01_00:35:00 on domain   1:    1.13958 elapsed seconds
Timing for main: time 2019-09-01_00:36:00 on domain   1:    1.14344 elapsed seconds
Timing for main: time 2019-09-01_00:37:00 on domain   1:    1.14957 elapsed seconds
Timing for main: time 2019-09-01_00:38:00 on domain   1:    1.27551 elapsed seconds
Timing for main: time 2019-09-01_00:39:00 on domain   1:    1.25798 elapsed seconds
Timing for main: time 2019-09-01_00:40:00 on domain   1:    1.25902 elapsed seconds
Timing for main: time 2019-09-01_00:41:00 on domain   1:    1.26951 elapsed seconds
Timing for main: time 2019-09-01_00:42:00 on domain   1:    1.26800 elapsed seconds
Timing for main: time 2019-09-01_00:43:00 on domain   1:    1.27193 elapsed seconds
Timing for main: time 2019-09-01_00:44:00 on domain   1:    1.27551 elapsed seconds
Timing for main: time 2019-09-01_00:45:00 on domain   1:    1.28397 elapsed seconds
Timing for main: time 2019-09-01_00:46:00 on domain   1:   11.02122 elapsed seconds
Timing for main: time 2019-09-01_00:47:00 on domain   1:    1.27023 elapsed seconds
Timing for main: time 2019-09-01_00:48:00 on domain   1:    1.27551 elapsed seconds
Timing for main: time 2019-09-01_00:49:00 on domain   1:    1.27852 elapsed seconds
Timing for main: time 2019-09-01_00:50:00 on domain   1:    1.28233 elapsed seconds
Timing for main: time 2019-09-01_00:51:00 on domain   1:    1.28715 elapsed seconds
Timing for main: time 2019-09-01_00:52:00 on domain   1:    1.28255 elapsed seconds
Timing for main: time 2019-09-01_00:53:00 on domain   1:    1.29147 elapsed seconds
Timing for main: time 2019-09-01_00:54:00 on domain   1:    1.29209 elapsed seconds
Timing for main: time 2019-09-01_00:55:00 on domain   1:    1.28485 elapsed seconds
Timing for main: time 2019-09-01_00:56:00 on domain   1:    1.28950 elapsed seconds
Timing for main: time 2019-09-01_00:57:00 on domain   1:    1.22633 elapsed seconds
Timing for main: time 2019-09-01_00:58:00 on domain   1:    1.18230 elapsed seconds
Timing for main: time 2019-09-01_00:59:00 on domain   1:    1.18209 elapsed seconds
Timing for main: time 2019-09-01_01:00:00 on domain   1:    1.22015 elapsed seconds
Timing for Writing wrfout_d01_2019-09-01_01:00:00 for domain        1:    0.30917 elapsed seconds
Timing for main: time 2019-09-01_01:01:00 on domain   1:   10.91046 elapsed seconds
Timing for main: time 2019-09-01_01:02:00 on domain   1:    1.16360 elapsed seconds
Timing for main: time 2019-09-01_01:03:00 on domain   1:    1.16502 elapsed seconds
Timing for main: time 2019-09-01_01:04:00 on domain   1:    1.16834 elapsed seconds
Timing for main: time 2019-09-01_01:05:00 on domain   1:    1.17009 elapsed seconds
Timing for main: time 2019-09-01_01:06:00 on domain   1:    1.16566 elapsed seconds
Timing for main: time 2019-09-01_01:07:00 on domain   1:    1.16701 elapsed seconds
Timing for main: time 2019-09-01_01:08:00 on domain   1:    1.16748 elapsed seconds
Timing for main: time 2019-09-01_01:09:00 on domain   1:    1.16884 elapsed seconds
Timing for main: time 2019-09-01_01:10:00 on domain   1:    1.16663 elapsed seconds
Timing for main: time 2019-09-01_01:11:00 on domain   1:    1.17054 elapsed seconds
Timing for main: time 2019-09-01_01:12:00 on domain   1:    1.17185 elapsed seconds
Timing for main: time 2019-09-01_01:13:00 on domain   1:    1.17250 elapsed seconds
Timing for main: time 2019-09-01_01:14:00 on domain   1:    1.17334 elapsed seconds
Timing for main: time 2019-09-01_01:15:00 on domain   1:    1.17249 elapsed seconds
Timing for main: time 2019-09-01_01:16:00 on domain   1:   10.60634 elapsed seconds
Timing for main: time 2019-09-01_01:17:00 on domain   1:    1.17593 elapsed seconds
Timing for main: time 2019-09-01_01:18:00 on domain   1:    1.17366 elapsed seconds
Timing for main: time 2019-09-01_01:19:00 on domain   1:    1.17816 elapsed seconds
Timing for main: time 2019-09-01_01:20:00 on domain   1:    1.17588 elapsed seconds
Timing for main: time 2019-09-01_01:21:00 on domain   1:    1.17973 elapsed seconds
Timing for main: time 2019-09-01_01:22:00 on domain   1:    1.17723 elapsed seconds
Timing for main: time 2019-09-01_01:23:00 on domain   1:    1.17647 elapsed seconds
Timing for main: time 2019-09-01_01:24:00 on domain   1:    1.18301 elapsed seconds
Timing for main: time 2019-09-01_01:25:00 on domain   1:    1.18147 elapsed seconds
Timing for main: time 2019-09-01_01:26:00 on domain   1:    1.17948 elapsed seconds
Timing for main: time 2019-09-01_01:27:00 on domain   1:    1.18113 elapsed seconds
Timing for main: time 2019-09-01_01:28:00 on domain   1:    1.18374 elapsed seconds
Timing for main: time 2019-09-01_01:29:00 on domain   1:    1.18213 elapsed seconds
Timing for main: time 2019-09-01_01:30:00 on domain   1:    1.23077 elapsed seconds
Timing for main: time 2019-09-01_01:31:00 on domain   1:   10.61306 elapsed seconds
Timing for main: time 2019-09-01_01:32:00 on domain   1:    1.18867 elapsed seconds
Timing for main: time 2019-09-01_01:33:00 on domain   1:    1.18642 elapsed seconds
Timing for main: time 2019-09-01_01:34:00 on domain   1:    1.18827 elapsed seconds
Timing for main: time 2019-09-01_01:35:00 on domain   1:    1.18740 elapsed seconds
Timing for main: time 2019-09-01_01:36:00 on domain   1:    1.18529 elapsed seconds
Timing for main: time 2019-09-01_01:37:00 on domain   1:    1.19899 elapsed seconds
Timing for main: time 2019-09-01_01:38:00 on domain   1:    1.20506 elapsed seconds
Timing for main: time 2019-09-01_01:39:00 on domain   1:    1.20008 elapsed seconds
Timing for main: time 2019-09-01_01:40:00 on domain   1:    1.21745 elapsed seconds
Timing for main: time 2019-09-01_01:41:00 on domain   1:    1.29799 elapsed seconds
Timing for main: time 2019-09-01_01:42:00 on domain   1:    1.30171 elapsed seconds
Timing for main: time 2019-09-01_01:43:00 on domain   1:    1.29963 elapsed seconds
Timing for main: time 2019-09-01_01:44:00 on domain   1:    1.32122 elapsed seconds
Timing for main: time 2019-09-01_01:45:00 on domain   1:    1.32148 elapsed seconds
Timing for main: time 2019-09-01_01:46:00 on domain   1:   10.96227 elapsed seconds
Timing for main: time 2019-09-01_01:47:00 on domain   1:    1.30514 elapsed seconds
Timing for main: time 2019-09-01_01:48:00 on domain   1:    1.30184 elapsed seconds
Timing for main: time 2019-09-01_01:49:00 on domain   1:    1.30325 elapsed seconds
Timing for main: time 2019-09-01_01:50:00 on domain   1:    1.30303 elapsed seconds
Timing for main: time 2019-09-01_01:51:00 on domain   1:    1.30823 elapsed seconds
Timing for main: time 2019-09-01_01:52:00 on domain   1:    1.30995 elapsed seconds
Timing for main: time 2019-09-01_01:53:00 on domain   1:    1.33068 elapsed seconds
Timing for main: time 2019-09-01_01:54:00 on domain   1:    1.31342 elapsed seconds
Timing for main: time 2019-09-01_01:55:00 on domain   1:    1.31491 elapsed seconds
Timing for main: time 2019-09-01_01:56:00 on domain   1:    1.31983 elapsed seconds
Timing for main: time 2019-09-01_01:57:00 on domain   1:    1.32208 elapsed seconds
Timing for main: time 2019-09-01_01:58:00 on domain   1:    1.32900 elapsed seconds
Timing for main: time 2019-09-01_01:59:00 on domain   1:    1.32734 elapsed seconds
Timing for main: time 2019-09-01_02:00:00 on domain   1:    1.36230 elapsed seconds
Timing for Writing wrfout_d01_2019-09-01_02:00:00 for domain        1:    0.31567 elapsed seconds
Timing for main: time 2019-09-01_02:01:00 on domain   1:   11.31574 elapsed seconds
Timing for main: time 2019-09-01_02:02:00 on domain   1:    1.32074 elapsed seconds
Timing for main: time 2019-09-01_02:03:00 on domain   1:    1.32928 elapsed seconds
Timing for main: time 2019-09-01_02:04:00 on domain   1:    1.32536 elapsed seconds
Timing for main: time 2019-09-01_02:05:00 on domain   1:    1.32726 elapsed seconds
Timing for main: time 2019-09-01_02:06:00 on domain   1:    1.32446 elapsed seconds
Timing for main: time 2019-09-01_02:07:00 on domain   1:    1.29825 elapsed seconds
Timing for main: time 2019-09-01_02:08:00 on domain   1:    1.25635 elapsed seconds
Timing for main: time 2019-09-01_02:09:00 on domain   1:    1.22087 elapsed seconds
Timing for main: time 2019-09-01_02:10:00 on domain   1:    1.22162 elapsed seconds
Timing for main: time 2019-09-01_02:11:00 on domain   1:    1.21900 elapsed seconds
Timing for main: time 2019-09-01_02:12:00 on domain   1:    1.22125 elapsed seconds
Timing for main: time 2019-09-01_02:13:00 on domain   1:    1.22183 elapsed seconds
Timing for main: time 2019-09-01_02:14:00 on domain   1:    1.22443 elapsed seconds
Timing for main: time 2019-09-01_02:15:00 on domain   1:    1.25006 elapsed seconds
Timing for main: time 2019-09-01_02:16:00 on domain   1:   10.69061 elapsed seconds
Timing for main: time 2019-09-01_02:17:00 on domain   1:    1.22428 elapsed seconds
Timing for main: time 2019-09-01_02:18:00 on domain   1:    1.22879 elapsed seconds
Timing for main: time 2019-09-01_02:19:00 on domain   1:    1.22716 elapsed seconds
Timing for main: time 2019-09-01_02:20:00 on domain   1:    1.23067 elapsed seconds
Timing for main: time 2019-09-01_02:21:00 on domain   1:    1.23143 elapsed seconds
Timing for main: time 2019-09-01_02:22:00 on domain   1:    1.22918 elapsed seconds
Timing for main: time 2019-09-01_02:23:00 on domain   1:    1.23200 elapsed seconds
Timing for main: time 2019-09-01_02:24:00 on domain   1:    1.23225 elapsed seconds
Timing for main: time 2019-09-01_02:25:00 on domain   1:    1.23180 elapsed seconds
Timing for main: time 2019-09-01_02:26:00 on domain   1:    1.23452 elapsed seconds
Timing for main: time 2019-09-01_02:27:00 on domain   1:    1.23333 elapsed seconds
Timing for main: time 2019-09-01_02:28:00 on domain   1:    1.23487 elapsed seconds
Timing for main: time 2019-09-01_02:29:00 on domain   1:    1.23170 elapsed seconds
Timing for main: time 2019-09-01_02:30:00 on domain   1:    1.28092 elapsed seconds
Timing for main: time 2019-09-01_02:31:00 on domain   1:   10.70489 elapsed seconds
Timing for main: time 2019-09-01_02:32:00 on domain   1:    1.23812 elapsed seconds
Timing for main: time 2019-09-01_02:33:00 on domain   1:    1.23837 elapsed seconds
Timing for main: time 2019-09-01_02:34:00 on domain   1:    1.24218 elapsed seconds
Timing for main: time 2019-09-01_02:35:00 on domain   1:    1.24013 elapsed seconds
Timing for main: time 2019-09-01_02:36:00 on domain   1:    1.24641 elapsed seconds
Timing for main: time 2019-09-01_02:37:00 on domain   1:    1.23694 elapsed seconds
Timing for main: time 2019-09-01_02:38:00 on domain   1:    1.24253 elapsed seconds
Timing for main: time 2019-09-01_02:39:00 on domain   1:    1.24703 elapsed seconds
Timing for main: time 2019-09-01_02:40:00 on domain   1:    1.24492 elapsed seconds
Timing for main: time 2019-09-01_02:41:00 on domain   1:    1.24504 elapsed seconds
Timing for main: time 2019-09-01_02:42:00 on domain   1:    1.24893 elapsed seconds
Timing for main: time 2019-09-01_02:43:00 on domain   1:    1.25067 elapsed seconds
Timing for main: time 2019-09-01_02:44:00 on domain   1:    1.25018 elapsed seconds
Timing for main: time 2019-09-01_02:45:00 on domain   1:    1.24893 elapsed seconds
Timing for main: time 2019-09-01_02:46:00 on domain   1:   10.73055 elapsed seconds
Timing for main: time 2019-09-01_02:47:00 on domain   1:    1.24412 elapsed seconds
Timing for main: time 2019-09-01_02:48:00 on domain   1:    1.25058 elapsed seconds
Timing for main: time 2019-09-01_02:49:00 on domain   1:    1.25061 elapsed seconds
Timing for main: time 2019-09-01_02:50:00 on domain   1:    1.25706 elapsed seconds
Timing for main: time 2019-09-01_02:51:00 on domain   1:    1.25132 elapsed seconds
Timing for main: time 2019-09-01_02:52:00 on domain   1:    1.25408 elapsed seconds
Timing for main: time 2019-09-01_02:53:00 on domain   1:    1.25327 elapsed seconds
Timing for main: time 2019-09-01_02:54:00 on domain   1:    1.25969 elapsed seconds
Timing for main: time 2019-09-01_02:55:00 on domain   1:    1.25727 elapsed seconds
Timing for main: time 2019-09-01_02:56:00 on domain   1:    1.26081 elapsed seconds
Timing for main: time 2019-09-01_02:57:00 on domain   1:    1.26228 elapsed seconds
Timing for main: time 2019-09-01_02:58:00 on domain   1:    1.26114 elapsed seconds
Timing for main: time 2019-09-01_02:59:00 on domain   1:    1.25900 elapsed seconds
Timing for main: time 2019-09-01_03:00:00 on domain   1:    1.30690 elapsed seconds
Timing for Writing wrfout_d01_2019-09-01_03:00:00 for domain        1:    0.30882 elapsed seconds
......
......
Timing for main: time 2019-09-10_23:55:00 on domain   1:    1.34217 elapsed seconds
Timing for main: time 2019-09-10_23:56:00 on domain   1:    1.33947 elapsed seconds
Timing for main: time 2019-09-10_23:57:00 on domain   1:    1.34104 elapsed seconds
Timing for main: time 2019-09-10_23:58:00 on domain   1:    1.34013 elapsed seconds
Timing for main: time 2019-09-10_23:59:00 on domain   1:    1.34154 elapsed seconds
Timing for main: time 2019-09-11_00:00:00 on domain   1:    1.38567 elapsed seconds
Timing for Writing wrfout_d01_2019-09-11_00:00:00 for domain        1:    0.30416 elapsed seconds
d01 2019-09-11_00:00:00 wrf: SUCCESS COMPLETE WRF

SolaProject

检查一下排放文件的字段是否齐全？这是我的文件的字段。FieldType（这个是整数就行，用ncdump看
会多个LL）之类的要有。

dimensions:
        Time = UNLIMITED ; // (12 currently)
        emissions_zdim = 1 ;
        south_north = 105 ;
        west_east = 105 ;
        DateStrLen = 19 ;
variables:
        float E_CO2(Time, emissions_zdim, south_north, west_east) ;
                E_CO2:FieldType = 104LL ;
                E_CO2:MemoryOrder = "XYZ" ;
                E_CO2:description = "EMISSIONS" ;
                E_CO2:units = "mol km^-2 hr^-1" ;
                E_CO2:stagger = "" ;
                E_CO2:coordinates = "XLONG XLAT" ;
                E_CO2:cell_methods = "Time: mean" ;
        char Times(Time, DateStrLen) ;
                Times:cell_methods = "Time: mean" ;

liwenjia

SolaProject 发表于 2024-4-16 10:38
检查一下排放文件的字段是否齐全？这是我的文件的字段。FieldType（这个是整数就行，用ncdump看
会多个LL ...

dimensions:
        Time = UNLIMITED ; // (241 currently)
        emissions_zdim = 8 ;
        south_north = 79 ;
        west_east = 79 ;
        DateStrLen = 19 ;
variables:
        char Times(Time, DateStrLen) ;
        float XLONG(south_north, west_east) ;
        float XLAT(south_north, west_east) ;
        float CO2_ANT(Time, emissions_zdim, south_north, west_east) ;
                CO2_ANT:FieldType = 104 ;
                CO2_ANT:MemoryOrder = "XYZ" ;
                CO2_ANT:description = "EMISSIONS" ;
                CO2_ANT:units = "mol km^-2 hr^-1" ;
                CO2_ANT:stagger = "Z" ;
                CO2_ANT:coordinates = "XLONG XLAT" ;
        float CO2_BIO(Time, emissions_zdim, south_north, west_east) ;
                CO2_BIO:FieldType = 104 ;
                CO2_BIO:MemoryOrder = "XYZ" ;
                CO2_BIO:description = "EMISSIONS" ;
                CO2_BIO:units = "mol km^-2 hr^-1" ;
                CO2_BIO:stagger = "Z" ;
                CO2_BIO:coordinates = "XLONG XLAT" ;
        float CO2_OCE(Time, emissions_zdim, south_north, west_east) ;
                CO2_OCE:FieldType = 104 ;
                CO2_OCE:MemoryOrder = "XYZ" ;
                CO2_OCE:description = "EMISSIONS" ;
                CO2_OCE:units = "mol km^-2 hr^-1" ;
                CO2_OCE:stagger = "Z" ;
                CO2_OCE:coordinates = "XLONG XLAT" ;
我模拟的是0901-0911，您看我的有问题吗

SolaProject

liwenjia 发表于 2024-4-16 14:55
dimensions:
        Time = UNLIMITED ; // (241 currently)
        emissions_zdim = 8 ;

看上去你给定了模拟时段的所有排放，一开始没注意到你的参数化方案中 io_style_emissions 选择的是 2。但是你排放文件的文件名好像是 wrfchemi_00z_d01，这是对应参数化方案 io_style_emissions 为 1 时的排放文件。

在 io_style_emissions 选择的是 2 时，排放文件的命名应当是 wrfchemi_d<domain_id>_<date/time>，例如 wrfchemi_d01_2019-09-01_00:00:00，每个文件里仅保留一个时次的排放。

在你的运行日志里可以看到，在每个小时开始时并没有读取排放数据。正常运行读取排放会有日志如下：

Timing for main: time 2019-05-26_12:59:54 on domain   1:    0.08309 elapsed seconds
Timing for Writing wrfout_d01_2019-05-26_13:00:00 for domain        1:    2.38604 elapsed seconds
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2019-05-26_13:00:00
mediation_integrate: med_read_wrf_chem_emissions: Open file wrfchemi_d01_2019-05-26_13:00:00
d01 2019-05-26_13:00:00  Error trying to read metadata
d01 2019-05-26_13:00:00  Input data is acceptable to use: wrfchemi_d01_2019-05-26_13:00:00
d01 2019-05-26_13:00:00  Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain        1:    2.29157 elapsed seconds
Timing for main: time 2019-05-26_13:00:06 on domain   1:    6.09932 elapsed seconds

liwenjia

SolaProject 发表于 2024-4-16 20:35
看上去你给定了模拟时段的所有排放，一开始没注意到你的参数化方案中 io_style_emissions 选择的是 2。但 ...

您好 是这样的，我之前是用的io_style_emission=2,我的排放文件也只有一个，命名为wrfchemi_d01_2019-09-01_00:00:00,这个日志是我用meic数据生成wrfchemi_00z_d01,wrfchemi_12z_d01这两个文件再运行生成的日志，可能是我上传错了，应该不是这个的问题，您还有别的建议吗，感谢您的回复

SolaProject

liwenjia 发表于 2024-4-16 22:42
您好 是这样的，我之前是用的io_style_emission=2,我的排放文件也只有一个，命名为wrfchemi_d01_2019-09- ...

在 io_style_emissions 选择的是 2 时，排放文件的命名应当是 wrfchemi_d<domain_id>_<date/time>，例如 wrfchemi_d01_2019-09-01_00:00:00，每个文件里仅保留一个时次的排放。

还是这个问题吧，有多少个小时就生成多少个小时的排放文件，然后都放在 wrf.exe 运行路径下。无法将所有时次的排放都放在一个文件里，WRF-Chem源码里这一块是写死的（或者可能也可以将一个包含所有时次的排放文件链接到各个小时，理论上可行，不过我没有这么试过），排放文件名、文件中 Times 对应的值都必须正确（程序会根据 Times 的时间检索数据印象）。

SolaProject

liwenjia 发表于 2024-4-16 14:55
dimensions:
        Time = UNLIMITED ; // (241 currently)
        emissions_zdim = 8 ;

补充说明，这里的排放变量应当是 E_CO2，而不是 CO2_ANT（这是模拟输出的结果）