Dear Rahimi,
Anyhow you will not be able to get aerosol optical depth as a direct output variable in any of WRF_Chem simulations. You have to calculate AOD from the output variable 'EXTCOF55' (Extinction coefficients for .55um).
########################################### Calculation Using MATLAB ####################################################
Height_in_Km = ((PH + PHB)/(9.81))*(1000) ##### This will calculate the height between WRF levels in Km where PH = Perturbation Geopotential and PHB = Base-State Geopotential.
Then if you have number of levels in your output is suppose K
for X=1:K
Distance_Between_Levels (:,:,X) = Height_in_Km(:,:,X+1) - Height_in_Km(:,:,X);
EXTCOF55_Converted(:,:,X) = (EXTCOF55(:,:,X).*Distance_Between_Levels(:,:,X));
end
AOD(:,:)= nansum(EXTCOF55_Converted,3); ############## Mean of the Level Dimension #########
#################################################### END ##############################################################
To calculate AOD, you just need to integrate the
extinction coefficient throughout the entire model column; i.e. multiple the
extinction coefficient by the thickness of the model box and sum all the boxes:
AOD = ∑ EXT(l) * dz(l)
In the user groups you will be able to find the details of this calculation.