예제 #1
0
	def __init__(self):
		catalog_url_head = "http://opendap.knmi.nl/knmi/thredds/catalog/"
		catalog_url_proj = "CLIPC/tier1_indicators/icclim_cerfacs"

		catalog_url      = catalog_url_head+catalog_url_proj+"/catalog.xml"

		wmsStart = 'https://climate4impact.eu/impactportal/ImpactService?source='
		#catalog = 'http://opendap.knmi.nl/knmi/thredds/catalog/'
		opendap = 'http://opendap.knmi.nl/knmi/thredds/dodsC/'
		wmsEnd = '&service=WMS&'



		ReadCatalog.__init__(self , catalog_url_head ,catalog_url_proj)


		print catalog_url

		links , ncs = parse_c4i_catalog.readcatalog4(catalog_url,opendap,wmsStart,wmsEnd )

		self.print_list_to_file('knmi.wms',links)
		self.print_list_to_file('knmi.ncs',ncs)

		print len(ncs)

		#class UniDict():

		print ""

		uni = {}
		uni['variable'] = []
		uni['experiment'] = []
		uni['model'] = []
		uni['ensemble'] = []
		uni['regional'] = []
		#uni['8'] = []

		for n in ncs:
			line = n.split('/')[11].split('_')
			print line
			if line[0] not in uni['variable']:
				uni['variable'].append(line[0])
			elif line[2] not in uni['model']:
				uni['model'].append(line[2])	
			elif line[3] not in uni['experiment']:
				uni['experiment'].append(line[3])	
			elif line[4] not in uni['ensemble']:
				uni['ensemble'].append(line[4])
			elif line[5] not in uni['regional']:
				uni['regional'].append(line[5])	
			#elif line[8] not in uni['8']:
			#	uni['8'].append(line[8])

		for k in uni.keys():
			print k+": "+str(uni[k])

		print ""
		self.printSelf()
예제 #2
0
	def __init__(self):
		catalog_url_head = "http://thredds.met.no/thredds/dodsC/"
		catalog_url_proj = "arcticdata/met.no/CLIPC"
		catalog_url      = catalog_url_head+catalog_url_proj+"/catalog.xml"

		wmsStart = 'http://thredds.met.no/thredds/wms/'#'https://climate4impact.eu/impactportal/ImpactService?source='
		#catalog = 'http://opendap.knmi.nl/knmi/thredds/catalog/'
		#opendap = 'http://opendap.knmi.nl/knmi/thredds/dodsC/'
		wmsEnd = '?service=WMS&version=1.3.0&request=GetCapabilities'



		ReadCatalog.__init__(self , catalog_url_head ,catalog_url_proj)


		layer = ''

		links , ncs = parse_c4i_catalog.readcatalog3B(catalog_url,catalog_url_head,wmsStart,wmsEnd , layer)

		self.print_list_to_file('met.no.wms',links)
		self.print_list_to_file('met.no.ncs',ncs)

		print len(ncs)

		#class UniDict():

		print ""

		uni = {}
		uni['variable'] = []
		uni['experiment'] = []
		uni['model'] = []
		uni['ensemble'] = []
		uni['regional'] = []
		#uni['8'] = []

		for n in ncs:
			line = n.split('_')
			#print line
			if line[1] not in uni['variable']:
				uni['variable'].append(line[1])
			elif line[3] not in uni['model']:
				uni['model'].append(line[3])	
			elif line[4] not in uni['experiment']:
				uni['experiment'].append(line[4])	
			elif line[5] not in uni['ensemble']:
				uni['ensemble'].append(line[5])
			elif line[6] not in uni['regional']:
				uni['regional'].append(line[6])	
			#elif line[8] not in uni['8']:
			#	uni['8'].append(line[8])

		for k in uni.keys():
			print k+": "+str(uni[k])

		print ""
		self.printSelf()