add some error message

Author: patrickscholz

templates_notebooks/template_hquiver.ipynb

@@ -55,7 +55,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 43,
+   "execution_count": 46,
    "metadata": {
     "slideshow": {
      "slide_type": ""
@@ -113,7 +113,7 @@
     "day            = None\n",
     "record         = None \n",
     "box            = None\n",
-    "depth          = 100\n",
+    "depth          = 'bottom'\n",
     "do_data_r2g    = False  # set to False if u,v data are already in geo-coordinates\n",
     "\n",
     "#___Define Reference Data, Year, Mon ...______________________________________________\n",
@@ -170,7 +170,7 @@
     "mesh_opt          = dict({'color':'k', 'linewidth':0.10})\n",
     "do_enum           = False  # do enumeration of panels\n",
     "do_reffig         = True   # plot reference fig when doing anomalies \n",
-    "do_clim           = True  # plot climatolgy values when doing absoluts\n",
+    "do_clim           = False  # plot climatolgy values when doing absoluts\n",
     "ax_title          = None\n",
     "cb_label          = None\n",
     "save_dpi          = 150\n",
@@ -1860,7 +1860,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 44,
+   "execution_count": 47,
    "metadata": {
     "slideshow": {
      "slide_type": ""
@@ -1872,14 +1872,12 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      " --> elasped time to load clim: 0.27 min.\n",
-      " --> clim uses 0.15 Gb:\n",
-      "100\n",
+      "bottom\n",
       "0 /albedo/work/user/pscholz/results/dart_linfs_pc0_ctrl_1/1/ dart test\n",
-      " --> elasped time to load data: 0.16 min.\n",
-      " --> data uses 0.15 Gb:\n",
+      " --> elasped time to load data: 0.97 min.\n",
+      " --> data uses 0.18 Gb:\n",
       "\n",
-      " --> total elasped time to process data: 25.58 sec.\n"
+      " --> total elasped time to process data: 58.34 sec.\n"
      ]
     }
    ],
@@ -1898,9 +1896,8 @@
     "    elif vname=='vec+unod+vnod' and  which_clim.lower()=='glorys12_v1': clim_vnameu,  clim_vnamev, onelem = 'uo', 'vo', False\n",
     "    else: raise ValueError('this variable is not supported for climatology') \n",
     "    # clim = xr.open_Dataset(clim_path).rename({clim_vnameu:'u', clim_vnamev:'v', })\n",
-    "    clim, depth = tpv.load_climatology_uv(mesh, clim_path, clim_vnameu, clim_vnamev, onelem=onelem, depth=depth, mon=mon)\n",
+    "    clim, depth = tpv.load_climatology_uv(mesh, clim_path, clim_vnameu, clim_vnamev, onelem=onelem, depth=depth, mon=mon, descript=which_clim)\n",
     "    clim = clim.rename({clim_vnameu:'u', clim_vnamev:'v'})\n",
-    "    clim.attrs['descript'] = which_clim\n",
     "    print(' --> elasped time to load clim: {:3.2f} min.'.format( (clock.time()-ts)/60  ))        \n",
     "    print(' --> clim uses {:3.2f} Gb:'.format(clim.nbytes/(1024**3)))\n",
     "    data_list.append(clim)    \n",

tripyview/sub_climatology.py

@@ -447,6 +447,8 @@ def load_climatology_uv(mesh, datapath, vnameu, vnamev, mon=None, depth=None, de
     #___________________________________________________________________________
     # do vertical interpolation
     if (depth) is not None:
+        if isinstance(depth, str):raise ValueError(' depth=\'bottom\' is not supported for the climatology')
+
         #_______________________________________________________________________
         # select depth level indices that are needed to interpolate the values
         # in depth list,array