Page 1 of 1

Error in convert_emiss.exe

PostPosted: Sat Jun 25, 2016 2:24 pm
by ozrr
Hi, I'm using WRF-Chem version 3.8 and PREP-CHEM-SRC-1.5.
I made antropogeinc emission input by prep_chem_sources_RADM_WRF_FIM_.exe.
The namelist is below.

Code: Select all
 $RP_INPUT
!#################################################################
!  CCATT-BRAMS/MCGA-CPTEC/WRF-Chem/FIM-Chem emission models CPTEC/INPE
!  version 1.5: Mar 2015
!  contact: gmai@cptec.inpe.br   - http://meioambiente.cptec.inpe.br
!#################################################################


!---------------- grid_type of the grid output
   grid_type= 'lambert',
   rams_anal_prefix = '../ANL/OPQUE',
!---------------- date of emission
    ihour=00,
    iday=13,
    imon=08,
    iyear=2015,

 !---------------- select the sources datasets to be used
   use_retro=1,  ! 1 = yes, 0 = not
   retro_data_dir='./datain/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',

   use_gocart=0,
   gocart_data_dir='./datain/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',

   use_bioge =1, ! 1 - geia, 2 - megan
   ! ######
   ! # BIOGENIC = 1
   bioge_data_dir ='./datain/Emission_data/biogenic_emissions',
   ! # MEGAN = 2
   ! ######
   !bioge_data_dir='./datain/Emission_data/MEGAN/2000',
   ! ######

   use_gfedv2=0,
   gfedv2_data_dir='./datain/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_data/WF_ABBA_v60_saulo/filt/',
   bbem_modis_data_dir    ='./datain/Emission_data/fires_data/MODIS/',
   bbem_inpe_data_dir     ='./datain/Emission_data/fires_data/DSA/',
   bbem_extra_data_dir    ='NONE',

!---------------- veg type data set (dir + prefix)
   veg_type_data_dir      ='./datain/surface_data/GL_IGBP_MODIS_INPE/MODIS',
!   veg_type_data_dir      ='./datain/surface_data/GL_IGBP_INPE_39classes/IGBP',


!---------------- 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',



!---------------- carbon density data set (dir + prefix)

   carbon_density_data_dir='./datain/surface_data/GL_OGE_INPE/OGE',

   fuel_data_dir      ='./datain/Emission_data/Carbon_density_Saatchi/amazon_biomass_final.gra',


!---------------- gocart background
   use_gocart_bg=0,
   gocart_bg_data_dir='./datain/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 ='./datain/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     = 215,274,100,        ! Number of x gridpoints
   NNYP     = 181,145,100,        ! Number of y gridpoints
   NXTNEST  = 0,1,2               ! Grid number which is the next coarser grid
   DELTAX   = 9000.,
   DELTAY   = 9000.,              ! 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,16,10,          ! Grid point on the next coarser
   NJNEST = 1,93,89,          !  nest where the lower southwest
                              !  corner of this nest will start.
                              !  If NINEST or NJNEST = 0, use CENTLAT/LON
   POLELAT  = 36.094,         ! If polar, latitude/longitude of pole point
   POLELON  = 127.169,         ! If lambert, lat/lon of grid origin (x=y=0.)

   STDLAT1  = 30.,           ! 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
   CENTLAT  = 36.094,
   CENTLON  = 127.169,



!---------------- 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 = 'new_emiss',
   chem_out_format = 'vfm',
!---------------- convert to WRF/CHEM (yes,not)
  special_output_to_wrf = 'YES',

$END



I made emission input files succesfully and I ran real.exe using the namelist.input below.

