Thanks for your response. If I’m understanding your problem correctly, there are at least two potential reasons for why your ISAM results for each region are zero - 1) the MASK_REGION file you generated with m3mask and m3merge does not provide the expected information, or 2) your CMAQ run script, EmissCtrl namelist flile, and isam_control.txt file have not been configured correctly to use the MASK_REGION file you generated.
Before going back to possibility 1) (which would be best discussed in the other thread you started), I’d like to rule out possibility 2).
Assuming the MASK_REGION file you generated is /home/pbwang/CMAQ5.3/data/maskdata1/files/cmaqmasks.nc and contains variables m01, m02, … m11 on the same grid as your emission files you’d have to make sure to do the following:
In your run script, setenv ISAM_REGIONS /home/pbwang/CMAQ5.3/data/maskdata1/files/cmaqmasks.nc
In the RegionsRegistry of your EmissCtrl file, reference the environment variable for the mask file set in your run script:
&RegionsRegistry
RGN_NML =
! | Region Label | File_Label | Variable on File
‘ALL’ ,‘ISAM_REGIONS’,‘ALL’,
/
[note that you can use any other name for the environment variable referencing your mask file, e.g. MASK_REGION, as long as the name is consistent between the run script and the EmissCtrl namelist file]
In your isam_control.txt file, define tags that reference one or more of the variables contained in your mask file (i.e. m01, m02, etc.), e.g.
!!! gridded emissions from region m01
TAG NAME |G01
REGION(S) |m01
EMIS STREAM(S) |GR_EMIS
where GR_EMIS has to be the label of an emissions stream defined in your run script (e.g. setenv GR_EMIS_LAB_001 GR_EMIS).
One way to confirm that CMAQ correctly reads the variables from your grid mask file is to check the following portion of the CTM_LOG files:
|> Check Emissions Mapping:
+============================
=============================================================================
|> SCALING EMISSIONS CONSISTENT WITH EMISSIONS CONTROL FILE SUPPLIED BY USER
=============================================================================
|> Regions Available for Scaling:
=================================
Number Region Label File Label Variable
------ ------------ ---------- --------
1 EVERYWHERE N/A N/A
2 m01 ISAM_REGIONS m01
3 m02 ISAM_REGIONS m02
etc.
If your setup matches what’s described above, it indeed would seem to suggest a problem with the mask file you generated with m3mask and m3merge. In that case, it would be best to follow up in the other thread.