Plotting of TKE budget terms series added

tools/plot_multtimser.py

@@ -19,8 +19,8 @@ font = {'family': 'serif',
 rc('font', **font)
 rc('text', usetex=True)
 rc('text.latex',unicode=True)
-rc('text.latex',preamble='\usepackage[utf8]{inputenc}')
-rc('text.latex',preamble='\usepackage[russian]{babel}')
+rc('text.latex',preamble=r'\usepackage[utf8]{inputenc}')
+rc('text.latex',preamble=r'\usepackage[russian]{babel}')
 
 nplots = int(input("Enter number of variables on a figure "))
 

tools/plot_timser.py

@@ -143,6 +143,7 @@ def plot_timeser(plotid,nvar,legendcommon) :
     uspeed_ = num; num +=1
     disch_ = num; num +=1
     epsilon_ = num; num +=1
+    TKE_budget_terms_ = num; num +=1
     ustar_ = num; num +=1
     windwork_ = num; num +=1
     mldmod_ = num; num +=1
@@ -196,6 +197,7 @@ def plot_timeser(plotid,nvar,legendcommon) :
      str(uspeed_) + " - horizontal component of speed \n" + \
      str(disch_) + " - x-component of discharge \n" + \
      str(epsilon_) + " - dissipation rate \n" + \
+     str(TKE_budget_terms_) + " - TKE budget terms \n" + \
      str(ustar_) + " - Friction velocity \n" + \
      str(windwork_) + " - Work by wind at the water surface \n" + \
      str(mldmod_) + " - Mixed-layer depth, m \n" + \
@@ -256,7 +258,7 @@ def plot_timeser(plotid,nvar,legendcommon) :
     #path_ = ['../results/IseoAugust2011_cdE-1/time_series/']#, \
         #         '../results/IseoAugust2011_cdE-1/time_series/']
     #path_ = ['../results/Uvs1979-2017/time_series/']#, \
-    path_ = ['../results/Mojai2016-2017_damw5/time_series/']
+    path_ = ['../results/river/time_series/']
     #path_ = ['../results/IseoJuneOctober2011_cdE-1/time_series/']
     #path_ = ['../results/IseoJuneOctober2011_botfric1backdiff1/time_series/']#, \
     #         '../results/IseoJuneOctober2011_botfric=1/time_series/'] #, \
@@ -475,6 +477,25 @@ def plot_timeser(plotid,nvar,legendcommon) :
         npaths = len(path_)
         sciformat = 1
         logscale = 1
+    if int(nplot) == TKE_budget_terms_: # dissipation rate, m^2/s^3
+        basename = 'TKE_budget_terms_series  1  1'
+        fileout  = 'TKE_budget_terms_wat'
+        nheader  = 10
+        ncol     = [6,7,8,9,10]
+        legend   = ['TKE tendency','TKE turbulent transport', \
+                    'TKE shear production','TKE production/sink by buoyancy', \
+                    'TKE dissipation rate, 0.4~m']
+        labelx   = 'Time, days'
+        labely   = 'TKE budget terms, m$^2$/s$^3$'
+        timescale = mdays # 365/12*24 - number of hours in a month
+        ctime = 5
+        #time0 = month_doy(5,ny)/mdays #Initial time
+        time0 = 0. #month_doy(7,ny)/mdays + 17./31. #Initial time, Bolshoi Vilui 2015
+        path = path_
+        nameexp = nameexp_
+        npaths = len(path_)
+        sciformat = 1
+        #logscale = 1
     if int(nplot) == ustar_: # Friction velocity
         basename = 'T_fluxes  1  1'
         fileout = 'Ustar_ser'