{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Square wave generator for MECSim\n", "=======\n", "\n", "The aim here is to mimic a square wave based on an input dc ramp from E_start to E_rev with a given scan rate. The ramp will be linear from E_start to E_rev and back to E_start with these parameters changed along with the scan rate.\n", "\n", "Output from script is the N_ac signals to add to the ``Master.inp`` file to perform this hack. MECSim is not designed for square waves as the function has no derivatives at the edges of the square wave, so a smooth version must be used. This notebook will output these terms as well as a PNG image file.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Parameters used for square wave(s)\n", "-----" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "E_start = 1.5 # Starting voltage (V)\n", "E_rev = 0.5 # reverse voltage (V)\n", "n_cycles = 1 # number of cycles of E_start to E_rev\n", "vscan = 1.0 # V/s\n", "N_squares = 2 # number of square waves\n", "N_ac = 100 # number of AC sinusoids to add to make ramp look like a square wave" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Import python packages as define some functions**" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def ConstSquareWave(Amp = 1, T = 1, nres = 10):\n", " a = np.zeros(nres)\n", " f = np.zeros(nres)\n", " for i in range(0,nres):\n", " n = i*2+1\n", " a[i] = 4.0/(np.pi * n)*1000 # output in mV\n", " f[i] = n / (2.0 * T)\n", " return a, f" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def EvaluateEAC(ACAmp, ACFreq, time):\n", " x = ACAmp*1.0e-3*np.sin(2.0*np.pi*ACFreq*time)\n", " return x\n", "def EvaluateRamp(Es, Er, vs, time):\n", " DE = Es - Er\n", " T = abs(DE) / vs\n", " if(time < T):\n", " x = Es - DE*time/T\n", " else:\n", " x = Er + DE*(time - T)/T\n", " return x" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def PrintSquareWave(a,f,nres=10):\n", " print \"Copy the following lines into Master.inp\"\n", " print nres,\" ! number of species (need n lines below)\"\n", " for i in range(0,nres):\n", " print a[i],\", \",f[i],\" ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\"\n", " return" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Re-calculate E_start (= E_end), E_rev and the scan rate**\n", "\n", "Minor change in these inputs is required to get a clean fit." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [], "source": [ "E_mid = 0.5*( E_start + E_rev )\n", "DE = E_start - E_rev\n", "T = n_cycles * DE / vscan\n", "amp, freq = ConstSquareWave(DE, T/N_squares, N_ac)\n", "DE_new = 1.0e-6\n", "E_start_new = E_mid + DE_new/2.0\n", "E_rev_new = E_mid - DE_new/2.0\n", "vscan_new = DE_new / T\n", "period_square = 2*T #/ N_squares" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Plot the resulting approximate square wave to file**\n", "\n", "More accurate the larger the number of ac sinusoids added (N_ac)." ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm0AAAHuCAYAAADNxztVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuUZWV95vHvj24IaqOIOg0ISBRI1CwjYhDFS8dbEBPM\nhUkgOjFOosZIMMQxGnUCWWOijllZxuswLjWOEY2jCUEDGjW2tzgochVBQEFBkSAIAo3aDb/5Y++C\nQ1OX962u857atb+ftWr1PnV2nfNW9VvPefZ+zzkVmYkkSZJWt51mPQBJkiQtzdImSZI0AJY2SZKk\nAbC0SZIkDYClTZIkaQAsbZIkSQMws9IWEftGxKcj4sKI+GpEHD/PPpsi4saIOKf/ePUsxipJkjRr\n62d431uBEzLz3IjYAHwlIj6RmRdtt99nMvOoGYxPkiRp1ZjZmbbM/F5mnttv3wxcBOw9z67RdGCS\nJEmr0Kp4TltE7A8cDJy53VUJPC4izouI0yPiYa3HJkmStBrMcnkUgH5p9EPAS/ozbpPOBvbNzC0R\n8QzgVOCgeW7Dv8UlSZIGIzOrVxJjln97NCJ2Bj4KnJGZbyzY/3LgkMy8frvP53K+eY1PRJyUmSfN\nehwaBueLSjlXVGO5vWWWrx4N4J3A1xYqbBGxsd+PiDiUrmReP9++kiRJa9ksl0cPB54DnB8R5/Sf\neyWwH0BmngwcDbwoIrYBW4BjZjFQSZKkWZvp8uhKcXlUpSJiU2ZunvU4NAzOF5VyrqjGcnuLpU2S\nJKmhwT2nTZIkSeUsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYCl\nTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAs\nbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAsbZIkSQNg\naZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAsbZIkSQNgaZMkSRoA\nS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIA\nWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTaMRwf+NYN2sx6FhiOCxEbxs1uPQ6me2qJXIzFmP\nYYdFRGZmzHocWt0iSGBDJrfMeixa/SL4Z+CoTMwWLcpsUa3l9hbPtGkUJo6Cb5/pQDQkW2c9AK1+\nZotasrRpLH6q/9ezJir1k1kPQINgtqgZS5vGYpf+X4NVpTzTphJmi5qxtGksDFbVsrSphNmiZixt\nGguXMFTL5VGVMFvUjKVNY+HRsGp5pk0lzBY1Y2nTWOyy9C7SXXimTSXMFjVjadNYuIShWpY2lTBb\n1IylTWPhEoZquTyqEmaLmrG0aSwMVtXyTJtKmC1qZs2Utgh+P4KHznocWn0i+B3gF+YuznIsGpSt\nYLZoYWaLWls/6wGsoF8HtgAXzXogWnWO4s5ANVhVam551GzRQswWNbWWStu9gHvMehBale7BnXPD\nYFWpueVRs0ULMVvU1JpZHgXuDRwewYGzHohWjwieAWwE7jvrsWhwbuv/NVt0N2aLZmEtlbb7As8D\nNkdwv1kPRrMVwa4RbABOBw7hzmD1aFil5uaK2aI7mC2apZmVtojYNyI+HREXRsRXI+L4BfZ7U0Rc\nGhHnRcTBi9zk3C/ObsD3I3huBLuu+MC16kVwAPBZ4FsTnzZYVWuytIHZMnpmi2ZtlmfatgInZObD\ngcOAF0fEXV6hFRFHAgdk5oHAC4C3L3J79+7/3a3/9++A10WQERwfwYYIfjbijuu1BkSwWwSPiuAh\nEfxVBNcCl9K9omuPiV3n5ofBqlJzc8VsGSGzRavRzF6IkJnfA77Xb98cERcBe3PXV2gdBbyn3+fM\niNg9IjZm5jWFd/NL/b9/C/wasAkgglOBI4APAbcANwIfAR4DXAt8kq5IbgP+H/AQupJ5Hd27X98M\n3BP4MfAfdM9r2AZcDTwASOAGuieoRn/d7f2/u/TbW/rbSOBWYNd+eyuwbvJHRffcmvXdj4GtEezc\n/Uzu2M7+NoO7Bsdtc7eVyba5r+vHsX7i9ufMlfit/Tgjkx/1ZxWi/3536feb+7rb+7ED3ATs3t/n\ndcD9gJ2BH9I9mXt9/7Xr+s/t03/dN4BH9j+v84CD+/2vAn62v90D+v+HTwHvpAvKVwOv6W/jIih6\nWwaDVaUWmytmy12zZf3EeMwWaVoyc+YfwP50p5s3bPf5jwCPm7j8SeCQeb4+IWf1cfPE9jcmtr87\nz743TGyfNbH9hYntK+b5ussmtjcv8HXfhtw2MZ7bIc+buP5TC9z3ddv9m5BnL3B/lyzxPZ0/sX3O\nMn+G0/7Ya9bz3Y9hfEAe13Bemi2Lf09mix9r6qOrX/VfN/O3/IiIDXRHpS/JzJvn22W7yzn/LZ00\nsb2p/5jX9+iO4PYGvkB3tHYo3Rm9h9Idzf4r8AS6o9tzgOcCpwLXA79Ht0y7G/AU4IN0pfO+/b6H\nAN/s7+cIuqK5le7I/F+ADcDDgS/RHdmtpzvquxi4jO7I7/H92AAe22/fC9gPuAT4d7q3I7gOuAD4\nOvCjft8v9rd5MPBluqPGjcDlwFfojvxvoDvqvKC/nccDn6c70v/5fr9/pjuy/Rbwtf77uQF4Yn8f\ntwFP739Wu/X398V+jLvSHeHOfV/fB47tb3Onfvu9/b6HAB+jO/twG3Au8Ox+7FcDzwf+qf+5PQ04\nrf9/OqD/vh/S/x9Ks2a2mC3SvCJiE4sUk2Izbpo7Ax8H/niB6/8XcMzE5YuBjfM11omjnVv7I8Gv\nQx4C+Y+Qb+r2yy9C/hJkQH4Y8h6Q94b8/f76jZCP6LfvD3mPfnu3O+8rd57Yjlm39SF+TP7c5n6e\nkLtC7tJv7w25rt9+av//tQHyBf3nngT58n77f0O+rN++CvIEyCMn5sO12x0N7z3r79+PYXxA/tF2\nc8dsWeUfZosfQ/lgmWfaov/i5iIi6I5Ar8vMExbY50jguMw8MiIOA96YmYfNs19OnIDbi+4dzI/N\n5AkRPBe4MZNTI3hAd3/cPo3vSe1FcG9gaya3RvA3wJvpnkt0DbAncDTwlokv2SeT77QfqYYmguPp\nnrM2x2wZEbNF0xQRmZnVz4Oc5fLo4cBzgPMj4pz+c6+kO61NZp6cmadHxJERcRndL8vzCm73FrrT\n7ld0t9O9kKHfvnblhq/VIJMfTmz/CUAE6+iWP66nmw+TfLKwSm0/V24Bvo3ZMgpmi1ajWb569PMU\nvOVIZh5XedNb6N708PTljEvDl8ltETwwk4xg8nmSN2Cwqtz2c2UL3XPH/mUGY9EqYLZo1tbSX0T4\nc2BLJrfNLf/OekCanYn//5vpzrw+u982WFUq6N6TzWzRHcwWzdJaKm3/wN1PV0u3ADdncgrdEx8N\nVpUKujMoZovmY7aoubVU2r4Ldz7HROpdQveWMrDg28VIC0rMFs3PbFFzM3v16Epa7qswNC4RfAt4\nUmb3RHJpMRG8FHjg3JPQpYWYLaq13N6yls60SUtxCUM1As+gqIzZoiYsbRoTg1U1LG0qZbaoCUub\nxsRgVQ1Lm0qZLWrC0qYxMVhVw9KmUmaLmrC0aUx8AFYt54xKOE/UhKVNY+PRsEo5V1TD+aKps7Rp\nTFzCUA2XR1XKbFETljaNicGqGpY2lTJb1ISlTWNisKqGpU2lzBY1YWnTmBisqmFpUymzRU1Y2jQm\nBqtqWNpUymxRE5Y2jYkPwKrlnFEJ54masLRpbDwaVinnimo4XzR1ljaNiUsYquHyqEqZLWrC0qYx\nMVhVw9KmUmaLmrC0aUwMVtWwtKmU2aImLG0aE4NVNSxtKmW2qAlLm8bEB2DVcs6ohPNETVjaNDYe\nDauUc0U1nC+aOkubxsQlDNVweVSlzBY1YWnTmBisqmFpUymzRU1Y2jQmBqtqWNpUymxRE5Y2jYnB\nqhqWNpUyW9SEpU1jYrCqhqVNpcwWNWFp05j4AKxazhmVcJ6oCUubxsajYZVyrqiG80VTZ2nTmLiE\noRouj6qU2aImLG0aE4NVNSxtKmW2qAlLm8bEYFUNS5tKmS1qwtKmMTFYVcPSplJmi5qwtGlMfABW\nLeeMSjhP1ISlTWPj0bBKOVdUw/miqbO0aUxcwlANl0dVymxRE5Y2jYnBqhqWNpUyW9SEpU1jYrCq\nhqVNpcwWNWFp05gYrKphaVMps0VNWNo0JgaraljaVMpsUROWNo2JD8Cq5ZxRCeeJmrC0aWw8GlYp\n54pqOF80dZY2jYlLGKrh8qhKmS1qwtKmMTFYVcPSplJmi5qwtGlMDFbVsLSplNmiJixtGhODVTUs\nbSpltqgJS5vGxGBVDUubSpktasLSJkmSNACWNo2JR8Oq4Zk2lTJb1ISlTWNisKqGpU2lzBY1YWnT\nmBisqmFpUymzRU1Y2jQmBqtqWNpUymxRE5Y2jYnBqhqWNpUyW9SEpU1j4gOwajlnVMJ5oiYsbRob\nj4ZVyrmiGs4XTZ2lTWPiEoZquDyqUmaLmrC0aUwMVtWwtKmU2aImLG0aE4NVNSxtKmW2qAlLm8bE\nYFUNS5tKmS1qwtKmMTFYVcPSplJmi5qwtGlMfABWLeeMSjhP1ISlTWPj0bBKOVdUw/miqbO0aUxc\nwlANl0dVymxRE5Y2jYnBqhqWNpUyW9TETEtbRLwrIq6JiAsWuH5TRNwYEef0H69uPUatKQaralja\nVMpsURPrZ3z/7wbeDPyfRfb5TGYe1Wg8WtsMVtWwtKmU2aImZnqmLTM/B/xgid38RdBK8QFYtZwz\nKuE8UROr/TltCTwuIs6LiNMj4mGzHpAGz4MAlXKuqIbzRVM36+XRpZwN7JuZWyLiGcCpwEHz7RgR\nJ01c3JyZm6c/PA2MSxiq4fKoSpktWlREbAI27ejtrOrSlpk3TWyfERFvi4g9MvP6efY9qengNEQG\nq2pY2lTKbNGi+hNJm+cuR8SJy7mdVb08GhEbIyL67UOBmK+wSYUMVtWwtKmU2aImZnqmLSLeDzwJ\nuH9EXAmcCOwMkJknA0cDL4qIbcAW4JhZjVVrgsGqGpY2lTJb1MRMS1tmHrvE9W8F3tpoOFr7DFbV\nsLSplNmiJlb18qi0wnwAVi3njEo4T9SEpU1j49GwSjlXVMP5oqmztGlMXMJQDZdHVcpsUROWNo2J\nwaoaljaVMlvUhKVNY2KwqoalTaXMFjVhadOYGKyqYWlTKbNFTVjaNCY+AKuWc0YlnCdqwtKmsfFo\nWKWcK6rhfNHUWdo0Ji5hqIbLoypltqgJS5vGxGBVDUubSpktasLSpjExWFXD0qZSZouasLRpTAxW\n1bC0qZTZoiYsbRoTg1U1LG0qZbaoCUubxsQHYNVyzqiE80RNWNo0Nh4Nq5RzRTWcL5o6S5vGxCUM\n1XB5VKXMFjVhadOYGKyqYWlTKbNFTVjaNCYGq2pY2lTKbFETljaNicGqGpY2lTJb1ISlTWPiA7Bq\nOWdUwnmiJixtGhuPhlXKuaIazhdNnaVNY+IShmq4PKpSZouasLRpTAxW1bC0qZTZoiYsbRoTg1U1\nLG0qZbaoCUubxsRgVQ1Lm0qZLWrC0qYxMVhVw9KmUmaLmrC0aUx8AFYt54xKOE/UhKVNY+PRsEo5\nV1TD+aKps7RpTFzCUA2XR1XKbFETljaNicGqGpY2lTJb1ISlTWNisKqGpU2lzBY1YWnTmBisqmFp\nUymzRU1Y2jQmPgCrlnNGJZwnasLSprHxaFilnCuq4XzR1FnaNCYuYaiGy6MqZbaoCUubxsRgVQ1L\nm0qZLWrC0qYxMVhVw9KmUmaLmrC0aUwMVtWwtKmU2aImLG0aE4NVNSxtKmW2qAlLm8bEB2DVcs6o\nhPNETVjaNDYeDauUc0U1nC+aOkubxsQlDNVweVSlzBY1YWnTmBisqmFpUymzRU1Y2jQmBqtqWNpU\nymxRE+uX2iEidgceC+xPNzGvAL6YmTdOdWTSyks8UFE5S5tKmS1qYsFJFhFPiIjTgM8CxwD70RW3\nY4HPRcRpEfH4JqOUVoZHw6phaVMps0VNLHam7deAl2bmpfNdGREHAX8AfH4aA5OmwAdg1XLOqITz\nRE0sVtrekJlXL3RlZl4C/MnKD0maKo+GVcq5ohrOF03dYmvw50TEJyPi9/rntUlD5xKGarg8qlJm\ni5pYrLTtA/w18ATg6xHxzxFxTETco83QpBVnsKqGpU2lzBY1sWBpy8xtmfmxzPxduhchvBt4FnB5\nRJzSaHzSSjJYVcPSplJmi5ooeolyZv4Y+BpwEXAT8NBpDkqaEoNVNSxtKmW2qIlFS1tE7BcRfxoR\nZwMfBdYBv5KZBzcZnbSyfABWLeeMSjhP1MSCrx6NiH+ne17bB4HnZ+ZXmo1Kmh6PhlXKuaIazhdN\n3WJv+fEK4POZeXurwUhT5hKGarg8qlJmi5pYbHn0ycADFroyIvaKiL9Y+SFJU2OwqoalTaXMFjWx\n2Jm2s4APRMQuwNnA1XSTck/gUcCP6d4SRBoKg1U1LG0qZbaoiQVLW2Z+FPhoROwLHE73th/Q/dmq\n12fmVQ3GJ60kg1U1LG0qZbaoicXOtAGQmVcCH2gwFmnaDFbVsLSplNmiJorep01aI3wAVi3njEo4\nT9SEpU1j49GwSjlXVMP5oqmztGlMXMJQDZdHVcpsURNLlraI+JmI+FREXNhffkREvHr6Q5NWnMGq\nGpY2lTJb1ETJmbZ3AK8EftJfvgA4dmojkqbHYFUNS5tKmS1qoqS03TMzz5y7kJkJbJ3ekKSpMVhV\nw9KmUmaLmigpbddGxAFzFyLiaLo32t0hEfGuiLgmIi5YZJ83RcSlEXFeRPhH6rWjfABWLeeMSjhP\n1ERJaTsOOBn42Yj4LnAC8KIVuO93A0csdGVEHAkckJkHAi8A3r4C9yl5NKxSzhXVcL5o6kreXPcb\nwFMi4l7ATpl500rccWZ+LiL2X2SXo4D39PueGRG7R8TGzLxmJe5fo+QShmq4PKpSZouaWLK0RcRL\nmQiuiAC4EfhKZp47vaHxQODKictXAfsAljYtl8GqGpY2lTJb1ETJ8ughwB/Qlah9gBcCzwDeEREv\nn+LY4O6/BAaodoTBqhqWNpUyW9TEkmfagH2BR2XmzQAR8efA6cCTgK8Ar5/S2L7T3/ecffrPzSsi\nTpq4uDkzN09nWBowg1U1LG0qZbZoURGxCdi0o7dTUtoewJ3v0Qbd231szMwtEfGjHR3AIk6jexHE\nByLiMOCGxZ7PlpknTXEsWhsMVtWwtKmU2aJF9SeSNs9djogTl3M7JaXtfcCZEXEq3aT8FeCU/oUJ\nX1vOnQJExPvpztbdPyKuBE4EdgbIzJMz8/SIODIiLgNuAZ633PuSej4Aq5ZzRiWcJ2qi5NWj/yMi\nPgYcTjcxX5iZZ/VXP3u5d5yZS/5Vhcw8brm3Ly3Ao2GVcq6ohvNFU1dypo3M/HJEfBvYFciI2C8z\nvz3doUkrziUM1XB5VKXMFjVR8gfjj4qIS4Fv0q3HXgGcMd1hSVNhsKqGpU2lzBY1UfKWH68BHgtc\nkpk/DTwFOHPxL5FWJYNVNSxtKmW2qImS0rY1M78P7BQR6zLz08CjpzwuaRoMVtWwtKmU2aImSp7T\n9oOI2A34HPC+iPgP4ObpDkuaCh+AVcs5oxLOEzVRcqbtWcAWuj8U/zHgMrq3/ZCGyKNhlXKuqIbz\nRVNXUtr+PDNvy8ytmfl3mfkm4E+nPTBpClzCUA2XR1XKbFETJaXt6fN87siVHojUgMGqGpY2lTJb\n1MSCz2mLiBcBfwg8JCIumLhqN+AL0x6YNAUGq2pY2lTKbFETi70Q4RS692N7HfBy7pyQN2XmddMe\nmDQFBqtqWNpUymxRE4uVtnXAD4EXs11wRcQemXn9NAcmTYHBqhqWNpUyW9TEYqXtbBYOrAQevPLD\nkabKB2DVcs6ohPNETSxY2jJz/4bjkFrxaFilnCuq4XzR1BX9wfiIeBbwRLqjic9k5kemOippOlzC\nUA2XR1XKbFETJX8w/nXA8cCFwEXA8RHx2mkPTJoCg1U1LG0qZbaoiZIzbc8EHpmZtwFExN8B5wJ/\nNsVxSdNgsKqGpU2lzBY1UfLmugnsPnF5dwwyDZPBqhqWNpUyW9REyZm21wJnR8Tm/vKTgFdMbUTS\n9PgArFrOGZVwnqiJxf4iwtuAUzLz/RHxGeAX6CbmKzLz6lYDlFaYR8Mq5VxRDeeLpm6xM22XAG+I\niL2BfwDen5nntBmWNBUuYaiGy6MqZbaoiQWf05aZb8zMx9Ith14PvCsivh4RJ0bEQc1GKK0cg1U1\nLG0qZbaoiSVfiJCZV2Tm6zLzYOAY4Nfo3vpDGhqDVbUsbSphtqiJkvdpWx8RR0XEKcDHgIuBX5/6\nyKSVZ7CqSEQ3TzItbSpitqiJxV6I8HS6M2vPBL4EvB94QWbe3Ghs0kozWFXKeaIaZouaWOyFCK+g\nK2r/LTOvbzQeaZo8a6IazheVcq6oicX+YPyTWw5EasSjYZVwnqiWc0ZTV/IXEaS1wiUMlfKVo6ph\ntqgJS5vGxGBVKUubapgtasLSpjExWFXK0qYaZouasLRpTAxWlbK0qYbZoiYsbRoTH4RVw/miUs4V\nNWFp09h4NKwSzhPVcs5o6ixtGhOXMFTK5VHVMFvUhKVNY2KwqpSlTTXMFjVhadOYGKwqZWlTDbNF\nTVjaNCYGq0pZ2lTDbFETljaNicGqUpY21TBb1ISlTWPig7BqOF9UyrmiJixtGhuPhlXCeaJazhlN\nnaVNY+IShkq5PKoaZouasLRpTAxWlbK0qYbZoiYsbRoTg1WlLG2qYbaoCUubxsRgVSlLm2qYLWrC\n0qYxMVhVytKmGmaLmrC0aUx8EFYN54tKOVfUhKVNY+PRsEo4T1TLOaOps7RpTFzCUCmXR1XDbFET\nljaNicGqUpY21TBb1ISlTWNisKqUpU01zBY1YWnTmBisKmVpUw2zRU1Y2jQmPgirhvNFpZwrasLS\nprHxaFglnCeq5ZzR1FnaNCYuYaiUy6OqYbaoCUubxsRgVSlLm2qYLWrC0qYxMVhVytKmGmaLmrC0\naUwMVpWytKmG2aImLG0aE4NVpSxtqmG2qAlLm8bEB2HVcL6olHNFTVjaNDYeDauE80S1nDOaOkub\nxsQlDJVyeVQ1zBY1YWnTmBisKmVpUw2zRU1Y2jQmBqtKWdpUw2xRE5Y2jYnBqlKWNtUwW9SEpU1j\n4oOwajhfVMq5oiYsbRobj4ZVwnmiWs4ZTZ2lTWPiEoZKuTyqGmaLmphpaYuIIyLi4oi4NCJePs/1\nmyLixog4p/949SzGqTXDYFUpS5tqmC1qYv2s7jgi1gFvAZ4KfAf4ckSclpkXbbfrZzLzqOYD1Fpk\nsKqUpU01zBY1McszbYcCl2XmFZm5FfgA8Kx59vMXQSvFYFUpS5tqmC1qYpal7YHAlROXr+o/NymB\nx0XEeRFxekQ8rNnotBYZrCplaVMNs0VNzGx5lLJAPBvYNzO3RMQzgFOBg+bbMSJOmri4OTM37/AI\ntdb4IKwazheVcq5oURGxCdi0o7czy9L2HWDficv70p1tu0Nm3jSxfUZEvC0i9sjM67e/scw8aVoD\n1Zri0bBKOE9UyzmjBfUnkjbPXY6IE5dzO7NcHj0LODAi9o+IXYDfAk6b3CEiNkZE9NuHAjFfYZMK\nuYShUi6PqobZoiZmdqYtM7dFxHHAx4F1wDsz86KIeGF//cnA0cCLImIbsAU4Zlbj1ZpgsKqUpU01\nzBY1McvlUTLzDOCM7T538sT2W4G3th6X1iyDVaUsbaphtqgJ/yKCxsRgVSlLm2qYLWrC0qYx8UFY\nNZwvKuVcUROWNo2NR8Mq4TxRLeeMps7SpjFxCUOlXB5VDbNFTVjaNCYGq0pZ2lTDbFETljaNicGq\nUpY21TBb1ISlTWNisKqUpU01zBY1YWnTmBisKmVpUw2zRU1Y2jQmPgirhvNFpZwrasLSprHxaFgl\nnCeq5ZzR1FnaNCYuYaiUy6OqYbaoCUubxsRgVSlLm2qYLWrC0qYxMVhVytKmGmaLmrC0aUwMVpWy\ntKmG2aImLG0aEx+EVcP5olLOFTVhadPYeDSsEs4T1XLOaOosbRoTlzBUyuVR1TBb1ISlTWNisKqU\npU01zBY1YWnTmBisKmVpUw2zRU1Y2jQmBqtKWdpUw2xRE5Y2jYnBqlKWNtUwW9SEpU1j4oOwajhf\nVMq5oiYsbRobj4ZVwnmiWs4ZTZ2lTWPiEoZKuTyqGmaLmrC0aUwMVpWytKmG2aImLG0aE4NVpSxt\nqmG2qAlLm8bEYFUpS5tqmC1qwtKmMfFBWDWcLyrlXFETljaNjUfDKuE8US3njKbO0qYxcQlDpVwe\nVQ2zRU1Y2jQmBqtKWdpUw2xRE5Y2jYnBqlKWNtUwW9SEpU1jYrCqlKVNNcwWNWFp05gYrCplaVMN\ns0VNWNo0Jj4Iq4bzRaWcK2rC0qax8WhYJZwnquWc0dRZ2jQmLmGolMujqmG2qAlLm8bEYFUpS5tq\nmC1qwtKmMTFYVcrSphpmi5qwtGlMDFaVsrSphtmiJixtGhODVaUsbaphtqgJS5vGxAdh1XC+qJRz\nRU1Y2jQ2Hg2rhPNEtZwzmjpLm8bEJQyVcnlUNcwWNWFp05j4IKxSljbVcK6oCUubRiOzC9YIj4i1\nJEubipktasXSpjEyWLUUS5uWw2zRVFnaNDY+EKuUc0U1nC+aOkubxsijYS3FOaLlcN5oqixtGhtf\n5aUSLo+qltmiqbO0aWwMVpWwtKmW2aKps7RpbAxWlbC0qZbZoqmztGlsDFaVsLSpltmiqbO0aWwM\nVpWwtKmW2aKps7RpbHwgVinnimo4XzR1ljaNkUfDWopzRMvhvNFUWdo0Ni5hqITLo6pltmjqLG0a\nG4NVJSxtqmW2aOosbRobg1UlLG2qZbZo6ixtGhuDVSUsbapltmjqLG0aGx+IVcq5ohrOF02dpU1j\nFBFsiuAmMNu2AAAOYElEQVQhsx6IVpcInhnBXnjGRMtjtmheE9myQyxtGpu5JYzfB54647Fo9XkJ\n8FhcHlU9s0WLmcuWHWJp09gk8DPA7sCGGY9Fq88G4EDgvljaVMds0WI2AAdGcN8duZH1KzQYaSgS\nOLPf/vIsB6JVaQPwun77XbMciAbHbNFi5rLlt4BHLfdGZnqmLSKOiIiLI+LSiHj5Avu8qb/+vIg4\nuPUYtebsOrF9L4AI3hXB3jMaj2Ysgt+O4Hf6i/eauMozbaphtuguFsiWHeoxMyttEbEOeAtwBPAw\n4NiIeOh2+xwJHJCZBwIvAN7efKBaa35qYntuCeMZdKet10XwshmMSY1FsE8Ez+kvHgYc2m9PLmtZ\n2lTDbFFptizbLM+0HQpclplXZOZW4APAs7bb5yjgPQCZeSawe0RsbDtMrWGHRfARYE9gI/Bg4H9G\nsCGCQyM4FSCCwyJ4Wr/9iAh277f3nLuhCNa1H77gzp99BPeM4D799s9FsEe/fUL/f7pHBF+NIIAn\nA38cwU8B/wn45QjeB9xv4qYtbVous2UNKMyW3ZaRLcs2y+e0PRC4cuLyVcBjCvbZB7hmukPTSBzM\nnaeq/2Hi85+leyL6nhG8GXgC8OMIHkf3fIQvRvBl4HURnAjcCvxZBC+lm5+PBz5IdwZ5K/A1ul/k\nLwDXAv8FeC/dQdMzgX+ie/LyTwMX9rdxI/ADuuc+nNXfx5HAvwA70x30fBm4J10YXE33APF9YAvw\nyP7rtgG/CHya7kzAQ4Gv9193j/4+9qH73foR3e/gF+kKyxOBz9Et+xwAfBO64AJ+CDwEuKz/uicD\n/9Z/T0/qv24D8CDgcmCvflw39GP/Sn/5N4EP9/fxdOBf6R7k9u1/FocA3wG+Czyv3/dW4ETgNf3Y\nnxPB2/v73RjBx4HfBS6O4HzgL4BN/TgfDpzbj+3B/efmPIi78sFSy2W2rK5s+SXg43TZsh/wVbps\n+S5dvjTNlh0p4pE5m4PJiPgN4IjMfH5/+TnAYzLzjyb2+Qjwusz8Qn/5k8CfZubZ291W0v3w5mzO\nzM1T/hY0QBGr/uzJDXQhC90v/a6L7PsTYJd++1rgAf32j7nrUs32fgjce56vu53Fz75fD90R5nbb\nS5ncd6nvaTV5b+Ydz0eRFmW2AGbLIjb3HwCv/0v40asys/r9IGd5pu07dEfTc/ala+SL7bNP/7m7\nycyTVnJwUu9jdM+7BHgZ8FrgOrqjsacC36A7Er0H3dHpI+iO8v6N7ghsK/ApuqM6gEvojgxvpjuT\nE8C36Y42f0R35Pzw/vMX0h2h7UQXTveiOxK8J3Ab8D1g//4+vkF3NL2+v7379193K93R8610R7Jb\n6c5U70d3xHs53e/Y+v7zc6E+F7I39bd1G93v54P76y+n+33cmS6c79N/zfp+7N/vx76tH9tD++su\n7se5c/9z3KP/ed2nv79v0p0BuA04m+45IRvoXpX3qP62b+jH/G3g5+geJD4FnEB3ZuBtwF/193MK\n8Nss7Aew4Evwd17k66QdZbZ0t7WN7nG9Jlv268c5oGzZ1H8AnPQGiFct8rULmuWZtvV0p1KfQneK\n8kvAsZl50cQ+RwLHZeaREXEY8MbMPGye28rlNFaNzxJHw9+n+yX9M7o5+Q26oNsAXJfJtumPUCsp\ngp3oAn03ugB9EF2Yvp7ugfFtwB/2u3+TOx84AD6cydHtRqshM1vGZQezZSPENYM605aZ2yLiOLp1\n5nXAOzPzooh4YX/9yZl5ekQcGRGXAbfQrTtLK+lCuqPP/0p35LQtk9u22+eW5qPSisjkdrqlnuv6\nj8v6q54WwT3pzg5cSXeW41vcNVh3QVo+s2UN24FsuYUdyJaZnWlbSZ5pU6n+aHgbdx6w7AnsnsnX\nZzcqzVoEhwO/QbcMMufjmXcsX0mLMls0n3my5SpgE8RlgzrTJs3QDXTPE/lkJtfgq5FHL5MvRHA9\n3avaHkl39t8zbapltugu5smWm9mBbPFvj2qMfgDckMmvznogWj0yuQh4MXcuWflCBNUyW3Q322XL\nj9mB0uaZNo3RDXSvFJK2dytdsN4bS5vqmS1ayFy2/ITF3zZlUZ5p0xhdTvdqZWl7VwP/3m+7PKpa\nZosWcjXdmyBPvgdeNUubxujyTJ4960Fo9cnk2om3+fBMm2qZLZpXny3/mR1cHrW0aYx+MusBaBAs\nbapltmgpLo9KlQxWlXB5VLXMFi3F5VGp0tZZD0CDYGlTLbNFS3F5VKrk0bBKWNpUy2zRUlwelSp5\nNKwSPqdNtcwWLcUzbVIlg1UlLG2qZbZoKZ5pkyq5hKESLo+qltmipfhCBKmSwaoSnmlTLbNFS3F5\nVKpksGop22Y9AA2S2aKluDwqVfJ5J1qKc0TL4bzRUnZoedQ/GK+x+SFw/qwHoVXvDcB9Zz0IDYrZ\nohLXAvde7hdHZq7gWGYjIjIzY9bjkCRJWspye4vLo5IkSQNgaZMkSRoAS5skSdIAWNokSZIGwNIm\nSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2\nSZKkAbC0SZIkDYClTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0\nSZIkDYClTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYCl\nTZIkaQAsbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAs\nbZIkSQNgaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAbC0SZIkDYClTZIkaQAsbZIkSQMw\nk9IWEXtExCci4pKI+NeI2H2B/a6IiPMj4pyI+FLrcWrtiYhNsx6DhsP5olLOFbUwqzNtrwA+kZkH\nAZ/qL88ngU2ZeXBmHtpsdFrLNs16ABqUTbMegAZj06wHoLVvVqXtKOA9/fZ7gF9dZN+Y/nAkSZJW\nt1mVto2ZeU2/fQ2wcYH9EvhkRJwVEc9vMzRJkqTVJzJzOjcc8Qlgz3muehXwnsy878S+12fmHvPc\nxl6ZeXVEPAD4BPBHmfm5efabzjchSZI0BZlZvZK4fhoDAcjMpy10XURcExF7Zub3ImIv4D8WuI2r\n+3+vjYh/Ag4F7lbalvONS5IkDcmslkdPA57bbz8XOHX7HSLinhGxW799L+DpwAXNRihJkrSKTG15\ndNE7jdgD+CCwH3AF8JuZeUNE7A28IzOfGREPBv6x/5L1wPsy87XNBytJkrQKzKS0SZIkqc5g/iJC\nRBwRERdHxKUR8fIF9nlTf/15EXFw6zFq9VhqvkTEpoi4sX/j5nMi4tWzGKdmLyLe1T/PdsGnX5gt\ngqXnirmiORGxb0R8OiIujIivRsTxC+xXlS2DKG0RsQ54C3AE8DDg2Ih46Hb7HAkckJkHAi8A3t58\noFoVSuZL7zP9GzcfnJmvaTpIrSbvppsr8zJbNGHRudIzVwSwFTghMx8OHAa8eCV6yyBKG92rRi/L\nzCsycyvwAeBZ2+1zxxv2ZuaZwO4RsdD7v2ltK5kv4Bs3C+jfRugHi+xitggomitgrgjIzO9l5rn9\n9s3ARcDe2+1WnS1DKW0PBK6cuHxV/7ml9tlnyuPS6lQyXxJ4XH9K+vSIeFiz0WlozBaVMld0NxGx\nP3AwcOZ2V1Vny9Tep22Flb5aYvsjHF9lMU4l/+9nA/tm5paIeAbd284cNN1hacDMFpUwV3QXEbEB\n+BDwkv6M29122e7yotkylDNt3wH2nbi8L10jXWyfffrPaXyWnC+ZeVNmbum3zwB27t+KRtqe2aIi\n5oomRcTOwIeBv8/Mu70fLcvIlqGUtrOAAyNi/4jYBfgtujfonXQa8DsAEXEYcMPE3zfVuCw5XyJi\nY0REv30o3dvfXN9+qBoAs0VFzBXN6efBO4GvZeYbF9itOlsGsTyamdsi4jjg48A64J2ZeVFEvLC/\n/uTMPD0ijoyIy4BbgOfNcMiaoZL5AhwNvCgitgFbgGNmNmDNVES8H3gScP+IuBI4EdgZzBbd1VJz\nBXNFdzoceA5wfkSc03/ulXR/VGDZ2eKb60qSJA3AUJZHJUmSRs3SJkmSNACWNkmSpAGwtEmSJA2A\npU2SJGkALG2SJEkDYGmTtGZExP0i4pz+4+qIuKrfviki3jKl+zwuIn53keuPioj/Po37ljQuvk+b\npDUpIk4EbsrMv5nifQTd35v8hczctsg+5/T7bJ3WWCStfZ5pk7SWzf1JoU0R8ZF++6SIeE9EfDYi\nroiIX4+Iv46I8yPijIhY3+93SERsjoizIuJjEbHnPLd/OHDxXGGLiOMj4sKIOK9/93yyOzL+IvD0\nFt+wpLXL0iZpjH4a+EXgKODvgU9k5iOAW4Fn9n/o+c3Ab2Tmo4F3A385z+08nu5v3c55OfDIzPx5\n4IUTn/8S8MQV/y4kjcog/vaoJK2gBM7IzNsi4qvATpn58f66C4D9gYOAhwOf7P/+9zrgu/Pc1n7A\n5ycunw+cEhGnAqdOfP67wBEr+U1IGh9Lm6Qx+glAZt4eEZPPM7udLhcDuDAzH1dwWzGx/Uy6M2q/\nArwqIn4uM2+nW9XwCcSSdojLo5LGJpbeha8DD4iIwwAiYueIeNg8+30L2LPfJ4D9MnMz8ArgPsCG\nfr+9+n0ladksbZLWspz4d75tuPsZsOxf5Xk08PqIOJfu1Z+Pnef2Pw88ut9eD7w3Is6ne0Xp32bm\nD/vrDgU+uyPfiCT5lh+StEwTb/nxmMz8yQL77NTv8+iF3hZEkkp4pk2Slql/O493AM9eZLdfBj5k\nYZO0ozzTJkmSNACeaZMkSRoAS5skSdIAWNokSZIGwNImSZI0AJY2SZKkAfj/y/Mu4CELoa8AAAAA\nSUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "time = np.arange(0, period_square, period_square/1000)\n", "Eapp = np.arange(0, period_square, period_square/1000)\n", "Es = E_start_new\n", "Er = E_rev_new\n", "vs = vscan_new\n", "i = -1\n", "for tt in time:\n", " i=i+1\n", " Eapp[i] = EvaluateRamp(Es, Er, vs, tt)\n", " EAC = 0.0\n", " for j in range(0,N_ac):\n", " EAC = EAC + EvaluateEAC(amp[j], freq[j], tt)\n", " Eapp[i] = Eapp[i] + EAC\n", "\n", "plt.figure(figsize=(10,8), dpi=100)\n", "plt.plot(time,Eapp)\n", "plt.xlabel('Time (s)')\n", "plt.ylabel('Voltage (V)')\n", "plt.savefig('SquareWave.png')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Output MECSim inputs\n", "------" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Change the following variables in Master.inp:\n", "1.0000005 ! E_start (V)\n", "0.9999995 ! E_rev (V)\n", "1e-06 ! scan speed (V/s)\n", "\n", "Copy the following lines into Master.inp\n", "100 ! number of species (need n lines below)\n", "1273.23954474 , 1.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "424.413181578 , 3.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "254.647908947 , 5.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "181.891363534 , 7.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "141.471060526 , 9.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "115.749049521 , 11.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "97.9415034412 , 13.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "84.8826363157 , 15.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "74.896443808 , 17.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "67.0126076176 , 19.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "60.6304545112 , 21.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "55.3582410754 , 23.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "50.9295817894 , 25.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "47.1570201754 , 27.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "43.9048118874 , 29.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "41.0722433786 , 31.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "38.5830165071 , 33.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "36.3782727067 , 35.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "34.4118795874 , 37.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "32.6471678137 , 39.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "31.0546230423 , 41.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "29.6102219706 , 43.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "28.2942121052 , 45.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "27.0902030795 , 47.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "25.9844805048 , 49.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "24.9654812693 , 51.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "24.0233876365 , 53.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "23.1498099043 , 55.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "22.3375358725 , 57.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "21.5803312667 , 59.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "20.8727794219 , 61.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "20.2101515037 , 63.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "19.5883006882 , 65.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "19.0035752946 , 67.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "18.4527470251 , 69.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "17.9329513343 , 71.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "17.4416375991 , 73.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "16.9765272631 , 75.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "16.5355785031 , 77.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "16.1169562625 , 79.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "15.7190067251 , 81.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "15.3402354787 , 83.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "14.9792887616 , 85.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "14.6349372958 , 87.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "14.3060623004 , 89.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "13.9916433487 , 91.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "13.6907477929 , 93.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "13.4025215235 , 95.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "13.1261808736 , 97.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "12.8610055024 , 99.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "12.6063321261 , 101.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "12.361548978 , 103.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "12.1260909022 , 105.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "11.8994349975 , 107.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "11.6810967407 , 109.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "11.4706265291 , 111.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "11.2676065906 , 113.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "11.0716482151 , 115.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.8823892712 , 117.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.6994919726 , 119.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.5226408656 , 121.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.3515410141 , 123.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.1859163579 , 125.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "10.0255082263 , 127.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.8700739902 , 129.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.71938583767 , 131.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.57322965966 , 133.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.43140403508 , 135.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.29371930464 , 137.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.15999672471 , 139.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "9.03006769316 , 141.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.90377304011 , 143.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.78096237748 , 145.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.6614935016 , 147.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.54523184386 , 149.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.43204996513 , 151.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.32182708977 , 153.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.21444867571 , 155.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.10980601742 , 157.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "8.00779587884 , 159.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.90832015363 , 161.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.81128555052 , 163.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.71660330143 , 165.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.62418889063 , 167.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.53396180317 , 169.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.44584529085 , 171.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.35976615454 , 173.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.27565454134 , 175.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.19344375557 , 177.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.11307008232 , 179.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "7.03447262285 , 181.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.95759314063 , 183.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.88237591749 , 185.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.8087676189 , 187.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.73671716791 , 189.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.66617562689 , 191.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.59709608671 , 193.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.52943356274 , 195.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.46314489713 , 197.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n", "6.39818866701 , 199.0 ! AC sin wave: amp (mV), freq(Hz) (REPEAT)\n" ] } ], "source": [ "print \"Change the following variables in Master.inp:\"\n", "print E_start_new, \" ! E_start (V)\"\n", "print E_rev_new, \" ! E_rev (V)\"\n", "print vscan_new, \" ! scan speed (V/s)\"\n", "print\n", "PrintSquareWave(amp,freq,N_ac)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }