Received: by 10.180.87.100 with SMTP id w4mr39117wiz.3.1349942191655; Thu, 11 Oct 2012 00:56:31 -0700 (PDT) Path: q10ni65082452wif.0!nntp.google.com!feeder3.cambriumusenet.nl!82.197.223.108.MISMATCH!feeder2.cambriumusenet.nl!feed.tweaknews.nl!194.109.133.84.MISMATCH!newsfeed.xs4all.nl!newsfeed5.news.xs4all.nl!xs4all!newspeer1.nac.net!border4.nntp.dca.giganews.com!border2.nntp.dca.giganews.com!news.glorb.com!npeer01.iad.highwinds-media.com!news.highwinds-media.com!feed-me.highwinds-media.com!border3.nntp.dca.giganews.com!border1.nntp.dca.giganews.com!nntp.giganews.com!local2.nntp.dca.giganews.com!news.giganews.com.POSTED!not-for-mail NNTP-Posting-Date: Thu, 04 Oct 2012 17:33:07 -0500 From: David Fanning Newsgroups: comp.lang.idl-pvwave Subject: Re: contour plots on map with Coyote Graphics System in postscript Date: Thu, 4 Oct 2012 16:33:00 -0600 Message-ID: References: <7b1bc81a-2543-4ec8-95dc-6f160d9aa840@googlegroups.com> Organization: Fanning Software Consulting MIME-Version: 1.0 User-Agent: MicroPlanet-Gravity/3.0.4 Lines: 45 X-Usenet-Provider: http://www.giganews.com X-Trace: sv3-EpVlpZXfHajvU2EzZNqfz/kPlIQhNMYU3djM4V84lC2+MiNgDEA8QOqbDFR7XiNTo1ZazKfP2t4KZLm!SzTre4rSYqwYUJzFmZay9a/n3rqyIB5EFsZobHqVjbeF1gj5TLhhiHkrmx3Pi8NL4c71D40p7Q== X-Complaints-To: abuse@giganews.com X-DMCA-Notifications: http://www.giganews.com/info/dmca.html X-Abuse-and-DMCA-Info: Please be sure to forward a copy of ALL headers X-Abuse-and-DMCA-Info: Otherwise we will be unable to process your complaint properly X-Postfilter: 1.3.40 X-Original-Bytes: 2508 X-Received-Bytes: 2649 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit sasha.mozg...@gmail.com writes: > Do you by any chance see what am I doing wrong? I would do it like this: PS_START, filename = work_path+'profiles_correlation_on_map.ps', $ /LANDSCAPE, /EUROPEAN levels = cgConLevels(Float(data), NLevels=nlevels+1, $ MinValue=Floor(Min(data)), STEP=step, Factor=1) LoadCT, 0 cgLoadCT, 2, /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, $ CENTER_LON=center_lon, $ Position=[0.1, 0.1, 0.9, 0.8], $ Limit=[Min(lat), Min(lon), Max(lat), Max(lon)]) map -> Draw cgContour, data, lon, lat, /Cell_Fill, /Overplot, /Outline, $ C_Colors=Indgen(nlevels)+1, Levels=levels, $ Map=map, OutColor='White' cgMap_Grid, Map=map, /Box cgMap_Continents, Map=map, Color='tomato' cgColorbar, Range=[Min(levels), Max(levels)-step], OOB_High=nlevels, $ NColors=nlevels-1, Bottom=1, /Discrete, $ Title='Temperature $\deg$') + 'K', $ Position=[0.1, 0.91, 0.9, 0.95] PS_END Cheers, David -- David Fanning, Ph.D. Fanning Software Consulting, Inc. Coyote's Guide to IDL Programming: http://www.dfanning.com/ Sepore ma de ni thue. ("Perhaps thos speakest truth.")