Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado
Grady Daub wrote: > I'm making orbit plots. There will be 2 plots, one each of the north and > south poles. > > MAP_SET,/cont,color=0,-90,0,0,pos=[.1,.75,.9,.90],scale=100e 6 > > would show you what I mean. (Using a window shaped like a sheet of > paper.) > > Is there an easy (or complicated, for that matter) way to determine the > SCALE= factor to set for the plot to include only up to a certain > latitude circle? > > What I need is a rectangular box (centered on the pole), about 20 > degrees tall and exactly 90 degrees wide (those in latitude). The data > will be plotted in a similar box, below the map plot, with the data > corresponding to the orbit position directly above it. > > This all could probably be done with LIMIT= , but, then I'd have to > worry about aspect ratio's. It seems simpler to make a box with > POSITION= , then scale the earth until it shows the +/- 45 degrees > latitude circle at the left/right edges. (top/bottom range doesn't > really matter.) When you use the SCALE keyword to MAP_SET, the only things defining the map coverage are the size of your graphics window, and the POSITION keyword if it is used. Consider the following examples: window, /free, xsize=300, ysize=300 map_set, 25, -81, scale=10e6, /cont window, /free, xsize=600, ysize=600 map_set, 25, -81, scale=10e6, /cont window, /free, xsize=600, ysize=600 map_set, 25, -81, scale=10e6, /cont, position=[0.1, 0.1, 0.5, 0.5] Notice how in each case, the scale of the map remains exactly the same (Florida is the same size), but the extent of the mapped area changes. So you just need to experiment with window size and POSITION, e.g. window, /free, xsize=714, ysize=924 map_set, -90, 0, scale=50e6, pos=[0.05,0.05,0.95,0.45], /cont map_grid, /label map_set, 90, 0, scale=50e6, pos=[0.05,0.55,0.95,0.95], /cont, /noerase map_grid, /label Cheers, Liam. -- Liam E. Gumley Space Science and Engineering Center, UW-Madison http://cimss.ssec.wisc.edu/~gumley
Grady Daub wrote: > > I'm making orbit plots. There will be 2 plots, one each of the north and > south poles. > > MAP_SET,/cont,color=0,-90,0,0,pos=[.1,.75,.9,.90],scale=100e 6 > > would show you what I mean. (Using a window shaped like a sheet of > paper.) > [...] > This all could probably be done with LIMIT= , but, then I'd have to > worry about aspect ratio's. isn't that what the /isotropic keyword in map_set is for? Martin -- |||||||||||||||\\\\\\\\\\\\\-------------------///////////// //||||||||||||||| Martin Schultz, DEAS, Harvard University, 29 Oxford St., Pierce 109, Cambridge, MA 02138 phone (617) 496 8318 fax (617) 495 4551 e-mail mgs@io.harvard.edu web http://www-as/people/staff/mgs/ ********* ADDRESS CHANGE : AFTER JULY 15, 1999 *************** Max-Planck-Institut fuer Meteorologie >>> N E W <<< Bundesallee 55 >>> N E W <<< 20147 Hamburg >>> N E W <<< Germany >>> N E W <<< phone (+49 40) 41173 - 0 >>> N E W <<< email martin.schultz@dkrz.de >>> N E W <<< (email to mgs@io.harvard.edu will be forwarded) ************************************************************ *******