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