| map_set in idl 5 [message #11529] |
Mon, 30 March 1998 00:00 |
Torben Schmith
Messages: 2
Registered: April 1998
|
Junior Member |
|
|
I have discovered an error in the idl-5.0.3 version of map_set when
using Lambert projection. The map is clipped in the horzontal direction
and the outermost latitude-circel interrupted. When '/noborder' is
omitted evrything is perfect. Also in idl 4 evrything works ok.
And here comes a test program illustrating the probling
**********
; define plotting environment
type = 'x'
if ( type eq 'x') then begin
set_plot,'x'
device, retain=2
window,0,xsize=800,ysize=800
endif
if ( type eq 'ps' ) then begin
set_plot,'ps'
device, filename="/home/tsc/tmp/plot.ps", $
/portrait, bits=24, xsize=20,ysize=26, xoffset=1.0,yoffset=1.0,
$
font_size=12
endif
loadct, 0
p0lat = 90.
p0lon = 0.
rot = 0.
; define sub-area
alat_s = 30.
alat_n = 90.
alon_w = -180.
alon_e = 180.
; plot
!p.multi = [0,1,1]
map_set, p0lat, p0lon, rot, /lamb, /isotropic, $
/noborder, /advance, $
limit = [alat_s, alon_w, alat_n, alon_e], $
xmargin=[3,3], ymargin=[3,3], $
title= 'Test'
map_continents, /coasts
map_grid, latdel=30, londel=90, $
glinestyle=0, glinethick=1
;exit
if ( type eq 'ps' ) then begin
device, /close
endif
end
**************************
--
************************************************************ *****
Torben Schmith
Danish Meteorological Institute
Meteorological and Oceanographic Research Division
Lyngbyvej 100
DK-2100 Copenhagen
DENMARK
Phone: +45 39 15 75 00
Direct: +45 39 15 74 44
Fax: +45 39 15 74 60
E-mail: tsc@dmi.dk
************************************************************ *****
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
