#### Written by Kim Wood
#### VERSIONS: python 3.7.3, matplotlib 3.0.3, numpy 1.16.3
####
import matplotlib.colors as mcolors
import numpy as np
colors = np.genfromtxt('rainbowIRsummer.txt', dtype=float) # adapated from SIFT's "Rainbow IR Summer" color scheme
cmap = mcolors.LinearSegmentedColormap.from_list('cmap', colors)
ccols = np.concatenate((np.linspace(0.09,0.8983,71),np.linspace(0.905,1.,61)))
cmap2 = cmap(ccols)
cmap = mcolors.LinearSegmentedColormap.from_list('cmap', cmap2)
levs = np.arange(-90.,41.,1.) # range of shading values in degrees Celsius
clevs = np.arange(-90.,40.1,10.) # range of values to label colorbar in degrees Celsius