[Thread Prev][Thread Next][Index]
Re: log contour intervals
Hello Alice,
unfortunately there is no direct solution. I had the same problem some
time ago. Here I forward the answer I got via the mailing list.
It contains an example that works also with the z-axis.
Hope it helps !
Best regards ...
Patrick
---------------- forwarded message -----------------------------------------
Hi Patrick,
Here is a more general technique that gives you full control over all
details of the axis labelling. I believe that it will work on Z axes, as
well.
- steve
===========================
! suppose you have some data on a grid x=1:100:1, y=1:100:1
define axis/x=1:100:1 xax
define axis/y=1:100:1 yax
let/title="test data" a = LOG(x[gx=xax]/8)+y[gy=yax]/100 ! some made up
data
! plot untransformed data
contour a
! define a new axis with the LOG transform of the X coordinates
define axis/from/name=xlog/x LOG(x[gx=xax])
! plot quick and dirty LOG-transformed data in the upper plot
! NOTE: X axis labels will be LOG values rather than exponents
set window/new
set view upper
contour a[gx=xlog@asn]
! *** here's how to do a nicely labelled plot ***
set view lower
! base plot with large tics, only, unlabelled
ppl axlabp 0
contour/xlimits=0:2:1/title=" " a[gx=xlog@asn]
! add 2 decades of minor tics
repeat/j=1:2:1 (repeat/i=1:10:1 ppl aline 1,`LOG(i)+j-1`,0,`LOG(i)+j-1`,-4)
! label the major tics with powers of 10
repeat/j=1:3:1 label `j-1`,-18,0,0,.1 "@SR10^`j-1`"
! label 2,3,5 of the minor tics
repeat/j=1:2:1 label `LOG(2)+j-1`,-15,0,0,.08 "@SR2"
repeat/j=1:2:1 label `LOG(3)+j-1`,-15,0,0,.08 "@SR3"
repeat/j=1:2:1 label `LOG(5)+j-1`,-15,0,0,.08 "@SR5"
============================
---- end of forwarded message -----------------------------------------------
---------------------------------------------------------------
Patrick Joeckel
Max Planck Institute for Chemistry
Joh.-J. Becher Weg 27
55128 Mainz phone:++49-6131-305452
Germany fax :++49-6131-305436
e-mail:joeckel@mpch-mainz.mpg.de
http://www.mpch-mainz.mpg.de/~joeckel
---------------------------------------------------------------
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement