prl900

from osgeo import gdal
from osgeo import ogr
from osgeo import gdalconst

#This raster is the model for our output (CRS, extent)
ndsm = 'geotif_output_extent_proj.tif'
#This shapefile contains the features we want to burn
shp = 'shapefile_features_fields.shp'

data = gdal.Open(ndsm, gdalconst.GA_ReadOnly)

prl900

prl900

import os
import requests # get the requsts library from https://github.com/requests/requests
import datetime
import urllib.request
from datetime import timedelta, datetime

def daterange(start_date, end_date):
    n = 0

    while True:

prl900

prl900

prl900

prl900

prl900

#!/g/data/xc0/software/python/miniconda2/bin/python

from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
server.retrieve({
    "class": "mc",
    "dataset": "cams_gfas",
    "date": "2019-01-01/to/2019-12-31",
    "expver": "0001",
    "levtype": "sfc",

prl900

# [START gae_python37_app]
from flask import Flask
from flask import send_file

from datetime import datetime
from datetime import timedelta
import pytz
import io
import os
import bz2

prl900