Waveguide Slot Antenna Calculation

wg-slot-antenna

{
 "metadata": {
  "name": "WG-Slot-Antenna"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "<h4>Design of a Waveguide Slot Antenna</h4>\n",
      "O.Nell, K. Solbach, J. Dreier: Rundumstrahlende Hohlleiter-Schlitzantenne f\u00fcr Horizontalpolarisation<br>\n",
      "UKW-Berichte 1/91 und 2/91\n"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "f0 = 5760 # [MHz] Frequency \n",
      "N = 10 # Number of slots\n",
      "\n",
      "Z_F0=120.0*pi\n",
      "c0=2.997831e8\n",
      "lambda_0= C0/f0/1e6\n",
      "\n",
      "print 'Freespace wavelength in: %6.2f mm' % (lambda_0*1e3)\n",
      "print 'Cutoff value for a: %6.2f mm' % (lambda_0*1e3/2)\n",
      "print 'Optimum width for a: %6.2f mm' % (lambda_0*sqrt(3.)*1e3/2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Freespace wavelength in:  52.05 mm\n",
        "Cutoff value for a:  26.02 mm\n",
        "Optimum width for a:  45.07 mm\n"
       ]
      }
     ],
     "prompt_number": 49
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "a = 43.8 # [mm] waveguide width\n",
      "b = 23.9 # [mm] waveguide height\n",
      "\n",
      "lambda_c=2.0*a/1e3\n",
      "lambda_h= lambda_0 / sqrt(1.0 - (lambda_0/lambda_c)**2)\n",
      "\n",
      "print 'lambda_h: %3.2f mm' % lambda_h\n",
      "\n",
      "g= 10*log10(N * lambda_h/lambda_0)\n",
      "print 'Approximate antenna gain: %3.1f dBi' % g\n",
      "\n",
      "dtheta=50.7 * lambda_0 / N / lambda_h * 2.0\n",
      "print 'Beamwidth: %3.1f deg' % dtheta\n",
      "\n",
      "# calculate offset from center\n",
      "x = 1e-3 * a/pi*arcsin(sqrt(lambda_0*b/(N*3.5*lambda_h*a*cos(pi*lambda_0/2.0/lambda_h)**2)))\n",
      "print 'slot offset from center: %6.2f mm' % (x*1e3)\n",
      "\n",
      "# slot width\n",
      "sb = lambda_h / 30\n",
      "print 'slot width: %6.2f mm' % (sb*1e3)\n",
      "\n",
      "# compute slot length\n",
      "# use figure 7 and interpolate\n",
      "x_23cm = f0/1270 * x\n",
      "sx=[ 0,    8.6,   10,    20,   29 ] # x in mm\n",
      "sl=[ .475, .491,  .494,  .519,  .55] # L/lambda_0\n",
      "\n",
      "L = interp(x_23cm*1e3,sx,sl)\n",
      "print 'slot length: %6.2f mm' % (L*1e3*lambda_0)\n",
      "\n",
      "# calculate slots\n",
      "print ''\n",
      "print '   Slot No    Distance'\n",
      "print '----------------------------------------'\n",
      "for i in range(N-1):\n",
      "    print '      %d       %6.2f mm' % ( i+1 ,1e3*(lambda_h/4+lambda_h/2*i) )\n",
      "\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "lambda_h: 0.07 mm\n",
        "Approximate antenna gain: 11.0 dBi\n",
        "Beamwidth: 8.1 deg\n",
        "slot offset from center:   5.08 mm\n",
        "slot width:   2.21 mm\n",
        "slot length:  27.50 mm\n",
        "\n",
        "   Slot No    Distance\n",
        "----------------------------------------\n",
        "      1        16.58 mm\n",
        "      2        49.73 mm\n",
        "      3        82.89 mm\n",
        "      4       116.05 mm\n",
        "      5       149.20 mm\n",
        "      6       182.36 mm\n",
        "      7       215.51 mm\n",
        "      8       248.67 mm\n",
        "      9       281.83 mm\n"
       ]
      }
     ],
     "prompt_number": 48
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}

I am interested in waveguide slot antenna and want to make simulate it in EM simulator or something. Are these codes for execution on EM simulators?