Accessing Argo data by float using Argopy

Import the usual suspects:

import xarray as xr
xr.set_options(display_style="html", display_expand_attrs=False);
from matplotlib import pyplot as plt
plt.style.use('default')

Import argopy and set-up a data fetcher:

from argopy import DataFetcher as ArgoDataFetcher
argo_loader = ArgoDataFetcher(src='erddap', cache=True, cachedir='tmp')

/Users/gmaze/miniconda3/envs/aos/lib/python3.10/site-packages/pyproj/__init__.py:89: UserWarning: pyproj unable to set database path.
  _pyproj_global_context_initialize()

You can load all profiles from a float:

apDS = argo_loader.float(6901254).load().data

or a single cycle:

argo_loader = argo_loader.profile(6901254, 1).load()
apDS2 = argo_loader.data

data = apDS2.argo.point2profile()
data

<xarray.Dataset>
Dimensions:                (N_PROF: 2, N_LEVELS: 97)
Coordinates:
  * N_PROF                 (N_PROF) int64 0 1
  * N_LEVELS               (N_LEVELS) int64 0 1 2 3 4 5 6 ... 91 92 93 94 95 96
    LATITUDE               (N_PROF) float64 29.16 29.18
    LONGITUDE              (N_PROF) float64 -15.49 -15.43
    TIME                   (N_PROF) datetime64[ns] 2018-10-23T20:54:00 2018-1...
Data variables: (12/13)
    CONFIG_MISSION_NUMBER  (N_PROF) int64 1 1
    CYCLE_NUMBER           (N_PROF) int64 1 1
    DATA_MODE              (N_PROF) <U1 'R' 'R'
    DIRECTION              (N_PROF) <U1 'D' 'A'
    PLATFORM_NUMBER        (N_PROF) int64 6901254 6901254
    POSITION_QC            (N_PROF) int64 1 1
    ...                     ...
    PRES_QC                (N_PROF) int64 1 1
    PSAL                   (N_PROF, N_LEVELS) float64 36.69 36.71 ... 35.14
    PSAL_QC                (N_PROF) int64 1 1
    TEMP                   (N_PROF, N_LEVELS) float64 23.4 23.35 ... 4.559 4.474
    TEMP_QC                (N_PROF) int64 1 1
    TIME_QC                (N_PROF) int64 1 1
Attributes: (8)

xarray.Dataset

• Dimensions:
  • N_PROF: 2
  • N_LEVELS: 97
• Coordinates: (5)
  • N_PROF
    (N_PROF)
    int64
    0 1

    casted :

    1

    array([0, 1])

  • N_LEVELS
    (N_LEVELS)
    int64
    0 1 2 3 4 5 6 ... 91 92 93 94 95 96

    casted :

    1

    array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17,
           18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
           36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
           54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,
           72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89,
           90, 91, 92, 93, 94, 95, 96])

  • LATITUDE
    (N_PROF)
    float64
    29.16 29.18

    _CoordinateAxisType :

    Lat

    actual_range :

    [29.162 29.182]

    axis :

    Y

    colorBarMaximum :

    90.0

    colorBarMinimum :

    -90.0

    ioos_category :

    Location

    long_name :

    Latitude of the station, best estimate

    standard_name :

    latitude

    units :

    degrees_north

    valid_max :

    90.0

    valid_min :

    -90.0

    casted :

    1

    array([29.162, 29.182])

  • LONGITUDE
    (N_PROF)
    float64
    -15.49 -15.43

    _CoordinateAxisType :

    Lon

    actual_range :

    [-15.488 -15.429]

    axis :

    X

    colorBarMaximum :

    180.0

    colorBarMinimum :

    -180.0

    ioos_category :

    Location

    long_name :

    Longitude of the station, best estimate

    standard_name :

    longitude

    units :

    degrees_east

    valid_max :

    180.0

    valid_min :

    -180.0

    casted :

    1

    array([-15.488, -15.429])

  • TIME
    (N_PROF)
    datetime64[ns]
    2018-10-23T20:54:00 2018-10-25T0...

    _CoordinateAxisType :

    Time

    actual_range :

    [1.54032804e+09 1.54044540e+09]

    axis :

    T

    ioos_category :

    Time

    long_name :

    Julian day (UTC) of the station relative to REFERENCE_DATE_TIME

    standard_name :

    time

    time_origin :

    01-JAN-1970 00:00:00

    casted :

    1

    array(['2018-10-23T20:54:00.000000000', '2018-10-25T05:30:00.000000000'],
          dtype='datetime64[ns]')

