<div>&nbsp; &nbsp;&nbsp;</div><div>Thanks Karin,</div><div>&nbsp; &nbsp;It's not the whole script, I choose the important part to show you. I just check VarSummary,&nbsp;</div><div><span style="line-height: 1.5;">&nbsp; &nbsp;</span><span style="line-height: 1.5;">level variable is listed below:</span></div><div>&nbsp; &nbsp;standard_name : &nbsp;air_pressure</div><div>&nbsp; &nbsp;long_name : &nbsp; vertical level</div><div>&nbsp; &nbsp;units : &nbsp; &nbsp; &nbsp; hPa</div><div>&nbsp; &nbsp;positive : &nbsp; &nbsp;down</div><div>&nbsp; &nbsp;axis : &nbsp; &nbsp; &nbsp; &nbsp;Z</div><div>(0) &nbsp; &nbsp; &nbsp;200</div><div>(1) &nbsp; &nbsp; &nbsp;150</div><div>(2) &nbsp; &nbsp; &nbsp;100</div><div>(3) &nbsp; &nbsp; &nbsp; 70</div><div>(4) &nbsp; &nbsp; &nbsp; 50</div><div>(5) &nbsp; &nbsp; &nbsp; 40</div><div>(6) &nbsp; &nbsp; &nbsp; 30</div><div>(7) &nbsp; &nbsp; &nbsp; 20</div><div>(8) &nbsp; &nbsp; &nbsp; 10</div><div>(9) &nbsp; &nbsp; &nbsp; &nbsp;7</div><div>(10) &nbsp; &nbsp; &nbsp; 5</div><div>(11) &nbsp; &nbsp; &nbsp; 4</div><div>(12) &nbsp; &nbsp; &nbsp; 3</div><div>(13) &nbsp; &nbsp; &nbsp; 2</div><div>(14) &nbsp; &nbsp; &nbsp; 1</div><div><br></div><div>&nbsp; &nbsp; and variable &nbsp;merra's VarSummary is&nbsp;<span style="line-height: 1.5;">listed below:</span></div><div>&nbsp; &nbsp; Dimensions and sizes: &nbsp; [lev | 15] x [lat | 73]</div><div>&nbsp; &nbsp; Coordinates:</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; lev: [ 200.. &nbsp; 1]</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; lat: [ &nbsp;90.. -90]</div><div><br></div><div>But I still don't know how to show level axis above 10hPa&nbsp;</div><div>via ¡°gsn_csm_pres_hgt¡±&nbsp;to show Pressure log coordinate system.</div><div><br></div><div>Yan</div><div><br>Hi,<br><br>if this is your full script comment the lines res@gsnDraw and <br>res@gsnFrame because you wouldn't get any plot.<br><br>Please, do a printVarSummary of your level variable (e.g. lev or level). <br>Is the unit attribute "hPa"? Do also a print of the level variable to <br>see which values are stored.<br><br>You don't have to add anything your variable should have all <br>information. Do also a printVarSummary(merra).<br><br>Bye,<br>Karin<br><br><br>Am 26.04.17 um 04:41 schrieb WIND:<br>&gt; Hi all,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; I'm trying to draw the pressure v.s. latitude contour via <br>&gt;&nbsp; "gsn_csm_pres_hgt", but the pressure axis just don't show pressure <br>&gt; data above 10hPa (from 10hPa to 1hPa). (Besides, I want to emphasize <br>&gt; the upper stratosphere data.)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; How should I add the whole pressure data?<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; Here is my script<br>&gt;<br>&gt; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>&gt; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>&gt; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>&gt; fq=addfile("/merra2/200-1/MERRA2.nc4","r")<br>&gt; merra=fq-&gt;T<br>&gt; wks&nbsp; = gsn_open_wks("eps" ,"t")<br>&gt; res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>&gt; res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>&gt; res@tmYRMode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "Automatic"<br>&gt; res@gsnDraw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>&gt; res@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>&gt; res@tmXBMinorOn&nbsp;&nbsp;&nbsp;&nbsp; = False<br>&gt; res@tiYAxisString="Pressure"<br>&gt; plot = gsn_csm_pres_hgt(wks,merra,res)<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; Any help would be grateful!<br>&gt;<br>&gt; Yan<br><br><br></div>