Difference between revisions of "Building Tables Of Aerosol Optical Properties"
| Line 5: | Line 5: | ||
You first need to find the necessary indices for your aerosols of choice. The HITRAN database is an excellent place to start: https://hitran.org/aerosols/ | You first need to find the necessary indices for your aerosols of choice. The HITRAN database is an excellent place to start: https://hitran.org/aerosols/ | ||
| − | + | You will end up with a table with the real and imaginary parts of the optical constants according to wavelength. If they are not there, you will have to dig into the literature. There are also several “Handbooks of Optical Constants”. | |
=== Mie code === | === Mie code === | ||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
2. Once you have the optical indices for your aerosol, divide the values into two separate files, one for the visible wavelengths and one for the infrared wavelengths. Run them through the Mie code. | 2. Once you have the optical indices for your aerosol, divide the values into two separate files, one for the visible wavelengths and one for the infrared wavelengths. Run them through the Mie code. | ||
If you have never used the Mie code before, read this advice from J.B. Madeleine: | If you have never used the Mie code before, read this advice from J.B. Madeleine: | ||
“In the Mie code, nsize is the number of bins used to integrate Qext, omega, and g over the particle radius. Usually 1E4 is a good number. You can try to increase this number to see if it changes the results, and if it does, keep increasing it until you reach constant and consistent results. Rmin and rmax are the minimum and maximum particle radius used for the integration. In this case, 1.E-10 and 1E-4 usually work because the reference radius is often around 1.E-6 (this may change depending on your aerosol). Nsun is the number of wavelengths in the input file optind_*.dat. So you will run the Mie code twice, once to create the properties for the visible spectrum, once to make the properties for the infrared part of the spectrum.” | “In the Mie code, nsize is the number of bins used to integrate Qext, omega, and g over the particle radius. Usually 1E4 is a good number. You can try to increase this number to see if it changes the results, and if it does, keep increasing it until you reach constant and consistent results. Rmin and rmax are the minimum and maximum particle radius used for the integration. In this case, 1.E-10 and 1E-4 usually work because the reference radius is often around 1.E-6 (this may change depending on your aerosol). Nsun is the number of wavelengths in the input file optind_*.dat. So you will run the Mie code twice, once to create the properties for the visible spectrum, once to make the properties for the infrared part of the spectrum.” | ||
Save these as different files with names like “optprop_aerosolnamevis_n20.dat”, where n## indicates how many radii you used in the Mie code (only to help you remember). Put these files in datagcm. | Save these as different files with names like “optprop_aerosolnamevis_n20.dat”, where n## indicates how many radii you used in the Mie code (only to help you remember). Put these files in datagcm. | ||
Revision as of 10:21, 21 February 2024
This page describes the main steps required to build new aerosol optical property tables for the Generic PCM.
Find a database
You first need to find the necessary indices for your aerosols of choice. The HITRAN database is an excellent place to start: https://hitran.org/aerosols/
You will end up with a table with the real and imaginary parts of the optical constants according to wavelength. If they are not there, you will have to dig into the literature. There are also several “Handbooks of Optical Constants”.
Mie code
2. Once you have the optical indices for your aerosol, divide the values into two separate files, one for the visible wavelengths and one for the infrared wavelengths. Run them through the Mie code. If you have never used the Mie code before, read this advice from J.B. Madeleine: “In the Mie code, nsize is the number of bins used to integrate Qext, omega, and g over the particle radius. Usually 1E4 is a good number. You can try to increase this number to see if it changes the results, and if it does, keep increasing it until you reach constant and consistent results. Rmin and rmax are the minimum and maximum particle radius used for the integration. In this case, 1.E-10 and 1E-4 usually work because the reference radius is often around 1.E-6 (this may change depending on your aerosol). Nsun is the number of wavelengths in the input file optind_*.dat. So you will run the Mie code twice, once to create the properties for the visible spectrum, once to make the properties for the infrared part of the spectrum.” Save these as different files with names like “optprop_aerosolnamevis_n20.dat”, where n## indicates how many radii you used in the Mie code (only to help you remember). Put these files in datagcm.