• Data variables: (13)
  • CONFIG_MISSION_NUMBER
    (N_PROF)
    int64
    1 1

    actual_range :

    [1 1]

    colorBarMaximum :

    100.0

    colorBarMinimum :

    0.0

    conventions :

    1...N, 1 : first complete mission

    ioos_category :

    Statistics

    long_name :

    Unique number denoting the missions performed by the float

    casted :

    1

    array([1, 1])

  • CYCLE_NUMBER
    (N_PROF)
    int64
    1 1

    long_name :

    Float cycle number

    convention :

    0..N, 0 : launch cycle (if exists), 1 : first complete cycle

    casted :

    1

    array([1, 1])

  • DATA_MODE
    (N_PROF)
    <U1
    'R' 'R'

    long_name :

    Delayed mode or real time data

    convention :

    R : real time; D : delayed mode; A : real time with adjustment

    casted :

    1

    array(['R', 'R'], dtype='<U1')

  • DIRECTION
    (N_PROF)
    <U1
    'D' 'A'

    long_name :

    Direction of the station profiles

    convention :

    A: ascending profiles, D: descending profiles

    casted :

    1

    array(['D', 'A'], dtype='<U1')

  • PLATFORM_NUMBER
    (N_PROF)
    int64
    6901254 6901254

    long_name :

    Float unique identifier

    convention :

    WMO float identifier : A9IIIII

    casted :

    1

    array([6901254, 6901254])

  • POSITION_QC
    (N_PROF)
    int64
    1 1

    long_name :

    Global quality flag of POSITION_QC profile

    convention :

    Argo reference table 2a

    casted :

    1

    array([1, 1])

  • PRES
    (N_PROF, N_LEVELS)
    float64
    14.0 24.0 ... 1.988e+03 2.011e+03

    long_name :

    Sea Pressure

    standard_name :

    sea_water_pressure

    units :

    decibar

    valid_min :

    0.0

    valid_max :

    12000.0

    resolution :

    0.1

    axis :

    Z

    casted :

    1

    array([[  14.,   24.,   34.,   45.,   54.,   65.,   75.,   84.,   94.,
             104.,  114.,  125.,  135.,  145.,  154.,  164.,  175.,  185.,
             195.,  212.,  237.,  262.,  287.,  312.,  337.,  362.,  387.,
             412.,  437.,  462.,  486.,  512.,  537.,  563.,  587.,  612.,
             637.,  662.,  687.,  712.,  737.,  762.,  787.,  812.,  837.,
             862.,  887.,  912.,  937.,  962.,  978.,   nan,   nan,   nan,
              nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,
              nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,
              nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,
              nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,   nan,
              nan,   nan,   nan,   nan,   nan,   nan,   nan],
           [   6.,    7.,    8.,    9.,   10.,   15.,   25.,   36.,   45.,
              55.,   66.,   76.,   86.,   96.,  105.,  115.,  125.,  135.,
             146.,  155.,  165.,  175.,  185.,  195.,  213.,  238.,  263.,
             288.,  313.,  338.,  363.,  388.,  414.,  438.,  463.,  488.,
             513.,  538.,  563.,  588.,  613.,  638.,  663.,  688.,  713.,
             738.,  763.,  788.,  813.,  838.,  863.,  888.,  913.,  938.,
             963.,  988., 1013., 1038., 1063., 1088., 1113., 1138., 1163.,
            1188., 1213., 1238., 1263., 1288., 1313., 1338., 1363., 1388.,
            1414., 1438., 1463., 1488., 1513., 1538., 1564., 1588., 1613.,
            1638., 1663., 1688., 1713., 1738., 1763., 1788., 1813., 1838.,
            1863., 1888., 1913., 1938., 1963., 1988., 2011.]])

  • PRES_QC
    (N_PROF)
    int64
    1 1

    long_name :

    Global quality flag of PRES_QC profile

    convention :

    Argo reference table 2a

    casted :

    1

    array([1, 1])

  • PSAL
    (N_PROF, N_LEVELS)
    float64
    36.69 36.71 36.71 ... 35.15 35.14

    long_name :

    PRACTICAL SALINITY

    standard_name :

    sea_water_salinity

    units :

    psu

    valid_min :

    0.0

    valid_max :

    43.0

    resolution :

    0.001

    casted :

    1

    array([[36.69400024, 36.70500183, 36.70899963, 36.67499924, 36.52199936,
            36.51699829, 36.52000046, 36.51699829, 36.51200104, 36.50099945,
            36.50400162, 36.48199844, 36.42900085, 36.41400146, 36.38800049,
            36.36000061, 36.3429985 , 36.31999969, 36.28099823, 36.19800186,
            36.07400131, 35.97100067, 35.89899826, 35.85900116, 35.82500076,
            35.78499985, 35.73899841, 35.7179985 , 35.68999863, 35.67699814,
            35.64400101, 35.60599899, 35.58300018, 35.56399918, 35.53499985,
            35.51499939, 35.48699951, 35.45500183, 35.43899918, 35.41699982,
            35.39599991, 35.38100052, 35.37400055, 35.36999893, 35.35699844,
            35.34600067, 35.33599854, 35.34199905, 35.35599899, 35.37400055,
            35.39199829,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan],
           [36.6969986 , 36.6969986 , 36.6969986 , 36.69800186, 36.69800186,
            36.69900131, 36.70600128, 36.70199966, 36.55799866, 36.47100067,
            36.43700027, 36.46300125, 36.49300003, 36.45999908, 36.42200089,
            36.41400146, 36.39799881, 36.3429985 , 36.3030014 , 36.27199936,
            36.23500061, 36.20600128, 36.17200089, 36.12799835, 36.06700134,
            35.98600006, 35.92399979, 35.85900116, 35.81100082, 35.77099991,
            35.7480011 , 35.72200012, 35.68799973, 35.66199875, 35.63600159,
            35.59999847, 35.57699966, 35.55599976, 35.5320015 , 35.50999832,
            35.48600006, 35.46500015, 35.44100189, 35.42900085, 35.4129982 ,
            35.40200043, 35.39099884, 35.38600159, 35.37799835, 35.36100006,
            35.34799957, 35.34600067, 35.35499954, 35.3769989 , 35.39799881,
            35.41400146, 35.43000031, 35.43600082, 35.44200134, 35.45399857,
            35.47000122, 35.49499893, 35.52999878, 35.54399872, 35.53300095,
            35.50099945, 35.47600174, 35.45500183, 35.43500137, 35.41799927,
            35.4109993 , 35.40800095, 35.38700104, 35.38600159, 35.37900162,
            35.375     , 35.36299896, 35.35499954, 35.3409996 , 35.32899857,
            35.32400131, 35.30400085, 35.29000092, 35.28300095, 35.27299881,
            35.26200104, 35.2519989 , 35.2480011 , 35.23500061, 35.22499847,
            35.21500015, 35.20600128, 35.1969986 , 35.18600082, 35.16899872,
            35.15299988, 35.14199829]])

  • PSAL_QC
    (N_PROF)
    int64
    1 1

    long_name :

    Global quality flag of PSAL_QC profile

    convention :

    Argo reference table 2a

    casted :

    1

    array([1, 1])

  • TEMP
    (N_PROF, N_LEVELS)
    float64
    23.4 23.35 23.33 ... 4.559 4.474

    long_name :

    SEA TEMPERATURE IN SITU ITS-90 SCALE

    standard_name :

    sea_water_temperature

    units :

    degree_Celsius

    valid_min :

    -2.0

    valid_max :

    40.0

    resolution :

    0.001

    casted :

    1

    array([[23.40299988, 23.34799957, 23.33499908, 23.00300026, 20.88699913,
            19.76300049, 18.9279995 , 18.45599937, 18.12400055, 17.85099983,
            17.68899918, 17.48800087, 17.12100029, 16.99699974, 16.80699921,
            16.62299919, 16.48999977, 16.34399986, 16.11400032, 15.66100025,
            14.94799995, 14.34200001, 13.85599995, 13.57999992, 13.31599998,
            13.02000046, 12.63799953, 12.47700024, 12.22999954, 12.11299992,
            11.80900002, 11.44200039, 11.22099972, 11.04300022, 10.7670002 ,
            10.56900024, 10.30300045, 10.00399971,  9.82800007,  9.57499981,
             9.30200005,  9.12800026,  9.00500011,  8.91800022,  8.74100018,
             8.58300018,  8.43700027,  8.36900043,  8.36699963,  8.32999992,
             8.34799957,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan,         nan,         nan,         nan,
                    nan,         nan],
           [23.41200066, 23.41200066, 23.41399956, 23.41200066, 23.41300011,
            23.40999985, 23.39800072, 23.36199951, 21.52000046, 19.5510006 ,
            18.40800095, 17.9470005 , 17.68899918, 17.37800026, 17.0739994 ,
            17.        , 16.86899948, 16.57600021, 16.33399963, 16.11300087,
            15.88199997, 15.70199966, 15.48600006, 15.23400021, 14.8920002 ,
            14.42199993, 14.02400017, 13.58100033, 13.24400043, 12.93299961,
            12.75      , 12.52600002, 12.22200012, 11.99400043, 11.72599983,
            11.39599991, 11.17599964, 10.98200035, 10.74100018, 10.52799988,
            10.29300022, 10.09799957,  9.86800003,  9.73099995,  9.55300045,
             9.39400005,  9.21399975,  9.09899998,  8.96500015,  8.79500008,
             8.62699986,  8.50199986,  8.38799953,  8.3739996 ,  8.3739996 ,
             8.36600018,  8.32900047,  8.24100018,  8.12800026,  8.08899975,
             8.10599995,  8.0909996 ,  8.17099953,  8.17500019,  8.02900028,
             7.79400015,  7.56799984,  7.40999985,  7.24399996,  7.08099985,
             6.95800018,  6.88199997,  6.69000006,  6.63999987,  6.53800011,
             6.45200014,  6.33699989,  6.25600004,  6.13100004,  6.01100016,
             5.95300007,  5.78999996,  5.67700005,  5.60900021,  5.52899981,
             5.44099998,  5.36299992,  5.32999992,  5.21799994,  5.1329999 ,
             5.04699993,  4.97200012,  4.89300013,  4.8039999 ,  4.67999983,
             4.55900002,  4.47399998]])

  • TEMP_QC
    (N_PROF)
    int64
    1 1

    long_name :

    Global quality flag of TEMP_QC profile

    convention :

    Argo reference table 2a

    casted :

    1

    array([1, 1])

  • TIME_QC
    (N_PROF)
    int64
    1 1

    long_name :

    Global quality flag of TIME_QC profile

    convention :

    Argo reference table 2a

    casted :

    1

    array([1, 1])

