Nonspherical core-shell oceanic particle single-scattering database

Updated Dec 5, 2024.

For questions, please contact: Dr. Ping Yang (project leader), [e-mail]

----------------------------------------------------Nonspherical core-shell oceanic particle single-scattering database

RRI: relative refractive index of aquatic material

File naming explanation:

Datasets begin with "set" and end with ".dat""****isca.dat" contains single-scattering properties (left to right 6 columns in total):size parameter; projected area; volume; extinction efficiency; single-scattering albedo; asymmetry factorRepeat for 93 sizes;Repeat above for 18 relative refractive indices (complex).In total 1674 rows.

"****phm.dat" contains phase matrices (left to right 6 columns in total):P_11   P_12/P_11   P_22/P_11   P_33/P_11   P_43/P_11   P_44/P_11721 scattering angles from 0 to 180 degrees with 0.25 degree interval;Repeat for 93 sizes;Repeat above for 18 relative refractive indices (complex).In total 1206954 rows.

The number after "set" varies from 0 to 7, indicating shell thickness setting:'0': Homogeneous'1': 5%'2': 10%'3': 20%'4': Average of cases 0 to 3.'5': Weighted average of cases 0 to 3, weights = [0.05,0.35,0.55,0.05]'6': Weighted average of cases 0 to 3, weights = [0.05,0.6,0.25,0.1]'7': Weighted average of cases 0 to 3, weights = [0.05,0.15,0.3,0.5]

The alphabet after number can be "a,b,c,d,x"There are in total 10 shapes numbered from 0 to 9. The shapes are rearranged in the order of sphericity.shape number = [4,9,6,8,5,2,3,7,0,1] (Sphericity from smallest to largest)sphericity = [0.64298,0.690126,0.704115,0.708696,0.72652,0.732807,0.750726,0.75268,0.761722,0.773472]'x': Equal weight for all 10 shapes.'a': Weight = [0.02,0.02,0.03,0.04,0.06,0.08,0.11,0.15,0.2,0.3]'b': Weight = [0.02,0.04,0.06,0.08,0.1,0.11,0.12,0.13,0.14,0.16]'c': Weight = [0.16,0.14,0.13,0.12,0.11,0.1,0.08,0.06,0.04,0.02]'d': Weight = [0.3,0.2,0.15,0.11,0.08,0.06,0.04,0.03,0.02,0.02]

-----------------------------------------------------------FORTRAN90 codes for retrieving single-scattering properties

spline.f90: Subroutines (written by John Barkadt, 1999) for spline interpolation

ref_interp.f90: The interpolation of single-scattering properties

main.f90: Main program

opsd.retrieve: Executable

------------------Retrieve a dataset

Input file examples:testsize.txttestwav.txttestref.dat (according to number of wavelengths)

Command line example:mkdir test1  #make a folder./opsd.retrieve 2 2 testsize.txt testwav.txt testref.dat x 4 test1Command line format:{executable} {number_of_sizes} {number_of_wavelengths} {size_input} {wavelength_input} {refractive_index_input} {shape_setting} {shell_setting} {folder_name}

-------------------------Example of output datasettest1/

7 files: isca.dat P11.dat P12.dat P22.dat P33.dat P43.dat P44.dat

isca.dat has 7 columns:

{wavelength} {particle_size} {projected_area} {volume} {extinction_efficiency} {single_scattering_albedo} {asymmetry_factor}

All data from P12.dat to P44.dat are already normalized as P_ij/P_11.