Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

PREP_CHEM SRC namelist error

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

ganesh_Mk

New member
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.
 
Top