Resampling very high resolution data 4Mhz to 10 Hz

Function pops_10Hz_grabber(pops_time, size_2d, time_grid)

 
    wave pops_time, size_2d,time_grid // time grid is a perfact 10hz timeseries that I use to average high resolution data of pops_time and size_2d.
    make/o/d/n=(numpnts(pops_time)) POPSTStore=nan
    make/o/d/n=(numpnts(pops_time),11) POPSTStore_size=nan
    variable i,k,j    
    variable cnt
    variable POPSStart, POPSEnd
    j=0

 
    for(i=0;i<numpnts(time_grid)-1;i+=1)
       variable chckVal = 0
       variable startT= pops_time[i] // here diff betweeen startT and endT is 0.1 sec always.
         variable endT = pops_time[i+1]

 
        //print i
       k = j
        do

 
        if(pops_time[k] >= startT && chckVal == 0 )
                POPSStart=k
                chckVal = 1 // once it finds the match chacval is 1 so the if statement doesn't fullfill and moves on to the next point
               //print POPSStart
                elseIf(chckVal == 1 && pops_time[k] > endT) // same way finds the endtime here in the high res timewave
              POPSEnd =k-1 //
              j=k          // saving the position of where endT was found in J, so the next loop start from that position
              k=numpnts(pops_time) // forcing k to be max value of the loop here and later adding +1 to it so I can come out of the loop
              endif
              k+=1
              while(k<numpnts(pops_time))

 

 
        duplicate/d/o/R=[POPSSTart, POPSEnd] pops_time, tempT // extracting data that fall into time gap  of 0.1 sec.
       wavestats/q tempT
      //      print V_npnts
        POPSTStore[i] = V_avg
     duplicate/d/o/R=[POPSSTart, POPSEnd][] POPSTStore_size tempP // not sure how to extract 2d data but I need to that too.

 
      //matrixop/o test_data=POPSTStore

 
    endFor

 
End