bennyistanto

# -*- coding: utf-8 -*-
"""
NAME
    download_era5land_precipitation.py
DESCRIPTION
    Download ERA5-Land daily total precipitation and save as 1 year data in 1 netCDF file.
REQUIREMENT
    You must registered as CDS Copernicus user, and access your profile to get UID and Api Key
EXAMPLES
    python download_era5land_precipitation.py

bennyistanto

# -*- coding: utf-8 -*-
"""
NAME
    chirps_precip_accumulation.py
    Calculate running dekad of 1-,3-,6-,9-,12- and 24-month accumulation precipitation using dekad data
DESCRIPTION
    Input data for this script will use https://data.chc.ucsb.edu/products/CHIRPS-2.0/global_dekad/tifs/
REQUIREMENT
    ArcGIS must installed before using this script, as it required arcpy module.
EXAMPLES

bennyistanto

# -*- coding: utf-8 -*-
"""
NAME
    chirps_longterm_stats.py
    CHIRPS precipitation statistics data 1991-2020, long-term average, max, min and stdev
DESCRIPTION
    Input data for this script will use data generated from chirps_precip_accumulation.py
REQUIREMENT
    ArcGIS must installed before using this script, as it required arcpy module.
EXAMPLES

bennyistanto

# -*- coding: utf-8 -*-
"""
NAME
    spi_monthly_blend.py
    Global CHIRPS's SPI-based dry/wet indicators blend.
DESCRIPTION
    Input data for this script will use CHIRPS SPI monthly data in GeoTIFF format
    This experimental SPI blends integrate several SPI scales into a single product.
    The combines 3-, 6-, 9-, 12-, and 24-month SPI to estimate the overall dry/wet conditions.
METHOD

bennyistanto

# -*- coding: utf-8 -*-
"""
NAME
    chirps_dekad_update.py
    Update Global CHIRPS's dekad collection in a folder, gunzip, rename and compress
DESCRIPTION
    This script will do:
    1. Check the exisiting file in the Dekad folder, match the date `..._YYYY.MM.D.tif` 
      with collection files in BASE_URL
    2. If the file exist in the Dekad folder, then skipped, proceed to next date. 

bennyistanto

bennyistanto

"""
NAME
    download_era5_agromet.py
DESCRIPTION
    Download ERA5 Agrometeorological Indicators and save as 1 month data in 1 netCDF file.
REQUIREMENT
    You must registered as CDS Copernicus user, and access your profile to get UID and Api Key
    You need to install `cdsapi` python package
EXAMPLES
    python download_era5_agromet.py

bennyistanto

From below script:

ds = gdal.Open('../SPI/mar/orig_geotiff/mar_cli_chirps_precip_month1_198101.tif') # Data location

a = ds.ReadAsArray()
nlat,nlon = np.shape(a)

b = ds.GetGeoTransform() #bbox, interval
lon = np.arange(nlon)*b[1]+b[0]

bennyistanto

bennyistanto