# How to set CMAQ/HDDM sensitivity input parameters for EMIS in specific REGIONS? **URL:** \ **Category:** CMAQ-DDM-3D\ **Created:** [July 25, 2022, 4:45pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444 "2022-07-25T16:45:28Z")\ **Posts on this page:** 11\ **Page:** 1
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 25, 2022, 4:45pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/1 "2022-07-25T16:45:29Z")
Hi, I am running CMAQv5.2 with HDDM successfully using the sensitivity input parameters (see attached [sensinput.dat.txt](https://forum.cmascenter.org/uploads/short-url/a9WZg1y6AvjW3KBGLaPkPJzRTcf.txt) (759 Bytes) for domain wide. Now I am going to run the CMAQ/HDDM focusing on sensitivity of only **anthropogenic emissions** (not all emissions including biogenic, wildfire, and lightning etc.) to NOx and VOC **in specific 2 regions** (I have developed region files in netcdf format used for CMAQ/ISAM successfully) such as **8-hour ozone nonattainment area** and \*\*entire county or state \*\*. How should I set my sensitivity input parameters? I developed a new sensitivity input parameter file (see attached [sensinput.dat.txt](https://forum.cmascenter.org/uploads/short-url/7b3wnwzgOIP6vfSOhpSNTZEXbhg.txt) (769 Bytes)) according to " \* The term ’ REGIONS’ (all-cap, one leading space). If this term is used, the next line must give an integer (2-leading spaces) number of regions to be defined. The following lines are used to specify regions from a regions netcdf file, and must exactly match the region name (2-leading spaces)." according to CMAQ-HDDM document at [CMAQ/CMAQ\_DDM.md at 5.2\_DDM-3D · USEPA/CMAQ · GitHub](https://github.com/USEPA/CMAQ/blob/5.2_DDM-3D/DOCS/Instrumented_Docs/CMAQ_DDM.md). But it does NOT work. It is really appreciated if you could provide an example of sensinput.dat according to my purpose because I did not find the sample on this anywhere. Thank you so much for your time. ---
**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\ **Post date:** [July 25, 2022, 5:26pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/3 "2022-07-25T17:26:58Z")
Here is an example of mixed emissions streams for CMAQ5.2 DDM. By the way, this stuff a lot more simple in the newly upcoming 5.4 version. ANX EMIS GRID 1 USANTHRO PT3D 2 STK\_GRPS\_07 STK\_EMIS\_07 STK\_GRPS\_09 STK\_EMIS\_09 SPECIES NO2 NO REGIONS 2 RGN1 RGN2 ---
**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\ **Post date:** [July 25, 2022, 5:27pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/4 "2022-07-25T17:27:39Z")
I can’t get the spacing to show up correctly on the forum, so follow the directions on that. Sergey ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 25, 2022, 5:35pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/5 "2022-07-25T17:35:13Z")
Hi Sergey, Thank you so much for sharing the example. I will try accordingly. ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 25, 2022, 6:10pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/6 "2022-07-25T18:10:09Z")
Sergey, One more question, can RGN1 and RGN2 overlay (e.g. RGN1 is covered by RGN2)? Thanks. ---
**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\ **Post date:** [July 25, 2022, 6:18pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/7 "2022-07-25T18:18:46Z")
You can create the ioapi file however you want, obviously. The DDM-3D model will read it and be ok with overlapping regions, because the first-order sensitivity parameters are independent of each other. The ISAM model will hopefully crash, because you are trying to have exact mass balance, and overlapping regions will possibly double count some emissions. Sergey ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 25, 2022, 6:25pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/8 "2022-07-25T18:25:50Z")
Sergey, You said the first-order sensitivity parameters are independent of each other. How are the higher order parameters working on this? I see the important difference between HDDM and ISAM (at least for the first-order sensitivity parameters so far). Thank you very much for the explanation. Feng ---
**Author:** ![sergey](https://avatars.discourse-cdn.com/v4/letter/s/9fc348/32.png) [@sergey](https://forum.cmascenter.org/u/sergey)\ **Post date:** [July 25, 2022, 9:14pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/9 "2022-07-25T21:14:22Z")
Higher order sensitivities are dependent on the existence of correctly defined lower-order (first) sensitivities in the same control file. ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 25, 2022, 10:55pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/10 "2022-07-25T22:55:07Z")
Hi Sergey, Thank you again for your reply. When I use mixed emissions streams in EMIS and REGIONS for CMAQ5.2 HDDM, I developed the sensitivity input parameters (sensinput.data.new [sensinput.dat.new.txt](https://forum.cmascenter.org/uploads/short-url/nyW1rIdjwxR2j8lIKPP4zKqF8zW.txt) (730 Bytes) attached ) according to the example your provided. I also have the following settings in my run script setenv DDM3D\_RGN Y # allows for a file specifying regions [T | Y | F | N] (default is N/F) set it below if Y/T setenv REGIONS\_1 **/INPUTDIR/GRIDMASK\_REGIONS.nc** # this is the same region file used for ISAM. The netcdf file includes RGN1 and RGN2. Is that correct? In addition, I do not know if I need and how to set up the following environment parameters? setenv DDM3D\_BCRGN N # allows for a file specifying boundary regions [T | Y | F | N] (default is N/F) set it below if Y/T setenv BCREGIONS\_1 Thanks, ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 26, 2022, 5:51pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/11 "2022-07-26T17:51:31Z")
Hi Sergey, I get the following error when I refer to your example of sensitivity input parameters above. The error indicates the model consider the term REGIONS as one of SPECIES. I am wondering if we put the term REGIONS in a wrong place. ``` *** ERROR ABORT in subroutine SEN_NUM on PE 001 Sensitivity parameter REGIONS not found in species name array ``` PM3EXIT: date&time specified as 0 Thank you very much for taking the time to check. ---
**Author:** ![FengLiu](https://avatars.discourse-cdn.com/v4/letter/f/49beb7/32.png) [@FengLiu](https://forum.cmascenter.org/u/FengLiu)\ **Post date:** [July 26, 2022, 5:57pm UTC](https://forum.cmascenter.org/t/how-to-set-cmaq-hddm-sensitivity-input-parameters-for-emis-in-specific-regions/3444/12 "2022-07-26T17:57:11Z")
Hi Sergey, I figure it out that The term ’ REGIONS’ (all-cap, one leading space). I had 2 leading space. Thanks.