Ignore:
Timestamp:
06/06/07 09:52:44 (15 years ago)
Author:
domlowe
Message:

netcdf output of subsetted qxf point series

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TI02-CSML/trunk/csml/API/ops_AbstractFeature.py

    r2547 r2548  
    1414 
    1515 
    16 def getDataFromChunks(self, fileList, minIdx, maxIdx): 
    17     minIndex=minIdx 
    18     maxIndex=maxIdx 
    19     data=[] 
    20     totalarraylength = 0     
    21     for file in fileList: 
    22         DI=csml.csmllibs.csmldataiface.DataInterface() 
    23         DI=DI.getUnknownInterfaceType(file)  
    24         DI.openFile(file) 
    25         DI.setVariable(self.name.CONTENT) 
    26         datachunk=DI.getDataForVar()        
    27         chunklength=len(datachunk) 
    28         startpoint=totalarraylength 
    29         totalarraylength=totalarraylength +chunklength 
    30                  
    31         if minIndex >= totalarraylength: 
    32             #print 'before selection, skipping %s'%file 
    33             continue 
    34          
    35         elif minIndex < totalarraylength: 
    36             if minIndex >= startpoint: 
    37                 sliceMin=minIndex-startpoint 
    38                 #print 'selection starts in %s at index %s'%(file, sliceMin) 
    39                 if maxIndex < totalarraylength: 
    40                     #print 'selection self contained in %s'%file 
    41                     sliceMax=sliceMin + (maxIndex-minIndex) +1 
    42                     data=datachunk[sliceMin:sliceMax] 
    43                     continue 
    44                 else: 
    45                     #print 'selection goes beyond this file' 
    46                     data.append(datachunk[sliceMin: len(datachunk)]) 
    47                     continue 
    48          
    49         if maxIndex  < startpoint: 
    50             #print 'gone past end of selection, skipping %s'%file 
    51             continue 
    52         elif maxIndex >= totalarraylength: 
    53             #get the whole chunk 
    54             data.append(datachunk) 
    55             continue 
    56          
    57         elif maxIndex <= totalarraylength:      
    58             sliceMax=(maxIndex-startpoint) +1 
    59             data.append(datachunk[0:sliceMax]) 
    60             continue     
    61     try: 
    62         result=MA.concatenate(data) 
    63     except: 
    64         result =data 
    65     return result 
     16 
    6617     
    6718 
Note: See TracChangeset for help on using the changeset viewer.