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:
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
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