Example #1
0
# -*- coding: utf-8 -*-
"""
Created on Fri Dec  6 19:53:24 2013

@author: heijer
"""

from jarkus.transects import Transects
import numpy as np

Jk = Transects(url='http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/jarkus/profiles/transect_r2013_filled.nc')
ids = Jk.get_data('id')
idxs = np.floor(ids/1e6) == 8
for idx in np.nonzero(idxs)[0]:
    Jk.set_filter(alongshore=idx)
    with open('transect_%i.jrk'%ids[idx], 'w') as f:
        f.write(Jk.get_jrk())
with open("C:/Users/cijzendoornvan/Documents/GitHub/jarkus/jarkus/Settings.txt"
          ) as file:
    settings = json.load(file)

Dir_per_transect = settings['Dir_C3']
Dir_per_variable = settings['Dir_D1']
Dir_variable_plots = settings['Dir_D4']

years_requested = list(range(1965, 2020))
labels_y = [str(yr) for yr in years_requested][0::5]

##################################
####      INITIALISATION      ####
##################################
# Collect the JARKUS data from the server
Jk = Transects(url=settings['url'])
ids = Jk.get_data('id')  # ids of all available transects

remove1 = list(range(2000303, 2000304))
remove2 = list(range(2000720, 2000721))
remove3 = list(range(2002019, 2002020))
remove4 = list(range(4002001, 4005917, 1))
remove5 = list(range(5003300, 5004000, 1))
remove6 = list(range(6000416, 6000900, 1))
remove7 = list(range(6003070, 6003453, 1))
remove8 = list(range(8000000, 8005626, 1))
remove9 = list(range(9010193, 9010194, 1))
remove10 = list(range(10000000, 10001411, 1))
remove11 = list(range(11000660, 11000661, 1))
remove12 = list(range(11000680, 11000681, 1))
remove13 = list(range(11000700, 11000841, 1))
Example #3
0
 def setUp(self):
     from jarkus.transects import Transects
     self.tr = Transects()
Example #4
0
from pybeach.beach import Profile

get_ipython().run_line_magic('matplotlib', 'auto')

#%%
##################################
####       RETRIEVE DATA      ####
##################################
# load basic settings from .txt file
with open("C:/Users/cijzendoornvan/Documents/GitHub/jarkus/jarkus/Settings.txt"
          ) as file:
    settings = json.load(file)

# Collect the JARKUS data from the server
Jk = Transects(url=settings['url'])
ids = Jk.get_data('id')  # ids of all available transects

#%%
##################################
####    USER-DEFINED REQUEST  ####
##################################
# Select whether this is a new run, or whether a dimension is recalculated or added to the dataframe
new_run = False

# Select which dimensions should be calculated from the transects
dune_height_and_location = False
mean_sea_level = False
mean_low_water = False
mean_high_water = False
intertidal_width = False
Example #5
0
class UnitTests(unittest.TestCase):
    def setUp(self):
        from jarkus.transects import Transects
        self.tr = Transects()
    def test_initurl(self):
        from jarkus.transects import Transects
        url = 'http://opendap.deltares.nl/thredds/dodsC/opendap/rijkswaterstaat/jarkus/profiles/transect.nc'
        self.assertEqual(Transects(url=url).url, url)
    def test_set_filter_bool(self):
        shape = self.tr.ds.variables['id'].shape
        B = np.ones(shape) == 0 # all false
        idx = -1
        B[-1] = True
        self.tr.set_filter(alongshore=B)
        self.assertEqual(self.tr.get_data('id'), self.tr.ds.variables['id'][idx])
    def test_set_filter_id(self):
        id = 7e6
        self.tr.set_filter(id=id)
        self.assertEqual(self.tr.get_data('id'), id)
    def test_set_filter_idx(self):
        idx = -10
        self.tr.set_filter(alongshore=idx)
        self.assertEqual(self.tr.get_data('id'), self.tr.ds.variables['id'][idx])
    def test_set_filter_year(self):
        self.tr.set_filter(year=2006)
    def test_get_jrk(self):
        idx = 100
        self.tr.reset_filter()
        self.tr.set_filter(time=-1)
        ids = self.tr.get_data('id')
        idx = np.nonzero(ids==8006000)[0]
        self.tr.set_filter(alongshore=idx[0])
        self.tr.get_jrk()