SUBROUTINE iniaerosol()


      use radinc_h, only: naerkind
      use aerosol_mod

      IMPLICIT NONE
c=======================================================================
c   subject:
c   --------
c   Initialization related to aerosols 
c   (CO2 aerosols, dust, water, chemical species, ice...)   
c
c   author: Laura Kerber
c   ------
c        
c=======================================================================


#include "callkeys.h"
      integer ia

      ! Initialize aerosol indexes to zero:
      ! NB: iaero_* indexes are in module aerosol_mod

      iaero_co2   =0
      iaero_h2o   =0
      iaero_dust  =0
      iaero_h2so4 =0
      
      ia=0
      if (aeroco2) then
      ia=ia+1
      iaero_co2=ia
      write(*,*) 'CO2 aerosol = ', iaero_co2
      endif
 
      if (aeroh2o) then
      ia=ia+1
      iaero_h2o=ia
      write(*,*) 'H2O aerosol = ', iaero_h2o
      endif

      if (dusttau.gt.0) then
      ia=ia+1
      iaero_dust=ia
      write(*,*) 'Dust aerosol = ', iaero_dust
      endif

      if (aeroh2so4) then
      ia=ia+1
      iaero_h2so4=ia
      write(*,*) 'H2SO4 aerosol = ', iaero_h2so4
      endif
      
      write(*,*) 'Number of aerosols= ', ia
      
      if (ia.eq.0) then  !For the zero aerosol case
      ia = 1
      endif

      if (ia.ne.naerkind) then
          print*, 'Aerosols counted not equal to naerkind'
          print*, 'Compile with tag -s',ia,'to run'
          print*, 'or change options in callphys.def'
          print*, 'Abort in iniaerosol.F'
      call abort
      endif

      end