• Attributes: (8)

  DATA_ID :

  ARGO

  DOI :

  http://doi.org/10.17882/42182

  Fetched_from :

  https://www.ifremer.fr/erddap

  Fetched_by :

  gmaze

  Fetched_date :

  2022/04/11

  Fetched_constraints :

  phy;WMO6901254_CYC1

  Fetched_uri :

  ['https://www.ifremer.fr/erddap/tabledap/ArgoFloats.nc?data_mode,latitude,longitude,position_qc,time,time_qc,direction,platform_number,cycle_number,config_mission_number,vertical_sampling_scheme,pres,temp,psal,pres_qc,temp_qc,psal_qc,pres_adjusted,temp_adjusted,psal_adjusted,pres_adjusted_qc,temp_adjusted_qc,psal_adjusted_qc,pres_adjusted_error,temp_adjusted_error,psal_adjusted_error&platform_number=~"6901254"&cycle_number=~"1"&distinct()&orderBy("time,pres")']

  history :

  Variables filtered according to DATA_MODE; Variables selected according to QC; Transformed with point2profile

The core-Argo profile files contain the core parameters provided by a float: pressure, temperature, salinity, conductivity (PRES, TEMP, PSAL, CNDC).

fig, ax = plt.subplots(1, 2, figsize=(8,10), sharey=True)

