antonv wrote: > Hi Jeff, > > Here is the whole script: > http://old.nabble.com/file/p27826457/_test_shiftgrid.py _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 >>>> 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 >>>> >>>> >>>> >>>> >>> >>> >> -- >> 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 >> >> >> > > -- 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