You're right, the lons was getting overwritten! Thanks again for the help!
Anton
________________________________
From: Jeff Whitaker <jsw...@fastmail.fm>
To: Anton Vasilescu <vasilescu_an...@yahoo.com>
Cc: matplotlib-users@lists.sourceforge.net
Sent: Wed, March 10, 2010 8:57:16 AM
Subject: Re: [Matplotlib-users] shiftgrid 'cyclic point not included' error
Anton Vasilescu wrote:
> Hi Jeff,
>
> Thank you for looking at the code for me! It works perfect for the first
> image but when it tries to do the second it errors out saying: "lon0 outside
> of range of lonsin". Any idea why is that happening? Is it happening for you
> too?
>
> Thanks,
> Anton
Anton: Probably because the original value of lons has been overwritten with
the shifted grid longitudes. You need to save the original lons array to pass
to shiftgrid for the second grid.
-Jeff
>
> ------------------------------------------------------------------------
> *From:* Jeff Whitaker <jsw...@fastmail.fm>
> *To:* antonv <vasilescu_an...@yahoo.com>
> *Cc:* matplotlib-users@lists.sourceforge.net
> *Sent:* Mon, March 8, 2010 12:50:48 PM
> *Subject:* Re: [Matplotlib-users] shiftgrid 'cyclic point not included' error
>
> antonv wrote:
> > Hi Jeff,
> >
> > Here is the whole script:
> > http://old.nabble.com/file/p27826457/_test_shiftgrid.py
> > _test_shiftgrid.py<http://_test_shiftgrid.py> And here are the grib data
> > files:
> > http://socalscubainfo.com/downloads/GRIB2_files.zip
> >
> > Also, if there is not too much to ask, is there a sample anywhere that shows
> > how to use shiftgrid for a mercator projection? If you look at the script, I
> > was initially trying to make the images using the mercator projection (there
> > is a commented basemap function there) but it would display just the data
> > starting from the origin. I was wondering how would I be able to display
> > let's say the whole globe but shift the grid 180 so that great britain would
> > be in the center of the map rather than being cut at the edge.
> >
> > Thanks,
> > Anton
> Anton: You weren't updating the lons and lats to be consistent with the
> shifted grid. I changed the lines
>
> #z, lons[i] = shiftgrid(180.0, z, lons[i], start=False, cyclic=360.0)
> #zq, lons_a[i] = shiftgrid(180.0, zq, lons_a[i], start=False,
> cyclic=360.0)
>
> in your script to
>
> lons1 = lons[0,:]
> z, lonsnew = shiftgrid(180.0, z, lons1, start=False)
> zq, lonsnew = shiftgrid(180.0, zq, lons1, start=False)
> lons, lonsnew = shiftgrid(180.0, lons, lons1, start=False)
> lats, lonsnew = shiftgrid(180.0, lats, lons1, start=False)
>
> and it works fine.
>
> -Jeff
> >
> >
> >
> > Jeff Whitaker wrote:
> > >> antonv wrote:
> >> >>> Hi Jeff,
> >>>
> >>> Thanks for the quick reply! I've updated to the latest shiftgrid and now
> >>> I
> >>> get the chart :D
> >>> I still have a small issue as there is a half degree sliver missing from
> >>> the
> >>> mapping at the greenwich meridian. Any idea how to solve that issue?
> >>> Here is an image showing that (above england):
> >>> http://old.nabble.com/file/p27824785/Global-HTSGW_DIRPW-0.png
> >>> Thanks,
> >>> Anton
> >>> >> Anton: Looks like the values at the Greenwich meridian got
> >>> messed up somehow. Can you post your script somewhere so I can try it
> >>> out?
> >>
> >> -Jeff
> >> >>> Jeff Whitaker wrote:
> >>> >>>> antonv wrote:
> >>>> >>>>> Hi all,
> >>>>>
> >>>>> I am getting a 'cyclic point not included' error using shiftgrid on a
> >>>>> grib
> >>>>> file from NOAA. I've checked and it seems that the file has the lons
> >>>>> from
> >>>>> 0
> >>>>> to 359.5 every .5 degree but I get that error whenever I try to run
> >>>>> this.
> >>>>> What I need to achieve is to shift the grid from 0 / 360 to -180 / 180
> >>>>> so
> >>>>> that I can create a map over Great Britain. I am using a Cassini
> >>>>> projection.
> >>>>> >>>> Anton: You can add the wraparound, or cyclic
> >>>>> point using the addcyclic function. Or, you can update to the latest
> >>>>> svn version of basemap which does not have this restriction.
> >>>>
> >>>>
> >>>> -Jeff
> >>>>
> >>>>
> >>>> -- Jeffrey S. Whitaker Phone : (303)497-6313
> >>>> Meteorologist FAX : (303)497-6449
> >>>> NOAA/OAR/PSD R/PSD1 Email : jeffrey.s.whita...@noaa.gov
> >>>> <mailto:jeffrey.s.whita...@noaa.gov>
> >>>> 325 Broadway Office : Skaggs Research Cntr 1D-113
> >>>> Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg
> >>>>
> >>>>
> >>>> ------------------------------------------------------------------------------
> >>>> Download Intel® Parallel Studio Eval
> >>>> Try the new software tools for yourself. Speed compiling, find bugs
> >>>> proactively, and fine-tune applications for parallel performance.
> >>>> See why Intel Parallel Studio got high marks during beta.
> >>>> http://p.sf.net/sfu/intel-sw-dev
> >>>> _______________________________________________
> >>>> Matplotlib-users mailing list
> >>>> Matplotlib-users@lists.sourceforge.net
> >>>> <mailto:Matplotlib-users@lists.sourceforge.net>
> >>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> >>>>
> >>>>
> >>>> >>> >> -- Jeffrey S. Whitaker Phone :
> >>>> (303)497-6313
> >> Meteorologist FAX : (303)497-6449
> >> NOAA/OAR/PSD R/PSD1 Email : jeffrey.s.whita...@noaa.gov
> >> <mailto:jeffrey.s.whita...@noaa.gov>
> >> 325 Broadway Office : Skaggs Research Cntr 1D-113
> >> Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg
> >>
> >>
> >> ------------------------------------------------------------------------------
> >> Download Intel® Parallel Studio Eval
> >> Try the new software tools for yourself. Speed compiling, find bugs
> >> proactively, and fine-tune applications for parallel performance.
> >> See why Intel Parallel Studio got high marks during beta.
> >> http://p.sf.net/sfu/intel-sw-dev
> >> _______________________________________________
> >> Matplotlib-users mailing list
> >> Matplotlib-users@lists.sourceforge.net
> >> <mailto:Matplotlib-users@lists.sourceforge.net>
> >> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> >>
> >>
> >> >
> >
>
> -- Jeffrey S. Whitaker Phone : (303)497-6313
> Meteorologist FAX : (303)497-6449
> NOAA/OAR/PSD R/PSD1 Email : jeffrey.s.whita...@noaa.gov
> <mailto:jeffrey.s.whita...@noaa.gov>
> 325 Broadway Office : Skaggs Research Cntr 1D-113
> Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg
>
>
-- Jeffrey S. Whitaker Phone : (303)497-6313
Meteorologist FAX : (303)497-6449
NOAA/OAR/PSD R/PSD1 Email : jeffrey.s.whita...@noaa.gov
325 Broadway Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg
------------------------------------------------------------------------------
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
