! (C) Copyright 2022- ECMWF.
!
! This software is licensed under the terms of the Apache Licence Version 2.0
! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
!
! In applying this licence, ECMWF does not waive the privileges and immunities
! granted to it by virtue of its status as an intergovernmental organisation
! nor does it submit to any jurisdiction.
!
! Author:  Robin Hogan
! Email:   r.j.hogan@ecmwf.int
!

! Program to test the functionality of the
! radiation_aerosol_optics_description module. Usage is for example:
! "test_aerosol_optics_description aerosol_ifs_49R1.nc"
program test_aerosol_optics_description

  use radiation_aerosol_optics_description

  type(aerosol_optics_description_type) :: aer_desc

  character(len=512) :: file_name

  integer :: istatus

  call get_command_argument(1, file_name, status=istatus)
  if (istatus /= 0) then
    stop 'Usage: test_aerosol_optics_description <aerosol_file.nc>'
  end if

  call aer_desc%read(file_name=file_name)

  ! These two should issue a warning that the aerosol type is ambiguous
  print *, 'DD: ', aer_desc%get_index('DD',.false.)
  print *, 'DD (bin=2): ', aer_desc%get_index('DD',.false.,ibin=2)
  ! Indicate preferred aerosol optical model, after which further
  ! calls will prefer one particular model
  print *, 'preferred_optical_model(DD,Fouquart)'
  call aer_desc%preferred_optical_model('DD','Fouquart')
  print *, 'DD (bin=2): ', aer_desc%get_index('DD',.false.,ibin=2)
  print *, 'DD (bin=2,model=Woodward): ', aer_desc%get_index('DD',.false.,ibin=2,optical_model_str="Woodward")
  print *, 'DD (bin=2,model=Woodward2001): ', aer_desc%get_index('DD',.false.,ibin=2,optical_model_str="Woodward2001")
  ! This should fail to find a match, returning zero
  print *, 'DD (model=Nobody): ', aer_desc%get_index('DD',.false.,optical_model_str="Nobody")
  print *, 'SS (bin=3): ', aer_desc%get_index('SS',.true.,ibin=3)

end program test_aerosol_optics_description