Commit cf33d386 authored by Martin Schrön's avatar Martin Schrön
Browse files

Bug fix for extremely high air humidities

parent b7ab003d
......@@ -73,6 +73,8 @@ index_column_number =
timestamp = no
# Separation character of the columns, e.g. \t for tab.
sep =
# Decimal letter, by default . but sometimes it can be ,
decimal = .
# Skip first lines
# Units: number of lines
skip_lines = 0
......@@ -194,11 +196,11 @@ interpolate_coords = no
neutrons_unit =
pressure_unit =
# Replace missing atmospheric data of air pressure (mbar), abs. air humidity (g/cm^3), and temperature (°C) by a constant.
# Replace missing atmospheric data of air pressure (mbar), abs. air humidity (%), and temperature (°C) by a constant.
# Type "DWD" to consult an external source to look up for a nearby German Weather Service station. Cite as: DWD Climate Data Center (CDC): Aktuelle 10-minütige Stationsmessungen des Luftdrucks, der Lufttemperatur (in 5cm und 2m Höhe), der Luftfeuchte und des Taupunkts in Deutschland, Version recent, 2019.
# Units: mbar or DWD
missing_pressure = DWD
# Units: g/m³ or DWD
# Units: % or DWD
missing_humidity = DWD
# Units: °C or DWD
missing_temperature = DWD
......@@ -225,10 +227,10 @@ neutron_n_sigma = 3
# Units: sec
timeres_range = 8, 14
# Allowed ranges (a,b) for atmospheric data of air pressure (mbar), abs. air humidity (g/cm^3), and temperature (°C):
# Allowed ranges (a,b) for atmospheric data of air pressure (mbar), rel. air humidity (%), and temperature (°C):
# Units: mbar
pressure_range = 1, 1100
# Units: g/m³
# Units: %
humidity_range = 0, 100
# Units: °C
temperature_range = -60, 60
......@@ -267,6 +269,7 @@ new_neutron_column = N
# New column name for absolute air humidity:
new_humidity_column = ah
# Method to calculate the humidity correction factor C_h (select one):
# - None # no humidity correction
# - Rosolem et al. (2013) # linear correction, = 1 + alpha*(h-h_ref)
# - Koehli et al. (2021) # polynom as a function of SM, = 1 (please select also N2sm_method = Schmidt)
humidity_method = Rosolem et al. (2013)
......
......@@ -73,6 +73,8 @@ index_column_number =
timestamp = no
# Separation character of the columns, e.g. \t for tab.
sep =
# Decimal letter, by default . but sometimes it can be ,
decimal = .
# Skip first lines
# Units: number of lines
skip_lines = 0
......
......@@ -1611,7 +1611,7 @@ def main(configfile=None, cfg_prefix=None):
if np.isnan(sm) or np.isnan(h) or np.isnan(p): return(np.nan)
if sm < 0.01: sm = 0.01
if sm > 0.49: sm = 0.49
if h > 30 : h = 30
if h > 30 : h = 29
#print(sm, h, p, int(round(100*sm)), int(round(h)))
return(footprint_data[int(round(100*sm))][int(round(h))] * 0.4922/(0.86-np.exp(-p/1013.25)))
......
Supports Markdown
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