In your run script you have one of the emissions streams labeled as PT_EGU
setenv STK_EMIS_LAB_002 PT_EGU
In your CMAQ_Control_DESID.nml.txt file, you define the StreamFamilyName to be identical as the StreamFamilyMembers.
I don’t know what the impact of having the same name for the stream family as the stream family member.
StreamFamilyName(2) = ‘PT_EGU’
StreamFamilyMembers(2,1:1)= ‘PT_EGU’
I would recommend using a unique name for the StreamFamilyName(2). To distinuguish between the two. You could use PT_EGU_FAM for the StreamFamilyName. I think we had suggested PT_EGUS in the responses above, which was not clear.
StreamFamilyName(2) = ‘PT_EGU_FAM’
StreamFamilyMembers(2,1:1)= ‘PT_EGU’
However, you may not need a new StreamFamilyName, as the PT_EGU is the only member in the StreamFamilyName(2). So, you could use the PT_EGU emission stream label, without defining a new family.
It doesn’t appear that you need any new stream families to be defined, as you only want to multiply the PT_EGU stream by 10 for SO2, and all emission streams by 10 for NOX.
I am going to assume that the clearest method is to not define any Stream Families. Note, I have commented out the StreamFamilyName(1) and StreamFamilyName(2) by using the ! at the front of those lines.
&Desid_StreamFamVars
Desid_N_Stream_Fams = 0 ! Exact number of stream family definitions
Desid_Max_Stream_Fam_Members = 20 ! Larger than the number of streams on all
! family definitions
/
&Desid_StreamFam
! For emission streams available in several run scripts under CCTM/scripts
! Set Desid_N_Stream_Fams to 1 !!!!!!!
!StreamFamilyName(1) = 'PT_SOURCES'
!StreamFamilyMembers(1,1:3)= 'PT_NONEGU','PT_EGU','PT_OTHER'
!StreamFamilyName(2) = 'PT_EGU_FAM'
!StreamFamilyMembers(2,1:1)= 'PT_EGU'
In the diagnostic rules section, you need to activate the DESID reporting of SO2 from PT_EGU and NOX from all sources
In this example, I’ve changed the number of Diagnosic Rules to 2, and listed two rules below. The first rule will print out the column sum for the SO2 variable for the PT_EGU emissions stream, and the second rule will print out the column sum for NOX from all emission streams.
&Desid_DiagVars
Desid_N_Diag_Rules = 2 ! Exact Number of Diagnostic Rules Below
&Desid_Diag
Desid_Diag_Streams_Nml(1,:)= 'PT_EGU'
Desid_Diag_Fmt_Nml(1) = 'COLSUM'
Desid_Diag_Spec_Nml(1,:) = 'SO2'
Desid_Diag_Streams_Nml(2,:)= 'ALL'
Desid_Diag_Fmt_Nml(2) = 'COLSUM'
Desid_Diag_Spec_Nml(2,:) = 'NOX'
In the CMAQ_Control_DESID_cb6r5_ae7_aq.nml file, I think you have the scaling set correctly to increase the SO2 by 10x for the PT_EGU Stream Label, and to increase NOX by 10x for all emission streams.
&Desid_Scaling
Desid_Rules_nml =
! Region | Stream Label |Emission | CMAQ- |Phase/|Scale |Basis |Op
! Label | |Species | Species |Mode |Factor| |
! Gases
'EVERYWHERE', 'ALL' ,'NO2' ,'NO2' ,'GAS' ,1 ,'UNIT','a',
'EVERYWHERE', 'ALL' ,'NO' ,'NO' ,'GAS' ,1 ,'UNIT','a',
'EVERYWHERE', 'ALL' ,'NOX' ,'NOX' ,'GAS' ,10 ,'UNIT','m',
'EVERYWHERE', 'ALL' ,'HONO' ,'HONO' ,'GAS' ,1. ,'UNIT','a',
'EVERYWHERE', 'ALL' ,'SO2' ,'SO2' ,'GAS' ,1 ,'UNIT','a',
'EVERYWHERE', 'PT_EGU' ,'SO2' ,'SO2' ,'GAS' ,10 ,'UNIT','m',
What I am unclear about, is whether multiplying by a factor of 10 for NOX only applies to emission streams that are defined in the run script. For your case, this would be the Point Source Emissions. Or, is DESID also able to multiply by a factor of 10 the NOX from the Gridded Emissions File that is defined here:
#> Gridded Emissions files
setenv N_EMIS_GR 1 #2
set EMISfile = all/emis_mole_all_${YYYYMMDD}_${APPL}_nobeis_2016fh_16j.ncf
I found the answer, the Gridded Emissions are labeled in your run script as GRIDDED_EMIS, so they will be included in the definition of ALL.
#> Gridded Emissions files
setenv N_EMIS_GR 1 #2
set EMISfile = all/emis_mole_all_${YYYYMMDD}_${APPL}_nobeis_2016fh_16j.ncf
setenv GR_EMIS_001 ${EMISpath}/${EMISfile}
setenv GR_EMIS_LAB_001 GRIDDED_EMIS
setenv GR_EM_SYM_DATE_001 F # To change default behaviour please see Users Guide for EMIS_SYM_DATE