From ea0c4719fe4d9d70f9f8973a98fc2419e9366298 Mon Sep 17 00:00:00 2001 From: seogi Date: Thu, 27 Feb 2014 09:15:32 -0800 Subject: [PATCH] added play notebook! --- simpegPF/notebooks/Play_analytic_fun.ipynb | 97 +++++++++++++++++++++- 1 file changed, 93 insertions(+), 4 deletions(-) diff --git a/simpegPF/notebooks/Play_analytic_fun.ipynb b/simpegPF/notebooks/Play_analytic_fun.ipynb index 580213cd..90321e3a 100644 --- a/simpegPF/notebooks/Play_analytic_fun.ipynb +++ b/simpegPF/notebooks/Play_analytic_fun.ipynb @@ -183,7 +183,18 @@ "cell_type": "code", "collapsed": false, "input": [ - "bx.size" + "import simpegPF as PF" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 13 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "PF.Magnetics.MagneticsDiffSecondary" ], "language": "python", "metadata": {}, @@ -191,13 +202,91 @@ { "metadata": {}, "output_type": "pyout", - "prompt_number": 7, + "prompt_number": 14, "text": [ - "524880" + "simpegPF.Magnetics.MagneticsDiffSecondary" ] } ], - "prompt_number": 7 + "prompt_number": 14 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import simpegPF.Magnetics.MagneticsDiffSecondary" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "ename": "ImportError", + "evalue": "No module named MagneticsDiffSecondary", + "output_type": "pyerr", + "traceback": [ + "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mImportError\u001b[0m Traceback (most recent call last)", + "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[1;32mimport\u001b[0m \u001b[0msimpegPF\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mMagnetics\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mMagneticsDiffSecondary\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m", + "\u001b[1;31mImportError\u001b[0m: No module named MagneticsDiffSecondary" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Quick test for converting from (inc., dec., intensity) to (bx, by, bz)" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import numpy as np" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 5 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "Btot = 59013.8\n", + "Inc = 82.71\n", + "Dec = 14.4\n", + "Bxa = 1862.7\n", + "Bya = 7256.3\n", + "Bza = 58536.4" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 8 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "Bx = Btot*np.cos(Inc/180.*np.pi)*np.sin(Dec/180.*np.pi)\n", + "By = Btot*np.cos(Inc/180.*np.pi)*np.cos(Dec/180.*np.pi)\n", + "Bz = -Btot*np.sin(Inc/180.*np.pi)\n", + "print Bx, By, Bz" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "1862.27664792 7253.08866839 -58536.7681111\n" + ] + } + ], + "prompt_number": 10 }, { "cell_type": "code",