Changeset 1072


Ignore:
Timestamp:
31/05/06 13:30:14 (13 years ago)
Author:
domlowe
Message:

now using cdtime for time conversion #2

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TI02-CSML/trunk/Scanner/csmllibs/csmltime.py

    r1071 r1072  
    117117        return FormattedFileTimeList     
    118118 
    119 def OLDgetFileTimeList(filelist, timedim): 
    120          
    121         #**************************************************************************************** 
    122         #Given a list of files and a name for the time dimension this function will return a list containing 
    123         # filename/time pairs eg tuple [(ta,fa),(tb,fb),(tc,fc)] in the correct temporal order. 
    124         #if files contain more than one time, then the tuple will be of the form: 
    125         #[(ta1,ta2,ta3,fa),(tb1,tb2,tb3,fb),(tc1,tc2,tc3,fc)] 
    126         #***************************************************************************************** 
    127         FileTimeList = [] 
    128  
    129         print "len of filelist = " + str(len(filelist)) 
    130         DI=csmllibs.csmldataiface.DataInterface() 
    131         #grab any (first) file                   
    132          
    133         DI=DI.getUnknownInterfaceType(filelist[0]) 
    134         for i in range(0,len(filelist)): 
    135         #for i in range(0,20): 
    136                 #print "opening" + filelist[i] 
    137                 #print timedim 
    138                 DI.openFile(filelist[i]) 
    139                 #print DI.getListOfAxes() 
    140                 #print DI.getListofVariables() 
    141                  
    142                 try: 
    143                         DI.setAxis(timedim) 
    144                 except: 
    145                         print "Could not set time axis" 
    146                         print "check you have the correct time axis (-t) set" 
    147                         sys.exit() 
    148                 data =DI.getDataForAxis() 
    149                 #fileobj=cdms.open( filelist[i]) 
    150                 #var = fileobj[timedim] 
    151                 #data=var.getValue() 
    152                 times=[] 
    153                 for j in range(0,len(data)): 
    154                         times.append(data[j]) 
    155                 #ftTuple = ( data[j], filelist[i]) 
    156                 ftTuple = (times, filelist[i]) 
    157                 FileTimeList.append(ftTuple) 
    158                 DI.closeFile() 
    159                 #fileobj.close() 
    160                 if i % 20 == 0: 
    161                         print i 
    162         FileTimeList.sort() 
    163         #print "FTLIST" 
    164         #print FileTimeList 
    165         #convert times to CSML format times 
    166         #grab any (first) file                   
    167         #fileobj=cdms.open(filelist[0]) 
    168         DI.openFile(filelist[i]) 
    169         #get the reference date for the time reference system. 
    170         DI.setAxis(timedim) 
    171         #var = fileobj[timedim] 
    172 #TO DO ERROR: if this line errors it's probably that user hasn't specified time dimension (-t) 
    173         #refTime=DI.getVariableAttribute('units') 
    174  
    175         refTime=DI.getTimeUnits() 
    176         #refTime = getattr(var,'units') 
    177         #remove characters:  'days since'   
    178         # to do - make this work for hours.minutes since too! 
    179         refTime = refTime[11:len(refTime)] 
    180         #print "REFTIME2" 
    181         #print refTime 
    182         FormattedFileTimeList = [] 
    183         for j in range (0,len(FileTimeList)): 
    184                 timeList = FileTimeList[j][0] 
    185                 #print 'timelist_____________________' +str(timeList) 
    186                 for k in range(0, len(timeList)): 
    187                         uduTime = udunitsTime(timeList[k], refTime) 
    188                         ftTuple = ( uduTime, filelist[j]) 
    189                         FormattedFileTimeList.append(ftTuple) 
    190         return FormattedFileTimeList     
    191119                 
    192120def GetOriginPosition(timestr): 
Note: See TracChangeset for help on using the changeset viewer.