import numpy as np
import pandas as pd
import gsw
import RandomFxns as RF
import BCClassFxns as BCF
import time
import cmocean.cm as cmo
import glob
import cartopy as ct
from cartopy.mpl.ticker import (LongitudeFormatter, LatitudeFormatter)
import matplotlib.gridspec as gridspec
import scipy.integrate as integrate
OFiles = glob.glob(
'/Users/Ellen/Documents/GitHub/BGCArgo_BCGyre/CSVFiles/OverlapFiles_BC/*')
WMODacDict = RF.WMODacPair()
#iso = [27.70, 27.75, 27.80,27.85]
iso = [27.70, 27.725, 27.75, 27.775]
# Interpolate to same space and time steps
# Calculate MLD at each point
good_QC = [1, 2, 5, 8]
lat_N = 70.000
lat_S = 48.00
lon_E = -45.00
lon_W = -70.00
figx = 10
@author: Ellen """ import xarray as xr import matplotlib.pyplot as plt import RandomFxns as RF import numpy as np from datetime import datetime import time import gsw import RandomFxns2 as RF2 dac = 'meds' dictfloat = RF.WMODacPair() WMO = 4902384 #4900879 dac = dictfloat[WMO] # 4900494 # f = RF.ArgoDataLoader(DAC=dac, WMO=WMO) # plt.figure() # f.plot.scatter(x='JULD',y='PRES_ADJUSTED',hue='DOXY_ADJUSTED') # plt.xticks(rotation=45) # plt.gca().invert_yaxis() # plt.show() ## # Random parameters you can change