diff --git a/Examples/Testdata/ecoc.nas b/Examples/Testdata/ecoc.nas
new file mode 100644
index 0000000000000000000000000000000000000000..e7f7956c924a662603267d319634719c0cc221e1
--- /dev/null
+++ b/Examples/Testdata/ecoc.nas
@@ -0,0 +1,106 @@
+101 1001
+Yttri, Karl Espen; Someone, Else
+NO01L, Norwegian Institute for Air Research, NILU, , Instituttveien 18, , 2007, Kjeller, Norway
+Hjellbrekke, Anne
+EMEP GAW-WDCA
+1 1
+2018 01 01 2020 10 15
+0.166667
+days from file reference point
+17
+1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+9.999999 9999.99 9999.99 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 99.999 9.999 9.999
+end_time of measurement, days from the file reference point
+pressure, hPa
+temperature, K
+elemental_carbon, ug C/m3, Detection limit=0.4 ug C/m3
+elemental_carbon, ug C/m3, Statistics=uncertainty
+organic_carbon, ug C/m3, Detection limit=0.4 ug C/m3
+organic_carbon, ug C/m3, Statistics=uncertainty
+organic_carbon, ug C/m3, Artifact=positive
+organic_carbon, ug C/m3, Artifact=positive, Statistics=uncertainty
+organic_carbon, ug C/m3, Fraction=OC1
+organic_carbon, ug C/m3, Fraction=OC2
+organic_carbon, ug C/m3, Fraction=OC3
+organic_carbon, ug C/m3, Fraction=OC4
+organic_carbon, ug C/m3, Fraction=OCPyr
+total_carbon, ug C/m3, Detection limit=0.4 ug C/m3
+total_carbon, ug C/m3, Statistics=uncertainty
+numflag, no unit
+0
+70
+Data definition: EBAS_1.1
+Set type code: TU
+Timezone: UTC
+File name: NO0002R.20180101000000.20201015093000.online_thermal-optical_analysis.oc_ec_tc.pm25.1y.4h.NO01L_NILU_sunset_mod4_Bir.NO01L_thermal-optical_analysis_EUSAAR_2.lev2.nas
+Startdate: 20180101000000
+Revision date: 20201015093000
+Version: 1
+Version description: initial revision, manually inspected
+Statistics: arithmetic mean
+Data level: 2
+Period code: 1y
+Resolution code: 4h
+Sample duration: 225mn
+Orig. time res.: 4h
+Station code: NO0002R
+Platform code: NO0002S
+Station name: Birkenes II
+Station WDCA-ID: GAWANO__BIR
+Station GAW-ID: BIR
+Station GAW-Name: Birkenes Atmospheric Observatory
+Station land use: Forest
+Station setting: Rural
+Station GAW type: R
+Station WMO region: 6
+Station latitude: 58.38853
+Station longitude: 8.252
+Station altitude: 219.0 m
+Measurement latitude: 58.2318
+Measurement longitude: 8.1507
+Measurement altitude: 212.0 m
+Measurement height: 4.0 m
+Regime: IMG
+Component: oc_ec_tc
+Unit: ug C/m3
+Matrix: pm25
+Laboratory code: NO01L
+Instrument type: online_thermal-optical_analysis
+Instrument name: NILU_sunset_mod4_Bir
+Instrument manufacturer: Sunset Laboratory
+Instrument model: Model-4 Semi-Continuous Field Analyzer
+Instrument serial number: 9635
+Method ref: NO01L_thermal-optical_analysis_EUSAAR_2
+Standard method: TEMP=EUSAAR_2
+Inlet type: Cyclone
+Inlet description: with rain hat
+Flow rate: 8.0 l/min
+Filter face velocity: 108.4 cm/s
+Exposed filter area: 1.23 cm2
+Filter description: Circular 16.5 mm diameter punches cut from WhatmanQM-A47mm
+Medium: Quartz
+Filter prefiring: Prefired by previous analysis
+Filter conditioning: None
+Sample preparation: None
+Blank correction: Not blank corrected
+Artifact correction: Positive only
+Artifact correction description: The positive artifact of OC was accounted for by placing a denuder upstream of the instrument. The denuder collection efficiency is 91 %. The remaining positive artifact due to the inefficiency is reported in a separate variable (Artifact=positive).
+Charring correction: By laser transmission
+Humidity/temperature control: None
+Volume std. temperature: ambient
+Volume std. pressure: ambient
+Detection limit expl.: Determined by taking 3 times standard deviation of ensemble of field blanks
+Measurement uncertainty expl.: Measurement uncertainty includes quantifiable uncertainty
+Zero/negative values code: Zero/negative possible
+Zero/negative values: Zero/negative values reported as measured if within detection limit, flagged with 780
+Originator: Yttri, Karl Espen, key@nilu.no, Norwegian Institute for Air Research, NILU, Atmosphere and Climate Department, Instituttveien 18, , 2007, Kjeller, Norway
+Originator: Someone, Else, Someone@somewhere.no, Some nice Institute, WOW, Super interesting division, Street 18, , X-9999, Paradise, Noway
+Submitter: Hjellbrekke, Anne, agh@nilu.no, Norwegian Institute for Air Research, NILU, Atmosphere and Climate Department, Instituttveien 18, , 2007, Kjeller, Norway
+Acknowledgement: Request acknowledgement details from data originator
+Comment:
+starttime endtime press temp EC EC_unc OC OC_unc OC_pos OC_unc_pos OC_peak1 OC_peak2 OC_peak3 OC_peak4 OC_Pyr TC TC_unc flag
+0.000000 0.156250 997.21 251.71 0.066 0.006 0.306 0.076 0.028 0.003 0.038 0.059 0.068 0.049 0.092 0.372 0.079 0.000
+0.166667 0.322917 997.73 250.23 0.150 0.014 0.923 0.126 0.083 0.008 0.097 0.132 0.076 0.080 0.535 1.073 0.136 0.000
+0.333333 0.489583 998.43 249.88 0.086 0.008 0.430 0.100 0.039 0.004 0.063 0.080 0.073 0.059 0.156 0.516 0.103 0.000
+0.500000 0.656250 998.62 249.21 0.175 0.017 0.419 0.099 0.038 0.004 0.042 0.059 0.076 0.066 0.174 0.594 0.107 0.000
+0.666667 0.822917 9999.99 9999.99 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 9.999 99.999 9.999 0.999
diff --git a/Examples/ebas_read_example.py b/Examples/ebas_read_example.py
index 3d09051f56edadd1462a98fb54832f2dca5e4a87..4cf076de077048597adf221a67555164cc201f6f 100755
--- a/Examples/ebas_read_example.py
+++ b/Examples/ebas_read_example.py
@@ -133,6 +133,29 @@ def use_file_example(nas, time):
print(('\n found variable ozone, ug/m3: '
'var_index={}'.format(i)))
+ # ec/oc find variable witth Fraction=OC2
+ for i in nas.find_variables({
+ 'comp_name': 'organic_carbon',
+ 'characteristics': {'Fraction': 'OC2'}}):
+ print(('\n found variable organic_carbon, Fraction=OC2: '
+ 'var_index={}'.format(i)))
+ # ec/oc find all variable where characteristics Fraction is set
+ # (no matter which value)
+ for i in nas.find_variables({
+ 'comp_name': 'organic_carbon',
+ 'characteristics': {'Fraction': True}}):
+ print(('\n found variable organic_carbon, with any Fraction: '
+ 'var_index={}'.format(i)))
+ # ec/oc find variable organic_carbon, arithmetic mean, without a characteristic
+ # Fraction or Artifact specified
+ for i in nas.find_variables({
+ 'comp_name': 'organic_carbon',
+ 'statistics': 'arithmetic mean',
+ 'characteristics': {'Fraction': False, 'Artifact': False}}):
+ print(('\n found variable organic_carbon, arithmetic mean, without '
+ 'any Fraction or Artifact: var_index={}'.format(i)))
+
+
EbasCommandline(
ebas_read_example,
custom_args=['CONFIG', 'LOGGING', 'TIME_CRIT', 'NASA_READ'],