polar stereographic plot

[Jennifer M. Adams]

On Jan 16, 2007, at 10:10 AM, Henrique Barbosa wrote:

> Hi,
> I am trying to plot geopotential height anomalies in
> a polar stereographic projection.... but when I do:
>
> 'set mproj nps'
> 'draw hgeo'
>
> I get the correct plot, but without axis and
> labels. Moreover, although the plot has a circular
> shape, grads draws a rectangular plot boundary.
>
> Is it possible to draw the longitude labels and
> a circular boundary around a polar plot?
The circular boundary is done with 'set frame circle'. To draw grid
lines on a non-lat-lon projection, you have to contour the internal
variables 'lat' and 'lon' as if you were drawing data variables. All
the usual contour control settings apply. Here's something I use for
COLA's weather maps:

* inside the main script you call it like this:
gridlines(10,5)


* * * * * * * * * * * * * * * * * * * * * * * * * *
* Draws grid lines at specified intervals
function gridlines(lonincr,latincr)
'set rgb 97 150 150 150'
'set gxout contour'
'set ccolor 97'
'set cint 'lonincr
'set cmin -175'
'set cmax 180'
'set clab off'
'set cstyle 5'
'set cthick 1'
'd lon'
'set ccolor 97'
'set clab off'
'set cmin -75'
'set cmax 75'
'set cint 'latincr
'set cstyle 5'
'set cthick 1'
'd lat'

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Beltsville, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070116/c9326123/attachment.html