-
Notifications
You must be signed in to change notification settings - Fork 0
/
SCS_search_settings.py
executable file
·80 lines (62 loc) · 3.47 KB
/
SCS_search_settings.py
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
#!/usr/bin/env python
'''File name: Denver-Water_XWT.py
Author: Andreas Prein
E-mail: [email protected]
Date created: 16.04.2018
Date last modified: 16.04.2018
##############################################################
Purpos:
Contains the setup for extreme weather typing (XWT) for
Denver Water watersheds
'''
from pdb import set_trace as stop
import numpy as np
import os
import pandas as pd
import datetime
# ###################################################
dStartDayPR=datetime.datetime(1990, 1, 1, 0) # earliest is 1990
dStopDayPR=datetime.datetime(2016, 12, 31, 23) # latest is 2016 #(2016, 6, 30,23)
rgdTime = pd.date_range(dStartDayPR, end=dStopDayPR, freq='d')
iMonths=[3,4,5] # [1,2,3,10,11,12] # [4,5,6,7,8,9]
# ---------
# Setup clustering algorithm
ClusterMeth='hdbscan' # current options are ['HandK','hdbscan']
ClusterBreakup = 1 # breakes up clusters that are unproportionally large (only for hdbscan)
# ---------
#sPlotDir='/glade/scratch/prein/projects/2020_SCS_XWT/plots/'
#sDataDir='/glade/scratch/prein/projects/2020_SCS_XWT/data/'
sPlotDir='/glade/work/fossell/preevents/XWT_Kates_code/plot_test/'
sDataDir='/glade/work/fossell/preevents/XWT_Kates_code/data_test/'
# Severe convective subregions
DW_Regions=['SPL.poly']
sRegion=0
Region=DW_Regions[sRegion]
sSubregionSCS='/glade/u/home/fossell/MET/met/data/poly/'+DW_Regions[sRegion]
rgsWTvars= ['var151', 'u', 'v', 'UV', 'tcw', 'FLX', 'q', 'q', 'z', 'cape', '06Shear', 'cape-shear', 'cape']
VarsFullName=['PSL', 'U850', 'V850', 'UV850', 'PW', 'FLX850', 'Q850', 'Q500', 'ZG500', 'CAPE', 'VS06', 'CAPE-Shear', 'MAX_CAPE']
rgsWTfolders=['/glade/scratch/prein/ERA-Interim/PSL/fin_PSL-sfc_ERA-Interim_12-0_',\
'/glade/scratch/prein/ERA-Interim/UV850/fin_UV850-sfc_ERA-Interim__',\
'/glade/scratch/prein/ERA-Interim/UV850/fin_UV850-sfc_ERA-Interim__',\
'/glade/scratch/prein/ERA-Interim/UV850/fin_UV850-sfc_ERA-Interim__',\
'/glade/scratch/prein/ERA-Interim/TCW/fin_TCW-sfc_ERA-Interim_12-0_',\
'/glade/scratch/prein/ERA-Interim/UV850xTCW/fin_FLX-pl_ERA-Interim_',\
'/glade/scratch/prein/ERA-Interim/Q850/Q850_daymean_',\
'/glade/scratch/prein/ERA-Interim/Q500/Q500_daymean_',\
'/glade/scratch/prein/ERA-Interim/Z500/Z500_daymean_',\
'/glade/scratch/prein/ERA-Interim/CAPE_ECMWF/fin_CAPE-ECMWF-sfc_ERA-Interim_12-9_',\
'/glade/scratch/prein/ERA-Interim/Shear06/Shear06/fin_0-6km-Shear_ERA-Interim_12-0_',\
'/glade/scratch/fossell/ERA-Interim/CAPE-Shear/MAX/fin_CAPE-Shear_ERA-Interim_12-9_',\
'/glade/scratch/fossell/ERA-Interim/CAPE_ECMWF/MAX_CAPE/fin_CAPE-MAX-sfc_ERA-Interim_12-9_']
# rgsWTvars= ['var151','u', 'v']
# VarsFullName=['PSL','U850','V850']
# rgsWTfolders=['/glade/scratch/prein/ERA-Interim/PSL/fin_PSL-sfc_ERA-Interim_12-0_',\
# '/glade/scratch/prein/ERA-Interim/UV850/fin_UV850-sfc_ERA-Interim__',\
# '/glade/scratch/prein/ERA-Interim/UV850/fin_UV850-sfc_ERA-Interim__']
rgrNrOfExtremes=[10,20,30] #[6,10,15,30]
WT_Domains=['S'] #['S','M','L'] # ['S','M','L','XXL']
DomDegreeAdd=[2] #[2,5,10] # [2,5,10,20]
Annual_Cycle=['1'] # '1' means that the annual cycle gets removed before clustering; '0' nothing is done
SpatialSmoothing=[0.5]
Metrics=['PSS', 'MRD', 'MRR', 'APR', 'PEX']
Dimensions=['Variables', 'Extreme Nr.', 'Domain Size', 'Annual Cycle', 'Smoothing', 'Split Sample', 'Metrics']