[gradsusr] labeling latitude/longitude on polar stereograph map projection

[Mikovitz, Colleen (LARC-E302)[SCIENCE SYSTEMS AND APPLICATIONS, INC]]

Kyle,

I have recently done it this way for a particular area.  You’ll probably need to modify the placement of xpos and ypos in the polargrid script.

'map nps '
'set lat 58 66'
'set lon 116 136’
(do your mapping here)

'set gxout contour'
'interval 58 66 2'
'set ccolor 1'
'set cstyle 5'
'set cthick 3'
'set clab off'
'd lat'
'interval 116 136 4'
'set ccolor 1'
'set cstyle 5'
'set cthick 3'
'set clab off'
'd lon'
'polargrid 116 136 4 58 66 2’

polargrid.gs is
* * * * * * * * * * * * * * * * * * * * * * * * * *
* Draws grid lines at specified intervals
* call is
* polargrid lonmin lonmax lonincr latmin latmax latincr
function polargrid(args)
* parse arguments
lonmin = subwrd(args,1)
lonmax = subwrd(args,2)
lonincr = subwrd(args,3)
latmin = subwrd(args,4)
latmax = subwrd(args,5)
latincr = subwrd(args,6)

midlat = (latmax+latmin)/2
quarterlat = (latmax+midlat)/2
'q dims'
xline = sublin(result,2)
yline = sublin(result,3)
lonstart = subwrd(xline,6)
lonend = subwrd(xline,8)
midlon = (lonstart+lonend)/2
*midlon = (lonmax+lonmin)/2
istart=latmin
while (istart <=latmax)
  'q w2xy 'midlon' ' istart
  xpos=subwrd(result,3)
  ypos=subwrd(result,6)
  'draw string 'xpos' 'ypos' 'istart
  istart=istart+latincr
endwhile
istart=lonmin
while (istart <=lonmax)
  'q w2xy 'istart' ' quarterlat+0.5
  xpos=subwrd(result,3)
  ypos=subwrd(result,6)
  'draw string 'xpos' 'ypos' 'istart
  istart=istart+lonincr
endwhile
return
*********************************
interval.gs is
*********************************

function interval(arg)
*limitation that you must know for sure the start & end will be correct with the interval
if (arg = '')
   say 'Enter Min Max and interval'
   pull cmdline
else
   cmdline=arg
endif
min=subwrd(cmdline,1)
max=subwrd(cmdline,2)
int=subwrd(cmdline,3)

number=(max-min)/int+1
lev = min

i=1
while (i <= number-1)
   next = min+int*i
   lev = lev%' '%next
   i=i+1
endwhile

'set clevs ' lev
*********************************

Colleen

On Sep 11, 2014, at 6:29 PM, Kyle Clem <Kyle.Clem at vuw.ac.nz<mailto:Kyle.Clem at vuw.ac.nz>> wrote:

Hi Everyone,

I am trying to display lat/lon labels on a polar stereographic map projection. Does anyone know how to do this?

Thanks,
Kyle

----------------------
Kyle Clem, M.Sc.
PhD Student
School of Geography, Environment and Earth Sciences
Victoria University of Wellington
Cotton Building, Room 222
kyle.clem at vuw.ac.nz<mailto:kyle.clem at vuw.ac.nz>
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140912/2d6c8dca/attachment.html