# Initial condition for Chemical constituents for CMAQ

**URL:** <https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018>\
**Category:** CMAQ-PREP\
**Created:** [October 8, 2019, 12:18am UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018 "2019-10-08T00:18:24Z")\
**Posts on this page:** 11\
**Page:** 1

<div class="post-metadata">

**Author:** ![CMAQ\_user](https://avatars.discourse-cdn.com/v4/letter/c/ba8739/32.png) [@CMAQ\_user](https://forum.cmascenter.org/u/CMAQ_user)\
**Post date:** [October 8, 2019, 12:18am UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/1 "2019-10-08T00:18:24Z")

</div>

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

---

<div class="post-metadata">

**Author:** ![huangj1311](https://avatars.discourse-cdn.com/v4/letter/h/8edcca/32.png) [@huangj1311](https://forum.cmascenter.org/u/huangj1311)\
**Post date:** [October 8, 2019, 2:17pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/2 "2019-10-08T14:17:51Z")

</div>

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  
[MOZART Download](https://www.acom.ucar.edu/wrf-chem/mozart.shtml)

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

> [@Time Stamp Issue to use MOZART output as CMAQ5.3 ICBC](https://forum.cmascenter.org/t/time-stamp-issue-to-use-mozart-output-as-cmaq5-3-icbc/922):
>
> CMAQ can use boundary conditions derived from global chemistry models (GCMs) like MOZART. I use mozart2camx to process mozart output and get CMAQ v5.3 IC, BC. However, CMAQ complains 1) IC is not time-dependent, and 2) BC has wrong time stamp. It is easy to fix the first issue by updating IC netcdf file TFLAG using NCO tool. For the 2nd, CMAQ requires hourly BC, but mozart2camx provides 6-hourly outputs only based on MOZART 6-hourly output time step. I wonder if anyone has the same issue or al…

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

> **[GitHub - barronh/geos2cmaq: GEOS2CMAQ maps GEOS-Chem outputs to CMAQ boundary...](https://github.com/barronh/geos2cmaq)**
>
> GEOS2CMAQ maps GEOS-Chem outputs to CMAQ boundary conditions for any chemical mechanism. - GitHub - barronh/geos2cmaq: GEOS2CMAQ maps GEOS-Chem outputs to CMAQ boundary conditions for any chemical ...

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/](https://views.cira.colostate.edu/iwdw/)

Hope this helps

Joey

---

<div class="post-metadata">

**Author:** ![CMAQ\_user](https://avatars.discourse-cdn.com/v4/letter/c/ba8739/32.png) [@CMAQ\_user](https://forum.cmascenter.org/u/CMAQ_user)\
**Post date:** [October 9, 2019, 12:46am UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/3 "2019-10-09T00:46:55Z")

</div>

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

---

<div class="post-metadata">

**Author:** ![esm](https://avatars.discourse-cdn.com/v4/letter/e/f08c70/32.png) [@esm](https://forum.cmascenter.org/u/esm)\
**Post date:** [May 4, 2020, 2:29pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/5 "2020-05-04T14:29:11Z")

</div>

Hi

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

Best regards,  
Ensiyeh

---

<div class="post-metadata">

**Author:** ![hogrefe.christian](https://avatars.discourse-cdn.com/v4/letter/h/74df32/32.png) [@hogrefe.christian](https://forum.cmascenter.org/u/hogrefe.christian)\
**Post date:** [May 4, 2020, 6:03pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/6 "2020-05-04T18:03:47Z")

</div>

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](http://www.camx.com/download/support-software.aspx):

"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](http://www.camx.com/contact.aspx)."

---

<div class="post-metadata">

**Author:** ![esm](https://avatars.discourse-cdn.com/v4/letter/e/f08c70/32.png) [@esm](https://forum.cmascenter.org/u/esm)\
**Post date:** [May 4, 2020, 6:11pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/7 "2020-05-04T18:11:23Z")

</div>

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

---

<div class="post-metadata">

**Author:** ![hogrefe.christian](https://avatars.discourse-cdn.com/v4/letter/h/74df32/32.png) [@hogrefe.christian](https://forum.cmascenter.org/u/hogrefe.christian)\
**Post date:** [May 4, 2020, 6:31pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/8 "2020-05-04T18:31:01Z")

</div>

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.

---

<div class="post-metadata">

**Author:** ![hogrefe.christian](https://avatars.discourse-cdn.com/v4/letter/h/74df32/32.png) [@hogrefe.christian](https://forum.cmascenter.org/u/hogrefe.christian)\
**Post date:** [May 4, 2020, 6:36pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/9 "2020-05-04T18:36:34Z")

</div>

I see that you also posted your question in a [third thread](https://forum.cmascenter.org/t/camx-v6-50-compilation/1470/3) 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.

---

<div class="post-metadata">

**Author:** ![amirdudu](https://avatars.discourse-cdn.com/v4/letter/a/bcef8e/32.png) [@amirdudu](https://forum.cmascenter.org/u/amirdudu)\
**Post date:** [June 26, 2021, 3:40pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/10 "2021-06-26T15:40:28Z")

</div>

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

---

<div class="post-metadata">

**Author:** ![cjcoats](https://avatars.discourse-cdn.com/v4/letter/c/ac8455/32.png) [@cjcoats](https://forum.cmascenter.org/u/cjcoats)\
**Post date:** [June 27, 2021, 5:54pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/11 "2021-06-27T17:54:50Z")

</div>

See [https://cjcoats.github.io/ioapi/GRIDS.html](https://cjcoats.github.io/ioapi/GRIDS.html): Equatorial Mercator `EQMGRD3 = 7`

Does CAMX support equatorial Mercator?

---

<div class="post-metadata">

**Author:** ![hogrefe.christian](https://avatars.discourse-cdn.com/v4/letter/h/74df32/32.png) [@hogrefe.christian](https://forum.cmascenter.org/u/hogrefe.christian)\
**Post date:** [June 30, 2021, 12:46pm UTC](https://forum.cmascenter.org/t/initial-condition-for-chemical-constituents-for-cmaq/1018/12 "2021-06-30T12:46:43Z")

</div>

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

```
