Initial condition for Chemical constituents for CMAQ

Hi,

Can anybody suggest me where to get actual initial condition data for the chemical constituents for CMAQv5.3? By that I mean a complete set of data including emissions etc. Is there such an inventory of data freely available?

If not, could you hint a little on what steps to follow to get the complete set of chemical data (except for meteorology which I can produce from WRF) to initialize the model.

I am planning to run CMAQ at about 1 km (a part of CONUS).

Thanks

Hi,

I am not sure which model year are you running the CMAQ. Following things are what I usually do:

  1. ICBC:
    a) from NOAA MOZART output
    you can download MOZART output here, which is a global air quality model
    WRF-Chem

you can then convert these data into CMAQ IC/BC by mozart2cams, please see following link

b) if you can obtain GEOS chem data, you can convert these data into CMAQ format by

c) if you can obtain hemi domain CMAQ output you can create you IC BC by icon and bcon in PREP

  1. Inventory data:
    EPA released all inventory data
    ftp://newftp.epa.gov/air/emismod/
    for example, this one includes all 2011 information that are SMOKE input ready
    ftp://newftp.epa.gov/air/emismod/2011/v3platform/

However, most of these data are in 12US1 (12km) or larger (36km) resolution (some cases EPA released 4km surrogate files), if you would like to run for 1km you have to generate 1km surrogate files

Also EPA will release 2016 EMP CMAQ input ready data for 12US1 here
https://views.cira.colostate.edu/iwdw/

Hope this helps

Joey

Many hanks for these links! I will try to work with them.

Hi

I want to compile and run mozart2camx. Could you guide me on how to do that, please?
Thank you

Best regards,
Ensiyeh

Hello Ensiyeh,

the mozart2camx program is not part of the CMAQ suite of tools and is not supported by the CMAS center or the CMAQ development team. For questions regarding the program, we suggest following the instructions on the page where you downloaded the tar file:

"Support Software

This page contains links to CAMx pre- and post-processors and other tools. These source codes are distributed by Ramboll as free software, and they contain notices like this:

Copyright © 2016-2019 Ramboll
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

Ramboll occasionally posts updates for certain widely-used programs when necessary, but does not actively support or maintain every one. Users can e-mail questions, comments, suggestions or improvements to CAMx."

Hi Christian

I know, but they didn’t answer me and I didn’t find another way to solve it. So I decided to ask it from people who compiled and run it. Could you help me, please?
Thank you

Hi Ensiyeh,

unless your question is specifically related to errors you are encountering when compiling the program and linking in the required netCDF and I/O API libraries, it is unlikely that the CMAS center or CMAQ development team would be able to help, although it is certainly possible that other users (such as @huangj1311) have used this tool in their work and might have additional insights.

Your post does not contain sufficient details to judge what you have already tried, what problems you are encountering, and what guidance you are seeking. General guidance on how to compile and run mozart2camx is beyond what either the CMAS center team or the CMAQ development team can provide because, to the best of my knowledge, no one from either team is using this tool for their work and because we are not the developers or maintainers of this tool.

I see that you also posted your question in a third thread and that @cjcoats provided a response in that thread. His response fits the exception I noted above - you can expect support on this forum if you are encountering problems with compiling code that requires the I/O API library.

Hello everyone,
Sorry to bother you, I know that this forum does not support mozart2camx but I did not find any forum related to this stuff so I am asking my question here.
I downloaded Mozart data and run ncf2ioapi_mozart successfully. But when I want to run mozart2camx_CB05_AE5__GEOS5 it gets stuck at this point:

amir-N53SM:~/mozart2camx.26feb19_1/mozart2camx_v3.2.1/src> ./mozart2camx_CB05_AE5__GEOS5

 This program uses the EPA-AREAL/MCNC-EnvPgms/BAMS Models-3      
 I/O Applications Programming Interface, [I/O API] which is      
 built on top of the netCDF I/O library (Copyright 1993, 1996    
 University Corporation for Atmospheric Research/Unidata         
 Program) and the PVM parallel-programming library (from         
 Oak Ridge National Laboratory).                                 
 Copyright (C) 1992-2002 MCNC,                                   
 (C) 1992-2018 Carlie J. Coats, Jr.,                             
 (C) 2003-2012 Baron Advanced Meteorological Systems, LLC, and   
 (C) 2014-2021 UNC Institute for the Environment.                
 Released under the GNU LGPL  License, version 2.1.  See URL     
                                                                 
     https://www.gnu.org/licenses/old-licenses/lgpl-2.1.html     
                                                                 
 for conditions of use.                                          
                                                                 
 ioapi-3.2: $Id: init3.F90 200 2021-05-10 14:06:20Z coats $
 Version with PARMS3.EXT/PARAMETER::MXVARS3= 2048
 netCDF version 4.7.4 of Jun 18 2021 03:29:43 $
  
  
 Missing environment variable EXECUTION_ID