#Temperature
ax[0].plot(data.TEMP[0],-data.PRES[0],'ro',label='N_PROF=0 D Fetcher')
ax[0].plot(data.TEMP[1],-data.PRES[1],'bo',label='N_PROF=1 A Fetcher')
ax[0].set_title(data['TEMP'].attrs['long_name'])
ax[0].set_ylabel(data['PRES'].attrs['long_name'])
ax[0].grid()
ax[0].legend();

ax[1].plot(data.PSAL[0],-data.PRES[0],'ro',label='N_PROF=0 D Fetcher')
ax[1].plot(data.PSAL[1],-data.PRES[1],'bo',label='N_PROF=1 A Fetcher')
ax[1].set_title(data['TEMP'].attrs['long_name'])
ax[1].grid()
ax[1].legend();

fig, ax = plt.subplots(figsize=(8,8))

ax.plot(data.LONGITUDE[0],data.LATITUDE[0],'ko',label='001D N_PROF=0 Descending')
ax.plot(data.LONGITUDE[1],data.LATITUDE[1],'bo',label='001 N_PROF=0 Ascending')
ax.grid()
ax.legend();

argo_loader.plot('trajectory', set_global=True)

(<Figure size 900x540 with 1 Axes>,
 <GeoAxesSubplot:xlabel='longitude', ylabel='latitude'>)