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