Then when I try to run it from scripts under the Job directory (mozart2camx.cmaq.cb05.job script) the following happens:

./mozart2camx.cmaq.cb05.job
rm: cannot remove ‘./output/bc.10km.cmaq.cb05.20150710.lst.ncf’: No such file or directory
rm: cannot remove ‘./output/ic.10km.cmaq.cb05.20150710.lst.hr0.ncf’: No such file or directory

 This program uses the EPA-AREAL/MCNC-EnvPgms/BAMS Models-3      
 I/O Applications Programming Interface, [I/O API] which is      
 built on top of the netCDF I/O library (Copyright 1993, 1996    
 University Corporation for Atmospheric Research/Unidata         
 Program) and the PVM parallel-programming library (from         
 Oak Ridge National Laboratory).                                 
 Copyright (C) 1992-2002 MCNC,                                   
 (C) 1992-2018 Carlie J. Coats, Jr.,                             
 (C) 2003-2012 Baron Advanced Meteorological Systems, LLC, and   
 (C) 2014-2021 UNC Institute for the Environment.                
 Released under the GNU LGPL  License, version 2.1.  See URL     
                                                                 
     https://www.gnu.org/licenses/old-licenses/lgpl-2.1.html     
                                                                 
 for conditions of use.                                          
                                                                 
 ioapi-3.2: $Id: init3.F90 200 2021-05-10 14:06:20Z coats $
 Version with PARMS3.EXT/PARAMETER::MXVARS3= 2048
 netCDF version 4.7.4 of Jun 18 2021 03:29:43 $
  
  
 Missing environment variable EXECUTION_ID

 *** Program expects to read in MCIP met files *** 
 Value for IOAPI_CHECK_HEADERS not defined;returning default:   FALSE
  
 "INFILE1" opened as OLD:READ-ONLY   
 File name "/home/amir/mozart2camx.26feb19_1/mozart2camx_v3.2.1/outputs/ncf/OUTFILE3D"
 File type GRDDED3 
 Execution ID "????????????????"
 Grid name "MOZART_GRID"
 Dimensions: 9 rows, 14 cols, 56 lays, 86 vbles
 NetCDF ID:     65536  opened as READONLY            
 Starting date and time  2015191:060000 (6:00:00   July 10, 2015)
 Timestep                          060000 (6:00:00 hh:mm:ss)
 Maximum current record number         4

 Mozart Origin   =    36.2500000       22.7368412    
 Mozart dx,dy    =    2.50000000       1.89473701    
 Mozart nx,ny,nz =          14           9          56


 Program attempting to open MCIP met files
  
 "INFILEMET3D" opened as OLD:READ-ONLY   
 File name "/home/amir/Build_CMAQ/CMAQ_Project/data/mcip/CMAQ-Iran/METCRO3D_Iran_10km.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "CMAQ-Iran_CROSS"
 Dimensions: 157 rows, 253 cols, 29 lays, 18 vbles
 NetCDF ID:    131072  opened as READONLY            
 Starting date and time  2015191:060000 (6:00:00   July 10, 2015)
 Timestep                          010000 (1:00:00 hh:mm:ss)
 Maximum current record number        19
  
 "INFILEMET2D" opened as OLD:READ-ONLY   
 File name "/home/amir/Build_CMAQ/CMAQ_Project/data/mcip/CMAQ-Iran/METCRO2D_Iran_10km.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "CMAQ-Iran_CROSS"
 Dimensions: 157 rows, 253 cols, 1 lays, 35 vbles
 NetCDF ID:    196608  opened as READONLY            
 Starting date and time  2015191:060000 (6:00:00   July 10, 2015)
 Timestep                          010000 (1:00:00 hh:mm:ss)
 Maximum current record number        19

Your GDTYP3D is           7
This cannot be processed.
end: Not in while/foreach.

Does anyone know what is wrong?
Thanks

See https://cjcoats.github.io/ioapi/GRIDS.html: Equatorial Mercator EQMGRD3 = 7

Does CAMX support equatorial Mercator?

Having a quick look at mozart2camx_v3.2.1/src/mozart2camx.f, Equatorial Mercator is not a currently supported projection in the mozart2camx processor:

    if (GDTYP3D .eq. 1) then
      print *, 'Output projection will be lat/lon'
    elseif (GDTYP3D .eq. 2) then
      print *, 'Output projection will be LCP'  
    elseif (GDTYP3D .eq. 5) then
      print *, 'Output projection will be UTM zone', P_ALP3D
    elseif (GDTYP3D .eq. 6) then
      print *, 'Output projection will be PSP'
    else
      print *, 'Your GDTYP3D is', GDTYP3D
      print *, 'This cannot be processed.'
      stop
    endif
1 Like