Dear all,
Good day, I am a beginner at using wrf chem, I am trying to run wrf chem 3.9.1 for the analysis of pm2.5. Now in order to get the proper emission data, I tried using prep chem src. The initial compilation of prep_chem src was successful and i was able to get the prep_chem executable. I edited the prep_chem Input file for mercator projection and tried executing the prep_chem src to get the binary files for wrf emission input. At this point, it showing that there is some kind of discrepancy or mismatch I guess in the prep chem namelist as the error is showing as
Opening prep_chem_sources.inp file
At line 99 of file prep_chem_sources.F90 (unit = 5, file = 'prep_chem_sources.inp')
Fortran runtime error: Cannot match namelist object name .
The prep_chem src namelist file for my domain configuration is as follows
$RP_INPUT
!#################################################################
! CCATT-BRAMS/MCGA-CPTEC/WRF-Chem/FIM-Chem emission models CPTEC/INPE
! version 1.5: Mar 2015
!#################################################################
!---------------- grid_type of the grid output
grid_type= 'mercator',
rams_anal_prefix = '../ANL/OPQUE',
!---------------- date of emission
ihour=00,
iday=15,
imon=10,
iyear=2008,
!---------------- select the sources datasets to be used
use_retro=0, ! 1 = yes, 0 = not
retro_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/RETRO/anthro',
use_edgar = 2, ! 0 - not,
! 1 - Version 3,
! 2 - Version 4 for some species
! 3 - Version HTAP
edgar_data_dir='/home/roms/BUILD_WRF/Global_emissions/EDGAR/anthro',
use_gocart=0,
gocart_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GOCART/emissions',
use_streets =0,
streets_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/STREETS',
use_seac4rs =0,
seac4rs_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/SEAC4RS',
use_fwbawb =0,
fwbawb_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/Emissions_Yevich_Logan',
use_bioge =2, ! 1 - geia, 2 - megan
! ######
! # BIOGENIC = 1
bioge_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/biogenic_emissions',
! # MEGAN = 2
! ######
!bioge_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/MEGAN/2000',
! ######
use_gfedv2=0,
gfedv2_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GFEDv2-8days',
use_bbem=1,
use_bbem_plumerise=0,
!--------------------------------------------------------------------------------------------------
!---------------- if the merging of gfedv2 with bbem is desired (=1, yes, 0 = no)
merge_GFEDv2_bbem =0,
!---------------- Fire product for BBBEM/BBBEM-plumerise emission models
bbem_wfabba_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/GOES/f',
bbem_modis_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/MODIS/Fires',
bbem_inpe_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/DSA/Focos',
bbem_extra_data_dir ='NONE',
!---------------- veg type data set (dir + prefix)
veg_type_data_dir ='/home/roms/BUILD_WRF/Global_emissions/surface_data/GL_IGBP_MODIS_INPE/MODIS',
!---------------- vcf type data set (dir + prefix)
use_vcf = 0,
vcf_type_data_dir ='/home/roms/BUILD_WRF/Global_emissions/surface_data/VCF/data_out/2005/VCF',
!---------------- olson data set (dir + prefix)
olson_data_dir ='/home/roms/BUILD_WRF/Global_emission/Emission_data/OLSON2/OLSON',
!---------------- carbon density data set (dir + prefix)
carbon_density_data_dir='/home/roms/BUILD_WRF/Global_emissions/surface_data/GL_OGE_INPE/OGE',
fuel_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/Carbon_density_Saatchi/amazon_biomass_final.gra',
!---------------- gocart background
use_gocart_bg=1,
gocart_bg_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GOCART',
!---------------- volcanoes emissions
use_volcanoes =0,
volcano_index =1143, !REDOUBT
use_these_values='NONE',
! define a text file for using external values for INJ_HEIGHT, DURATION,
! MASS ASH (units are meters - seconds - kilograms) and the format for
begin_eruption='201303280000', !begin time UTC of eruption YYYYMMDDhhmm
!---------------- degassing volcanoes emissions
use_degass_volcanoes =0,
degass_volc_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/VOLC_SO2',
!---------------- user specific emissions directory
!---------------- Update for South America megacities
user_data_dir='NONE',
!--------------------------------------------------------------------------------------------------
pond=1, ! mad/mfa 0 -> molar mass weighted
! 1 -> Reactivity weighted
!---------------- for grid type 'll' or 'gg' only
grid_resolucao_lon=1.0,
grid_resolucao_lat=1.0,
nlat=320, ! if gg (only global grid)
lon_beg = -180., ! (-180.:+180.) long-begin of the output file
lat_beg = -90., ! ( -90.:+90. ) lat -begin of the output file
delta_lon = 360, ! total long extension of the domain (360 for global)
delta_lat = 180, ! total lat extension of the domain (180 for global)
!---------------- For regional grids (polar or lambert)
NGRIDS = 3, ! Number of grids to run
NNXP = 100,235,163, ! Number of x gridpoints
NNYP = 79,193,145, ! Number of y gridpoints
NXTNEST = 1,1,2, ! Grid number which is the next coarser grid
DELTAX = 26700.,
DELTAY = 26700., ! X and Y grid spacing
! Nest ratios between this grid and the next coarser grid.
NSTRATX = 1,3,3, ! x-direction
NSTRATY = 1,3,3, ! y-direction
NINEST = 1,14,99, ! Grid point on the next coarser
NJNEST = 1,10,71, ! nest where the lower southwest
! corner of this nest will start.
! If NINEST or NJNEST = 0, use CENTLAT/LON
POLELAT = 28.184., !-89.99, ! If polar, latitude/longitude of pole point
POLELON = 73.93., ! If lambert, lat/lon of grid origin (x=y=0.)
STDLAT1 = 28.184., ! If polar for BRAMS, use 90.0 in STDLAT2
STDLAT2 = 0., ! If lambert, standard latitudes of projection
!(truelat2/truelat1 from namelist.wps, STDLAT1 < STDLAT2)
! If mercator STDLAT1 = 1st true latitude
CENTLAT = 28.184.,!-89.99, -23., 27.5, 27.5,
CENTLON = 73.93., -46.,-80.5, -80.5,
!---------------- model output domain for each grid (only set up for rams)
lati = -90., -90., -90.,
latf = +90., +90., +90.,
loni = -180., -180., -180.,
lonf = 180., 180., 180.,
!---------------- project rams grid (polar sterogr) to lat/lon: 'YES' or 'NOT'
proj_to_ll='YES',
!---------------- output file prefix (may include directory other than the current)
chem_out_prefix = 'emissions',
chem_out_format = 'vfm',
!---------------- convert to WRF/CHEM (yes,not)
special_output_to_wrf = 'YES',
$END
I am not able to get what exactly the problem is here, please help me with this issue.
Good day, I am a beginner at using wrf chem, I am trying to run wrf chem 3.9.1 for the analysis of pm2.5. Now in order to get the proper emission data, I tried using prep chem src. The initial compilation of prep_chem src was successful and i was able to get the prep_chem executable. I edited the prep_chem Input file for mercator projection and tried executing the prep_chem src to get the binary files for wrf emission input. At this point, it showing that there is some kind of discrepancy or mismatch I guess in the prep chem namelist as the error is showing as
Opening prep_chem_sources.inp file
At line 99 of file prep_chem_sources.F90 (unit = 5, file = 'prep_chem_sources.inp')
Fortran runtime error: Cannot match namelist object name .
The prep_chem src namelist file for my domain configuration is as follows
$RP_INPUT
!#################################################################
! CCATT-BRAMS/MCGA-CPTEC/WRF-Chem/FIM-Chem emission models CPTEC/INPE
! version 1.5: Mar 2015
!#################################################################
!---------------- grid_type of the grid output
grid_type= 'mercator',
rams_anal_prefix = '../ANL/OPQUE',
!---------------- date of emission
ihour=00,
iday=15,
imon=10,
iyear=2008,
!---------------- select the sources datasets to be used
use_retro=0, ! 1 = yes, 0 = not
retro_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/RETRO/anthro',
use_edgar = 2, ! 0 - not,
! 1 - Version 3,
! 2 - Version 4 for some species
! 3 - Version HTAP
edgar_data_dir='/home/roms/BUILD_WRF/Global_emissions/EDGAR/anthro',
use_gocart=0,
gocart_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GOCART/emissions',
use_streets =0,
streets_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/STREETS',
use_seac4rs =0,
seac4rs_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/SEAC4RS',
use_fwbawb =0,
fwbawb_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/Emissions_Yevich_Logan',
use_bioge =2, ! 1 - geia, 2 - megan
! ######
! # BIOGENIC = 1
bioge_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/biogenic_emissions',
! # MEGAN = 2
! ######
!bioge_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/MEGAN/2000',
! ######
use_gfedv2=0,
gfedv2_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GFEDv2-8days',
use_bbem=1,
use_bbem_plumerise=0,
!--------------------------------------------------------------------------------------------------
!---------------- if the merging of gfedv2 with bbem is desired (=1, yes, 0 = no)
merge_GFEDv2_bbem =0,
!---------------- Fire product for BBBEM/BBBEM-plumerise emission models
bbem_wfabba_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/GOES/f',
bbem_modis_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/MODIS/Fires',
bbem_inpe_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/fires_data/DSA/Focos',
bbem_extra_data_dir ='NONE',
!---------------- veg type data set (dir + prefix)
veg_type_data_dir ='/home/roms/BUILD_WRF/Global_emissions/surface_data/GL_IGBP_MODIS_INPE/MODIS',
!---------------- vcf type data set (dir + prefix)
use_vcf = 0,
vcf_type_data_dir ='/home/roms/BUILD_WRF/Global_emissions/surface_data/VCF/data_out/2005/VCF',
!---------------- olson data set (dir + prefix)
olson_data_dir ='/home/roms/BUILD_WRF/Global_emission/Emission_data/OLSON2/OLSON',
!---------------- carbon density data set (dir + prefix)
carbon_density_data_dir='/home/roms/BUILD_WRF/Global_emissions/surface_data/GL_OGE_INPE/OGE',
fuel_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/Carbon_density_Saatchi/amazon_biomass_final.gra',
!---------------- gocart background
use_gocart_bg=1,
gocart_bg_data_dir='/home/roms/BUILD_WRF/Global_emissions/Emission_data/GOCART',
!---------------- volcanoes emissions
use_volcanoes =0,
volcano_index =1143, !REDOUBT
use_these_values='NONE',
! define a text file for using external values for INJ_HEIGHT, DURATION,
! MASS ASH (units are meters - seconds - kilograms) and the format for
begin_eruption='201303280000', !begin time UTC of eruption YYYYMMDDhhmm
!---------------- degassing volcanoes emissions
use_degass_volcanoes =0,
degass_volc_data_dir ='/home/roms/BUILD_WRF/Global_emissions/Emission_data/VOLC_SO2',
!---------------- user specific emissions directory
!---------------- Update for South America megacities
user_data_dir='NONE',
!--------------------------------------------------------------------------------------------------
pond=1, ! mad/mfa 0 -> molar mass weighted
! 1 -> Reactivity weighted
!---------------- for grid type 'll' or 'gg' only
grid_resolucao_lon=1.0,
grid_resolucao_lat=1.0,
nlat=320, ! if gg (only global grid)
lon_beg = -180., ! (-180.:+180.) long-begin of the output file
lat_beg = -90., ! ( -90.:+90. ) lat -begin of the output file
delta_lon = 360, ! total long extension of the domain (360 for global)
delta_lat = 180, ! total lat extension of the domain (180 for global)
!---------------- For regional grids (polar or lambert)
NGRIDS = 3, ! Number of grids to run
NNXP = 100,235,163, ! Number of x gridpoints
NNYP = 79,193,145, ! Number of y gridpoints
NXTNEST = 1,1,2, ! Grid number which is the next coarser grid
DELTAX = 26700.,
DELTAY = 26700., ! X and Y grid spacing
! Nest ratios between this grid and the next coarser grid.
NSTRATX = 1,3,3, ! x-direction
NSTRATY = 1,3,3, ! y-direction
NINEST = 1,14,99, ! Grid point on the next coarser
NJNEST = 1,10,71, ! nest where the lower southwest
! corner of this nest will start.
! If NINEST or NJNEST = 0, use CENTLAT/LON
POLELAT = 28.184., !-89.99, ! If polar, latitude/longitude of pole point
POLELON = 73.93., ! If lambert, lat/lon of grid origin (x=y=0.)
STDLAT1 = 28.184., ! If polar for BRAMS, use 90.0 in STDLAT2
STDLAT2 = 0., ! If lambert, standard latitudes of projection
!(truelat2/truelat1 from namelist.wps, STDLAT1 < STDLAT2)
! If mercator STDLAT1 = 1st true latitude
CENTLAT = 28.184.,!-89.99, -23., 27.5, 27.5,
CENTLON = 73.93., -46.,-80.5, -80.5,
!---------------- model output domain for each grid (only set up for rams)
lati = -90., -90., -90.,
latf = +90., +90., +90.,
loni = -180., -180., -180.,
lonf = 180., 180., 180.,
!---------------- project rams grid (polar sterogr) to lat/lon: 'YES' or 'NOT'
proj_to_ll='YES',
!---------------- output file prefix (may include directory other than the current)
chem_out_prefix = 'emissions',
chem_out_format = 'vfm',
!---------------- convert to WRF/CHEM (yes,not)
special_output_to_wrf = 'YES',
$END
I am not able to get what exactly the problem is here, please help me with this issue.