radar_processing_scripts/ERA5_download_pressurelvl_by_year-month.py at main · JulianGiles/radar_processing_scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
import cdsapi
import sys
from optparse import OptionParser
parser = OptionParser(usage="usage: %prog [year] [variable] ",\
                      version='%prog v1.0')
# general options
parser.add_option("-q", "--quiet",\
                   action="store_false", dest="verbose", default=True,\
                  help="don't print status messages to stdout")
(options, args) = parser.parse_args()
c = cdsapi.Client()
try:
    iyear = str(args[0])
except IndexError:
    sys.stderr.write("____ERROR____: Year not defined!!")
    sys.exit(200)
try:
    imonth = str(args[1])
except IndexError:
    sys.stderr.write("____ERROR____: Month not defined!!")
    sys.exit(200)
try:
    ivar=str(args[2])
except IndexError:
    sys.stderr.write("____ERROR____: Variable not defined!!")
    sys.exit(200)
try:
    lvl=str(args[3])
except IndexError:
    sys.stderr.write("____ERROR____: Level not defined!!")
    sys.exit(200)
try:
    north=float(args[4])
except IndexError:
    sys.stderr.write("____ERROR____: North not defined!!")
    sys.exit(200)
try:
    south=float(args[5])
except IndexError:
    sys.stderr.write("____ERROR____: South not defined!!")
    sys.exit(200)
try:
    west=float(args[6])
except IndexError:
    sys.stderr.write("____ERROR____: West not defined!!")
    sys.exit(200)
try:
    east=float(args[7])
except IndexError:
    sys.stderr.write("____ERROR____: East not defined!!")
    sys.exit(200)


if lvl == "all":
    lvls = ['1', '2', '3', '5', '7', '10',
            '20', '30', '50', '70', '100', '125',
            '150', '175', '200', '225', '250', '300',
            '350', '400', '450', '500', '550', '600',
            '650', '700', '750', '775', '800', '825',
            '850', '875', '900', '925', '950', '975', '1000'
            ]
else:
    lvls=[lvl]

if imonth == "all":
    months = [
            '01','02','03','04','05','06','07','08','09','10','11','12'
            ]
else:
    months=[imonth]

c.retrieve(
    'reanalysis-era5-pressure-levels',
    {
        'product_type':'reanalysis',
        'data_format':'netcdf',
        'variable':ivar,
	'pressure_level': lvls,
        'year':[
            iyear,
        ],
        'month':months,
        'day':[
            '01','02','03',
            '04','05','06',
            '07','08','09',
            '10','11','12',
            '13','14','15',
            '16','17','18',
            '19','20','21',
            '22','23','24',
            '25','26','27',
            '28','29','30',
            '31'
        ],
        'time':[
            '00:00','01:00','02:00',
            '03:00','04:00','05:00',
            '06:00','07:00','08:00',
            '09:00','10:00','11:00',
            '12:00','13:00','14:00',
            '15:00','16:00','17:00',
            '18:00','19:00','20:00',
            '21:00','22:00','23:00'
        ],
        'area':[
            north, west, south, east,
        ],
    },
    f'{ivar}_{lvl}_{iyear}-{imonth}.nc')