I have a gridded emission file with 23 levels, but when I run CMAQv533 model, it prints “CTM_EMLAYS | 0 (default)”.
The log files show “Number of Emissions Layers is 23”, which is causing confusion. I wondering if I should set CTM_EMLAYS=23 in the runscript? weather my emission file is valid to read?
Thank you for your help.
In my recollection, the environment variable CTM_EMLAYS can be used to instruct the CMAQ inline plume rise calculations to set a maximum layer to which plume rise should be restricted. The default value of 0 indicates that there is no restriction, i.e. the maximum layer to which plumes are allowed to rise is the topmost model layer. The actual plume rise of course is determined by the inline stack parameters and ambient meteorological conditions and would typically result in plume rise much below the top model layers - I believe the option to impose a maximum layer by setting CTM_EMLAYS was meant to allow users to guard against cases where erroneous stack parameters would lead to unrealistically high plume rise.
If CTM_EMLAYS is specified (in which case the default value of 0 is overridden), it must be set to a value between 4 and the number of layers in the meteorology files, as noted in one of the comments in PT3D_DEFN.F.
I should also note that starting in CMAQv5.4, the functionality discussed above is no longer controlled by environment variable CTM_EMLAYS, but instead by variable
Desid_MaxLays specified in CMAQ_Control_DESID.nml.
Oh, I got it. Thanks for the answer!