source: TI02-CSML/trunk/csml/testfiles/apicalls.py @ 2207

'''apicalls.py  - test the api'''

import csml,sys

f='./testfiles/gridseries/testout.xml'

#Initialise the CSML container
csmlContainer=csml.API.csmlContainer.Container(csmlpath='myfile.xml',docID="mycsmldoc_1", metadata="mymetadatastring")

#Initialise and parse the dataset
csmldoc = csml.parser.Dataset()
csmldoc.parse(f)


#get list of features in the dataset
flist= csmldoc.getFeatureList() 
print '\n Here are all the features in %s:' %f
print flist

#select a feature by name (gml:id)
print '\n Selecting feature with gml:id = %s' %flist[4]
feature=csmldoc.getFeature(flist[4])

#feature 'bestname'
print feature.id
#print feature.description

#get the domain of the feature

print '\n The feature has domain:' 
print feature.getDomain()

#get list of allowed subsettings
print '\n the following feature subsetting operations are allowed:'
print feature.getAllowedSubsettings()


#define a selection (based on the domain)
subsetDictionary={}
subsetDictionary['t']=('2794-12-1T0:0:0.0','2795-12-1T0:0:0.0')
#subsetDictionary['t']=()
subsetDictionary['latitude']=(-90,60)
subsetDictionary['longitude']=(-31.265,30.258)

#request subsetted data from features (can set paths here) and add to container
#subset a feature
csmlContainer.add(feature.subsetToGridSeries(ncpath='myfile.nc',**subsetDictionary))
#subset another feature (in fact it's the same feature again, but it proves the concept)
subsetDictionary['t']=('2794-12-1T0:0:0.0','2795-12-1T0:0:0.0')
subsetDictionary['latitude']=(75)
subsetDictionary['longitude']=(41.25,45)
csmlContainer.add(feature.subsetToProfileSeries(ncpath='myfile2.nc',**subsetDictionary))

#when you have finished, get the contents of the container
csmlbundle=csmlContainer.getContents()
print 'CSML %s'%csmlbundle[0] # csml document
paths=csmlbundle[1:]
print paths #paths to netcdf files