# -*- coding: utf-8 -*-
"""
Created on Mon Jun 7 09:30:26 2021
@author: Daniel
"""
import hydrogeosines as hgs
import numpy as np
import pandas as pd
#%% Testing MVC principal
fowlers = hgs.Site('Fowlers Gap', geoloc=[141.73099, -31.2934, 160])
fowlers.import_csv('tests/data/fowlers_gap/acworth_all.csv',
input_category=['BP', 'GW', 'GW', 'GW', 'ET'],
utc_offset = 10,
unit=['m', 'm', 'm', 'm', 'm**2/s**2'],
loc_names = ["Baro", "FG822-1", "FG822-2", "Smith", "ET"],
how="add", check_duplicates=True)
#%%
process = hgs.Processing(fowlers) #.decimate(2).by_dates(start='2015-11-01', stop='2016-02-01').by_gwloc("FG822-2")
#%%
process.info()
#%% test gw_correct
gw_correct_results = process.GW_correct(lag_h=24, et_method='hals')
#%%
# -*- coding: utf-8 -*-
"""
Created on Mon Jun 7 09:30:26 2021
@author: Daniel
"""
import hydrogeosines as hgs
import numpy as np
import pandas as pd
#%% Testing MVC principal
death_valley = hgs.Site('death valley', geoloc=[-116.471360, 36.408130, 688])
death_valley.import_csv('tests/data/death_valley/Rau_et_al_2021.csv',
input_category=["GW", "BP", "ET"],
utc_offset=0,
unit=["m", "m", "nstr"],
how="add",
check_duplicates=True)
#%%
# tmp = death_valley.data.hgs.get_loc_unit('ET')
#%% Processing
# create Instance of Processing with csiro_site
process = hgs.Processing(death_valley)
process.info()
#%%
# test hals method
hals_results = process.hals()
import hydrogeosines as hgs
#%% Testing MVC principal
## MODEL
acworth_site = hgs.Site('acworth', geoloc=[141.762065, -31.065781, 160])
print(acworth_site.data)
#%%
acworth_site.import_csv('tests/data/fowlers_gap/acworth_gw.csv', input_category=["GW","BP","GW"], utc_offset=10, unit=["Cm","mm","M"], how="add", check_duplicates=True)
acworth_site.import_csv('tests/data/fowlers_gap/acworth_bp.csv', input_category='BP', utc_offset=10, unit=["Hpa"], how="add", check_duplicates=False)
acworth_site.import_csv('tests/data/fowlers_gap/acworth_bp.csv', input_category='BP', utc_offset=10, unit="Cm", how="add", check_duplicates=True)
#acworth_site.import_csv('tests/data/fowlers_gap/acworth_et.csv', input_category='ET', utc_offset=10, unit='nm/s^2', how="add", check_duplicates=True)
#acworth_site.import_csv('test_data/fowlers_gap/acworth_short_gaps.csv', utc_offset=10, input_type=["BP", 'GW', 'GW', 'GW', 'ET'], unit=["m", 'm', 'm', 'm', 'nm/s^2'], method="add", check_duplicates=True) #, dt_fmt='%d/%m/%Y %H:%M'
#%%
## Model
acworth_site = hgs.Site('acworth', geoloc=[141.762065, -31.065781, 160])
# read
acworth_site.import_csv('tests/data/fowlers_gap/acworth_gw.csv',
input_category=["GW","BP","GW"], utc_offset=10, unit=["Cm","mm","M"],
header = ["Site_A","Site_B","Site_C"], how="add", check_duplicates=True)
# data
data = acworth_site.data
# hgs methods
data.hgs.resample(freq=5)
# datetime methods
data.hgs.dt.to_num
#%%
print(data.hgs.dt.to_utc)
import hydrogeosines as hgs
import numpy as np
import pandas as pd
from copy import deepcopy
#%% Csiro Data
## Model
csiro_site = hgs.Site('csiro', geoloc=[141.762065, -31.065781, 160])
# read data
csiro_site.import_csv('tests/data/csiro/test_sample/CSIRO_GW_short.csv',
input_category=["GW"] * 3,
utc_offset=10,
unit=["m"] * 3,
loc_names=["Loc_A", "Loc_B", "Loc_C"],
how="add",
check_duplicates=True)
csiro_site.import_csv('tests/data/csiro/test_sample/CSIRO_BP_short.csv',
input_category="BP",
utc_offset=10,
unit="mbar",
loc_names="Baro",
how="add",
check_duplicates=True)
#%%
csiro_site.add_ET(et_comp='g', et_cat=4)
#%% Processing
# create Instance of Processing with csiro_site
process_csiro = hgs.Processing(csiro_site)
Created on Fri Mar 19 12:22:07 2021
@author: Daniel
"""
#%% Load HGS and packages for testing
import numpy as np
import pandas as pd
import string
import collections.abc
import hydrogeosines as hgs
#%% Load Site Data
# Site Model
csiro_site = hgs.Site('csiro', geo=[141.762065, -31.065781, 160])
# Import
csiro_site.import_csv('tests/data/csiro/test_sample/CSIRO_GW_short.csv',
input_category=["GW"] * 3,
utc_offset=10,
unit=["m"] * 3,
loc_names=["Site_A", "Site_B", "Site_C"],
how="add",
check_duplicates=True)
csiro_site.import_csv('tests/data/csiro/test_sample/CSIRO_BP_short.csv',
input_category="BP",
utc_offset=10,
unit="mbar",
loc_names="Baro",
how="add",
# -*- coding: utf-8 -*-
"""
Created on Mon Jun 7 09:30:26 2021
@author: Daniel
"""
import hydrogeosines as hgs
import numpy as np
import pandas as pd
#%% Testing MVC principal
fowlers = hgs.Site('csiro', geoloc=[141.73099, -31.2934, 160])
fowlers.import_csv('tests/data/fowlers_gap/acworth_all.csv',
input_category=['BP', 'GW', 'GW', 'GW', 'ET'],
utc_offset=10,
unit=['m', 'm', 'm', 'm', 'm**2/s**2'],
loc_names=["Baro", "FG822-1", "FG822-2", "Smith", "ET"],
how="add",
check_duplicates=True)
#%% Processing
# create Instance of Processing with csiro_site
process = hgs.Processing(fowlers).by_gwloc(['FG822-2', 'Smith'])
# test hals method
fft_results = process.fft(update=True)
#%% Output
csiro_output = hgs.Output(fft_results)
# -*- coding: utf-8 -*-
"""
Created on Mon Jun 7 09:30:26 2021
@author: Daniel
"""
import hydrogeosines as hgs
import numpy as np
import pandas as pd
#%% Testing MVC principal
site = hgs.Site('Porto Alegre', geoloc=[-30.03779, -51.19535, 40])
site.import_csv('tests/data/brasil/Brito.csv',
input_category=["GW", "BP", "ET"],
utc_offset=-3,
unit=["m", "m", "nm/s**2"],
loc_names = ["GW", "Baro", "ET"],
how="add", check_duplicates=True)
#%% Processing
# create Instance of Processing with csiro_site
process = hgs.Processing(site)
process.info()
#%% test gw_correct
gw_correct_results = process.GW_correct(lag_h=8, et_method='hals')
correct_output = hgs.Output(gw_correct_results)
#%%