Code: Select all
&time_control
 run_days                            = 0,
 run_hours                           = 120,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2015, 2015, 2015,
 start_month                         = 08,   08,   08,
 start_day                           = 13,   13,   13,
 start_hour                          = 00,   00,   00,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2015, 2015, 2015,
 end_month                           = 08,   08,   08,
 end_day                             = 17,   17,   17,
 end_hour                            = 18,   18,   18,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 interval_seconds                    = 10800
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 720,  60,   60,
 frames_per_outfile                  = 1,    1,    1,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 auxinput1_inname                    = 'met_em.d<domain>.<date>',
 auxinput5_inname                    = 'wrfchemi_d<domain>_<date>',
 auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
 io_form_auxinput2                    = 2
 io_form_auxinput5                    = 2
 io_form_auxinput8                    = 2
 auxinput5_interval_m                = 720
 auxinput8_interval_m                = 720
 frames_per_auxinput5                = 1
 frames_per_auxinput8                = 1
 debug_level                         = 0
 /

 &domains
 time_step                           = 54,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 3,
 e_we                                = 215,   274,   100,
 e_sn                                = 181,   145,   100,
 e_vert                              = 30,    30,    30,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 9000, 3000,  1000,
 dy                                  = 9000, 3000,  1000,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     16,    10,
 j_parent_start                      = 1,     93,    89,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 mp_physics                          = 6,     6,     6,
 ra_lw_physics                       = 1,     1,     1,
 ra_sw_physics                       = 1,     1,     1,
 radt                                = 9,     3,     1,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 2,     2,     2,
 bl_pbl_physics                      = 1,     1,     1,
 bldt                                = 0,     0,     0,
 cu_physics                          = 0,     0,     0,
 cudt                                = 5,     5,     5,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 num_land_cat                        = 21,
 sf_urban_physics                    = 0,     0,     0,
 cu_diag                             = 0,     0,     0,
 /

 &fdda
 /

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

 &bdy_control
 spec_bdy_width                      = 5,
 spec_zone                           = 1,
 relax_zone                          = 4,
 specified                           = .true., .false.,.false.,
 nested                              = .false., .true., .true.,
 /

 &grib2
 /

 &chem
 chem_opt                            = 0,
 emiss_opt                           = 2,
 chem_in_opt                         = 1,
 bio_emiss_opt                       = 1,
 chem_conv_tr                        = 0,
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /



Also I ran real.exe succesfully.

I linked the emission files in test/em_real directory.

ln -sf ../../chem/PREP-CHEM-SRC-1.5/bin/new_emiss-T-2015-08-13-000000-g1-ab.bin wrfem_00to12z_d01
ln -sf ../../chem/PREP-CHEM-SRC-1.5/bin/new_emiss-T-2015-08-13-000000-g1-bb.bin emissfire_d01



But when I run convert_emiss.exe, I've got error.
The namelist.input for convert_emiss.exe and the rsl.error.0000 file are below.

Code: Select all
&time_control
 run_days                            = 0,
 run_hours                           = 12,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2015, 2015, 2015,
 start_month                         = 08,   08,   08,
 start_day                           = 13,   13,   13,
 start_hour                          = 00,   00,   00,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2015, 2015, 2015,
 end_month                           = 08,   08,   08,
 end_day                             = 13,   17,   17,
 end_hour                            = 12,   18,   18,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 interval_seconds                    = 10800
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 720,  60,   60,
 frames_per_outfile                  = 1,    1,    1,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 auxinput1_inname                    = 'met_em.d<domain>.<date>',
 auxinput5_inname                    = 'wrfchemi_d<domain>_<date>',
 auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
 io_form_auxinput2                    = 2
 io_form_auxinput5                    = 2
 io_form_auxinput8                    = 2
 auxinput5_interval_m                = 720
 auxinput8_interval_m                = 720
 frames_per_auxinput5                = 1
 frames_per_auxinput8                = 1
 debug_level                         = 0
 /

 &domains
 time_step                           = 54,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 215,   274,   100,
 e_sn                                = 181,   145,   100,
 e_vert                              = 30,    30,    30,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 9000, 3000,  1000,
 dy                                  = 9000, 3000,  1000,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     16,    10,
 j_parent_start                      = 1,     93,    89,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 mp_physics                          = 6,     6,     6,
 ra_lw_physics                       = 1,     1,     1,
 ra_sw_physics                       = 1,     1,     1,
 radt                                = 9,     3,     1,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 2,     2,     2,
 bl_pbl_physics                      = 1,     1,     1,
 bldt                                = 0,     0,     0,
 cu_physics                          = 0,     0,     0,
 cudt                                = 5,     5,     5,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 num_land_cat                        = 21,
 sf_urban_physics                    = 0,     0,     0,
 cu_diag                             = 0,     0,     0,
 /

 &fdda
 /

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

 &bdy_control
 spec_bdy_width                      = 5,
 spec_zone                           = 1,
 relax_zone                          = 4,
 specified                           = .true., .false.,.false.,
 nested                              = .false., .true., .true.,
 /

 &grib2
 /

 &chem
 chem_opt                            = 1,
 emiss_opt                           = 2,
 chem_in_opt                         = 1,
 bio_emiss_opt                       = 1,
 chem_conv_tr                        = 0,
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /



