Changeset 2088 in lmdz_wrf

Timestamp:

Aug 22, 2018, 9:07:39 PM (7 years ago)

Author:

lfita

Message:

Adding histogram

Location:

trunk/tools

Files:

• drawing.py (modified) (1 diff)
• drawing_tools.py (modified) (17 diffs)

trunk/tools/drawing.py

@@ -10557 +10557 @@
 
     return
+#
 #foldn='/media/lluis/ExtDiskC_ext3/DATA/estudios/FPS_Alps/additional/IOP/sims/wrfinput_select/'
-#fils=foldn + '120lev_cdxwrf2/simin_vars.nc,' + foldn + '80lev_cdxwrf2/simin_vars.nc,' + foldn +      \
-#  '50lev_cdxwrf2/simin_vars.nc,' + foldn + '50lev_assigned/simin_vars.nc,' + foldn +         \
+#fils=foldn + '120lev_cdxwrf2/simin_vars.nc,' + foldn + '120lev_assigned/simin_vars.nc,' + foldn +      \
+#  '80lev_cdxwrf2/simin_vars.nc,' + foldn + '50lev_cdxwrf2/simin_vars.nc,' + foldn +         \
 #  '38lev_cdxwrf2/simin_vars.nc'
-#vals='120lev,80lev,50lev,50leva,38lev:auto:auto:0|4|2:FPS!Alps!vertical!levels:pdf:no:true'
+#vals='120lev,120leva,80lev,50lev,38lev:auto:auto:0|4|2:FPS!Alps!vertical!levels:pdf:no:true'
 #print fils
-draw_WRFeta_levels(fils, vals)
+#draw_WRFeta_levels(fils, vals)
 
 #quit()

trunk/tools/drawing_tools.py

@@ -12775 +12775 @@
     plt.rc('text', usetex=True)
 
-    fig, axs = plt.subplots(8)
-
-    plt.subplot(4,2,1)
+    fig, axs = plt.subplots(10)
+
+    plt.subplot(5,2,1)
     # Absolute heights
     Nlabs = len(labels)
@@ -12785 +12785 @@
         znw = aetaw[ilab]
 
