*embarrassed*
Turns out I was looking at a "compressed" version of the image that didn't
have all the grid lines on. Examining the full version the image looks
correct. I apologize for asking what is not a really dumb question.
However, the second question about the colorbar is still puzzling.
-Patrick
On Sat, Aug 23, 2008 at 9:49 AM, Patrick Marsh <
[EMAIL PROTECTED] <[EMAIL PROTECTED]>>wrote:
> 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) # (12*12*20 for PLAINS).....(29*43*20
> for EUR)
> freq2 = freq2 / (12 * 12 * 20) # (12*12*20 for PLAINS).....(29*43*20
> 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)
>
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
