{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "
currently SET data sets:\n", " 1> http://vesg.ipsl.upmc.fr/thredds/dodsC/work_thredds/p86caub/IPSLCM6/DEVT/pdControl/CM608-LR-pdCtrl-01/MONITORING/files/ATM_t2m_global_ave.nc (default)\n", " name title I J K L\n", " T2M_GLOBAL\n", " ... ... ... 1:1680\n", " \n", "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%ferret\n", "\n", "use \"http://vesg.ipsl.upmc.fr/thredds/dodsC/work_thredds/p86caub/IPSLCM6/DEVT/pdControl/CM608-LR-pdCtrl-01/MONITORING/files/ATM_t2m_global_ave.nc\"\n", "\n", "show data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%ferret\n", "\n", "let var = T2M_GLOBAL[l=@sbx:120]\n", "plot var" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
Message: vardict is now available in python as a dictionary containing the variable's metadata and data array." ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "{'title': 'var', 'data': array([[[[[[-1.e+34]],\n", " \n", " [[-1.e+34]],\n", " \n", " [[-1.e+34]],\n", " \n", " ...,\n", " \n", " [[-1.e+34]],\n", " \n", " [[-1.e+34]],\n", " \n", " [[-1.e+34]]]]]]), 'missing_value': array([-1.e+34]), 'data_unit': '', 'axis_types': (7,\n", " 7,\n", " 7,\n", " 4,\n", " 7,\n", " 7), 'axis_names': ('', '', '', 'TIME_COUNTER', '', ''), 'axis_units': ('',\n", " '',\n", " '',\n", " 'CALTYPE_NOLEAP',\n", " '',\n", " ''), 'axis_coords': (None,\n", " None,\n", " None,\n", " array([[ 16, 1, 1950, 12, 0, 0],\n", " [ 15, 2, 1950, 0, 0, 0],\n", " [ 16, 3, 1950, 12, 0, 0],\n", " ...,\n", " [ 16, 10, 2089, 12, 0, 0],\n", " [ 16, 11, 2089, 0, 0, 0],\n", " [ 16, 12, 2089, 12, 0, 0]], dtype=int32),\n", " None,\n", " None)}" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "%ferret_getdata vardict = var\n", "vardict" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Prepare data extracted from pyferret" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import datetime" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([nan, nan, nan, ..., nan, nan, nan])" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "values = vardict['data'].squeeze()\n", "values[values == -1E+34] = np.nan\n", "values" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 16 1 1950 12 0 0]\n" ] }, { "data": { "text/plain": [ "'1950-01-16 12:00:00'" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "datesNum = vardict['axis_coords'][3]\n", "(day, month, year, hour, minute, second) = datesNum[0]\n", "print(datesNum[0])\n", "d = datetime.datetime(year, month, day, hour, minute, second)\n", "d.strftime('%Y-%m-%d %H:%M:%S')" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[datetime.datetime(1950, 1, 16, 12, 0),\n", " datetime.datetime(1950, 2, 15, 0, 0),\n", " datetime.datetime(1950, 3, 16, 12, 0),\n", " datetime.datetime(1950, 4, 16, 0, 0)]" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dates = []\n", "datesString = []\n", "for d in datesNum:\n", " date = datetime.datetime(d[2], d[1], d[0], d[3], d[4], d[5])\n", " dates.append(date)\n", " datesString.append(date.strftime('%Y-%m-%d %H:%M:%S'))\n", "dates[0:4]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot using bokeh" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "colors = [\"#8c564b\",\"#1f77b4\",\"#2ca02c\",\"#d62728\",\"#9467bd\",\"#e377c2\",\"#7f7f7f\",\"#bcbd22\",\"#17bec\"]" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "\n", "(function(root) {\n", " function now() {\n", " return new Date();\n", " }\n", "\n", " var force = true;\n", "\n", " if (typeof root._bokeh_onload_callbacks === \"undefined\" || force === true) {\n", " root._bokeh_onload_callbacks = [];\n", " root._bokeh_is_loading = undefined;\n", " }\n", "\n", " var JS_MIME_TYPE = 'application/javascript';\n", " var HTML_MIME_TYPE = 'text/html';\n", " var EXEC_MIME_TYPE = 'application/vnd.bokehjs_exec.v0+json';\n", " var CLASS_NAME = 'output_bokeh rendered_html';\n", "\n", " /**\n", " * Render data to the DOM node\n", " */\n", " function render(props, node) {\n", " var script = document.createElement(\"script\");\n", " node.appendChild(script);\n", " }\n", "\n", " /**\n", " * Handle when an output is cleared or removed\n", " */\n", " function handleClearOutput(event, handle) {\n", " var cell = handle.cell;\n", "\n", " var id = cell.output_area._bokeh_element_id;\n", " var server_id = cell.output_area._bokeh_server_id;\n", " // Clean up Bokeh references\n", " if (id != null && id in Bokeh.index) {\n", " Bokeh.index[id].model.document.clear();\n", " delete Bokeh.index[id];\n", " }\n", "\n", " if (server_id !== undefined) {\n", " // Clean up Bokeh references\n", " var cmd = \"from bokeh.io.state import curstate; print(curstate().uuid_to_server['\" + server_id + \"'].get_sessions()[0].document.roots[0]._id)\";\n", " cell.notebook.kernel.execute(cmd, {\n", " iopub: {\n", " output: function(msg) {\n", " var id = msg.content.text.trim();\n", " if (id in Bokeh.index) {\n", " Bokeh.index[id].model.document.clear();\n", " delete Bokeh.index[id];\n", " }\n", " }\n", " }\n", " });\n", " // Destroy server and session\n", " var cmd = \"import bokeh.io.notebook as ion; ion.destroy_server('\" + server_id + \"')\";\n", " cell.notebook.kernel.execute(cmd);\n", " }\n", " }\n", "\n", " /**\n", " * Handle when a new output is added\n", " */\n", " function handleAddOutput(event, handle) {\n", " var output_area = handle.output_area;\n", " var output = handle.output;\n", "\n", " // limit handleAddOutput to display_data with EXEC_MIME_TYPE content only\n", " if ((output.output_type != \"display_data\") || (!output.data.hasOwnProperty(EXEC_MIME_TYPE))) {\n", " return\n", " }\n", "\n", " var toinsert = output_area.element.find(\".\" + CLASS_NAME.split(' ')[0]);\n", "\n", " if (output.metadata[EXEC_MIME_TYPE][\"id\"] !== undefined) {\n", " toinsert[toinsert.length - 1].firstChild.textContent = output.data[JS_MIME_TYPE];\n", " // store reference to embed id on output_area\n", " output_area._bokeh_element_id = output.metadata[EXEC_MIME_TYPE][\"id\"];\n", " }\n", " if (output.metadata[EXEC_MIME_TYPE][\"server_id\"] !== undefined) {\n", " var bk_div = document.createElement(\"div\");\n", " bk_div.innerHTML = output.data[HTML_MIME_TYPE];\n", " var script_attrs = bk_div.children[0].attributes;\n", " for (var i = 0; i < script_attrs.length; i++) {\n", " toinsert[toinsert.length - 1].firstChild.setAttribute(script_attrs[i].name, script_attrs[i].value);\n", " }\n", " // store reference to server id on output_area\n", " output_area._bokeh_server_id = output.metadata[EXEC_MIME_TYPE][\"server_id\"];\n", " }\n", " }\n", "\n", " function register_renderer(events, OutputArea) {\n", "\n", " function append_mime(data, metadata, element) {\n", " // create a DOM node to render to\n", " var toinsert = this.create_output_subarea(\n", " metadata,\n", " CLASS_NAME,\n", " EXEC_MIME_TYPE\n", " );\n", " this.keyboard_manager.register_events(toinsert);\n", " // Render to node\n", " var props = {data: data, metadata: metadata[EXEC_MIME_TYPE]};\n", " render(props, toinsert[toinsert.length - 1]);\n", " element.append(toinsert);\n", " return toinsert\n", " }\n", "\n", " /* Handle when an output is cleared or removed */\n", " events.on('clear_output.CodeCell', handleClearOutput);\n", " events.on('delete.Cell', handleClearOutput);\n", "\n", " /* Handle when a new output is added */\n", " events.on('output_added.OutputArea', handleAddOutput);\n", "\n", " /**\n", " * Register the mime type and append_mime function with output_area\n", " */\n", " OutputArea.prototype.register_mime_type(EXEC_MIME_TYPE, append_mime, {\n", " /* Is output safe? */\n", " safe: true,\n", " /* Index of renderer in `output_area.display_order` */\n", " index: 0\n", " });\n", " }\n", "\n", " // register the mime type if in Jupyter Notebook environment and previously unregistered\n", " if (root.Jupyter !== undefined) {\n", " var events = require('base/js/events');\n", " var OutputArea = require('notebook/js/outputarea').OutputArea;\n", "\n", " if (OutputArea.prototype.mime_types().indexOf(EXEC_MIME_TYPE) == -1) {\n", " register_renderer(events, OutputArea);\n", " }\n", " }\n", "\n", " \n", " if (typeof (root._bokeh_timeout) === \"undefined\" || force === true) {\n", " root._bokeh_timeout = Date.now() + 5000;\n", " root._bokeh_failed_load = false;\n", " }\n", "\n", " var NB_LOAD_WARNING = {'data': {'text/html':\n", " \"
\\n\"+\n", " \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n", " \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n", " \"
\\n\"+\n", " \"\\n\"+\n",
" \"from bokeh.resources import INLINE\\n\"+\n",
" \"output_notebook(resources=INLINE)\\n\"+\n",
" \"
\\n\"+\n",
" \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"
\\n\"+\n \"