Download notebook.

Download Swot from AVISO

This notebook explains how to explore and retrieve available SWOT LR data from AVISO’s THREDDS Data Server and FTP Server, presented here.

L2 data can be explored at: https://tds.aviso.altimetry.fr/thredds/L2/L2-SWOT-DATA/L2-SWOT.html

L3 data can be explored at: https://tds.aviso.altimetry.fr/thredds/catalog/dataset-l3-swot-karin-nadir-validated/l3_lr_ssh/catalog.html

Tutorial Objectives

Discover Aviso’s Thredds Data Server catalogue
Download files through HTTP and FTP
Download data subset through Opendap

Import + code

[1]:

import os
import re
from getpass import getpass
import requests as rq
import ftplib
import xarray as xr
from xarray.backends import PydapDataStore

[2]:

def _get_filename_from_cd(cd):
    """ Get filename from content-disposition """
    if not cd:
        return None
    fname = re.findall('filename="(.+)"', cd)
    if len(fname) == 0:
        return None
    return fname[0]

def http_download_file(dataset_url:str, output_path:str):
    response = rq.get(dataset_url, auth=(username, password))

    filename = _get_filename_from_cd(response.headers['content-disposition'])
    file_path = os.path.join(output_path, filename)

    open(file_path, 'wb').write(response.content)
    return file_path

def _download_file(ftp:str, filename:str, target_directory:str):
    try:
        local_filepath = os.path.join(target_directory, filename)
        with open(local_filepath, 'wb') as file:
            ftp.retrbinary('RETR %s' % filename, file.write)
            print(f"Downloaded {filename} to {target_directory}")
            return local_filepath
    except Exception as e:
        print(f"Error downloading {filename}: {e}")

def ftp_download_file(ftp_path:str, filename:str, output_path:str):
    # Set up FTP server details
    ftpAVISO = 'ftp-access.aviso.altimetry.fr'
    try:
        # Logging into FTP server using provided credentials
        with ftplib.FTP(ftpAVISO) as ftp:
            ftp.login(username, password)
            ftp.cwd(ftp_path)
            print(f"Connection Established {ftp.getwelcome()}")

            # Check if the file exists in the directory
            if filename in ftp.nlst():
                return _download_file(ftp, filename, output_path)
            else:
                print(f"File {filename} does not exist in the directory {ftp}.")
    except ftplib.error_perm as e:
        print(f"FTP error: {e}")
    except Exception as e:
        print(f"Error: {e}")


def open_opendap_dataset(dataset_url):
    session = rq.Session()
    session.auth = (username, password)
    store = PydapDataStore.open(dataset_url, session=session, timeout=300, user_charset='UTF-8')
    return xr.open_dataset(store)

Parameters

Define output folder to save results

[3]:

output_dir = "downloads"

Authentication parameters

Enter your AVISO+ credentials

[ ]:

password = getpass(f"Enter password for {username}:")

[ ]:

username = input("Enter username:")

Download files through HTTP

To find the data you want to download, you can explore the Swot LR L2 or Swot LR L3 catalogues.

Copy the FileServer url to download data through HTTP :

Swot LR L2 data

[5]:

l2_dataset_url = "https://tds.aviso.altimetry.fr/thredds/fileServer/dataset-l2-swot-karin-lr-ssh-pre-validated/PIC0/Basic/cycle_014/SWOT_L2_LR_SSH_Basic_014_001_20240417T111950_20240417T121118_PIC0_01.nc"

[6]:

l2_file = http_download_file(l2_dataset_url, output_dir)
l2_file

[6]:

'downloads/SWOT_L2_LR_SSH_Basic_014_001_20240417T111950_20240417T121118_PIC0_01.nc'

Swot LR L3 data

[7]:

l3_dataset_url = "https://tds.aviso.altimetry.fr/thredds/fileServer/dataset-l3-swot-karin-nadir-pre-validated/l3_lr_ssh/v1_0/Basic/cycle_014/SWOT_L3_LR_SSH_Basic_014_001_20240417T111950_20240417T121116_v1.0.nc"

[8]:

l3_file = http_download_file(l3_dataset_url, output_dir)
l3_file

[8]:

'downloads/SWOT_L3_LR_SSH_Basic_014_001_20240417T111950_20240417T121116_v1.0.nc'

Download files through FTP

Files available through FTP follow the same file tree as in the THREDDS Data Server.

Data	FTP starting point	Browsing files URL
Swot LR L2	/swot_products/l2_karin/l2_lr_ssh	https://tds.aviso.altimetry.fr/thredds/catalog/dataset-l2-swot-karin-lr-ssh-pre-validated/catalog.html
Swot LR L3	/swot_products/l3_karin_nadir/l3_lr_ssh	https://tds.aviso.altimetry.fr/thredds/catalog/dataset-l3-swot-karin-nadir-validated/l3_lr_ssh/catalog.html

Swot LR L2 data

[9]:

l2_ftp_path = '/swot_products/l2_karin/l2_lr_ssh/PGC0/WindWave/cycle_001/'
l2_ftp_file = 'SWOT_L2_LR_SSH_WindWave_001_155_20230726T173636_20230726T182805_PGC0_01.nc'

[10]:

l2_file = ftp_download_file(l2_ftp_path, l2_ftp_file, output_dir)
l2_file

Connection Established 220 192.168.10.119 FTP server ready
Downloaded SWOT_L2_LR_SSH_WindWave_001_155_20230726T173636_20230726T182805_PGC0_01.nc to downloads

[10]:

'downloads/SWOT_L2_LR_SSH_WindWave_001_155_20230726T173636_20230726T182805_PGC0_01.nc'

Swot LR L3 data

[11]:

l3_ftp_path = '/swot_products/l3_karin_nadir/l3_lr_ssh/v1_0/Basic/cycle_011/'
l3_ftp_file = 'SWOT_L3_LR_SSH_Basic_011_550_20240305T114902_20240305T124029_v1.0.nc'

[12]:

l3_file = ftp_download_file(l3_ftp_path, l3_ftp_file, output_dir)
l3_file

Connection Established 220 192.168.10.119 FTP server ready
Downloaded SWOT_L3_LR_SSH_Basic_011_550_20240305T114902_20240305T124029_v1.0.nc to downloads

[12]:

'downloads/SWOT_L3_LR_SSH_Basic_011_550_20240305T114902_20240305T124029_v1.0.nc'

Load subset via OPeNDAP

To find the data you want to subset, you can explore the Swot LR L2 or Swot LR L3 catalogues.

Copy the OPeNDAP url to subset data through Opendap :

Store subset to NetCDF file

[ ]:

ds.to_netcdf(os.path.join(output_dir, 'dataset.nc'))

Download notebook.

Download Swot from AVISO

Tutorial Objectives

Import + code

Parameters

Authentication parameters

Download files through HTTP

Swot LR L2 data

Swot LR L3 data

Download files through FTP

Swot LR L2 data

Swot LR L3 data

Load subset via OPeNDAP

Swot LR L2 data

Select variables

Select lines

Swot LR L3 data

Select variables

Select lines

Store subset to NetCDF file