[gradsusr] plotting over 0E from netCDF

James T. Potemra jimp at hawaii.edu
Wed Jun 29 13:49:52 EDT 2016


Kathy,

I think the problem is you don't have a continuous grid across 0/360.
If your data were, for example, evenly spaced and you had a single value
at 360 then you could do 340 to 415.  In this case you don't have any
longitudes greater than 360 defined.  You might need to either regrid
or use an external tool (e.g., nco) to redefine longitude.  Of course you
could also use "parea" to plot the western side and eastern sides 
individually.

Jim

On 6/29/16 4:57 AM, Cathy Smith (NOAA Affiliate) wrote:
> Hi all
> This should be an easy question but I'm stumped. I'm trying to plot
> Africa (20W to 55E). the file I'm plotting from is netCDF and the
> longitudes are 0 to 359.27
>
> e.g.
>    lon = 0, 0.7031252, 1.40625, 2.109376, 2.812501, 3.515626, 4.218751,
>       4.921877, 5.625002, 6.328127, 7.031252, 7.734378, 8.437503, 9.140628,
>       9.843754, 10.54688, 11.25, 11.95313, 12.65625, 13.35938, 14.0625,
> ......
>       354.3751, 355.0782, 355.7814, 356.4845, 357.1876, 357.8907, 358.5939,
>       359.297 ;
>
> my grads code has
>
>    if( proj = 'africa' )
>       'set mproj latlon'
>       'set mpdset mres'
>       'set lat -40 45'
>       'set lon -20 55'
>       'set grid on'
>       'set frame on'
>       'set ylint 5'
>        istate=0
>     endif
>
> which plots the data from 0 to 55 only and not the data between 340 and 360.
>
> I also tried 340 to 415 instead of -20 to 55 and that plotted the data
> from 340 to 0 and not 0 to 55.
>
> What am I doing wrong? I know this can be done in GrADS.
>
> Thanks
> Cathy
>



More information about the gradsusr mailing list