VolumeRender.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "5e2689d5",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning: sysroot \"/Applications/Xcode_12.4.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX.sdk\" not found (ignoring for now)."
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Vapor 3.8.2\n",
      "Python 3.9.16 (/Users/zongrunli/opt/anaconda3/envs/vapor)\n",
      "OpenGL 4.1 INTEL-16.5.14\n"
     ]
    }
   ],
   "source": [
    "from vapor import session, renderer, dataset, camera, utils, transferfunction\n",
    "from vapor.animation import Animation\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "b8e28ec4",
   "metadata": {},
   "outputs": [],
   "source": [
    "filenames = [\"/Users/zongrunli/Desktop/WRF_Fire_Visualization/wrfout_d01_2022-04-25_00:00:00\"]\n",
    "ses = session.Session()\n",
    "data = ses.OpenDataset(dataset.WRF, filenames)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "708d2715",
   "metadata": {},
   "source": [
    "### TimeStep?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "042f073e",
   "metadata": {},
   "outputs": [],
   "source": [
    "ses.SetTimestep(60)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "66426feb",
   "metadata": {},
   "source": [
    "## Smoke 3D Render"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "86227371",
   "metadata": {},
   "source": [
    "## Read Variable"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "fe998579",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "smoke_tracer_name = \"tr17_2\"\n",
    "smoke_tracer_idx = data.GetDataVarNames(3).index(smoke_tracer_name)\n",
    "# render\n",
    "ren = data.NewRenderer(renderer.VolumeRenderer)\n",
    "ren.SetVariableName(data.GetDataVarNames(3)[smoke_tracer_idx])\n",
    "render_min = data.GetDataRange(smoke_tracer_name)[0]\n",
    "render_max = data.GetDataRange(smoke_tracer_name)[1]\n",
    "data.IsTimeVarying(smoke_tracer_name)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5e8f5577",
   "metadata": {},
   "source": [
    "## Transfer Function"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "4f73c9fa",
   "metadata": {},
   "outputs": [],
   "source": [
    "smoke_max = 100\n",
    "smoke_min = 35\n",
    "opacities = [(render_min, 0), (smoke_min, 0), (smoke_max, 1), (render_max, 1)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "cd5bfebe",
   "metadata": {},
   "outputs": [],
   "source": [
    "ses.GetCamera().ViewAll()\n",
    "tf = ren.GetPrimaryTransferFunction()\n",
    "tf.SetOpacityControlPoints(opacities)\n",
    "tf.SetMaxMapValue(100)\n",
    "tf.SetMinMapValue(35)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "985357bb",
   "metadata": {},
   "source": [
    "## Colorbar"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4fa97feb",
   "metadata": {},
   "outputs": [],
   "source": [
    "colorbar = ren.GetColorbarAnnotation()\n",
    "colorbar.SetEnabled(True)\n",
    "# Customized colorbar annotation\n",
    "colorbar.SetSize((0.1, 0.25))\n",
    "colorbar.SetTitle(\"PM2.5\")\n",
    "colorbar.SetCornerPosition((0.03, 0.03))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "16f595e2",
   "metadata": {},
   "source": [
    "## Scene"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "0bf3cbc0",
   "metadata": {},
   "outputs": [],
   "source": [
    "scene = ses.GetSceneAnnotations()\n",
    "scene.SetTimeType(scene.TimeAnnotationType.Formatted)\n",
    "scene.SetTimeLLX(0.6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "e7b54a1a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAHgCAIAAAD2dYQOAAAfqElEQVR4nO3df3QU1f3/8TvZJESCWZpVCVKw/Ba1gEK/CqZNGhNiRC0IBjkCHipWDz0qaFH8RT+oWAtKrQWDfulXRbRRUQiHNCaRWho1YimhCII21kIjLMJuNj/AELLczx/b7nfdX0wyS5Ll/Xyc/DF79869dyYkL+6d2YlSAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAccno6gHEAa11Vw8BADrCMPglj1jQWpOFAOIOv7tMSujqAQAA0JUSu3oA8cEwDK31osTvhn1XBywx628vN/tfhm4EVtAmFqnD1vT/Zy96I4HVwo4t6N2w9cM2GNh10F5hjz1wl9B9jW+3HMqIerZV5HaMb7+rInQUeiyhFVTUOqc8hEiMcOcn0rsq3Lc+9JRGOl1RGjczSPPf38ABh4ryPQotNEI6Cv3WR28h+lGEthDpBy2ociRhBxlpwNFPUccWNx9tq/P94urQ3rIQhGYZhrHI1i/8WxG2A1+GbgS+NPMPPWzNKF1HqhZlbGaaClsh+tjCdhRUGGUwUbqOMvLohxylIzPfkeh1OnxNxuShRaoc9pRGOedRGo8iyrfV5L/G6HWiFJr85ppp1sy/yeijau/PSJQK7f3JNYmrg+axNAoAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEAADRCEIAgGgEIQBANIIQACAaQQgAEI0gBACIRhACAEQjCAEAohGEp2YYRlcPAQBwuiS2dwet9ekYRzdHFgLAmYoZIQBANIIQACAaQQgAEI0gBACIZjUIe/funZycnJKSkpKScv755992220NDQ3+t0aOHBlYubW1tXfv3pmZmb6X77///uWXX/6d73xnyJAhL7/8cmjjZ599dkqAzz//3OJoAQAIEoMZ4caNG1taWlpaWj766KO9e/fOmzfP/5bL5dqzZ4//ZWVlZc+ePX3bHo9n4sSJP//5z91u95o1a+64447du3cHNnvy5MmjR49++eWXLf81bNgw66MFACBQLJdGBwwYsGDBgj/+8Y/+kpycnDfeeMP/8s033/zRj37k225tbX366adnzZplGMb48eNHjRq1devWwNaam5u11r169YrhCAEACBLja4Rerzc5Odn/8tprr/UH4YkTJ8rLyydMmOB7ed55582ZM8e3rbXet2/fgAEDAptqbGw0DCM1NTW2IwQAIFDMglBr/cUXX/z617+eOnWqv3DMmDHHjx/ftWuXUurdd9/9wQ9+kJaWFrrvU0891a9fv5ycnMDCpqYmm812/fXXOxyOCy64YOnSpbEaKgAAfu1+skyoyZMn22w2wzD69Olzww03LF68OPDdwsLCN95445JLLlm3bl1hYWHo7qtWrfrtb3+7ZcuWhIRvpXJqauqcOXNmzJgxfvz47du35+XlDRw48MYbb7Q+YAAA/GIwI1y/fn1zc3NTU1Ntbe3SpUvPOuuswHd9QdjW1vbOO+9cf/31QfsuX7586dKlf/nLXwYPHhz01oABA4qKiq688krDMMaMGTNjxozAq48AAMTEaf8c4ejRo5VSL7zwwtixY4PWRV9//fXf/e53W7ZsGTRoUOiObrf7008/9b88fvy4/45TAABipTM+UF9YWLhkyZKgVc2mpqa777577dq1/fv3D7vXp59+esUVV/huJf3oo4+Ki4snT57cCaMFAIjSGUE4bdo0l8sVtC76zjvvHDp06KqrrvJ/Xn7+/PlKqdGjR69YsUIplZmZuXz58hkzZvTu3Xv27NnPPvtsbm5uJ4wWACBKu/+6kOQ/w7TI1q+rBwIApjzq/UoppbXmr8idEs8aBQCIRhACAEQjCAEAohGEAADRCEIAgGhxGYTV1dXDhg1btWqV76XWet68eXl5eVlZWdXV1WFL/NpVGQBwxovBs0Y7WWtr69q1awMfW1pRUVFfX19ZWel0OgsKCmpqakJLOlYZAHDGi78ZYXJy8sqVK1NSUvwlVVVV+fn5SqmMjIzExESXyxVa0rHKAIAzntUZYWnShTEZh1Jq4om9vo1IH/+M9Fl+t9vtf4qp3W53u92hJQ6HowOVAQBnPMtLo6fhkQXtfXiNw+HweDy+bY/H43A4Qks6VhkAcMazujSacFZCrL46PIbs7OzS0lKlVF1dnWEY6enpoSUdqwwAOONZnRHaLARYx3zyySd33nnnvn37kpKSiouLS0pKcnJyysrKcnNzvV5vUVGRUiq0xOl0ZmZm1tbWmqkMAJDD6kO3N18wMlZDuWrfzlg1FXM8dBtAfOGh2+ZZnREmJMfffacAAPhZXholCAEA8czyjDCJIAQAxDHLM0KCEAAQz5gRAgBEszwjTCQIAQBxzPKMkCAEAMQzlkYBAKIxIwQAiGb9GqEtJuMAAKBLMCMEAIhmOQhtBCEAII5ZnxGyNAoAiGMsjQIARCMIAQCiWQ1Cw8bSKAAgjjEjBACIZv2uUWaEAIA4ZnlplBkhACCeMSMEAIjGB+oBAKJZXxrt7BlhS0vL7NmzvV5vU1PTnDlzpkyZ8sorr/z+97/3vfvXv/710KFDvXr18r3UWs+fP3/37t2tra1PPvnkuHHjQks6efwAgG7F+scnOntGuGHDBofDsWLFiqNHj1588cVTpkyZOXPmzJkzlVKbN2/esGGDPwWVUhUVFfX19ZWVlU6ns6CgoKamJrSkk8cPAOhW4u8a4bBhw1566aXGxsZdu3YNHjzYX97W1vbII4+UlZUFVq6qqsrPz1dKZWRkJCYmulyu0BKHw9HJhwAA6D6644zQMIyw5VprpdSIESN69Ogxffp0p9N5//33+98tLi7Oz8+32+2Bu7jd7rS0NN+23W53u92hJQQhAEjWHa8R+gIvkpUrV1555ZX33XdfW1vbJZdcMnHixNTUVKXU888//+KLLwZVdjgcHo/Ht+3xeBwOR2hJzMcPAIgjVudzhs0Wqy+TPXq9Xq/Xq5QyDKOlpcVX2NzcfODAgSFDhgRVzs7OLi0tVUrV1dUZhpGenh5aYvEMAADiWndcGo3ujjvumD17dmFhYVNT0wMPPOCbDu7bt69fv37+Ok6nMzMzs7a2Nicnp6ysLDc31+v1FhUVKaVCSwAAkoW/GhdF0Lpl4zP3xGooafOWx6qpmPNdtlxk63fKmgDQHTzq/UoppbWOdNcF/PjrEwAA0eJvaRQAgBhiRggAEM1qEKoEghAAEMeYEQIARLM8I+QaIQAgnjEjBACIxjVCAIBozAgBAKIxIwQAiMbNMgAA0ZgRAgBEsz4jJAgBAHGMGSEAQDRmhAAA0azPCLlZBgAQx6wGoWZpFAAQz7hGCAAQjSAEAIjG0igAQDRmhAAA0aze86kTEmL1ZbLHV155Jfu/UlNTm5ubfeVNTU39+/fftGnTt4an9bx58/Ly8rKysqqrq8OWAAAki7+l0ZkzZ86cOVMptXnz5g0bNvTq1ctX/tBDDw0fPjyockVFRX19fWVlpdPpLCgoqKmpCS3p5PEDALqV7hiEhmGE70tr/3ZbW9sjjzxSVlbme7l169bW1taxY8cG7VJVVZWfn6+UysjISExMdLlcoSUOhyPmhwAAiBeWPw6fYIvZ13/pCAK7LS4uzs/Pt9vtSqkTJ04sXLhwyZIloaNzu91paWm+bbvd7na7Q0usngEAQDyzPCM0uuZmmeeff/7FF1/0bS9btmzWrFlhJ3YOh8Pj8fi2PR6Pw+EILemM4QIAuqv4u1lGKdXc3HzgwIEhQ4b4X5aUlEyaNOntt99esmTJzp07/TWzs7NLS0uVUnV1dYZhpKenh5ZYPAMAgLhmdUZ4sitmhPv27evXr5//5RNPPOHbWLhwYWZm5siRI51OZ2ZmZm1tbU5OTllZWW5urtfrLSoqUkqFlgAAJAt/W0oUQdfq3DurYjWU9JE/jFVTMee7f2eRrd8pawJAd/Co9yullNY60u2H8IvXa4QAAMSE9Y9P8GeYAABxLC6vEQIAECssjQIARLM8I2RpFAAQz5gRAgBEs36NkBkhACCOcbMMAEA0ghAAIJrlILT+9ysAAOg6XCMEAIhm+a5RTRACAOKY9aVRrhECAOIY1wgBAKJZDkLNH/gAAMQxZoQAANGszwgJQgBAHCMIAQCiWV8a5RohACCOMSMEAIhmNcZO6oRYfZnvdM2aNRMmTLjiiiueffZZX0l1dfWwYcNWrVoVVFNrPW/evLy8vKysrOrq6rAlAADJ4u/jE06nc/369eXl5W1tbTU1NUqp1tbWtWvXFhYWhlauqKior6+vrKx0Op0FBQU1NTWhJZ08fgBAt2J9RmjE6stkj5WVlTab7ZZbbsnKyvryyy+VUsnJyStXrkxJSQmtXFVVlZ+fr5TKyMhITEx0uVyhJRbPAAAgrlmdEXpPwzVCwwgfilprpdShQ4eSk5PXrFnT1NT0/e9/f+rUqTZbxMe8ud3utLQ037bdbne73aElDocj1kcAAIgb3XFp1Bd4kdjt9qFDhyqlzj777HPOOcflcp133nmRKjscDo/H49v2eDwOhyO0JFbDBgDEo/hbGs3MzKysrDx58mRDQ0NDQ8O5554bpXJ2dnZpaalSqq6uzjCM9PT00BKLZwAAENe644wwuhEjRkyePDkvL6+1tfWZZ54xDOOTTz6588479+3bl5SUVFxcXFJS8s0332RmZtbW1ubk5JSVleXm5nq93qKiIqVUaAkAQLJ2x1jQumXF31tjNZQJo5Jj1VTM+S5bLrL16+qBAIApj3q/UkpprSPddQE/yzfLnOQUAwDiWPwtjQIAEEOWg5AZIQAgnln/HGFMhgEAQNdgRggAEM36jJAgBADEMeszwpgMAwCArkEQAgBEY2kUACAaM0IAgGjWnywTk2EAANA1mBECAERjRggAEI0ZIQBANGaEAADRmBECAESzHoQxfOo2H0kEAHQ2y0uj3pgMAwCArsGMEAAgGjfLAABEszwj5C/zAgDiWYLF/b0nY/bVLk1NTf3799+0aZNSqqWlZfr06YWFhQUFBW+99VZgNa31vHnz8vLysrKyqqurw5YAACTrVtcI2+Ghhx4aPny4b3vDhg0Oh2PFihVHjx69+OKLp0yZ4q9WUVFRX19fWVnpdDoLCgpqampCS7pk/ACAbsLyjNCrY/Xlb9OIwF9h69atra2tY8eO9b0cNmxYbW1tY2Pj3//+98GDBwcOr6qqKj8/XymVkZGRmJjocrlCSyyeAQBAXOuOM0Kto7V54sSJhQsXrlu3btmyZb6SESNG9OjRY/r06U6n8/777w+s7Ha709LSfNt2u93tdoeWOByOmB8CACBeWP8cYWcvjS5btmzWrFmB6bVy5corr7zyvvvua2tru+SSSyZOnJiamup7y+FweDwe37bH43E4HKElnTp6AEA3Y3Vp9KRXx+rLZI/Nzc0lJSWTJk16++23lyxZsnPnTq/X6/V6lVKGYbS0tARWzs7OLi0tVUrV1dUZhpGenh5aYvEMAADimvXPEXb2jPCJJ57wbSxcuDAzM3PkyJEXXHDB7NmzCwsLm5qaHnjggdTUVKfTmZmZWVtbm5OTU1ZWlpub6/V6i4qKlFKhJQAAydr9MJegC3j3rjwaq6E8/fPUWDUVc75bdRbZ+nX1QADAlEe9XymltNaBdxoirPibEQIAEEM8WQYAIJr1IORhowCAOBZ/H58AACCGmBECAESzfrMMQQgAiGPcLAMAEI2lUQCAaNZvliEIAQBxjBkhAEA0ghAAIBpBCAAQzfo1Qm9MxgEAQJdgRggAEI0gBACIRhACAETjGiEAQDRmhAAA0awHITNCQK6MURePv+f2gVnjUvuc2/bNN0c+++KT10s+Xvmi90RbYDV7//N/uPDOofnZvfr2aW1qPvC3nR8+83+/qNzSsdZMVovE5GACXfDDy2dvftNISHh92u2fvlUakzOTmNLjkeYvIu2+4dZ7al5+w0xHgQIP7c6mxm3btoWt1r9//wcffPDqq6/u27dvU1PTtm3bfvOb31RUVJjpwuS+3byLIEZ7d9D6W0/Zzrt5ewd6Davy1cti1VTMGYahlFpk69fVAwG6kVEzpkxa/XRCYuI39Q2uz7/o6UhPH/I9pdT+Dz5+OX96W8txX7WMkRfNfu+tFPvZx1z1X+/+LPWc9HMvGqaUKr/vsQ+XP9/e1kxWi8TkYAIl90qdW1P5nYEDlFImg9DMIM/ue94v/r39ZFubc+enoS28t3j556XvnrKjKId27Jy0iy66SCn1i1/84umnn/7/Yxs1asuWLXa73eVy7dq169xzzw1bLfxxmdu3m3cRymoQXnXTX9vbQiSbi39gplpLS8utt9565MiRxsbGe+6558Ybb/z4449/8pOfDB8+XCk1ZsyYwLOgtZ4/f/7u3btbW1uffPLJcePGhZaY6dQwDK31Xba+HTs04MxzzuDvPfjJlsQeyZufeq500a/bjrcqpYb9OPNnG9ck9zyr5P7HNj/1nFIqwWZ78JMt5w0f/P7za96a97C39YRS6sK8rNs2vJyYnLTs/+TX1ewy35rJapGYHEyQm55/avycmxsOOO3nZ/y/wtt2vLUpJmemz4ihD+36i3tf3f8MMvWrL7rQQ3vWezA/P7+8vPzkyZNjx46tqalRStlstt27dw8fPnzVqlV33313a2urUmrChAklJSXJycn+amGZ3LebdxEb+tuyp1bH6kubs379+qVLl2qtDx8+PHDgQK11WVnZXXfdFbbyO++8M2vWLK31wYMHR48eHbbEjNN7ToE4tGTJEq31hx9+GFT++OOPa63ff/9938spU6Zorffu3Wuz2QKrPfbYY1rr1157rV2tmawWicnBBLrmmmu01hs3bly3bp3WeurUqdG7MD9I3//Ld+zYccoGzQh7aP7fYP5D68AZiN5F6L7dvIuwEtq7QxBv7JjscdKkSQsWLPB6vdXV1RdeeKFSyuVyVVZWXnvttWPHjt24cWNg5aqqqvz8fKVURkZGYmKiy+UKLbF4BgCZysvL77nnnl/+8pdB5b7/jPft+5/lk+uuu04p9dprrwX9jK9du1YpNXHixMTERPOtmawWicnB+KWnp69evdrtdt9+++3RWw5kcpC9e/dWSnk8HvMtRxHp0Hz8h2b+DOzYsUNr/fDDD5+yi6B9u1UXJlm+WaYt9neN+i7IhQqcmeXk5Bw6dKioqEgplZeXN27cuEGDBh0+fHjs2LHXXHON/yy43e60tDTftt1ud7vdoSUOh8PMqJgXAuYNGjQo8Edm8eLFixcvDq2WlpZ24sSJ9rZmsVoHBnPgwAHfxptvvnnK9qMLGmRWVlYMf7eEPbTGxsa0tLShQ4fu2bNnzJgxSqmPP/44qM5nn30WWC1s4yb37eZdxIbJtcRO0NDQMGjQoObm5sDCyy677MCBA/6XDz/88CuvvOLbHjNmjMvlCi0x09dpPKHAGcQwjA8++EBrPX/+fF/JsWPHtNa+GxmC7Nq1S2t97bXXmm/NSrX2DmbatGla6z/84Q++l+aXRk0Ocu7cuVrr4uLihQsX7tix49ixY0ePHt2xY8cDDzxw1llntbd9k4dm/gyETtfisQuTrC6Ndr5nnnlm+fLlSqmkpKSWlhbDMLZu3VpVVaWUcjqdDQ0NGRkZ/srZ2dmlpaVKqbq6OsMw0tPTQ0u66DiAM9CCBQvGjx//+eef+1ZrkpKSfL/TGxoaQiv7Cn0rhGZas1itXYPp27fvc889d+DAgblz50Zp07zQQfr6mjZt2q9+9auePXtu3br18OHDo0aNeuKJJz744AO73W6+cZOH1q4zMHr0aMMwHn/88bjuwiSrS6Od75Zbbpk1a9amTZuOHj26aNGinj179unT56c//anWuq2t7YUXXjAMw+l0ZmZm1tbW5uTklJWV5ebmer1e37+/0BIAMbFgwYInn3zS7XZPmjSppaVFKZWSkuJ7y3dfX5Djx48H1jlla1aqtXcwq1evTk9PLygoqK+vj9KmSWEHmZCQcOjQoYMHD86dO7e6utpXmJeXt3bt2ksvvXTFihUzZ8402b7JQ7Py7TgzuoiZdq1enjFO7zkF4pzNZisqKtJaHzx4cOTIkf7ypKQk309Qv35hPoPrWyq8+eabTbbWsWodGMxtt92mtQ76j3LHlkbbO0il1A033KC19nq955xzjsleTB5ax74dZ1IXkcTf0iiAbiUtLa20tPSOO+7YvXv3FVdcsXPnTv9bJ06cOHbsmIqwVOUrDJpyRWnNfDWPxxP4H9kjR46YH8zAgQOXL1++Z8+ee++9N/qBh+2lA8cSpKSkpKWlJSEh4bLLzD5jxOShdeDbcYZ1AQCxZ7fbt2/frrUuLy/3348dyHdDROjNCwkJCb473YYMGWK+NZPVdu7c+a8A27dvNz+Ye++995RLRL7b9CP10q5jCevf//631nry5MnmdzF5ntv17TgjuwCAWOrZs6dvJerVV19NSkoKW+eFF17QWj/22GNB5ZdeeqnW+uuvv05ISDDfmvlqHR7Mddddty6curo6rXV1dfW6devuuuuu6B2ZGWRCQsL5558f+m6PHj28Xq/Wevz48bE9NPPVzuAuACCWXn31Va31+vXrgx7wEaigoEBr/c9//jPoN/6KFSu01qtWrWpXa+arWRxMqHZdIzQzyP3792utb7311qDyn/3sZ1rrxsbGHj16mOnLx+ShWTkDZ0YXABAzV199tdZ6//79vXr1ilLNMIy//e1vWuvVq1cnJyf7Cm+88ca2trbjx48PHDiwXa2ZrGZxMGGZD0KTg/Q9ia2+vv7HP/6xvzArK8v34eZHH33U3DH9h8lDM38GQj/kF49dAMBpVFFR4VuG2hbBd7/7XV/NoUOHHjp0SGt9+PDh9957b+/evVrrtra2GTNmtLc1851GYmYwYZkPQpODTElJ+fOf/+y76PiPf/zjT3/60549e3wv161b196HhJk/NJPVQlMqHrsAgNNo27ZtEW4l+Y/AGxb69OmzYsWKf/3rX8ePH//666/Xr19/+eWXd6C1dnUaySkHE5b5IDQ/SJvNdtNNN7377rv79+9vbW09cuRIeXn5TTfdFOkxk7E6NDPVwqZU3HUBAAAAAAAAAAAAxJMZM2ZEv9bo99JLL3X1YOPb/wItllMXSjEiFwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=600x480>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ses.Show()\n",
    "# ses.Save(\"/Volumes/Shield/Coupling/10_min_timestep/tutorial.vs3\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "c1e1d457",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "a69a12a825bc47e09a0964e9b6fd5ecd",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "VBox(children=(VaporVisualizerWidget(imageData='/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDB…"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "ename": "TypeError",
     "evalue": "void Session::SetTimestep(int ts) =>\n    TypeError: could not convert argument 1 (int conversion expects an integer object)",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/ipywidgets/widgets/widget.py:773\u001b[0m, in \u001b[0;36mWidget._handle_msg\u001b[0;34m(self, msg)\u001b[0m\n\u001b[1;32m    771\u001b[0m         \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mbuffer_paths\u001b[39m\u001b[38;5;124m'\u001b[39m \u001b[38;5;129;01min\u001b[39;00m data:\n\u001b[1;32m    772\u001b[0m             _put_buffers(state, data[\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mbuffer_paths\u001b[39m\u001b[38;5;124m'\u001b[39m], msg[\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mbuffers\u001b[39m\u001b[38;5;124m'\u001b[39m])\n\u001b[0;32m--> 773\u001b[0m         \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mset_state\u001b[49m\u001b[43m(\u001b[49m\u001b[43mstate\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    775\u001b[0m \u001b[38;5;66;03m# Handle a state request.\u001b[39;00m\n\u001b[1;32m    776\u001b[0m \u001b[38;5;28;01melif\u001b[39;00m method \u001b[38;5;241m==\u001b[39m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mrequest_state\u001b[39m\u001b[38;5;124m'\u001b[39m:\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/ipywidgets/widgets/widget.py:655\u001b[0m, in \u001b[0;36mWidget.set_state\u001b[0;34m(self, sync_data)\u001b[0m\n\u001b[1;32m    652\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m name \u001b[38;5;129;01min\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mkeys:\n\u001b[1;32m    653\u001b[0m     from_json \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mtrait_metadata(name, \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mfrom_json\u001b[39m\u001b[38;5;124m'\u001b[39m,\n\u001b[1;32m    654\u001b[0m                                     \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_trait_from_json)\n\u001b[0;32m--> 655\u001b[0m     \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mset_trait(name, from_json(sync_data[name], \u001b[38;5;28mself\u001b[39m))\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/contextlib.py:126\u001b[0m, in \u001b[0;36m_GeneratorContextManager.__exit__\u001b[0;34m(self, typ, value, traceback)\u001b[0m\n\u001b[1;32m    124\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m typ \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[1;32m    125\u001b[0m     \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[0;32m--> 126\u001b[0m         \u001b[38;5;28;43mnext\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mgen\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    127\u001b[0m     \u001b[38;5;28;01mexcept\u001b[39;00m \u001b[38;5;167;01mStopIteration\u001b[39;00m:\n\u001b[1;32m    128\u001b[0m         \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;01mFalse\u001b[39;00m\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/traitlets/traitlets.py:1502\u001b[0m, in \u001b[0;36mHasTraits.hold_trait_notifications\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m   1500\u001b[0m \u001b[38;5;28;01mfor\u001b[39;00m changes \u001b[38;5;129;01min\u001b[39;00m cache\u001b[38;5;241m.\u001b[39mvalues():\n\u001b[1;32m   1501\u001b[0m     \u001b[38;5;28;01mfor\u001b[39;00m change \u001b[38;5;129;01min\u001b[39;00m changes:\n\u001b[0;32m-> 1502\u001b[0m         \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mnotify_change\u001b[49m\u001b[43m(\u001b[49m\u001b[43mchange\u001b[49m\u001b[43m)\u001b[49m\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/ipywidgets/widgets/widget.py:701\u001b[0m, in \u001b[0;36mWidget.notify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m    698\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m name \u001b[38;5;129;01min\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mkeys \u001b[38;5;129;01mand\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_should_send_property(name, \u001b[38;5;28mgetattr\u001b[39m(\u001b[38;5;28mself\u001b[39m, name)):\n\u001b[1;32m    699\u001b[0m         \u001b[38;5;66;03m# Send new state to front-end\u001b[39;00m\n\u001b[1;32m    700\u001b[0m         \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39msend_state(key\u001b[38;5;241m=\u001b[39mname)\n\u001b[0;32m--> 701\u001b[0m \u001b[38;5;28;43msuper\u001b[39;49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mnotify_change\u001b[49m\u001b[43m(\u001b[49m\u001b[43mchange\u001b[49m\u001b[43m)\u001b[49m\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/traitlets/traitlets.py:1517\u001b[0m, in \u001b[0;36mHasTraits.notify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m   1515\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mnotify_change\u001b[39m(\u001b[38;5;28mself\u001b[39m, change):\n\u001b[1;32m   1516\u001b[0m \u001b[38;5;250m    \u001b[39m\u001b[38;5;124;03m\"\"\"Notify observers of a change event\"\"\"\u001b[39;00m\n\u001b[0;32m-> 1517\u001b[0m     \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_notify_observers\u001b[49m\u001b[43m(\u001b[49m\u001b[43mchange\u001b[49m\u001b[43m)\u001b[49m\n",
      "File \u001b[0;32m~/opt/anaconda3/envs/vapor/lib/python3.9/site-packages/traitlets/traitlets.py:1564\u001b[0m, in \u001b[0;36mHasTraits._notify_observers\u001b[0;34m(self, event)\u001b[0m\n\u001b[1;32m   1561\u001b[0m \u001b[38;5;28;01melif\u001b[39;00m \u001b[38;5;28misinstance\u001b[39m(c, EventHandler) \u001b[38;5;129;01mand\u001b[39;00m c\u001b[38;5;241m.\u001b[39mname \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[1;32m   1562\u001b[0m     c \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mgetattr\u001b[39m(\u001b[38;5;28mself\u001b[39m, c\u001b[38;5;241m.\u001b[39mname)\n\u001b[0;32m-> 1564\u001b[0m \u001b[43mc\u001b[49m\u001b[43m(\u001b[49m\u001b[43mevent\u001b[49m\u001b[43m)\u001b[49m\n",
      "Cell \u001b[0;32mIn[12], line 4\u001b[0m, in \u001b[0;36msliderChanged\u001b[0;34m(change)\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21msliderChanged\u001b[39m(change):\n\u001b[0;32m----> 4\u001b[0m     \u001b[43mses\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mSetTimestep\u001b[49m\u001b[43m(\u001b[49m\u001b[43mchange\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m      5\u001b[0m     viz\u001b[38;5;241m.\u001b[39mRender(fast\u001b[38;5;241m=\u001b[39m\u001b[38;5;28;01mTrue\u001b[39;00m)\n",
      "\u001b[0;31mTypeError\u001b[0m: void Session::SetTimestep(int ts) =>\n    TypeError: could not convert argument 1 (int conversion expects an integer object)"
     ]
    }
   ],
   "source": [
    "from jupyter_vapor_widget import *\n",
    "viz = VaporVisualizerWidget(ses)\n",
    "def sliderChanged(change):\n",
    "    ses.SetTimestep(change)\n",
    "    viz.Render(fast=True)\n",
    "\n",
    "slider = widgets.IntSlider(\n",
    "    value=105, # initial value\n",
    "    min=100,\n",
    "    max=110,\n",
    "    step=1\n",
    ")\n",
    "slider.observe(sliderChanged, names='value')\n",
    "widgets.VBox([\n",
    "    viz,\n",
    "    widgets.HBox([widgets.Label(\"TimeStep:\"), slider])\n",
    "])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "814819f1",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}