maxMinTOD.f90

A Simple M3 example in Fortran of TOD processing which derives the maximum and minumum value from all of the timestreams in a Run Config.

!first try, Columbia, Nov. 18th
!Nicolas Ponthieu

program minmaxtod
  
  use m3_typesf
  use m3_f
  implicit none

  type(m3_f77_runconfig) :: runconfig
  type(m3_f77_dataset) :: dataset
  type(m3_f77_interval) :: intervalnode !intervallist

  character(len=200) :: runconfig_file

  integer :: m3err, m3err1, m3err2, n_intervals, i_int, ios
  INTEGER(KIND=int8) ::firstsample, lastsample, n_samples, i

  !! m3err = 1 means "fine"

  real(kind=double), allocatable, dimension(:) :: tod
  real(kind=double) :: min, max

  !--------------------------------------------------------------------------
  max = -100000
  min =  100000

  runconfig_file = "madmap_test_littleendian_small.xml"
  

  !! Reading the runconfig_file
  m3err = m3_f77_runconfig_read( runconfig, runconfig_file, m3_springtide_run_type)
  
  !Get the first dataset (= detector)
  m3err = m3_f77_runconfig_getdatasetroot( runconfig, dataset)

  !!Loop over the data sets (descripted in the xml file)
  do while (m3err == 1)
     
     !get 1st interval in the list
     m3err2 = m3_f77_dataset_getcoveredintervalroot( dataset, intervalnode)
     
     do while (m3err2==1)
        
        !get this interval
        m3err2 = m3_f77_interval_getinterval( intervalnode, firstsample, lastsample )
        n_samples = lastsample - firstsample + 1

        !allocate memory
        allocate( tod(0:n_samples-1), stat=ios)
        if (ios /=0) then
           print*, "allocation error"
           stop
        endif
        
        !read the data in that interval
        m3err2 = m3_f77_dataset_gettod( dataset, firstsample, lastsample, tod )
        
        do i=0, n_samples-1
           if (tod(i) < min) min = tod(i)
           if (tod(i) > max) max = tod(i)
        enddo
        
        deallocate( tod)
        
        m3err2 = m3_f77_interval_getnextnodeinlist( intervalnode )
     enddo

     !to go the next dataset in the list
     m3err = m3_f77_dataset_getnextnodeinlist( dataset)
  enddo

  print*, "min = ", min
  print*, "max = ", max

end program minmaxtod

---