Changeset 2082 in lmdz_wrf

Timestamp:

Aug 12, 2018, 12:05:20 AM (7 years ago)

Author:

lfita

Message:

Final version of the `draw_WRFeta_levels' function:

• Adding file with data
• pressure values
• number of columns into legend

Location:

trunk/tools

Files:

• drawing.py (modified) (5 diffs)
• drawing_tools.py (modified) (9 diffs)

trunk/tools/drawing.py

@@ -122 +122 @@
 
 errormsg = 'ERROR -- error -- ERROR -- error'
+infmsg = 'INFORMATION -- information -- INFORMATION -- information'
 warnmsg = 'WARNING -- waring -- WARNING -- warning'
 fillValue=1.e20
@@ -10304 +10305 @@
     """ Function to plot vertical levels from WRF wrfinput_d[nn] file
       filenames= ',' list of name of files to use
-      values = [labs]:[colors]:[markers]:[legvals]:[imgtit]:[imgkind]:[close]
+      values = [labs]:[colors]:[markers]:[legvals]:[imgtit]:[imgkind]:[newfile]:[close]
         [labs]= ',' list of labels for the plot (LaTeX like)
         [colors]= ',' list of colors for the lines ('auto' for automatic)
         [markers]= '@' list of markers for the lines ('auto' for automatic)
-        [legvals]=[loclegend]|[fonstsize] values for the legend
+        [legvals]=[loclegend]|[fonstsize]|[ncol] values for the legend
           [locleg]: location of the legend (0, automatic)
             1: 'upper right', 2: 'upper left', 3: 'lower left', 4: 'lower right',
@@ -10314 +10315 @@
             9: 'upper center', 10: 'center'
           [fontsize]: font size for the legend ('auto' for 12)
+          [ncol]: number of columns ('auto' for 1)
         [imgtit]= title of the image ('!' for spaces)
         [imgkind]= kind of file output of the image (ps, pns, pdf, ...)
+        [newfile]= whether the possible pre-existing file with the required data be 
+          should removed or not ('yes/no')
         [close]= whether figure should be closed or not
     """
+    import subprocess as sub
+
     fname = 'draw_WRFeta_levels'
 
@@ -10325 +10331 @@
         quit()
 
-    expectargs = '[labs]:[colors]:[markers]:[legvals]:[imgtit]:[imgkind]:[close]'
+    expectargs = '[labs]:[colors]:[markers]:[legvals]:[imgtit]:[imgkind]:' +         \
+      '[newfile]:[close]'
     drw.check_arguments(fname,values,expectargs,':')
 
@@ -10336 +10343 @@
     imgtit = values.split(':')[4].replace('!', ' ')
     imgkind = values.split(':')[5]
-    close = gen.Str_Bool(values.split(':')[6])
+    newfile = gen.Str_Bool(values.split(':')[6])
+    close = gen.Str_Bool(values.split(':')[7])
     
     files = filenames.split(',')
+    Nfiles = len(files)
+
     allhgtsea = []
     allxhgt = []
+    allpsea = []
+    allpxhgt = []
     allhgtxhgt = []
     alldhgtsea = []
     alldhgtxhgt = []
+    alldpsea = []
+    alldpxhgt = []
     alletaw = []
     alletau = []
-
-    for filen in files:
-        if not os.path.isfile(filen):
-            print errormsg
-            print '  ' + fname + ": file '" + filen + "' does not exist !!"
-            quit(-1)
-
-        onc = NetCDFFile(filen, 'r')
-        ophb = onc.variables['PHB']
-        olandseamask = onc.variables['LANDMASK']
-        ohgt = onc.variables['HGT']
-        oznw = onc.variables['ZNW']
-        oznu = onc.variables['ZNU']
-
-        dimx = ophb.shape[3]
-        dimy = ophb.shape[2]
-        dimz = ophb.shape[1]
-        dimt = ophb.shape[0]
-
-        for j in range(10,dimy):
-            for i in range(10,dimx):
-                if olandseamask[0,j,i] == 0:
-                    hgtssea = ophb[0,:,j,i]/9.8
-                    break
-
-        maxhgt = np.max(ohgt[0,:,:])
-        ijxhgt = gen.index_mat(ohgt[0,:,:], maxhgt)
-        hgtsxhgt = ophb[0,:, ijxhgt[0], ijxhgt[1]]/9.8
-
-        allhgtsea.append(hgtssea)
-        allxhgt.append(maxhgt)
-        allhgtxhgt.append(hgtsxhgt)
-        alldhgtsea.append(hgtssea[1:dimz]-hgtssea[0:dimz-1])
-        alldhgtxhgt.append(hgtsxhgt[1:dimz]-hgtsxhgt[0:dimz-1])
-        alletaw.append(oznw[0,:])
-        alletau.append(oznu[0,:])
-
+    alldimz = []
+
+    # Creation of a file with the data
+    titS = imgtit.replace(' ','-')
+    ofilen = 'WRFeta_levels_' + titS + '.nc'
+
+    # Removing pre-existing file
+    if newfile:
+        if os.path.isfile(ofilen):
+            print '  ' + fname + ": removing pre-existing file with required '" +    \
+              "data '" + ofilen + "' !!"
+            sub.call('rm ' + ofilen, shell=True)
+
+    if not os.path.isfile(ofilen):
+
+        maxdz = -10
+        for filen in files:
+            if not os.path.isfile(filen):
+                print errormsg
+                print '  ' + fname + ": file '" + filen + "' does not exist !!"
+                quit(-1)
+
+            onc = NetCDFFile(filen, 'r')
+            ophb = onc.variables['PHB']
+            opb = onc.variables['PB']
+            olandseamask = onc.variables['LANDMASK']
+            ohgt = onc.variables['HGT']
+            oznw = onc.variables['ZNW']
+            oznu = onc.variables['ZNU']
+    
+            dimx = ophb.shape[3]
+            dimy = ophb.shape[2]
+            dimz = ophb.shape[1]
+            dimt = ophb.shape[0]
+            dimzp = opb.shape[1]
+
+            if maxdz < dimz: maxdz = dimz
+
+            for j in range(10,dimy):
+                for i in range(10,dimx):
+                    if olandseamask[0,j,i] == 0:
+                        hgtssea = ophb[0,:,j,i]/9.8
+                        pssea = opb[0,:,j,i]
+                        break
+
+            maxhgt = np.max(ohgt[0,:,:])
+            ijxhgt = gen.index_mat(ohgt[0,:,:], maxhgt)
+            hgtsxhgt = ophb[0,:, ijxhgt[0], ijxhgt[1]]/9.8
+            psxhgt = opb[0,:, ijxhgt[0], ijxhgt[1]]
+
+            allhgtsea.append(hgtssea)
+            allxhgt.append(maxhgt)
+            allhgtxhgt.append(hgtsxhgt)
+            allpsea.append(pssea)
+            allpxhgt.append(psxhgt)
+            alldhgtsea.append(hgtssea[1:dimz]-hgtssea[0:dimz-1])
+            alldhgtxhgt.append(hgtsxhgt[1:dimz]-hgtsxhgt[0:dimz-1])
+            alldpsea.append(pssea[0:dimzp-1]-pssea[1:dimzp])
+            alldpxhgt.append(psxhgt[0:dimzp-1]-psxhgt[1:dimzp])
+            alletaw.append(oznw[0,:])
+            alletau.append(oznu[0,:])
+            alldimz.append(dimz)
+
+            onc.close()
+
+        print infmsg
+        print '  ' + fname + ": creation of file '" + ofilen + "' to keep required"+ \
+          "values to plot !!"
+        onewnc = NetCDFFile(ofilen, 'w')
+
+        # Dimensions
+        newdim = onewnc.createDimension('Nfiles', Nfiles)
+        newdim = onewnc.createDimension('Lstring', 250)
+        newdim = onewnc.createDimension('evert', maxdz)
+
+        # Variabledimensions
+        newvar = onewnc.createVariable('Nfiles', 'c', ('Nfiles', 'Lstring'))
+        ncvar.writing_str_nc(newvar, files, 250)
+        ncvar.basicvardef(newvar,'files','Files used to retrieve data','-')
+
+        newvar = onewnc.createVariable('evert_file', 'i', ('Nfiles'))
+        ncvar.basicvardef(newvar,'evert_file','number vertical layers per file','-')
+        newvar[:] = alldimz
+
+        # Variables
+        newvar = onewnc.createVariable('hgtsea', 'f', ('Nfiles', 'evert'),           \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'hgtsea','heights above sea point', 'm')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]] = allhgtsea[iff]
+
+        newvar = onewnc.createVariable('pressea', 'f', ('Nfiles', 'evert'),          \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'pressea','pressures above sea point', 'Pa')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-1] = allpsea[iff]
+
+        newvar = onewnc.createVariable('hgtmax', 'f', ('Nfiles'))
+        ncvar.basicvardef(newvar,'hgtmax','maximum height', 'm')
+        newvar[:] = allxhgt
+
+        newvar = onewnc.createVariable('hgtxhgt', 'f', ('Nfiles', 'evert'),          \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'hgtxhgt','heights above maximum height', 'm')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]] = allhgtxhgt[iff]
+
+        newvar = onewnc.createVariable('presxhgt', 'f', ('Nfiles', 'evert'),         \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'presxhgt','pressures above maximum height', 'Pa')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-1] = allpxhgt[iff]
+
+        onewnc.sync()
+
+        newvar = onewnc.createVariable('dhgtsea', 'f', ('Nfiles', 'evert'),          \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'dhgtsea','delta heights above sea point', 'm')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-1] = alldhgtsea[iff]
+
+        newvar = onewnc.createVariable('dhgtxhgt', 'f', ('Nfiles', 'evert'),         \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'dhgtxhgt','delta heights above maximum height', 'm')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-1] = alldhgtxhgt[iff]
+
+        newvar = onewnc.createVariable('dpressea', 'f', ('Nfiles', 'evert'),         \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'dpressea','delta pressures above sea point', 'Pa')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-2] = alldpsea[iff]
+
+        newvar = onewnc.createVariable('dpresxhgt', 'f', ('Nfiles', 'evert'),        \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'dpresxhgt','delta pressures above maximum height', \
+         'Pa')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-2] = alldpxhgt[iff]
+
+        newvar = onewnc.createVariable('eta_full', 'f', ('Nfiles', 'evert'),         \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'eta_full','full eta layers', '-')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]] = alletaw[iff]
+
+        newvar = onewnc.createVariable('eta_half', 'f', ('Nfiles', 'evert'),         \
+          fill_value=gen.fillValueF)
+        ncvar.basicvardef(newvar,'eta_half','half eta layers', '-')
+        for iff in range(Nfiles): newvar[iff,0:alldimz[iff]-1] = alletau[iff]
+
+        onewnc.sync()
+        onewnc.close()
+        print fname + ": succesfull written of file '" + ofilen + "' with all data !!"
+
+    else:
+        print fname + ": readding existing file '" + ofilen + "' with all data !!"
+        print "    if a new file is required set variable 'newfile' to 'yes'" 
+        onc = NetCDFFile(ofilen, 'r')
+
+        ovar = onc.variables['hgtsea']
+        allhgtseav = ovar[:]
+        ovar = onc.variables['pressea']
+        allpseav = ovar[:]
+        ovar = onc.variables['hgtmax']
+        allxhgtv = ovar[:]
+        ovar = onc.variables['hgtxhgt']
+        allhgtxhgtv = ovar[:]
+        ovar = onc.variables['presxhgt']
+        allpxhgtv = ovar[:]
+        ovar = onc.variables['dhgtsea']
+        alldhgtseav = ovar[:]
+        ovar = onc.variables['dhgtxhgt']
+        alldhgtxhgtv = ovar[:]
+        ovar = onc.variables['dpressea']
+        alldpseav = ovar[:]
+        ovar = onc.variables['dpresxhgt']
+        alldpxhgtv = ovar[:]
+        ovar = onc.variables['eta_full']
+        alletawv = ovar[:]
+        ovar = onc.variables['eta_half']
+        alletauv = ovar[:]
+ 
         onc.close()
 
+        for iff in range(Nfiles):
+            allhgtsea.append(allhgtseav[iff,:].compressed())
+            allxhgt.append(allxhgtv[iff])
+            allhgtxhgt.append(allhgtxhgtv[iff,:].compressed())
+            allpsea.append(allpseav[iff,:].compressed())
+            allpxhgt.append(allpxhgtv[iff,:].compressed())
+            alldhgtsea.append(alldhgtseav[iff,:].compressed())
+            alldhgtxhgt.append(alldhgtxhgtv[iff,:].compressed())
+            alldpsea.append(alldpseav[iff,:].compressed())
+            alldpxhgt.append(alldpxhgtv[iff,:].compressed())
+            alletaw.append(alletawv[iff,:].compressed())
+            alletau.append(alletauv[iff,:].compressed())
+
     # Legend
-    locleg, legfontsize = drw.legend_values(legvals,'|')
-
-    drw.plot_WRFeta_levels(allhgtsea, allxhgt, allhgtxhgt, alldhgtsea, alldhgtxhgt,  \
-      alletaw, alletau, labs, colors, markers, [locleg, legfontsize], imgtit,        \
-      imgkind, close)
+    locleg, legfontsize, legncol = drw.legend_values(legvals,'|')
+
+    drw.plot_WRFeta_levels(allhgtsea, allpsea, allxhgt, allhgtxhgt, allpxhgt,        \
+      alldhgtsea, alldhgtxhgt, alldpsea, alldpxhgt, alletaw, alletau, labs, colors,  \
+      markers, [locleg, legfontsize, legncol], imgtit, imgkind, close)
 
     return
-#foldn='/media/lluis/ExtDiskC_ext3/DATA/estudios/FPS_Alps/additional/IOP/sims/'
-#fils=foldn + '120lev/simin_vars.nc,' + foldn + '80lev/simin_vars.nc,' + foldn +      \
-#  '50lev/simin_vars.nc,' + foldn + '50lev_assigned/simin_vars.nc,' + foldn +         \
-#  '38lev/simin_vars.nc'
-#vals='120lev,80lev,50lev,50leva,38lev:auto:auto:0|8:FPS!Alps!vertical!levels:pdf:true'
+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 +         \
+  '38lev_cdxwrf2/simin_vars.nc'
+vals='120lev,80lev,50lev,50leva,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

@@ -1823 +1823 @@
             9: 'upper center', 10: 'center'
           [fontsize]: font size for the legend ('auto' for 12)
+          [ncol]: number of columns for the legend ('auto' for 1)
       char= separation character
     """
@@ -1833 +1834 @@
     else:
         fontsize = int(fontsize0)
-
-    return locleg, fontsize
+    if len(legstring.split(char)) == 2:
+        return locleg, fontsize
+    else:
+        ncol = int(legstring.split(char)[2])
+    return locleg, fontsize, ncol
 
 #######    #######    #######    #######    #######    #######    #######    #######    #######    #######
@@ -12733 +12737 @@
         plt.plot(stlon, stlat, stmrk, color=stcol, markersize=stsiz)
         if stlab != 'None':
-            plt.annotate(gen.latex_text(stlab), xy=(stlon+ddlon,stlat),         \
+            plt.annotate(gen.latex_text(stlab), xy=(stlon+ddlon,stlat),              \
               xycoords='data', color=stcol, fontsize=fontsiz)
 
@@ -12743 +12747 @@
     return
 
-def plot_WRFeta_levels(ahgtsea, axhgt, ahgtxhgt, adhgtsea, adhgtxhgt, aetaw, aetau,  \
-     labels, cols, marks, legvs, figtitle, kfig, close):
+def plot_WRFeta_levels(ahgtsea, apsea, axhgt, ahgtxhgt, apxhgt, adhgtsea, adhgtxhgt, \
+  adpsea, adpxhgt, aetaw, aetau, labels, cols, marks, legvs, figtitle, kfig, close):
     """ Function to plot different sets of WRF eta-levels
       ahgtsea: all hegights of a point above sea
+      apsea: all pressures of a point above sea
       axhgt: all maximum hegights
       ahgtxhgt: all hegights of the point above mximum height
+      apxhgt: all pressures of the point above mximum height
       adhgtsea: all relative heights of a point above sea
       adhgtxhgt: all relative heights of the point above maximum topo
+      adpsea: all relative pressures of a point above sea
+      adpxhgt: all relative pressures of the point above maximum topo
       aetaw: all eta-full values
       aetau: all eta-half values
@@ -12756 +12764 @@
       cols: colors
       marks: markers
-      legvs: [legloc, legfs] legend values
+      legvs: [legloc, legfs, legncol] legend values
       figtitle: title of figure
       kfig: kind of file of graphical output
@@ -12767 +12775 @@
     plt.rc('text', usetex=True)
 
-    fig, ax = plt.subplots(4)
-
-    plt.subplot(2,2,1)
+    fig, axs = plt.subplots(8)
+
+    plt.subplot(4,2,1)
     # Absolute heights
     Nlabs = len(labels)
@@ -12787 +12795 @@
               markersize=4, label=None, color=cols[ilab])
 
-    plt.xlabel('$\eta$-level', color='k')
-    plt.ylabel('absolute height ($m$)', color='black')
+    plt.xlabel('$\eta$-level', color='k', fontsize=6)
+    plt.ylabel('absolute height ($m$)', color='black', fontsize=6)
+    plt.xlim(1.,0.)
     plt.yscale('log')
 
-    plt.title('Absolute Heights',position=titleloc)
-    plt.legend(loc=legvs[0], prop={'size':legvs[1]})
-
-    plt.subplot(2,2,2)
+    plt.title('Heights',position=titleloc, fontsize=8)
+    plt.legend(loc=legvs[0], prop={'size':legvs[1]}, ncol=legvs[2])
+
+    plt.subplot(4,2,2)
+    # Absolute pressures
+    Nlabs = len(labels)
+    nx = 100000.
+    xx = -100000.
+    for ilab in range(Nlabs):
+        hgtssea = apsea[ilab]/100.
+        hgtsxhgt = apxhgt[ilab]/100.
+        znw = aetau[ilab]
+
+        nv = np.min(hgtssea)
+        if nx > nv: nx = nv
+        nv = np.min(hgtsxhgt)
+        if nx > nv: nx = nv
+        xv = np.max(hgtssea)
+        if xx < nv: xx = xv
+        xv = np.max(hgtsxhgt)
+        if xx < nv: xx = xv
+
+        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+          label='$p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
+        if ilab == 0:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label='$p_{hgtmax}^{'+labels[ilab]+'}$', color=cols[ilab])
+        else:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label=None, color=cols[ilab])
+
+    plt.xlabel('$\eta$-level', color='k', fontsize=6)
+    plt.ylabel('absolute pressure ($hPa$)', color='black', fontsize=6)
+    plt.xlim(1.,0.)
+    plt.ylim(xx,nx)
+    #plt.yscale('log')
+
+    plt.title('pressures',position=titleloc, fontsize=8)
+    #plt.legend(loc=legvs[0], prop={'size':legvs[1]})
+
+    plt.subplot(4,2,3)
+    # Relative heights
+    Nlabs = len(labels)
+    for ilab in range(Nlabs):
+        hgtssea = adhgtsea[ilab]
+        hgtsxhgt = adhgtxhgt[ilab]
+        znw = aetau[ilab]
+
+        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+          label='$\delta hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
+        if ilab == 0:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label='$\delta hgt_{hgtmax}^{'+labels[ilab]+'}$',        \
+              color=cols[ilab])
+        else:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label=None, color=cols[ilab])
+    plt.xlabel('$\eta$-level', color='k', fontsize=6)
+    plt.ylabel('$\delta height$ ($m$)', color='black', fontsize=6)
+    plt.xlim(1.,0.)
+    plt.yscale('log')
+
+    #plt.title('Relative Heights',position=titleloc, fontsize=8)
+    #plt.legend(loc=0, prop={'size':10})
+
+    plt.subplot(4,2,4)
+    # Relative heights
+    Nlabs = len(labels)
+    for ilab in range(Nlabs):
+        hgtssea = adpsea[ilab]
+        hgtsxhgt = adpxhgt[ilab]
+        dz = len(aetaw[ilab])
+        znw = aetaw[ilab][1:dz-1]
+
+        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
+          label='$\delta p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
+        if ilab == 0:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label='$\delta p_{hgtmax}^{'+labels[ilab]+'}$',        \
+              color=cols[ilab])
+        else:
+            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
+              markersize=4, label=None, color=cols[ilab])
+    plt.xlabel('$\eta$-level', color='k', fontsize=6)
+    plt.ylabel('$\delta pressure$ ($Pa$)', color='black', fontsize=6)
+    plt.xlim(1.,0.)
+    plt.yscale('log')
+
+    #plt.title('Relative pressures',position=titleloc, fontsize=8)
+    #plt.legend(loc=0, prop={'size':10})
+
+    plt.subplot(4,2,5)
     # Absolute heights & relative
     Nlabs = len(labels)
+    nx = 10000.
+    xx = -10000.
     for ilab in range(Nlabs):
         hgtssea = ahgtsea[ilab]
@@ -12814 +12913 @@
 #              linewidth=1.5, markersize=4, label=None, color=cols[ilab])
 
-    plt.xlabel('aboslute height ($m$)', color='k')
-    plt.ylabel('$\delta height$ ($m$)', color='black')
+    plt.xlabel('absolute height ($m$)', color='k', fontsize=6)
+    plt.ylabel('$\delta height$ ($m$)', color='black', fontsize=6)
     plt.xscale('log')
     plt.yscale('log')
 
-    plt.title('Absolute \& relative heights',position=titleloc)
+    #plt.title('Absolute \& relative heights',position=titleloc, fontsize=8)
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(2,2,3)
-    # Relative heights
+    plt.subplot(4,2,6)
+    # Absolute pressures & relative
     Nlabs = len(labels)
+    nx = 100000.
+    xx = -100000.
     for ilab in range(Nlabs):
-        hgtssea = adhgtsea[ilab]
-        hgtsxhgt = adhgtxhgt[ilab]
-        znw = aetau[ilab]
-
-        plt.plot(znw, hgtssea, '-', marker=marks[ilab], linewidth=1.5, markersize=4, \
-          label='$\delta hgt_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
-        if ilab == 0:
-            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
-              markersize=4, label='$\delta hgt_{hgtmax}^{'+labels[ilab]+'}$',        \
-              color=cols[ilab])
-        else:
-            plt.plot(znw, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,        \
-              markersize=4, label=None, color=cols[ilab])
-    plt.xlabel('$\eta$-level', color='k')
-    plt.ylabel('$\delta height$ ($m$)', color='black')
+        hgtssea = apsea[ilab]/100.
+        hgtsxhgt = apxhgt[ilab]/100.
+        dhgtssea = adpsea[ilab]
+        dhgtsxhgt = adpxhgt[ilab]
+        dz = len(hgtssea)
+
+        nv = np.min(hgtssea[1:dz])
+        if nx > nv: nx = nv
+        nv = np.min(hgtsxhgt[1:dz])
+        if nx > nv: nx = nv
+        xv = np.max(hgtssea[1:dz])
+        if xx < nv: xx = xv
+        xv = np.max(hgtsxhgt[1:dz])
+        if xx < nv: xx = xv
+
+        plt.plot(hgtssea[1:dz], dhgtssea, '-', marker=marks[ilab], linewidth=1.5,    \
+          markersize=4, label='$p_{sea}^{'+labels[ilab]+'}$', color=cols[ilab])
+#        if ilab == 0:
+#            plt.plot(hgtsxhgt[1:dz], dhgtsxhgt, '-.', marker=marks[ilab],            \
+#              linewidth=1.5, markersize=4, label='$hgt_{hgtmax}^{'+labels[ilab]+'}$',\
+#              color=cols[ilab])
+#        else:
+#            plt.plot(hgtsxhgt[1:dz], dhgtsxhgt, '-.', marker=marks[ilab],            \
+#              linewidth=1.5, markersize=4, label=None, color=cols[ilab])
+
+    plt.xlabel('absolute pressure ($hPa$)', color='k', fontsize=6)
+    plt.ylabel('$\delta pressure$ ($Pa$)', color='black', fontsize=6)
+    plt.xlim(xx,nx)
+
+    #plt.xscale('log')
     plt.yscale('log')
 
-    plt.title('Relative Heights',position=titleloc)
+    #plt.title('Absolute \& relative pressures',position=titleloc, fontsize=8)
     #plt.legend(loc=0, prop={'size':10})
 
-    plt.subplot(2,2,4)
+    plt.subplot(4,2,7)
     # Absoulte heights & grid-level
     Nlabs = len(labels)
@@ -12866 +12982 @@
               markersize=0, label=None, color=cols[ilab])
     plt.grid()
-    plt.xlabel('\% z-grid point', color='k')
-    plt.ylabel('height ($m$)', color='black')
+    plt.xlabel('\% z-grid point', color='k', fontsize=6)
+    plt.ylabel('height ($m$)', color='black', fontsize=6)
     #plt.xscale('log')
     plt.yscale('log')
 
-    plt.title('grid point',position=titleloc)
+    #plt.title('grid point',position=titleloc, fontsize=8)
     #plt.legend(loc=0, prop={'size':10})
 
+    plt.subplot(4,2,8)
+    # Absoulte pressures & grid-level
+    Nlabs = len(labels)
+    nx = 100000.
+    xx = -100000.
+    for ilab in range(Nlabs):
+        hgtssea = apsea[ilab]/100.
+        hgtsxhgt = apxhgt[ilab]/100.
+        znw = aetau[ilab]
+        dz = len(hgtssea)
+        zi = np.arange(1.,dz+1.)*100./dz
+
+        nv = np.min(hgtssea[1:dz])
+        if nx > nv: nx = nv
+        nv = np.min(hgtsxhgt[1:dz])
+        if nx > nv: nx = nv
+        xv = np.max(hgtssea[1:dz])
+        if xx < nv: xx = xv
+        xv = np.max(hgtsxhgt[1:dz])
+        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])
+        if ilab == 0:
+            plt.plot(zi, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,         \
+              markersize=0, label='$\delta p_{hgtmax}^{' + labels[ilab] + '}$',    \
+              color=cols[ilab])
+        else:
+            plt.plot(zi, hgtsxhgt,  '-.', marker=marks[ilab], linewidth=1.5,         \
+              markersize=0, label=None, color=cols[ilab])
+    plt.grid()
+    plt.xlabel('\% z-grid point', color='k', fontsize=6)
+    plt.ylabel('pressure ($hPa$)', color='black', fontsize=6)
+
+    plt.ylim(xx,nx)
+    #plt.xscale('log')
+    #plt.yscale('log')
+
+    #plt.title('grid point',position=titleloc, fontsize=8)
+    #plt.legend(loc=0, prop={'size':10})
+
     fig.suptitle(gen.latex_text(figtitle))
-    plt.subplots_adjust(wspace=0.3, hspace=0.4)
+    plt.subplots_adjust(wspace=0.4, hspace=0.6)
 
     figname = 'WRFeta_levels'