bianchi2_lut_gen.f90

Go to the documentation of this file.

!----------------------------------------------------------------------------
! bianchi2_lut_gen -- BIANCHI2 writting LUT program
!
!----------------------------------------------------------------------------

program bianchi2_lut_gen

  use s2_types_mod
  use bianchi2_lut_mod

  implicit none

  type(bianchi2_lut) :: lut
  integer :: lmax, ntheta
  character(len=S2_STRING_LEN) :: filename_lut

  !Get the arguments.
  call parse_options

  ! Initialize the lut.
  lut = bianchi2_lut_init(lmax,ntheta)

  ! Write the lut into a file.
  call bianchi2_lut_write(lut,filename_lut)
  write(*,'(a)') 'LUT written in '//trim(filename_lut)


  contains

  !------------------------------------------------------------------------
  ! Parse options
  !
  !------------------------------------------------------------------------

    subroutine parse_options

      use extension, only: getargument, narguments
     
      implicit none

      integer :: n, i
      character(len=S2_STRING_LEN) :: arg
      character(len=S2_STRING_LEN) :: opt
      
      n = narguments()
     
      do i = 1,n,2
        
         call getargument(i,opt)
     
         if (i == n .and. trim(opt) /= '-help') then
            write(*,'(a,a,a)') 'Option ', trim(opt), ' has no argument'
            stop
         end if
     
         if(trim(opt) /= '-help') call getargument(i+1,arg)
         
         ! Read each argument in turn
         select case (trim(opt))
            
         case ('-help')
            write(*,'(a)') 'Usage: bianchi2_lut_gen [-lmax lmax]'
            write(*,'(a)') '                        [-Ntheta Ntheta]'
            write(*,'(a)') '                        [-out filename_LUT]'
            stop
            
         case ('-lmax')
            read(arg,*) lmax

         case ('-Ntheta')  
            read(arg,*) ntheta
            
         case ('-out')
            filename_lut = trim(arg)

         case default
            print '("Unknown option ",a," ignored")', trim(opt)            
            
         end select
      end do

    end subroutine parse_options
    

end program bianchi2_lut_gen