Page 1 of 1

PREP_CHEM_SOURCES RUNNING PROBLEM

PostPosted: Wed Dec 28, 2016 3:25 am
by sunnykant
Hello,
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