Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_RCP45.nc =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_hist.nc =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_ERA5.p =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_hist.p =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/~$st day 1.docx =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py =================================================================== diff -u --- trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py (revision 0) +++ trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py (revision 103) @@ -0,0 +1,46 @@ +# -*- coding: utf-8 -*- +""" +Created on Mon Feb 3 14:12:50 2020 + +@author: alvareza +""" + +import os +import numpy as np +import pandas as pd +import xarray as xr +from netCDF4 import Dataset +from matplotlib import pyplot as plt + + +#import sys +#sys.path.insert(0,r'C:\Users\alvareza\checkouts\SD_Toolbox\trunk') +from SDToolBox import extract_data as sd_eda + +#%% Paths +class Paths: + def __init__(self,pth): + self.base = pth + self.data = os.path.join(self.base,'data') + self.out = os.path.join(self.base,'out') + self.plots = os.path.join(self.out,'plots') + + def _makedirs(self): + for path in vars(self).values(): + if not os.path.exists(path): + os.mkdir(path) + +pth = Paths(os.getcwd());pth._makedirs() + +#%% Get the coordinates +with Dataset(os.path.join(pth.data,'NS_WAM_ERA5.nc'),'r') as ds: + lat = ds['lat'][:] + lon = ds['lon'][:] +#era5 = pd.read_pickle(os.path.join(pth.data,'NS_WAM_ERA5.p')) +plt.plot(lat,lon,'or');plt.axis('equal') +#%% Extract data +# 1st initialize the class +inda = sd_eda.InputData('input_coordinates',{'lon':lon,'lat':lat})#, + input_variables=['swh'],input_scenarios=['era5'], + input_years:[1986]) +sdt.extract_data.InputData() \ No newline at end of file Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_RCP45.p =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_RCP85.nc =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_ERA5.nc =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/data/NS_WAM_RCP85.p =================================================================== diff -u Binary files differ Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/Test day 1.docx =================================================================== diff -u Binary files differ