-        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+        plt.plot(znw, hgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=4, \
           label='$hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
         if ilab == 0:
@@ -12803 +12803 @@
     plt.legend(loc=legvs[0], prop={'size':legvs[1]}, ncol=legvs[2])
 
-    plt.subplot(4,2,2)
+    plt.subplot(5,2,2)
     # Absolute pressures
     Nlabs = len(labels)
@@ -12822 +12822 @@
         if xx < nv: xx = xv
 
-        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+        plt.plot(znw, hgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=4, \
           label='$p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
         if ilab == 0:
@@ -12840 +12840 @@
     #plt.legend(loc=legvs[0], prop={'size':legvs[1]})
 
-    plt.subplot(4,2,3)
+    plt.subplot(5,2,3)
     # Relative heights
     Nlabs = len(labels)
@@ -12848 +12848 @@
         znw = aetau[ilab]
 
-        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+        plt.plot(znw, hgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=4, \
           label='$\delta hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
         if ilab == 0:
@@ -12865 +12865 @@
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(4,2,4)
+    plt.subplot(5,2,4)
     # Relative heights
     Nlabs = len(labels)
@@ -12874 +12874 @@
         znw = aetaw[ilab][1:dz-1]
 
-        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+        plt.plot(znw, hgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=4, \
           label='$\delta p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
         if ilab == 0:
@@ -12891 +12891 @@
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(4,2,5)
+    plt.subplot(5,2,5)
     # Absolute heights & relative
     Nlabs = len(labels)
@@ -12903 +12903 @@
         dz = len(hgtssea)
 
-        plt.plot(hgtssea[1:dz], dhgtssea, '-', marker=marks[ilab], linewidth=1.5,    \
+        plt.plot(hgtssea[1:dz], dhgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5,    \
           markersize=4, label='$hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
 #        if ilab == 0:
@@ -12921 +12921 @@
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(4,2,6)
+    plt.subplot(5,2,6)
     # Absolute pressures & relative
     Nlabs = len(labels)
@@ -12942 +12942 @@
         if xx < nv: xx = xv
 
-        plt.plot(hgtssea[1:dz], dhgtssea, '-', marker=marks[ilab], linewidth=1.5,    \
+        plt.plot(hgtssea[1:dz], dhgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5,    \
           markersize=4, label='$p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
 #        if ilab == 0:
@@ -12962 +12962 @@
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(4,2,7)
+    plt.subplot(5,2,7)
     # Absoulte heights & grid-level
     Nlabs = len(labels)
@@ -12972 +12972 @@
         zi = np.arange(1.,dz+1.)*100./dz
 
-        plt.plot(zi, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=0,  \
+        plt.plot(zi, hgtssea, linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=0,  \
           label='$\delta hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
         if ilab == 0:
@@ -12990 +12990 @@
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(4,2,8)
+    plt.subplot(5,2,8)
     # Absoulte pressures & grid-level
     Nlabs = len(labels)
@@ -13011 +13011 @@
         if xx < nv: xx = xv
 
-        plt.plot(zi, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=0,  \
-          label='$\delta p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
+        plt.plot(zi, hgtssea, color=cols[ilab], linestyle='solid', marker=marks[ilab], linewidth=1.5, markersize=0,  \
+          label='$\delta p_{sea}^{'+labels[ilab]+'}$')
         if ilab == 0:
-            plt.plot(zi, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,         \
+            plt.plot(zi, hgtsxhgt,  linestyle='-.', marker=marks[ilab], linewidth=1.5,         \
               markersize=0, label='$\delta p_{hgtmax}^{' + labels[ilab] + '}$',    \
               color=cols[ilab])
@@ -13031 +13031 @@
     #plt.legend(loc=0, prop={'size':10})
 
+    xhgtmax = np.max(ahgtsea[0])
+    hgtbins = np.arange(0.,xhgtmax,1000.)
+
+    plt.subplot(5,2,9)
+    # Histogram absoulte heights & grid-level
+    Nlabs = len(labels)
+    for ilab in range(Nlabs):
+        hgtssea = ahgtsea[ilab]
+        n, bins, patches = plt.hist(hgtssea, hgtbins, normed=True,                   \
+          color=cols[ilab], histtype='step')
+    plt.xlabel('height ($m$)', color='black', fontsize=6)
+    plt.ylabel('levels ($\%$)', color='black', fontsize=6)
+
+    ypos, ylabs = plt.yticks()
+    newlabs = []
+    for ylab in ypos: newlabs.append('{:.3g}'.format(np.float(ylab)*100.))
+    plt.yticks(ypos, newlabs)
+
+
+    plt.subplot(5,2,10)
+    # Histogram absoulte pressures & grid-level
+    pbins = np.arange(0.,105000., 5000.)
+    Nlabs = len(labels)
+    for ilab in range(Nlabs):
+        hgtssea = apsea[ilab]
+        n, bins, patches = plt.hist(hgtssea, pbins, normed=True,                     \
+          color=cols[ilab], histtype='step')
+
+    plt.xlim(100000.,0.)
+    plt.xlabel('pressure ($Pa$)', color='black', fontsize=6)
+    plt.ylabel('levels ($\%$)', color='black', fontsize=6)
+
+    ypos, ylabs = plt.yticks()
+    newlabs = []
+    for ylab in ypos: newlabs.append('{:.3g}'.format(np.float(ylab)*100.))
+    plt.yticks(ypos, newlabs)
+
     fig.suptitle(gen.latex_text(figtitle))
-    plt.subplots_adjust(wspace=0.4, hspace=0.6)
+    plt.subplots_adjust(wspace=0.4, hspace=0.8)
 
     figname = 'WRFeta_levels'