Ignore:
Timestamp:
08/02/07 16:22:22 (13 years ago)
Author:
domlowe
Message:

fixed problems with longitude subsetting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TI02-CSML/trunk/csml/csmllibs/csmldataiface.py

    r2121 r2128  
    423423                lonMax =kwargs[lonkey][1] 
    424424                kwargs[lonkey]=(0.0, lonMax) 
    425                 sel=cdms.selectors.Selector(**kwargs) 
    426                 subset1=self.file(self.varobj.id,sel) 
     425                subset1=self.file(self.varobj.id,**kwargs) 
    427426                kwargs[lonkey]=(lonMin,359.9999) 
    428                 sel=cdms.selectors.Selector(**kwargs) 
    429                 subset2=self.file(self.varobj.id,sel) 
    430                 #concatenate arrays along longitude               
    431                 longitudeAxis=subset1.getAxisIndex(lonkey) 
    432                 subset = cdms.MV.concatenate([subset1,subset2],axis=longitudeAxis) 
     427                print 'kwargs: %s'%kwargs 
     428                try: 
     429                    subset2=self.file(self.varobj.id,**kwargs) 
     430                    longitudeAxis=subset1.getAxisIndex(lonkey) 
     431                    #concatenate arrays along longitude               
     432                    subset = cdms.MV.concatenate([subset1,subset2],axis=longitudeAxis) 
     433                except: 
     434                    subset=subset1 
    433435        if type(subset) is not cdms.tvariable.TransientVariable: 
    434             sel=cdms.selectors.Selector(**kwargs) 
    435             subset=self.file(self.varobj.id,sel) 
     436            subset=self.file(self.varobj.id,**kwargs) 
    436437        return subset 
    437438 
Note: See TracChangeset for help on using the changeset viewer.