url_thredds = 'http://gliders.ioos.us/thredds/dodsC/deployments/aoml/SG668-20190819T1217/SG668-20190819T1217.nc3.nc'
var_name = 'temperature'
#var = 'salinity'
date_ini = '2019/09/01/00' # year/month/day/hour
date_end = '2019/09/10/00' # year/month/day/hour
scatter_plot = 'no'
kwargs = dict(date_ini=date_ini, date_end=date_end)
tempg, timeg, latg, long, depthg, dataset_id = \
read_glider_data_thredds_server(url_thredds,var_name,scatter_plot,**kwargs)
delta_z = 0.4 # bin size in the vertical when gridding the variable vertical profile
# default value is 0.3
contour_plot = 'yes' # default value is 'yes'
tempg_gridded, timegg, depthg_gridded = \
grid_glider_data(var_name,dataset_id,tempg,timeg,latg,long,depthg,delta_z,contour_plot)
#%% Cell #4: Search for glider data sets given
# a latitude and longitude box and time window
from read_glider_data import retrieve_dataset_id_erddap_server
# Server location
url_erddap = 'https://data.ioos.us/gliders/erddap'
'''
# MAB
lon_lim = [-75.0,-72.0]
lat_lim = [38.0,40.0]
# date limits
date_ini = '2018/09/01/00'
if not empty_dataset:
print(dataset_id)
kwargs = dict(date_ini=date_ini,date_end=date_end)
scatter_plot = 'no'
tempg, saltg, timeg, latg, long, depthg = read_glider_data_erddap_server(url_erddap,dataset_id,\
lat_lim,lon_lim,scatter_plot,**kwargs)
# Get glider density
densg = sw.dens(saltg,tempg,depthg)
# Grid glider variables according to depth
delta_z = 0.5
tempg_gridded, timegg, depthg_gridded = \
grid_glider_data('temp',dataset_id,tempg,timeg,depthg,delta_z,contour_plot='no')
saltg_gridded, _, _ = \
grid_glider_data('temp',dataset_id,saltg,timeg,depthg,delta_z,contour_plot='no')
densg_gridded, _, _ = \
grid_glider_data('temp',dataset_id,densg,timeg,depthg,delta_z,contour_plot='no')
# Changing times to timestamp
tstamp_glider = mdates.date2num(timeg)
#%% Get glider transect from GOFS 3.1
# Conversion from glider longitude and latitude to GOFS convention
target_lon, target_lat = glider_coor_to_GOFS_coord(long,latg)
# interpolating glider lon and lat to lat and lon on model time
print(gliders)
dataset_id = gliders[6]
kwargs = dict(date_ini=date_ini, date_end=date_end)
scatter_plot = 'no'
tempg, saltg, timeg, latg, long, depthg = read_glider_data_erddap_server(url_erddap,dataset_id,\
lat_lim,lon_lim,scatter_plot,**kwargs)
#%%
contour_plot = 'no' # default value is 'yes'
delta_z = 0.4 # default value is 0.3
tempg_gridded, timegg, depthg_gridded = \
grid_glider_data('temperature',dataset_id,tempg,timeg,depthg,delta_z,contour_plot)
#saltg_gridded, timegg, depthg_gridded = \
# grid_glider_data('salinity',dataset_id,saltg,timeg,depthg,delta_z,contour_plot)
#%% RU33
# variable to retrieve
var_name = 'temperature'
varg = tempg
color_map = cmocean.cm.thermal
clabel = '($^oC$)'
timeg_matrix = np.tile(timeg.T, (depthg.shape[0], 1))
ttg = np.ravel(timeg_matrix)
dg = np.ravel(depthg)
teg = np.ravel(varg)
import numpy as np
ti = mdates.date2num(datetime(2020, 10, 15))
te = mdates.date2num(datetime(2020, 10, 27))
okt = np.where(np.logical_and(mdates.date2num(timeg) >= ti,\
mdates.date2num(timeg) <= te))[0]
tempgg = tempg[:, okt]
timegg = timeg[okt]
latgg = latg[okt]
longg = long[okt]
depthgg = depthg[:, okt]
tempg_gridded, timeg_gridded, depthg_gridded = \
grid_glider_data(var_name_glider,dataset_id,tempgg,timegg,latgg,longg,depthgg,delta_z,contour_plot='yes')
#%% Plot trayectory
import matplotlib.pyplot as plt
import cartopy
import cartopy.feature as cfeature
#import matplotlib.ticker as mticker
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
fig, ax = plt.subplots(figsize=(7, 3),
subplot_kw=dict(projection=cartopy.crs.PlateCarree()))
coast = cfeature.NaturalEarthFeature('physical', 'coastline', '10m')
ax.add_feature(coast, edgecolor='black', facecolor='none')
ax.add_feature(cfeature.BORDERS) # adds country borders
ax.add_feature(cfeature.STATES) # adds statet borders
#%% Sam
dataset_id = [id for id in gliders if id.split('-')[0] == 'sam'][0]
kwargs = dict(date_ini=date_ini, date_end=date_end)
scatter_plot = 'no'
tempg, saltg, timeg, latg, long, depthg = read_glider_data_erddap_server(url_erddap,dataset_id,\
lat_lim,lon_lim,scatter_plot,**kwargs)
#%%
contour_plot = 'no' # default value is 'yes'
delta_z = 0.4 # default value is 0.3
tempg_gridded, timegg, depthg_gridded = \
grid_glider_data('temperature',dataset_id,tempg,timeg,depthg,delta_z,contour_plot)
saltg_gridded, timegg, depthg_gridded = \
grid_glider_data('salinity',dataset_id,saltg,timeg,depthg,delta_z,contour_plot)
#%%
# variable to retrieve
var_name = 'temperature'
varg = tempg
if var_name == 'temperature':
color_map = cmocean.cm.thermal
clabel = '($^oC$)'
else:
if var_name == 'salinity':
color_map = cmocean.cm.haline
