Ignore:
Timestamp:
19/03/07 09:51:59 (13 years ago)
Author:
domlowe
Message:

subset ProfileSeries to Profile, netcdf working

File:
1 edited

Legend:

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

    r2288 r2290  
    6363        self.addAxis('time',floatTimes,isTime=1,units=calunits,calendar=caltype) 
    6464        axesdone.append('time') 
    65         for ord in enumerate(ordinates): 
    66             vals=[] 
    67             lon,lat=None,None 
    68             if ord[1].coordAxisLabel.CONTENT=='time': 
    69                 continue 
    70             else: 
    71                 for val in ord[1].coordAxisValues.coordinateList.CONTENT.split(','): 
    72                     if val != ' ': 
    73                         vals.append(float(val))   
    74             if ord[1].coordAxisLabel.CONTENT=='longitude': 
    75                 lon=1 
    76                 name='longitude' 
    77             elif ord[1].coordAxisLabel.CONTENT=='latitude': 
    78                 lat=1 
    79                 name='latitude' 
    80             else: 
    81                 name=ord[1].coordAxisLabel.CONTENT 
    82             for ax in enumerate(axisorder): 
    83                 if ax[1]==name: 
    84                     position=ax[0] 
    85             axesdone.append(name) 
    86             self.addAxis(name,vals,isLon=lon,isLat=lat,units=unitlist[position])#to do, units attribute for CF compliance 
     65        if ordinates is not None: 
     66            for ord in enumerate(ordinates): 
     67                vals=[] 
     68                lon,lat=None,None 
     69                if ord[1].coordAxisLabel.CONTENT=='time': 
     70                    continue 
     71                else: 
     72                    for val in ord[1].coordAxisValues.coordinateList.CONTENT.split(','): 
     73                        if val != ' ': 
     74                            vals.append(float(val))   
     75                if ord[1].coordAxisLabel.CONTENT=='longitude': 
     76                    lon=1 
     77                    name='longitude' 
     78                elif ord[1].coordAxisLabel.CONTENT=='latitude': 
     79                    lat=1 
     80                    name='latitude' 
     81                else: 
     82                    name=ord[1].coordAxisLabel.CONTENT 
     83                for ax in enumerate(axisorder): 
     84                    if ax[1]==name: 
     85                        position=ax[0] 
     86                axesdone.append(name) 
     87                self.addAxis(name,vals,isLon=lon,isLat=lat,units=unitlist[position])#to do, units attribute for CF compliance 
    8788         
    8889        if kwargs is not None: 
Note: See TracChangeset for help on using the changeset viewer.