diff --git a/notebooks/SeismicFWD.ipynb b/notebooks/SeismicFWD.ipynb index 21cc7bf..7133c75 100644 --- a/notebooks/SeismicFWD.ipynb +++ b/notebooks/SeismicFWD.ipynb @@ -1,7 +1,7 @@ { "metadata": { "name": "", - "signature": "sha256:43cef8feaa23d5020ffbd5874e54f80865dd6632477d24adbab1e45f5867176c" + "signature": "sha256:b8a51c40a152664559ff5f36e935013a6d1019641ab852ea89f76ef72cce6b44" }, "nbformat": 3, "nbformat_minor": 0, @@ -27,24 +27,24 @@ ] } ], - "prompt_number": 23 + "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "# Step1: set time\n", - "time = np.linspace(0, 0.04, 2**9)\n", + "time = np.linspace(0, 0.03, 2**9)\n", "dt = time[1]-time[0]\n", "# Step2: set Tx and Rx\n", "options={'tlag':0.0025, 'fmain':400} # You need to set waveform to set Tx\n", - "rx = AcousticRx(np.vstack((np.r_[0, 1], np.r_[0, 1])))\n", + "rx = AcousticRx(np.r_[10., 1])\n", "tx = AcousticTx(np.r_[0, 1], time, [rx], **options)\n", "# Step3: set survey (pass txlist)\n", "survey = SurveyAcoustic([tx])\n", "wave = tx.RickerWavelet()\n", "# Step4: set mesh\n", - "cs = 0.5\n", + "cs = 0.25\n", "hx = np.ones(150)*cs\n", "hy = np.ones(150)*cs\n", "mesh = Mesh.TensorMesh([hx, hy], 'CC')\n", @@ -67,14 +67,14 @@ "text": [ "You are good to go:)\n", ">> Stability information\n", - " dt: 7.83e-05 s\n", - " Optimal dt: 1.25e-04 s\n", - " Cell per wavelength (G): 4.00e+01\n", + " dt: 5.87e-05 s\n", + " Optimal dt: 6.25e-05 s\n", + " Cell per wavelength (G): 8.00e+01\n", " Optimal G: 1.60e+01\n" ] } ], - "prompt_number": 32 + "prompt_number": 2 }, { "cell_type": "code", @@ -91,12 +91,12 @@ "stream": "stdout", "text": [ ">> Start Computing Acoustic Wave\n", - ">> dt: 7.83e-05 s\n", - ">> Optimal dt: 1.25e-04 s\n", + ">> dt: 5.87e-05 s\n", + ">> Optimal dt: 6.25e-05 s\n", ">> Main frequency, fmain: 1.00e+02 Hz\n", - ">> Cell per wavelength (G): 4.00e+01\n", + ">> Cell per wavelength (G): 8.00e+01\n", " Tx at ( 0.00, 0.00): 1/ 1\n", - ">>Elapsed time: 2.55e+00 s" + ">>Elapsed time: 2.48e+00 s" ] }, { @@ -107,7 +107,7 @@ ] } ], - "prompt_number": 25 + "prompt_number": 3 }, { "cell_type": "code", @@ -119,26 +119,27 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 26 + "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ - "print data[0].shape" + "20./(0.0119-0.0025)" ], "language": "python", "metadata": {}, "outputs": [ { - "output_type": "stream", - "stream": "stdout", + "metadata": {}, + "output_type": "pyout", + "prompt_number": 5, "text": [ - "(2L, 512L)\n" + "2127.659574468085" ] } ], - "prompt_number": 27 + "prompt_number": 5 }, { "cell_type": "code", @@ -303,12 +304,12 @@ "\n", "\n", "