Changeset 1968 for TI02-CSML/trunk


Ignore:
Timestamp:
08/01/07 14:49:18 (13 years ago)
Author:
domlowe
Message:

crs idenfication tests real model data

File:
1 edited

Legend:

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

    r1966 r1968  
    11import string 
     2import cdms # just while testing 
    23 
    34class CRSystem(object): 
     
    8182    print 'srsDimension = %s' %crs.srsDimension 
    8283    print 'axisLabels = %s' %crs.axisLabels 
     84    
     85    
     86    
     87    crsdictionary={} 
     88    f=cdms.open('/home/dom/Desktop/SVN/trunk/csml/testfiles/gridseries/xaaqda@pxs19c1.nc') 
     89    for var in f.variables: 
     90        v=f.variables[var] 
     91        axes=v.getAxisIds() 
     92 #   print axes 
     93        units =[] 
     94        for axis in axes: 
     95            ax=f.getAxis(axis) 
     96            units.append(ax.attributes['units']) 
     97#    print units 
     98        crsdictionary[str(axes)]=units 
     99    print crsdictionary 
     100    n=2 
     101    #test each definition in dictionary 
     102    for key in crsdictionary: 
     103        n=n+1 
     104        print 'TEST %d'%n 
     105        print '*************************************' 
     106        axs =key 
     107        uns=crsdictionary[key] 
     108        print 'AXES: %s'%axs 
     109        print 'UNITS: %s'%uns 
     110        crs= cat.getCRS(axes=axs, units=uns) 
     111        print 'srsName = %s' %crs.srsName 
     112        print 'srsDimension = %s' %crs.srsDimension 
     113        print 'axisLabels =%s' %crs.axisLabels 
     114        print '**************************************' 
     115    
     116    
     117    
     118     
     119     
    83120 
    84121if __name__=="__main__": 
Note: See TracChangeset for help on using the changeset viewer.