Ignore:
Timestamp:
01/02/07 10:43:07 (13 years ago)
Author:
domlowe
Message:

making coordinate reference systems and netcdf output CF compliant

Location:
TI02-CSML/trunk/csml/API
Files:
2 edited

Legend:

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

    r2083 r2091  
    66 
    77 
    8 def _subsetDomain(time,times, domain,**kwargs): 
     8def _subsetDomain(timeaxis,times, domain,**kwargs): 
    99    '''takes the domain and returns a subset according to the keyword selection criteria 
    1010    time = name of time dimension 
     
    1717        straxisValues='' 
    1818        if key in kwargs: 
    19             if key ==time: 
     19            if key ==timeaxis: 
    2020                straxisValues=strTimes 
    2121                arraySize=len(strTimes.split()) 
     
    4444                                    
    4545        subsetDomain[key]=straxisValues 
    46         if key != time: 
     46        if key != timeaxis: 
    4747            arraySize=len(subsetDomain[key].split(',')) 
    4848        else: 
  • TI02-CSML/trunk/csml/API/ops_GridSeriesFeature.py

    r2084 r2091  
    136136        time=csml.API.ops_AbstractFeature.__getCDtime(time).torel(calunits) 
    137137        floatTimes.append(time.value) 
    138     nc.addAxis('t',floatTimes,isTime=1,units=calunits,calendar=caltype) 
     138    nc.addAxis('time',floatTimes,isTime=1,units=calunits,calendar=caltype) 
     139     
    139140    for ord in ords: 
    140         print ord.coordAxisLabel.CONTENT 
    141141        vals=[] 
    142142        lon,lat=None,None 
    143         if ord.coordAxisLabel.CONTENT=='Time': 
     143        if ord.coordAxisLabel.CONTENT=='time': 
    144144            continue 
    145145        else: 
     
    147147                if val != ' ': 
    148148                    vals.append(float(val))   
    149                     print vals 
    150         if ord.coordAxisLabel.CONTENT=='Lon': 
     149        if ord.coordAxisLabel.CONTENT=='longitude': 
    151150            lon=1 
    152151            name='longitude' 
    153         elif ord.coordAxisLabel.CONTENT=='Lat': 
     152        elif ord.coordAxisLabel.CONTENT=='latitude': 
    154153            lat=1 
    155154            name='latitude' 
     
    159158        nc.addAxis(name,vals,isLon=lon,isLat=lat,units='')#to do, units attribute for CF compliance 
    160159    if len(ords)==3: 
    161         axes=['t',axisorder[1],axisorder[2]] 
     160        axes=['time',axisorder[1],axisorder[2]] 
    162161    elif len(ords)==2: 
    163         axes=['t',axisorder[1]] 
    164     #print fulldata 
    165     print axes 
     162        axes=['time',axisorder[1]] 
    166163    nc.addVariable(fulldata,self.id, axes,units='') #to do, units attribute for CF compliance 
    167     print 'added' 
    168164    nc.closeFinishedFile() 
    169165    return subsettedFeature, pathToSubsetNetCDF 
Note: See TracChangeset for help on using the changeset viewer.