2013/2/16

This directory contains files to calculate ion-mediated nucleation (IMN)
rates as a function of [H2SO4], temperature (T), relative humidity (RH),
ionization rate (Q), and surface area (S). IMN reduces to binary 
homogeneuous nucleation (BHN) when Q=0.

A scheme to calculate ionization rates as a function of latitude,
lognitude, and altitude is also provided.

References: 
 Ion-mediated nucleation (IMN): 
   Yu, F., Ion-mediated nucleation in the atmosphere: Key controlling parameters, 
        implications, and look-up table, J. Geophy. Res., 115, D03206, 
        doi:10.1029/2009JD012630, 2010.
   Yu, F., From molecular clusters to nanoparticles: Second-generation ion-mediated 
        nucleation model, Atmos. Chem. Phys.,6, 5193-5211, 2006.
   Yu, F. and R. P. Turco, From molecular clusters to nanoparticles: The role of 
        ambient ionization in tropospheric aerosol formation, J. Geophy. Res., 
        106, 4797-4814, 2001.

 Binary homogeneuous nucleation (BHN):
   Yu, F., Updated H2SO4-H2O binary homogeneous nucleation rate look-up tables, 
       J. Geophy. Res.,113, D24201, doi:10.1029/2008JD010527, 2008.
   Yu, F., Improved quasi-unary nucleation model for binary H2SO4-H2O homogeneous 
       nucleation, J. Chem. Phys., 127, 054301, 2007.

 Scheme to calculate ionization rates:
    Yu, F., G. Luo , T. Bates , B. Anderson , A. Clarke , V. Kapustin , R. Yantosca , 
       Y. Wang , S. Wu, Spatial distributions of particle number concentrations in 
       the global troposphere: Simulations, observations, and implications for nucleation 
       mechanisms, J. Geophys. Res., 115, D17205, doi:10.1029/2009JD013473, 2010.


Files included in the package:

IMNBHNtable/: IMN and BHN lookup table,
              2-D annual mean galactic cosmic ray ionization rates.

apm_nucl_mod.f: Stand alone module for IMN and BHN rate calculation.
test.f:         Test program to demonstrate how the lookup table is used.

inputGridInfo.txt: Information of grid box location for calculating
                   ionization rate for the test program.
inputGridVars.txt: Sample variables for nucleation rate calculation
                   for the test program.

Use "f90 *.f" to compile and get executable program a.out

To run the test program:

>a.out

For your reference, screen output of running the program on my
workstation has been saved to outputTest.txt.


I hope that the comments provided in the subroutines make
the code self-explanatory. Contact me if you have questions.

Enjoy and have fun!

Fangqun Yu
SUNY-Albany
(fyu@albany.edu)