I\'ve used xarray to read in a set of netCDF files that contain soil temperature data.
filin = xr.open_mfdataset(\'TSOI/*.001*.nc\')
Because it is soi