Code: Select all
taskid: 0 hostname: nid00040
 module_io_quilt_old.F        2931 T
 Ntasks in X            4, ntasks in Y            4
WRF-CHEM V3.8 EMISSIONS PREPROCESSOR
 *************************************
 Parent domain
 ids,ide,jds,jde            1         215           1         181
 ims,ime,jms,jme           -4          61          -4          52
 ips,ipe,jps,jpe            1          54           1          45
 *************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1,              521620796 bytes allocated
 open file wrfinput_d01
  med_initialdata_input: calling input_input
Timing for processing wrfinput file (stream 0) for domain        1:    1.06993 elapsed seconds
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
 LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33 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 V3.8 MODEL
*                                                                   *
*    PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at                     *
*                                                                   *
*              wrfchemhelp.gsd@noaa.gov                             *
*                                                                   *
*********************************************************************
   past opening wrfinput_d01
Time period #   1 to process = 2015-08-13_00:00:00
Time period #   2 to process = 2015-08-13_12:00:00
Total analysis times to input =    2
Total number of times to input =    2    2
 Number of emissions:         29
 EMISSIONS INPUT FILE TIME PERIOD (GMT):       -1 TO        0
forrtl: severe (67): input statement requires too much data, unit 91, file /home/WRFV3.8-chem/WRFV3/test/em_real/wrfem_00to12z_d01
Image              PC                Routine            Line        Source
convert_emiss.exe  0000000003C78746  Unknown               Unknown  Unknown
convert_emiss.exe  0000000003C9F148  Unknown               Unknown  Unknown
convert_emiss.exe  0000000003C9CE56  Unknown               Unknown  Unknown
convert_emiss.exe  000000000040D334  Unknown               Unknown  Unknown
convert_emiss.exe  000000000040133E  Unknown               Unknown  Unknown
convert_emiss.exe  0000000003D65231  Unknown               Unknown  Unknown
convert_emiss.exe  00000000004011B5  Unknown               Unknown  Unknown




In this case the wrf-chem tutorial recommends to check these :
* have a bad link,

* a zero-size file, or

* a missing the GOCART backgroud file.

But I cannot find any problems with these.


I think the namelist.input for convert_emiss.exe needs more information in &chem section,
but I have no idea how I add or modify the namelist.input.


Please answer if you have any idea or comments.

Error in running prep chem sources

PostPosted: Fri Dec 30, 2016 8:33 am
by sunnykant
Hello,
I am try to running prep_chem_sources.exe. The namelist.inp is given below:

$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 =1, ! 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/Global_emissions_v3/Emission_data/biogenic_emissions',

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 = 3, ! 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 = 88,115,62, ! Number of x gridpoints for e_we
NNYP = 98,121,61, ! Number of y gridpoints for e_sn
NXTNEST = 0,1,2, ! Grid number which is the next coarser grid
DELTAX = 100000.,
DELTAY = 100000., ! 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,3, ! x-direction
NSTRATY = 1,3,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,20,50, ! Grid point on the next coarser
NJNEST = 1,25,62, ! nest where the lower southwest

!POLELAT = 85.0, !-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 = 25.0., ! If polar for BRAMS, use 90.0 in STDLAT2
STDLAT2 = 10.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., 20.18., !ref_lat/ref_lon from namelist.wps
CENTLON = 85.50., 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

After running of 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 .

Actually i don't know what type of error is this. Anyone can help me to solve this problem.

Regards
Sunny

Re: Error in convert_emiss.exe

PostPosted: Tue Feb 28, 2017 9:18 am
by arianna.v
Did anyone find the solution?

If yes please post it because it's quite hard to find it across internet.. and the guides does not help at all.

Thank you!

Re: Error in convert_emiss.exe

PostPosted: Sat Jul 27, 2019 3:24 pm
by sunnykant
I am using a different method for producing emission data which can give as input for WRF-Chem. For generating the emission data is already mention in the user guide of WRF-Chem or WRF-Chem website.