!
  
  REAL(kind=dp), parameter :: rho_sulf = 1.7e3_dp
  REAL(kind=dp), parameter :: rho_orgc = 1.8e3_dp
  REAL(kind=dp), parameter :: rho_bc   = 1.0e3_dp

  REAL(kind=dp) :: es, qvs, rh
  REAL(kind=dp) :: aer_so4, aer_oc2, aer_bc2
