MathGroup Archive 2002

[Date Index] [Thread Index] [Author Index]

Search the Archive

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






  • Prev by Date: Re: Optional
  • Next by Date: Re: Re: Particular structure 2
  • Previous by thread: RE: Question - Karmarkar's Algorithm
  • Next by thread: RE: combination of two ContourPlots - impossible?