combination of two ContourPlots - impossible?
- To: mathgroup at smc.vnet.net
- Subject: [mg33714] combination of two ContourPlots - impossible?
- From: Harald von der Osten-Woldenburg <hvdosten at lb.netic.de>
- Date: Mon, 8 Apr 2002 03:04:55 -0400 (EDT)
- Organization: LF.net GmbH, Stuttgart, Germany
- Sender: owner-wri-mathgroup at wolfram.com
Hi,
first of all: Thanks a lot for your help. It was really a stupid
question ( ** --> ^), sorry...
But I have still a problem and whatever I tried in the last hours, I had
no success. I want two combine two ContourPlots (plot1: negative values,
plot2: positive values) and I did the following:
-----------START----------------------
Needs["Graphics`Colors`"];
Needs["Graphics`Graphics`"];
majorTicks =
Table[{x, x, {0.015, 0}, {GrayLevel[0.], Thickness[0.003]}}, {x,
-1.5,
2.5, .5}];
minorTicks =
Table[{x, "", {0.012, 0}, {GrayLevel[0.], Thickness[0.0025]}}, {x,
-1.5,
2.5, .25}];
obenunten = Join[majorTicks, minorTicks];
majorTicks =
Table[{x, x, {0.015, 0}, {GrayLevel[0.], Thickness[0.003]}}, {x, -1,
0, .2}];
minorTicks =
Table[{x, "", {0.012, 0}, {GrayLevel[0.], Thickness[0.0025]}}, {x,
-1,
0, .1}];
rechtslinks = Join[majorTicks, minorTicks];
a = 0.5; b = 0.5; c = 0.5; y = 0;
negpart =
ContourPlot[
Evaluate[(x*(x - a) + y*y + z*z)/(Sqrt[(x*x + y*y + z*z)^3]*
Sqrt[((x - a)^2 + y*y + z*z)^3]) - ((x - a)*(x - a - b -
c) +
y*y + z*z)/(Sqrt[((x - a)^2 + y*y + z*z)^3]*
Sqrt[((x - a - b - c)^2 + y*y + z*z)^3]) - (x*(x - a - b)
+
y*y + z*z)/(Sqrt[(x*x + y*y + z*z)^3]*
Sqrt[((x - a - b)^2 + y*y + z*z)^3]) + ((x - a - b)*(x - a
- b -
c) + y*y + z*z)/(Sqrt[((x - a - b)^2 + y*y +
z*z)^3]*
Sqrt[((x - a - b - c)^2 + y*y + z*z)^3])], {x, -1, 2.5},
{z,
0.001, -1}, ContourShading -> False, PlotRange -> {-500, 00},
Contours -> 12, PlotPoints -> 150,
ContourStyle -> Dashing[{0.005, 0.006}], AspectRatio -> Automatic,
ImageSize -> {800, 450}, FrameStyle -> {Thickness[0.003]}, Frame ->
True,
FrameTicks -> {obenunten, rechtslinks, obenunten, rechtslinks},
FrameLabel -> {"\n Profile [m]", "\n depth [m]"}, RotateLabel ->
True,
PlotRange -> {{-1, 2.5}, {0, -1}}, AspectRatio -> .5,
DefaultFont -> {"Times-Bold", 14}, FormatType -> OutputForm]
pospart =
ContourPlot[
Evaluate[(x*(x - a) + y*y + z*z)/(Sqrt[(x*x + y*y + z*z)^3]*
Sqrt[((x - a)^2 + y*y + z*z)^3]) - ((x - a)*(x - a - b -
c) +
y*y + z*z)/(Sqrt[((x - a)^2 + y*y + z*z)^3]*
Sqrt[((x - a - b - c)^2 + y*y + z*z)^3]) - (x*(x - a - b)
+
y*y + z*z)/(Sqrt[(x*x + y*y + z*z)^3]*
Sqrt[((x - a - b)^2 + y*y + z*z)^3]) + ((x - a - b)*(x - a
- b -
c) + y*y + z*z)/(Sqrt[((x - a - b)^2 + y*y +
z*z)^3]*
Sqrt[((x - a - b - c)^2 + y*y + z*z)^3])], {x, -1, 2.5},
{z,
0.001, -1}, ContourShading -> False, PlotRange -> {00, 500},
Contours -> 12, PlotPoints -> 150, AspectRatio -> Automatic,
ImageSize -> {800, 450}, FrameStyle -> {Thickness[0.003]}, Frame ->
True,
FrameTicks -> {obenunten, rechtslinks, obenunten, rechtslinks},
FrameLabel -> {"\n Profile [m]", "\n depth [m]"}, RotateLabel ->
True,
PlotRange -> {{-1, 2.5}, {0, -1}}, AspectRatio -> .5,
DefaultFont -> {"Times-Bold", 14}, FormatType -> OutputForm]
ois = Show[negpart, pospart]
-------THE END----------------------
why is this possible, or is there another way to do this?
Again thanks a lot,
Harry
--
Harald von der Osten-Woldenburg
Geophysical Prospection of Archaeological Sites
National Heritage Department of Baden-Wuerttemberg
Silberburgstrasse 193, D-70178 Stuttgart
Fax Office: +49-(0)711-1694-707
http://www.lb.netic.de/hvdosten : Geomagnetics, Geoelectrics, Radar, EMI