Greetings,
I’m having considerable difficulties in getting the
axes grid to display they way I would like it to. Basically, I’m
wanting to grid to be shown for all of the x/y-ticks displayed, but
only an irregular number show up. What I’m wanting to know is if my
code is wrong or am I’m completely overlooking something. My code is
below and an image can be found at http://forwarn.org/images/USDIFFtest.png.
Also, I might as well ask…is there a way to manually specify the
labels that are displayed on the colorbar? For instance, in the
example below, my values range from -5 to 4…how would I go about
making sure the last 5 was included in the colorbar? I’ll admit I
haven’t looked into this much as of now since the first question was /
is more pressing.
Thanks!
-Patrick Marsh
from future import division
import numpy as np
import pylab as P
import copy
def open_text(REGION,FIELD):
···
global DATAPATH1
global DATAPATH2
xaxis1,yaxis1,freq1 = np.loadtxt(DATAPATH1+‘bin_CUR_%s_%s.txt’ % (REGION, FIELD), unpack=True)
xaxis2,yaxis2,freq2 = np.loadtxt(DATAPATH2+‘bin_FUT_%s_%s.txt’ % (REGION, FIELD), unpack=True)
return(xaxis1,yaxis1,freq1,xaxis2,yaxis2,freq2)
if name == ‘main’:
#Main Program
TIME = ‘CUR’
REGION = ‘PLAINS’
FIELD = ‘CAPE-SHR6’
PATHS TO THE INPUT (DATA) AND OUTPUT (IMAGES) DIRECTORIES
DATAPATH1 = ‘E:/RESEARCH/CCSM_OUTPUT/CURRENT/BIN/’
DATAPATH2 = ‘E:/RESEARCH/CCSM_OUTPUT/FUTURE/BIN/’
IMAGEPATH = ‘E:/RESEARCH/IMAGES/%s/%s/%s/’ % (REGION, TIME, FIELD)
xaxis1,yaxis1,freq1,xaxis2,yaxis2,freq2 = open_text(REGION, FIELD)
TRANSFORM INTO THE CORRECT FORMAT
xaxis1 = xaxis1.reshape(51,-1).T
yaxis1 = yaxis1.reshape(51,-1).T
freq1 = freq1.reshape(51,-1).T
freq2 = freq1.reshape(51,-1).T
freq1 = freq1 / (12 * 12 * 20) # (121220 for PLAINS)…(294320 for EUR)
freq2 = freq2 / (12 * 12 * 20) # (121220 for PLAINS)…(294320 for EUR)
diff = freq1 - freq2
xaxis = np.log10(xaxis1)
yaxis = np.log10(yaxis1)
xaxis = xaxis1
yaxis = yaxis1
levs = np.arange(0,10,.1)
levs = np.arange(-5,5,.1)
xticks = [1,2,5,10,20,50,100,200,500,1000,2000,5000]
yticks = [1,2,5,10,20,50,100]
x = np.arange(1,10000,1)
P.figure(figsize=(8,8))
CS = P.contourf(xaxis,yaxis,diff,levs)
P.xlabel(’\nCAPE’)
P.ylabel(‘0-6 km Shear’)
P.xlim(0,3)
P.ylim(0,2)
P.xticks(np.log10(xticks),xticks,)
P.yticks(np.log10(yticks),yticks,)
P.plot(np.log10(x),(np.log10(10000/x)),color=‘black’, linewidth=2)
P.colorbar(drawedges=False, orientation=‘horizontal’)
P.grid(True)
P.gca().xaxis.grid(True, which=‘minor’,xdata=xticks,ydata=yticks)
P.title(‘Difference\nCentral U.S. Normalized Frequency Distribution’)
P.savefig(‘USDIFFtest.png’,dpi = 200)