scanMap.f90

  TYPE(m3_f77_runconfig) :: runConfig
  TYPE(m3_f77_dataset) :: dataSet
  TYPE(m3_f77_interval) :: intervalNode
  integer(i4b) :: storeSize
  integer(i4b), dimension(:), allocatable :: storeSpace
  real(dp), dimension(:,:), allocatable :: pointing

  M3_f77_RunConfig_Read_mpi( runConfig, "runConfig.xml", MPI_COMM_WORLD)

  errp = M3_f77_runconfig_getpixelclassroot( runconfig, pixelclass )
  i = 0
  do while( errp == 1 )
    M3_pixelclass_getnumpixelinmap( pixelclass, thisnumpixel )
    if( numpixel == -1 ) then
      numpixel = thisnumpixel
    else if( numpixel /= thisnumpixel ) then
      print*, "ERROR number of pixels in maps is not uniform"
      stop
    end if 
    errp = M3_f77_pixelclass_getnextnodeinlist( pixelclass )
    i = i + 1
  end do

  if( i /= 18 ) then
    print*, "ERROR number of pixel classes must be 18"
    stop
  end if

  allocate( mapcmb( 0 : numpixel - 1, 0 : 17 ) )
!  allocate( mapfgd150( 0 : numpixel - 1, 0 : 17 ) )
!  allocate( mapfgd250( 0 : numpixel - 1, 0 : 17 ) )
!  allocate( mapfgd410( 0 : numpixel - 1, 0 : 17 ) )

  errp = M3_f77_runconfig_getpixelclassroot( runconfig, pixelclass )
  i = 0
  do while( errp == 1 )
    M3_f77_pixelclass_getfullmap( pixelclassnode,  pixel, mapcmb(:,i) )
    errp = M3_f77_pixelclass_getnextnodeinlist( pixelclass )
  end do

  M3_f77_RunConfig_GetGCPointingGroupRoot( runConfig, gcpointinggroup )

  numnz = -1
  maxIntervalLength = 0
  ! loop over datasets (detectors)
  errd = M3_f77_RunConfig_getdatasetroot( runConfig, dataset )
  do while( errd == 1 )
    ! loop over intervals 
    err = M3_f77_dataset_getcoveredintervalroot( dataset, interval )
    do while( err == 1 )
      M3_f77_datasetll_sampleintervaltotimeinterval( dataset, interval, timeinterval )
      M3_f77_interval_getinterval( interval, firstsample, lastsample )
      if( lastsample - firstsample + 1 > maxIntervalLength ) then
        maxIntervalLnegth = lastsample - firstsample + 1
      end if

      M3_f77_dataset_getnumnz( dataset, thisnumnz )
      if( numnz == -1 ) then
        numnz = thisnumnz
      else if( thisnumnz /= numnz ) then
        print*, "ERROR: number of non-zeros is not uniform over the datasets."
        stop
      end if
      M3_f77_timeinterval_insertinterval( timeintervalList, timeinterval )
      err = M3_f77_Interval_getnextnodeinlist( interval )
    end do
    errd = M3_f77_datasetll_getnextnodeinlist( dataset )
  end do
  
  M3_f77_gcpointinggroup_getstoresize( gcpointinggroup, timeIntervalList, storeSize )
  storeSize = int(storeSize/4) + 1
  allocate( storeSpace( 0 : storeSize - 1), stat=ios )
  M3_f77_gcpointinggroup_initializestore( gcpointinggroup, timeintervallist, storeSpace )

  allocate( pointing( 0 : numnz - 1, 0 : maxintervallength - 1), stat = ios )
  allocate( tod( 0 : maxintervallength - 1), stat = ios )
  allocate( pixel( 0 : 17, 0 : maxintervallength - 1), stat = ios )
  allocate( weight( 0 : 17, 0 : maxintervallength - 1), stat = ios )

  ! loop over datasets
  errd = M3_f77_RunConfig_getdatasetroot( runConfig, dataset )
  firstrun = 1
  do while( errd == 1 )
    ! loop over intervals 
    err = M3_f77_dataset_getcoveredintervalroot( dataset, interval )
    do while( err == 1 )
     tod = 0.d0
     M3_dataset_getpointing( dataset, interval, pixel, weight )
     m3_interval_getinterval( interval, firstsample, lastsample )
     do i = firstsample, lastsample
       do j = 0, 17
         tod(i) = tod(i) + weight(j,i)*map(pixel(j,i), j )
       end do
    end do
    if ( firstrun == 1) then
       file = "tod.bin"
       call output_tl( lastsample-firstsample+1, tod, file)
    end if
    firstrun = 0

    err = M3_f77_Interval_getnextnodeinlist( interval )
    end do
    errd = M3_f77_datasetll_getnextnodeinlist( dataset )
  end do

---