main.F90

program main
   use data_types
   use misc
   implicit none
   character(len=*), parameter :: args_error = "Usage:"//new_line('A')//&
      '  ./Paramo params_file output_file option1 option2 option3'//new_line('A')//&
      'Options:'//new_line('A')//&
      '  params_file     Parameters file'//new_line('A')//&
      '  output_file     name of the output file'//new_line('A')//&
#ifdef TEST
#elif BLAZ
#elif AGLOW
      '  with-wind       constant or wind-like external medium:'//new_line('A')//&
#elif TURB
#elif MEZCAL
      '  with IC         with inverse Compton? T/F (True/False)'//new_line('A')//&
#elif INTERJETS
#elif MEZCAL
#endif
      '  KN cooling      Klein-Nishina cooling: T/F (True/False)'//new_line('A')//&
      '  blob            is the emission region a blob?: T/F (True/False)'//new_line('A')
   integer :: numArgs, i
   character(len=256) :: program_name
   character(len=256), allocatable, dimension(:) :: args
   logical, allocatable, dimension(:) :: with_arg

   numArgs = command_argument_count()
   !!!!! TODO: Reduce number of argmuents to 2: params_file and output_file. The rest of the arguments should be set into the parameters file
   if ( numArgs /= 5 ) call an_error(args_error)
   call get_command_argument(0, program_name)
   allocate(args(numArgs), with_arg(numArgs - 2))
   do i=1, numArgs
      call get_command_argument(i, args(i))
      write(*, *) i, args(i), numArgs, size(with_arg)
      if ( i > 2 ) then
         if ( args(i) == 'T' ) then
            with_arg(i - 2) = .true.
         else if ( args(i) == 'F' ) then
            with_arg(i - 2) = .false.
         else
            call an_error(args_error)
         end if
      end if
   end do

#ifdef TEST
#elif BLAZ
   call blazMag(trim(args(1)), args(2), with_arg(1), with_arg(2))
#elif MEZCAL
   call mezcal(trim(args(1)), args(2), with_arg(1), with_arg(2), with_arg(3))
#elif AGLOW
   call bw1D_afterglow(trim(args(1)), args(2), with_arg(1), with_arg(2), &
         with_arg(3))
#elif INTERJETS
   call interJets(trim(args(1)), args(2))
#elif TURB
   call turBlaz(trim(args(1)), args(2), with_arg(1), with_arg(2))
#endif

end program main