ESMGFZ Product Repository - extractlatlon_bilinintp_remote Login, user settings, help
URL:Download
script to extract station timeseries

The script downloads the time series for the 4 gridcells surrounding the station location and performs bilinear interpolation. You have to download only this script and run it locally on your machine. In order to use the option -m you also have to set the NCO library path in the script.

In addition to single products (NTAL, NTOL, HYDL, SLEL), you can also download sums of loading products:
NTAL+NTOL
NTAL+HYDL
NTOL+SLEL
NTAL+NTOL+HYDL
NTAL+NTOL+HYDL+SLEL
The individual time series will be downloaded, resampled, and summed up before the bilinear interpolation.

ATTENTION:
You can use the -m option to avoid interpolation over the coastline as the deformation response might be completely different between ocean and land due to a high contrast in the mass load along the coastline (e.g. IB correction over the ocean). The appropriate land/ocean mask file is lwmask_AOD_0.5.nc. With option -m the script uses nearest-neighbour interpolation, if the station is located along the coastline to avoid bicubic interpolation over the coastline between ocean and land gridcells. For time series of total loading (NTAL+NTOL+HYDL) it is recommended to add the gridded NTAL, NTOL, and HYDL contributions in advance and extract the station time series afterwards, the -m option might become redundant.

NOTE:
Requests of loading time series can take a bit of time to fulfill, especially for larger aggregations like NTAL time series over 40 years, so please be patient. You can approximately calculate with:

  • 40-60 sec for 3-hourly data (e.g. NTAL, NTOL, NTAL24h, NTAL+HYDL, ...)
  • 3 sec for 24-hourly (e.g. HYDL, SLEL)
per station, per year.

USAGE:
to extract single station time series
> bash extractlatlon_bilinintp_remote dataset frame startdate enddate latitude longitude [-m] [-v variable1 [variable2 ...]] [-o outfilename]

to extract multi station time series, saved as single files for each station
>bash extractlatlon_bilinintp_remote dataset frame startdate enddate -s stationlistfile [-c columnorder] [-m] [-v variable1 [variable2 ...]] [-o outfilename]

arguments: dataset one of NTAL, NTOL, HYDL or
NTAL24h, NTOL24h or
NTAL+NTOL, NTOL+SLEL, NTAL+HYDL, NTAL+NTOL+HYDL, NTAL+NTOL+HYDL+SLEL
frame CF or CM
startdate 1st epoch (format: yyyy-mm-dd)
enddate last epoch (format: yyyy-mm-dd)
latitude latitude of station [degree]
longitude longitude of station [degree]

options: -s stationlistfile name of stationlistfile containing station coordinates (no header, 2 or 3 columns: [stationname/id] longitude latitude, UNIX format)
-c columnorder columnorder "latlon" or "lonlat", columnorder in stationlistname file (default=latlon)
"xyz" for cartesian coordinates, coverted by WGS84 parameters
-m uses land/ocean maskfile to prevent interpolation over the coastline (0=land/1=ocean)
-v list of vars list of vars specify variables from netCDF files (default: duV, duN, duEW)
-o outfilename specify file names for output (->outfilename.stationID.txt)
default is >extractlatlon<.

EXAMPLES:
> extractlatlon_bilinintp_remote HYDL CF 2000-01-01 2015-01-01 -3.023 -60.055

This will extract the variables duV, duEW, duNS from hydrological loading (CF frame) at the location Manaus City (3.023°S, 60.055°W) and writes the time series into the default output file extractlatlon.duV.

> extractlatlon_bilinintp_remote HYDL CF 2000-01-01 2000-12-31T21:00:00Z -3.023 -60.055 -v duV -o manaus

This will extract only the variable duV for the year 2000 from HYDL (CF) at the location Manaus City (3.1333°S, 60.0167°W) and writes the time series to then file manaus.txt

> extractlatlon_bilinintp_remote NTAL+HYDL CM 2000-01-01 2015-01-01 -s stationlist.txt -m -v duV duEW duNS -o esmgfz_ntalhydl

This will extract the variables duV, duEW, duNS from the sum of NTAL and HYDL (CM) at the locations listed in the file stationlist.txt and writes the time series into single files named esmgfz_ntalhydl.stationID.txt. The file stationlist.txt consist of 2 or 3 columns with station name or ID (optional), latitude, longitude. If latitude and longitude are switched, use option -c lonlat. The -m option avoids interpolation over the coastline defined in land-water maskfile lwmask_AOD_0.5.nc (mask that was used for IB correction).

> extractlatlon_bilinintp_remote NTAL CF 2000-01-01 2000-12-31T21:00:00Z -s stationlist.txt -c xyz -m -v duV duEW duNS -o ESMGFZ_NTAL

this will extract the variables duV, duEW, duNS from atmospheric loading (NTAL, CF frame) at the locations listed in the file stationlist.txt and writes the time series into single files named ESMGFZ_NTAL.stationID.txt. The file stationlist.txt consist of 3 or 4 columns with station name or ID (optional), x-coordinate, y-coordinate, and z-coordinate (-c xyz). The -m option avoids interpolation over the coastline defined in land-water maskfile lwmask_AOD_0.5.nc (mask that was used for IB correction).
File: extractlatlon_bilinintp_remote  Download  30 KB
Kind: Software Tool
Software Use:Analysis
Software Type:Desktop
Science Domain:Earth
Platform:Unix
License:Open Source
Development Status:beta
Capabilities:extract station time series
from netCDF loading data files

needs NCO for option -m

Powered by Geode Systems and RAMADDA