------------------------------------------------------------------------
*From:* Jeff Whitaker <jswhit@...146...>
*To:* antonv <vasilescu_anton@...9...>
*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
>>> 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.Whitaker@...259... <mailto:Jeffrey.S.Whitaker@…259…>
>>>> 325 Broadway Office : Skaggs Research Cntr 1D-113
>>>> Boulder, CO, USA 80303-3328 Web : Jeffrey S. Whitaker: NOAA Physical Sciences Laboratory
>>>>
>>>> ------------------------------------------------------------------------------
>>>> 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>
>>>> matplotlib-users List Signup and Options
>>>>
>>>> >>> >> -- Jeffrey S. Whitaker Phone : (303)497-6313
>> Meteorologist FAX : (303)497-6449
>> NOAA/OAR/PSD R/PSD1 Email : Jeffrey.S.Whitaker@...259... <mailto:Jeffrey.S.Whitaker@…259…>
>> 325 Broadway Office : Skaggs Research Cntr 1D-113
>> Boulder, CO, USA 80303-3328 Web : Jeffrey S. Whitaker: NOAA Physical Sciences Laboratory
>>
>> ------------------------------------------------------------------------------
>> 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>
>> matplotlib-users List Signup and Options
>>
>> >
>
-- Jeffrey S. Whitaker Phone : (303)497-6313
Meteorologist FAX : (303)497-6449
NOAA/OAR/PSD R/PSD1 Email : Jeffrey.S.Whitaker@...259... <mailto:Jeffrey.S.Whitaker@…259…>
325 Broadway Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web : Jeffrey S. Whitaker: NOAA Physical Sciences Laboratory