Could not extract FINN emission file in CMAQv5.4

Dear All

I use finn2cmaq to output FINN emission and add it to CMAQv5.4 (with cb6r3_ae7 mechanism) but get an error.


I check this post and add the missing species (the value is 0) with gspro of finn2cmaq but get the same error.

Here is my log
CTM_LOG_001.v54_intel_tibet27km_20220727.txt (22.1 KB)

the FINN emission file

the CMAQ_Control_DESID_cb6r3_ae7_aq.nml file
CMAQ_Control_DESID_cb6r3_ae7_aq.nml.txt (14.8 KB)

the gspro of finn2cmaq
gc12_to_cb6r3_ae7_v2.txt (3.1 KB)

Any advice would be appreciated.
Li

You say “the missing species.” Which species is missing? I don’t see a species name in the log.

Hi barronh
Thank you very much for your reply.

I compared the var-list between the output of finn2cmaq and SMOKE which ran successfully and found these species are missing: ALDX, ISOP, MEOH, NVOL, OLE, PH2O, SULF, TERP, UNK, UNR, VOC_INC, XYL.

However, I found the number of species in my emission file of MEGANv2.1 is 25. I also marked the species that were missing from the output of finn2cmaq: ISOP, TERP, XYL, OLE, NR, MEOH, CH4, ETOH, ALDX.

The same missing species are ALDX, ISOP, MEOH, OLE, TERP, XYL. I don’t know if the error was caused by missing these species or if there were other causes.

And also I didn’t use the original txt2daily.py scripts I changed it to use FINNv2.5 data. Could this cause the problem?

Thanks a lot.
Li

I am not familiar with @barronh 's finn2cmaq tool, but it might help if you could post a text file containing the output of ncdump -h /public5/home/sch6618/soft/model/CMAQv5.4/data/emis/tibet27km/FINNv2.5_2022.GEOSCHEM.nc . Something in the structure or content of this file causes CMAQ’s read function to fail.

Hi hogrefe.christian

Thank you very much for your reply.

Here is the txt file.
FINN.txt (27.1 KB)

Thanks for posting the full header of your file.

Overall, nothing really inconsistent stood out to me, but I did notice that the HISTORY attribute shows that ncrcat was used to string together about 40 individual day-specific FINN files into this final file with 975 hourly time steps.

Was this done as part of running finn2cmaq, or afterwards? Generally, using NCO or CDO operators on files conforming to I/O API conventions will generate output files that violate these conventions and cannot be used as input to CMAQ or other tools that use the I/O API library for file input. In particular, using such operators will often mess up TFLAG which then renders the file not usable by I/O API functions. A simple test to see if file FINNv2.5_2022.GEOSCHEM.nc conforms to conventions is to use it as input to m3stat or other ‘m3tools’ programs. If this fails, it will also fail when used as input to CMAQ.

If using ncrcat outside of finn2cmaq was indeed the problem while the original day-specific files generated by finn2cmaq were fine, you could instead use ‘m3tools’ like m3xtract to concatenate these files into the longer time period for which you want to create your FINN emissions file.

1 Like

Hi hogrefe.christian

Thank you very much.

Yes, I use ncrcat outside of finn2cmaq to combine multiple days as one. And this causes the error. I used m3cple instead and it finally worked successfully.

Thanks again, It helps me a lot.

1 Like

This is not a correctly formed I/O API file. In particular, it is missing many of the required header-attributes, including particularly VAR-LIST (and also SDATE, STIME, TSTEP, NCOLS, NROWS, NTHIK, NVARS, GDTYP, P_ALP, P_BET, p_GAM, XCENT, YCENT, XCELL, YCELL, XORIG, YORIG)

Hi cjcoats

Thank you very much for your reply. I will remember this.