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:
The individual time series will be downloaded, resampled, and summed up before the bilinear interpolation.
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.
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:
per station, per year.
- 40-60 sec for 3-hourly data (e.g. NTAL, NTOL, NTAL24h, NTAL+HYDL, ...)
- 3 sec for 24-hourly (e.g. HYDL, SLEL)
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]
||one of NTAL, NTOL, HYDL or
NTAL24h, NTOL24h or
NTAL+NTOL, NTOL+SLEL, NTAL+HYDL, NTAL+NTOL+HYDL, NTAL+NTOL+HYDL+SLEL
||CF or CM
||1st epoch (format: yyyy-mm-dd)
||last epoch (format: yyyy-mm-dd)
||latitude of station [degree]
||longitude of station [degree]
||name of stationlistfile containing station coordinates (no header, 2 or 3 columns: [stationname/id] longitude latitude, UNIX format)
||columnorder "latlon" or "lonlat", columnorder in stationlistname file (default=latlon)
"xyz" for cartesian coordinates, coverted by WGS84 parameters
||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)
||specify file names for output (->outfilename.stationID.txt)
default is >extractlatlon<.
> 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 30 KB|
|Kind: ||Software Tool|
|Capabilities:||extract station time series|
from netCDF loading data files
needs NCO for option -m