# Questions about sensinput file

**URL:** <https://forum.cmascenter.org/t/questions-about-sensinput-file/1658>\
**Category:** CMAQ-DDM-3D\
**Created:** [July 19, 2020, 1:26pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658 "2020-07-19T13:26:00Z")\
**Posts on this page:** 20\
**Page:** 1

<div class="post-metadata">

**Author:** ![franky](https://avatars.discourse-cdn.com/v4/letter/f/958977/32.png) [@franky](https://forum.cmascenter.org/u/franky)\
**Post date:** [July 19, 2020, 1:26pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/1 "2020-07-19T13:26:00Z")

</div>

Hi, I’m new to CAMQ-HDDM. Although I’ve already run it succesfully, I had some questions. Looking forward to your reply, thanks!

1. EMIS term. I have some problems using “EMIS\_1”，and fixed it with replacing EMIS\_1 with TOTA. But i don’t know how to use more than 1 emision file, something like “EMIS\_2”?
2. AMOUNT term. If i want to cut 20% NOx in the region\_1 using positive number, the next line after AMOUNT should be 20, 0.2, 80 or 0.8?
3. Layer term. For now, CMAQ-5.2-DDM only support 1 layer, will it be the plan for the next version.

Looking forward to your reply, many thanks!

---

<div class="post-metadata">

**Author:** ![franky](https://avatars.discourse-cdn.com/v4/letter/f/958977/32.png) [@franky](https://forum.cmascenter.org/u/franky)\
**Post date:** [July 20, 2020, 2:46pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/2 "2020-07-20T14:46:31Z")

</div>

Although the above questions remains unsolved, i have new understanding about it, but i don’t know whether it’s right or wrong.

1. EMIS： Suppose i have 2 gridded emssion files A and B. I have to add B to A and make a new gridded file C to run CMAQv5.2 if i turn off HDDM. If i hope to get the sensitvity of B, is it true that setting EMIS\_1=A and put B in the seninput file. And I am confused about the result file. Does it mean the difference between the result of C and A?
2. Amount: Looking forward to someone’s response. It means emission rate or the ratio of emission. Many thanks to you.

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [July 20, 2020, 8:15pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/3 "2020-07-20T20:15:02Z")

</div>

Hi Franky,

Its very important to note that my answer to your question is dependent on using DDM in CMAQv5.2. The interface will be substantially changed in the upcoming release, because CMAQ emissions treatment significantly changed and the following approaches will not work!

For the EMIS term, using the keyword “TOTA” will use all of the emissions. You can specify a subset through using environment variables of the emissions streams. EMIS\_1 is probably not the correct environment variable names. I don’t remember the release scripts exactly, but you it should be something like “GRIDDED\_EMIS” for the gridded file and something like “PT\_EGU” for the point source EGUs (for example). You also need to specify the type of file you are using (this may the problem you ran into).

Here is an example of a sensitivity definition I ran a while ago (again, this will not work in 5.3+):

ANX  
EMIS  
GRID  
1  
USANTHRO  
PT3D  
5  
STK\_GRPS\_07  
STK\_EMIS\_07  
STK\_GRPS\_06  
STK\_EMIS\_06  
STK\_GRPS\_08  
STK\_EMIS\_08  
STK\_GRPS\_02  
STK\_EMIS\_02  
STK\_GRPS\_01  
STK\_EMIS\_01  
SPECIES  
NO2  
NO

This example uses 1 (one) gridded file labeled “USANTHRO” and 5 (five) pairs of stack group/stack emissions files. Species NO2 and NO are pulled out of each of these 6 emissions inputs for sensitivity.

The “AMOUNT” can not be used to cut emissions by percentages. It is an actual emissions rate amount. This option is basically never used and there will be much easier ways to accomplish this functionality in the next release, so its probably going away anyways.

The “LAYER” term is also similarly never used, so it will not be upgraded. There will be easier ways to handle these things through the new CMAQ emissions streams upgrade.

I don’t exactly understand what you are trying to accomplish with your first question in the second set. You can, however, calculate sensitivities to emissions not used in the base CMAQ simulations. So, you can specify CMAQ to use FILE\_C to calculate concentrations as it normally does. At the same time, you can calculate sensitivity to species in FILE\_B. You just need to make sure you define everything in your runscript.

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [July 20, 2020, 8:16pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/4 "2020-07-20T20:16:26Z")

</div>

For some reason, the forum collapsed the spaces out of my sensinput.dat example. You still need to preserve the correct formatting…

---

<div class="post-metadata">

**Author:** ![franky](https://avatars.discourse-cdn.com/v4/letter/f/958977/32.png) [@franky](https://forum.cmascenter.org/u/franky)\
**Post date:** [July 21, 2020, 1:11am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/5 "2020-07-21T01:11:36Z")

</div>

Hi, Sergey,

Thanks so much for your prompt reply.

For the EMIS term, the corrct environment variables is EGRIDFILE1-n.

I am doing some work about the senstivity of cutting emissions by different percentages. As the set 2 said in the first question, i want to know the effect of cutting 20% NOx of Emission File\_B in different regions.

Another question, i am confused about the result of HDDM and ISAM. Does the result of ISAM mean the difference of pulling out the emission like HDDM? In my knowledge, ISAM should include the effcect of reactions on local species(primary & secondary) and species (primary & secondary) from other regions.

Thanks so much and best regards

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [July 22, 2020, 1:11pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/6 "2020-07-22T13:11:03Z")

</div>

Franky,

To calculate the impact of a 20% cut of NOx from Emission File\_B in a region, calculate the sensitivity to NOx from that region and multiply the resulting field by 0.20.

For the difference question… ISAM calculates source apportionment. DDM calculates sensitivity. I am not exactly sure what you are asking, but it may be helpful to review some of the literature on this to get a better understanding. There are papers from our group and also from the CAMx group and some some of the European models as well.

---

<div class="post-metadata">

**Author:** ![franky](https://avatars.discourse-cdn.com/v4/letter/f/958977/32.png) [@franky](https://forum.cmascenter.org/u/franky)\
**Post date:** [July 22, 2020, 1:37pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/8 "2020-07-22T13:37:17Z")

</div>

Hi, Sergey,  
Since the sensitivity analysis cannot be used for studying the sources of secondary pollutants, the source receptor relationships for secondary pollutants should be taken into consideration. To be specific, the result of ISAM should caculate the contribution of every reaction.

I will search for more literature about ISAM. Thank you for taking the time.  
Best regrads

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [July 23, 2020, 6:50pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/9 "2020-07-23T18:50:00Z")

</div>

Franky,

To clarify, sensitivity certainly can be used to study the sources of secondary pollutants. And ISAM does calculate contributions from every reaction relevant to the thing you are tracking.

Sergey

---

<div class="post-metadata">

**Author:** ![franky](https://avatars.discourse-cdn.com/v4/letter/f/958977/32.png) [@franky](https://forum.cmascenter.org/u/franky)\
**Post date:** [August 2, 2020, 12:16pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/10 "2020-08-02T12:16:31Z")

</div>

Thanks for your reply. I’ll check it later.

---

<div class="post-metadata">

**Author:** ![Yaluk](https://avatars.discourse-cdn.com/v4/letter/y/eb9ed0/32.png) [@Yaluk](https://forum.cmascenter.org/u/Yaluk)\
**Post date:** [August 20, 2020, 9:27am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/11 "2020-08-20T09:27:21Z")

</div>

Hello Professor @Sergey, and @franky

I am still a beginner in using CMAQ-DDM-3D, fortunately I have run successfully the CMAQ\_v52\_DDM-3D (with the resulting ASENS\*, SENGRID\*, SENWDEP\* and SENDDEP\* outputs). I need help on the POST PROCESSSING of these sensitivity outputs; which tools should I use?

I have been able to use combine and VERDI to analyse the ACONC\* output successfully.

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [August 24, 2020, 4:56pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/12 "2020-08-24T16:56:42Z")

</div>

Yaluk,

CMAQ-DDM-3D outputs are in the same file format as your ACONC\* files and will open with VERDI and other tools and post processing software you use for standard CMAQ.

Sergey

---

<div class="post-metadata">

**Author:** ![Yaluk](https://avatars.discourse-cdn.com/v4/letter/y/eb9ed0/32.png) [@Yaluk](https://forum.cmascenter.org/u/Yaluk)\
**Post date:** [October 15, 2020, 4:35am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/13 "2020-10-15T04:35:01Z")

</div>

Hello @sergey,

The user guide describes clearly the use of IREGIONS term and the relevant associated specifications as well as the correct setting in the SEN\_INPUT control file. But what i am struggling with is how to create such (netcdf) file. Kindly, shed some light on how to prepare this REGIONS file… and (or) is it okay to use the ISAM regions file?

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [December 8, 2020, 10:37pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/14 "2020-12-08T22:37:10Z")

</div>

Yaluk,

The DDM regions files and the ISAM regions files use identical formatting and are interchangeable.

Sergey

---

<div class="post-metadata">

**Author:** ![Yaluk](https://avatars.discourse-cdn.com/v4/letter/y/eb9ed0/32.png) [@Yaluk](https://forum.cmascenter.org/u/Yaluk)\
**Post date:** [December 14, 2020, 2:32pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/15 "2020-12-14T14:32:38Z")

</div>

Thanks @sergey… this is very useful information

---

<div class="post-metadata">

**Author:** ![wfy](https://avatars.discourse-cdn.com/v4/letter/w/4af34b/32.png) [@wfy](https://forum.cmascenter.org/u/wfy)\
**Post date:** [February 11, 2021, 1:04am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/16 "2021-02-11T01:04:12Z")

</div>

> [@sergey](#):
>
> again, this will not work in 5.3+

Hi Sergey, I am running DDM with CMAQv5.3.2. Is there anything needed to change about this version. I had set DDM3D\_ES to Y and set the environment variable for EGRIDFILE1 to point to the emissions file in run\_cctm.csh script. But an error about “User specified DDM Emissions Stream - EGEIRFILE1 - not found in available emissions streams” occurred.

---

<div class="post-metadata">

**Author:** ![wfy](https://avatars.discourse-cdn.com/v4/letter/w/4af34b/32.png) [@wfy](https://forum.cmascenter.org/u/wfy)\
**Post date:** [March 13, 2021, 9:45am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/17 "2021-03-13T09:45:33Z")

</div>

Hi, Sergey,  
Could you help me if i want to calculate the sensitivity of PM2.5 to its precursor emissions, how should I set the species in sensinput.cmas.dat?  
I set the emissions species related to PM2.5 but an error is showed “parameter not found in species name array”.  
Thank you for taking the time.

 ![image](https://canada1.discourse-cdn.com/flex027/uploads/cmas/original/2X/f/fb36c27b0397874ab2defbe57e65dbed8d05ba70.png)

---

<div class="post-metadata">

**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\
**Post date:** [March 25, 2021, 9:23pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/18 "2021-03-25T21:23:58Z")

</div>

You need to use the CMAQ CGRID species for PM sensitivities and not the emissions surrogates.

Also, DDM for PM species is not yet implemented in the latest CMAQ release.

Sergey

---

<div class="post-metadata">

**Author:** ![wfy](https://avatars.discourse-cdn.com/v4/letter/w/4af34b/32.png) [@wfy](https://forum.cmascenter.org/u/wfy)\
**Post date:** [March 26, 2021, 11:06am UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/19 "2021-03-26T11:06:03Z")

</div>

Thank you for your reply.

---

<div class="post-metadata">

**Author:** ![erica](https://avatars.discourse-cdn.com/v4/letter/e/9dc877/32.png) [@erica](https://forum.cmascenter.org/u/erica)\
**Post date:** [May 24, 2022, 3:01pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/20 "2022-05-24T15:01:54Z")

</div>

Hello Professor @sergey ,  
Sorry for bothering you, your answer to this question confuses me. I thought that while CMAQ 5.2 version dose not supports second-order sensitivity calculation for PM species, it does support first-order sensitivity calculation. Am I right?  
Looking forward to your reply, thanks!

---

<div class="post-metadata">

**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\
**Post date:** [May 26, 2022, 10:01pm UTC](https://forum.cmascenter.org/t/questions-about-sensinput-file/1658/21 "2022-05-26T22:01:25Z")

</div>

Hi Sergey,  
I would like to make sure DDM does not work in CMAQ v5.3+ because I noticed that the compilation option of ‘set DDM3D\_CCTM’ is available in bldit\_cctm.csh and can be built CCTM\_v533\_DDM3D.exe successfully. However, I should go back to CMAQ-5.2\_DDM-3D if I need to use HDDM, am I right?

Thanks.

Feng
