I am try to run prep_chem_sources...exe with RETRO and mercator projection. The prep_chem_sources.inp is given as below:
- Code: Select all
$RP_INPUT
!#################################################################
! CCATT-BRAMS/MCGA-CPTEC/WRF-Chem/FIM-Chem emission models CPTEC/INPE
! version 1.3: May 2012
! contact: gmai@cptec.inpe.br - http://meioambiente.cptec.inpe.br
!#################################################################
!---------------- grid_type of the grid output
!grid_type= 'll',
grid_type= 'mercator',
rams_anal_prefix = '../ANL/OPQUE',
!---------------- date of emission
ihour=00,
iday=16,
imon=02,
iyear=2015,
!---------------- select the sources datasets to be used
use_retro=1, ! 1 = yes, 0 = not
!retro_data_dir='./datain/EMISSION_DATA/RETRO/anthro',
retro_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/RETRO/anthro',
use_edgar =0, ! 0 - not,
! 1 - Version 3,
! 2 - Version 4 for some species
! 3 - Version HTAP
!edgar_data_dir='./datain/EMISSION_DATA/EDGAR/anthro/hdf',
edgar_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/EDGARV4',
use_gocart=1,
!gocart_data_dir='./datain/EMISSION_DATA/GOCART/emissions',
gocart_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/GOCART/emissions',
use_streets =0,
streets_data_dir='./datain/EMISSION_DATA/STREETS',
use_seac4rs =0,
seac4rs_data_dir='./datain/EMISSION_DATA/SEAC4RS',
use_fwbawb =0,
!fwbawb_data_dir ='./datain/EMISSION_DATA/Emissions_Yevich_Logan',
fwbawb_data_dir ='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/Emissions_Yevich_Logan',
use_bioge =2, ! 1 - geia, 2 - megan
! ######
! # BIOGENIC = 1
!bioge_data_dir ='./datain/EMISSION_DATA/biogenic_emissions',
! # MEGAN = 2
! ######
!bioge_data_dir='./datain/EMISSION_DATA/MEGAN/2000',
! ######
bioge_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/MEGAN',
use_gfedv2=0,
!gfedv2_data_dir='./datain/EMISSION_DATA/GFEDv2-8days',
gfedv2_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/GFEDv2-8days',
use_bbem=0,
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 ='./datain/EMISSION_DATA/FIRES/GOES/f',
!bbem_modis_data_dir ='./datain/EMISSION_DATA/FIRES/MODIS/Fires',
!bbem_inpe_data_dir ='./datain/EMISSION_DATA/FIRES/DSA/Focos',
!bbem_extra_data_dir ='NONE',
bbem_wfabba_data_dir ='/home/maslab-1/data/wrfchemData/Data_WRFChem/Fire_data/WF_ABBA',
bbem_modis_data_dir ='/home/maslab-1/data/wrfchemData/Data_WRFChem/Fire_data/MODIS/Fires.',
bbem_inpe_data_dir ='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/fires_data/DSA/Focos',
bbem_extra_data_dir ='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/fires_data/BLMALASKA/current.dat',
!---------------- veg type data set (dir + prefix)
!veg_type_data_dir='./datain/SURFACE_DATA/GL_IGBP_INPE_39classes/IGBP',
veg_type_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/surface_data/GL_IGBP_MODIS_INPE/MODIS',
!---------------- vcf type data set (dir + prefix)
use_vcf = 0,
vcf_type_data_dir ='./datain/SURFACE_DATA/VCF/data_out/2005/VCF',
!---------------- olson data set (dir + prefix)
!olson_data_dir='./datain/EMISSION_DATA/OLSON2/OLSON',
olson_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/OLSON2/OLSON',
!---------------- carbon density data set (dir + prefix)
!carbon_density_data_dir='./datain/SURFACE_DATA/GL_OGE_INPE/OGE',
carbon_density_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/surface_data/GL_OGE_INPE/OGE',
!fuel_data_dir='./datain/EMISSION_DATA/Carbon_density_Saatchi/amazon_biomass_final.gra',
fuel_data_dir='',
!---------------- gocart background
use_gocart_bg=1,
!gocart_bg_data_dir='./datain/EMISSION_DATA/GOCART',
gocart_bg_data_dir='/home/maslab-1/data/wrfchemData/Data_WRFChem/Global_emissions_v3/Emission_data/GOCART',
!---------------- volcanoes emissions
use_volcanoes =0,
!volcano_index =1143, !REDOUBT
volcano_index =0, !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 ='./datain/EMISSION_DATA/VOLC_SO2',
degass_volc_data_dir ='',
!---------------- 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 = 1, ! Number of grids to run
!NNXP = 150,50,86,46, ! Number of x gridpoints
!NNYP = 90,50,74,46, ! Number of y gridpoints
!NXTNEST = 0,1,1,1, ! Grid number which is the next coarser grid
!DELTAX = 111487.,
!DELTAY = 111487., ! X and Y grid spacing
NNXP = 82,103, ! Number of x gridpoints
NNYP = 96,118, ! Number of y gridpoints
NXTNEST = 0,1,1, ! Grid number which is the next coarser grid
DELTAX = 10000.,
DELTAY = 10000., ! X and Y grid spacing
! Nest ratios between this grid and the next coarser grid.
!NSTRATX = 1,2,3,4, ! x-direction
!NSTRATY = 1,2,3,4, ! y-direction
NSTRATX = 1,3, ! x-direction
NSTRATY = 1,3, ! y-direction
!NINEST = 1,10,0,0, ! Grid point on the next coarser
!NJNEST = 1,10,0,0, ! nest where the lower southwest
! corner of this nest will start.
! If NINEST or NJNEST = 0, use CENTLAT/LON
NINEST = 1,0, ! Grid point on the next coarser
NJNEST = 1,0, ! nest where the lower southwest
!POLELAT = 85., !-89.99, ! If polar, latitude/longitude of pole point
!POLELON = 0., ! If lambert, lat/lon of grid origin (x=y=0.)
POLELAT = 20.18., !-89.99, ! If polar, latitude/longitude of pole point
POLELON = 85.50., ! If lambert, lat/lon of grid origin (x=y=0.)
!STDLAT1 = 0., ! If polar for BRAMS, use 90.0 in STDLAT2
!STDLAT2 = -60., ! If lambert, standard latitudes of projection
!(truelat2/truelat1 from namelist.wps, STDLAT1 < STDLAT2)
! If mercator STDLAT1 = 1st true latitude
STDLAT1 = 20.18., ! If polar for BRAMS, use 90.0 in STDLAT2
STDLAT2 = 00.0., ! If lambert, standard latitudes of projection
!CENTLAT = 85.,!-89.99, -23., 27.5, 27.5,
!CENTLON = 0., -46.,-80.5, -80.5,
CENTLAT = 20.18., 20.18., !ref_lat/ref_lon from namelist.wps
CENTLON = 85.50., 85.50.,
!---------------- 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',
!proj_to_ll='NO',
!---------------- output file prefix (may include directory other than the current)
!chem_out_prefix = 'matrixfire',
chem_out_prefix = 'Emission',
chem_out_format = 'vfm',
!---------------- convert to WRF/CHEM (yes,not)
special_output_to_wrf = 'YES',
$END
While running prep_chem_sources.exe. I am getting this error
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 .
i have no clue what will cause this problem.
Any suggestion would help a lot to solve this problem.
Best Regards
Sunny