Commit 6741169f authored by David Schäfer's avatar David Schäfer
Browse files

make new config work

parent e17a0173
......@@ -7,21 +7,21 @@ varname ; test
# quality control
# ---------------
N1 ; outliers.flagRange(min=100, max=3000)
N1ET ; outliers.flagRange(min=1, max=10000)
p1 ; outliers.flagRange(min=1, max=1100)
p4 ; outliers.flagRange(min=1, max=1100)
RH1 ; outliers.flagRange(min=0, max=100)
N1RH ; outliers.flagRange(min=0, max=100)
RHX ; outliers.flagRange(min=0, max=100)
T1 ; outliers.flagRange(min=-60, max=60)
TX ; outliers.flagRange(min=-60, max=60)
N1T ; outliers.flagRange(min=-60, max=60)
Vbat ; outliers.flagRange(min=10, max=20)
N1Cts ; outliers.flagRange(min=100, max=3000)
N1ET_sec ; outliers.flagRange(min=1, max=10000)
p1 ; outliers.flagRange(min=1, max=1100)
p4_mb ; outliers.flagRange(min=1, max=1100)
RH1 ; outliers.flagRange(min=0, max=100)
N1RH ; outliers.flagRange(min=0, max=100)
RHX ; outliers.flagRange(min=0, max=100)
T1 ; outliers.flagRange(min=-60, max=60)
T_CS215 ; outliers.flagRange(min=-60, max=60)
N1T ; outliers.flagRange(min=-60, max=60)
Vbat ; outliers.flagRange(min=10, max=20)
# convert units
# -------------
N ; generic.process(func= N1 / N1ET * 3600)
N ; generic.process(func= N1Cts / N1ET_sec * 3600)
N_cleaned ; generic.process(func=N)
# quality control - new variables
......@@ -31,26 +31,22 @@ N_cleaned ; outliers.flagRange(min=500, max=5000)
# data corrections/processing
# ---------------------------
ah ; generic.process(func=6.112 * exp(17.67 * TX/(243.5 + TX))/(273.15 + TX) * 2.1674 * RHX)
correct_h ; generic.process(func=1 + 0.0054 * (ah - 0))
correct_p ; generic.process(func=exp((p4 - 1013.25) / 136))
correct_inc ; generic.process(func=1/(1 - 1 * (1 - NM / 150)))
N_corrected ; generic.process(func=N * correct_p * correct_h * correct_inc)
ah ; generic.process(func=6.112 * exp(17.67 * T_CS215/(243.5 + T_CS215))/(273.15 + T_CS215) * 2.1674 * RHX)
correct_h ; generic.process(func=1 + 0.0054 * (ah - 0))
correct_p ; generic.process(func=exp((p4_mb - 1013.25) / 136))
correct_inc ; generic.process(func=1/(1 - 1 * (1 - NM / 150)))
N_corrected ; generic.process(func=N * correct_p * correct_h * correct_inc)
# soil moisture calculation
# -------------------------
SM ; calcSoilMositure(neutrons_field="N_corrected", n0=5100, lattice_water=0.0043, bulk_density=1.6, soil_org_carbon=0.0050)
SM ; outliers.flagRange(min=0, max=0.7)
SM ; calcSoilMositure(neutrons_field="N_corrected", n0=5100, lattice_water=0.0043, bulk_density=1.6, soil_org_carbon=0.0050)
SM ; outliers.flagRange(min=0, max=0.7)
# rename variables to match the DMP names
# ---------------------------------------
N1 ; tools.copy(new_field="N1Cts")
N1ET ; tools.copy(new_field="N1ET_sec")
p4 ; tools.copy(new_field="P4_mb")
RH1 ; tools.copy(new_field="h1")
RHX ; tools.copy(new_field="RH_CS215")
T1 ; tools.copy(new_field="t1")
TX ; tools.copy(new_field="T_CS215")
Vbat ; tools.copy(new_field="bat")
# flag everything
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment