MERRA2 aerosol

I try to run the pyCRTM using MERRA2 aerosol data.
MERRA2 3d data has 5 types of dust and 5 types of sea salt.
In pyCRTM, based on the test sample, just 1 aerosol data is used here.

profiles.aerosols[i,:,0,0] = np.asarray(h5[‘aerosolConcentration’])
profiles.aerosols[i,:,0,1] = np.asarray(h5[‘aerosolEffectiveRadius’])

How can I select one value of these data?
And also, how can I calculate the aerosol effective radius?
Or can you recommend some aerosol data for CRTM?

I found a way to change the number of input aerosols.
Now I’m wondering how I can combine the 5 dust bins into one dust variable.

Hello, thanks for reaching out!
You don’t need to combine the 5 dust bins. You may use 5 dust species, and assign the mode/average radii of the bins to the corresponding 5 dust species. With the default setup, CRTM will adjust the optical properties of dust based on input dust effective radius.

Thanks for answering! :smiley:

The term “assign mode” refers to specifying the aerosol type, is it correct?
I was planning to input the aerosol concentrations and use the effective radius values (Reff) from the “AerosolCoeff.nc4” file to run CRTM.

If I want to use 5 dust bins, do I have to manually calculate the effective radius of each layer and put it as a variable?

Hello

The term “assign mode” refers to specifying the aerosol type, is it correct?

Yes. With the default table “AerosolCoeff.nc4” file, you need to specify the aerosol type, and you also need to input/assign the profiles of effective radius.

If I want to use 5 dust bins, do I have to manually calculate the effective radius of each layer and put it as a variable?

Yes, that would be the case. You may use the effective radius of MERRA2 aerosols as input for CRTM.