I don't know of a general code. THe following
loops will do it
-------------------------------------
1088 format(i2,x,'N')
1087 format(i2,x,'S')
do 15,lat=-60,60,20
do 15 lon=27,360,20
flon=lon
flat=lat
call maptra(flat,flon,uu,vv)
if(lat .lt. 0) then
write(ctemp,1087)-lat
else
write(ctemp,1088) lat
endif
call gslwsc(2.0)
call plchmq( uu,vv,ctemp,1.,0.,0.)
call gslwsc(1.0)
15 continue
do 16,lat=-33,60,20
do 16 lon=0,360,20
flon=lon
flat=lat
call maptra(flat,flon,uu,vv)
if( lon .lt. 180) then
write(ctemp,1089)lon
1089 format(i3,'E')
1091 format(i3,'W')
else
write(ctemp,1091)(360-lon)
endif
call gslwsc(2.0)
call plchmq( uu,vv,ctemp,1.,0.,0.)
call gslwsc(1.0)
16 continue
-----------------------------------------------------------
You may want to make the numbers bigger.. these ar really small.
The fourth argument in plchmq is the number size.
I inserted this just after the mapdrw call in mpex05.f
for the stereographic projection and got reasonable labels.
(I also called mapsti with increment of 10 rather than 15
in mpex05.f))
(you also need to create character*4 ctemp)
This archive was generated by hypermail 2b29 : Wed Jun 28 2000 - 09:40:34 MDT