Skip to content
Snippets Groups Projects
Commit 48b6bd1a authored by weihe's avatar weihe
Browse files

co_filter keys

parent 42dd26fc
Branches
No related tags found
No related merge requests found
...@@ -338,7 +338,7 @@ class ObsPackObservations(Observations): ...@@ -338,7 +338,7 @@ class ObsPackObservations(Observations):
def add_model_data_mismatch(self, filename, filename_cofilter): def add_model_data_mismatch(self, filename):
""" """
Get the model-data mismatch values for this cycle. Get the model-data mismatch values for this cycle.
...@@ -383,6 +383,7 @@ class ObsPackObservations(Observations): ...@@ -383,6 +383,7 @@ class ObsPackObservations(Observations):
site_hourly = {} # option added to include only certain hours of the day (for e.g. PAL) IvdL site_hourly = {} # option added to include only certain hours of the day (for e.g. PAL) IvdL
site_foot = {} # option added to check for available footprints per observation site_foot = {} # option added to check for available footprints per observation
site_incalt = {} # option to increase sampling altitude for sites specified in sites and weights file site_incalt = {} # option to increase sampling altitude for sites specified in sites and weights file
co_filter = {}
for key, value in sites_weights.iteritems(): for key, value in sites_weights.iteritems():
if 'co2_' in key or 'sf6' in key: # to be fixed later, do not yet know how to parse valid keys from rc-files yet.... WP if 'co2_' in key or 'sf6' in key: # to be fixed later, do not yet know how to parse valid keys from rc-files yet.... WP
sitename, sitecategory = key, value sitename, sitecategory = key, value
...@@ -401,13 +402,16 @@ class ObsPackObservations(Observations): ...@@ -401,13 +402,16 @@ class ObsPackObservations(Observations):
if 'site.incalt' in key: if 'site.incalt' in key:
identifier, incalt = value.split(';') identifier, incalt = value.split(';')
site_incalt[identifier.strip()] = (int(incalt)) site_incalt[identifier.strip()] = (int(incalt))
if 'co.filter' in key:
identifier_cofilter, cofilter = value.split(';')
co_filter[identifier_cofilter.strip()] = (str(cofilter))
#for obs in self.datalist: #for obs in self.datalist:
do_not_simulate=[] do_not_simulate=[]
do_simulate=[] do_simulate=[]
eventids=[] eventids=[]
for line in open(filename_cofilter): for line in open(str(cofilter)):
columns = line.split() columns = line.split()
if columns[0] == "eventid": if columns[0] == "eventid":
pass pass
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment