Ignore:
Timestamp:
03/12/07 15:10:24 (13 years ago)
Author:
domlowe
Message:

subsetting of profileseries, pointseries retested and fixed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TI02-CSML/trunk/csml/testfiles/gridseries/subset.py

    r2974 r3096  
    5757#csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',latitude=(-90,90), longitude=(-180,180),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0'))) 
    5858 
    59 csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',latitude=(60), longitude=(180),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0'))) 
     59#csmlContainer.add(feature.subsetToGridSeries(outputdir,ncname='mygridseries.nc',latitude=(60), longitude=(180),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0'))) 
    6060 
    61 #subset another feature (in fact it's the same feature again, but this time subsetting to a profileseries) 
    62 subsetDictionary={} 
    63 subsetDictionary['time']=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0') 
    64 subsetDictionary['latitude']=(75) 
    65 subsetDictionary['longitude']=(45) 
     61##subset another feature (in fact it's the same feature again, but this time subsetting to a profileseries) 
     62#subsetDictionary={} 
     63#subsetDictionary['time']=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0') 
     64#subsetDictionary['latitude']=(75) 
     65#subsetDictionary['longitude']=(45) 
    6666 
    67 ##assigning profileseries feature to variable 
    68 #print 'subset grid to profileseries' 
    69 psfeature, file, descriptor =feature.subsetToProfileSeries(outputdir, ncname='myprofileseries.nc',**subsetDictionary) 
    70 csmlContainer.add([psfeature,file, descriptor]) 
     67###assigning profileseries feature to variable 
     68##print 'subset grid to profileseries' 
     69#psfeature, file, descriptor =feature.subsetToProfileSeries(outputdir, ncname='myprofileseries.nc',**subsetDictionary) 
     70#csmlContainer.add([psfeature,file, descriptor]) 
    7171 
    7272#csmlContainer.add(feature.subsetToProfileSeries(outputdir,ncname='myprofileseries.nc',latitude=(75), longitude=(45),time=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0'))) 
    7373 
    74 #Now subset this New ProfileSeries feature into a Profile feature 
     74##Now subset this New ProfileSeries feature into a Profile feature 
    7575#print 'subset  profileseries to profile' 
    7676#subsetDictionary={} 
    77 #subsetDictionary['time']=('2794-12-1T0:0:0.0') 
    78 #csmlContainer.add(psfeature.subsetToProfile(outputdir, ncname='myprofile1.nc',**subsetDictionary)) 
     77#subsetDictionary['time']=('2804-12-01T00:00:00.0') 
     78#csmlContainer.add(psfeature.subsetToProfile(outputdir, ncname='myprofile.nc',**subsetDictionary)) 
    7979 
    8080 
    8181#Now try the subset to profile method of the grid series feature -BROKEN 
    82 #subsetDictionary={} 
    83 #subsetDictionary['time']=('2804-12-1T0:0:0.0') 
    84 #subsetDictionary['latitude']=(80) 
    85 #subsetDictionary['longitude']=(55) 
    86 #csmlContainer.add(feature.subsetToProfile(outputdir, ncname='myprofile.nc',**subsetDictionary)) 
     82print 'subsetting gridseries to profile' 
     83subsetDictionary={} 
     84subsetDictionary['time']=('2804-12-01T00:00:00.0') 
     85subsetDictionary['latitude']=(80) 
     86subsetDictionary['longitude']=(55) 
     87csmlContainer.add(feature.subsetToProfile(outputdir, ncname='myprofile.nc',**subsetDictionary)) 
    8788 
    8889#Now try the subset to PointSeries method of the grid series feature 
    89 #subsetDictionary={} 
    90 #subsetDictionary['time']=('2794-12-1T0:0:0.0','2904-12-1T0:0:0.0') 
    91 #subsetDictionary['latitude']=(80) 
    92 #subsetDictionary['longitude']=(55) 
    93 ##subsetDictionary['pressure']=(0)  #BROKEN 
    94 #csmlContainer.add(feature.subsetToPointSeries(outputdir,ncname='mypointseries.nc',**subsetDictionary)) 
     90print 'subsetting gridseries to pointseries' 
     91subsetDictionary={} 
     92subsetDictionary['time']=('2804-12-01T00:00:00.0','2904-12-01T00:00:00.0') 
     93subsetDictionary['latitude']=80 
     94subsetDictionary['longitude']=55 
     95#subsetDictionary['pressure']=(0)  #BROKEN 
     96csmlContainer.add(feature.subsetToPointSeries(outputdir,ncname='mypointseries.nc',**subsetDictionary)) 
    9597 
    9698 
Note: See TracChangeset for help on using the changeset viewer.