How to extract pm2.5 species in WRF-Chem with MOSAIC 4BIN?

[gc13141112]

No description provided.

[douglowe]

Hi,

The scripts should be able to extract PM2.5 species for MOSAIC 4 bin aerosol in the same way they do for 8 bin aerosol. Can you give me more details on what you've tried, and I'll see if I can help you get the scripts working.

cheers,
Doug

[gc13141112]

Thanks for your reply.
I use python scripts to read wrfout file and extract concentrations of PM2.5 species, and the source code is as follows:
data = xr.open_dataset('/data2/phd/wrfchem/'+feedback+'/2017' + mons[ttt] +'/wrfout_d01_'+str(date)[0:10]+'_00:00:00')
PM25_SO4 = (data.data_vars['so4_a01'].values+data.data_vars['so4_a02'].values+data.data_vars['so4_a03'].values)/data.data_vars['ALT'].values
PM25_NO3 = (data.data_vars['no3_a01'].values+data.data_vars['no3_a02'].values+data.data_vars['no3_a03'].values)/data.data_vars['ALT'].values
PM25_CL = (data.data_vars['cl_a01'].values+data.data_vars['cl_a02'].values+data.data_vars['cl_a03'].values)/data.data_vars['ALT'].values
PM25_NH4 = (data.data_vars['nh4_a01'].values+data.data_vars['nh4_a02'].values+data.data_vars['nh4_a03'].values)/data.data_vars['ALT'].values
PM25_NA = (data.data_vars['na_a01'].values+data.data_vars['na_a02'].values+data.data_vars['na_a03'].values)/data.data_vars['ALT'].values
PM25_UNSPEC1 = (data.data_vars['oin_a01'].values+data.data_vars['oin_a02'].values+data.data_vars['oin_a03'].values)/data.data_vars['ALT'].values
PM25_EC = (data.data_vars['bc_a01'].values+data.data_vars['bc_a02'].values+data.data_vars['bc_a03'].values)/data.data_vars['ALT'].values
PM25_OC = (data.data_vars['oc_a01'].values+data.data_vars['oc_a02'].values+data.data_vars['oc_a03'].values)/data.data_vars['ALT'].values
Can you help check it ?
Chao

[douglowe]

hi Chao,

Is that the altitude value you are dividing the mass concentrations by? Are you trying to get the total column amount of aerosol?

cheers,
Doug