Ooups, forgot the code producing all that fron the previous post.
=====
set_plot, 'PS'
PageParam = PSWINDOW(/LANDSCAPE, /CM, PAGESIZE='A4')
PS_START, filename = 'name.ps', /LANDSCAPE, /EUROPEAN
!P.MULTI=[0, 2, 2, 0, 0]
step = 1.0
levels = cgConLevels(NLevels=nlevels, MinValue=minimum, MaxValue=maximum)
LoadCT, 0
cgLoadCT, 25, /Reverse, /Brewer, NColors=nlevels-1, Bottom=1
TVLCT, cgColor('white', /Triple), nlevels
map = Obj_New('cgMap', 'Equirectangular', Ellipsoid=19, $
XRange=xrange, YRange=yrange, /LatLon_Ranges, $
Position=[0.1, 0.1, 0.9, 0.8], $
Limit=[Min(lat), Min(lon), Max(lat), Max(lon)])
;map -> Erase
map -> Draw
cgContour, data, lon, lat, /Fill, /Overplot, /Outline, $
missingvalue=!Values.F_NAN, C_Colors=Indgen(nlevels)+1, label =0, Levels=levels, Title='title.ps', Map=map, OutColor='white'
cgMap_Grid, Map=map, /Box
cgMap_Continents, Map=map, Color='black', thick=5
cgColorbar, Range=[Min(levels), Max(levels)+step], OOB_High=nlevels, $
NColors=nlevels-1, Bottom=1, /Discrete, /Top, $
Tlocation='TOP', Position=[0.1, 0.96, 0.9, 0.99]
PS_END
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
