Quick start

This tutorial shows how to open and plot your data.

import xarray_well_ice.xw_load as xwl
import xarray_well_ice.xw_tilt
import matplotlib.pyplot as plt

Piezometer

Load data from file

adr_P='data/Planpincieux_2_Press.dat'
ds_p=xwl.load_xw_pressure(adr_P)

ds_p

<xarray.DataArray (time: 4291)> Size: 34kB
array([5.761, 5.76 , 5.762, ..., 5.882, 5.88 , 5.88 ])
Coordinates:
  * time     (time) datetime64[ns] 34kB 2024-06-13T12:55:00 ... 2024-08-06T05...
Attributes:
    unit:     bar

xarray.DataArray

• time: 4291

• 5.761 5.76 5.762 5.764 5.763 5.762 ... 5.878 5.883 5.882 5.88 5.88

  array([5.761, 5.76 , 5.762, ..., 5.882, 5.88 , 5.88 ])

• Coordinates: (1)
  • time
    (time)
    datetime64[ns]
    2024-06-13T12:55:00 ... 2024-08-...

    array(['2024-06-13T12:55:00.000000000', '2024-06-13T12:56:00.000000000',
           '2024-06-13T12:57:00.000000000', ..., '2024-08-06T05:10:00.000000000',
           '2024-08-06T05:20:00.000000000', '2024-08-06T05:30:00.000000000'],
          dtype='datetime64[ns]')

• Indexes: (1)
  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2024-06-13 12:55:00', '2024-06-13 12:56:00',
                   '2024-06-13 12:57:00', '2024-06-13 12:58:00',
                   '2024-06-13 12:59:00', '2024-06-13 13:00:00',
                   '2024-06-13 13:01:00', '2024-06-13 13:02:00',
                   '2024-06-13 13:03:00', '2024-06-13 13:04:00',
                   ...
                   '2024-08-06 04:00:00', '2024-08-06 04:10:00',
                   '2024-08-06 04:20:00', '2024-08-06 04:30:00',
                   '2024-08-06 04:40:00', '2024-08-06 04:50:00',
                   '2024-08-06 05:00:00', '2024-08-06 05:10:00',
                   '2024-08-06 05:20:00', '2024-08-06 05:30:00'],
                  dtype='datetime64[ns]', name='time', length=4291, freq=None))

• Attributes: (1)

  unit :

  bar

Plot pressure evolution with time

ds_p.sel(time=slice("2024-07-13", "2024-08-06")).plot(label='presuure')
plt.grid()

Inclinometer

Load data from file

adr_i='data/Planpincieux_2_Tilt.dat'
ds_i=xwl.load_xw_tilt(adr_i)

ds_i.sel(time=slice("2024-07-13", "2024-08-06")).xwt.plot_raw()
plt.grid()
plt.legend()

<matplotlib.legend.Legend at 0x74b67adf74d0>