<div dir="ltr"><div>Looking at a dump of the '<a href="http://temp.nc">temp.nc</a>'  file, I see that the source file is a "Monthly Mean of Air temperature" file from NOAA's ESRL/PSD archive.<br><br>   <a href="http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.htm">http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.htm</a><br><br>====<br></div><div>I suggest downloading the '<a href="http://hgt.mon.mean.nc">hgt.mon.mean.nc</a>" and creating a file analogous to "<a href="http://temp.nc">temp.nc</a>". Then ..<br><br></div><div>==== untested ======<br><br></div><div>   fh    = addfile("<a href="http://hgt.nc">hgt.nc</a>" , "r")<br></div><div>   ft     = addfile("<a href="http://temp.nc">temp.nc</a>","r")<br></div><div>   hgt  = fh->hgt                      ; (time,level,lat,lon)<br></div><div>   tmp = ft->air                        ;            "<br><br></div><div>   dimh = dimsizes(hgt)<br></div><div>   dimt  = dimsizes(tmp)<br></div><div>   if (.not.all(dimh.eq.dimt)) then<br></div><div>       print("variable size mismatch")<br></div><div>       print(dimh)<br></div><div>       print(dimt)<br></div><div>       exit<br></div><div>   end if<br></div><div>       <br></div><div>   linlog = 0      ; user choice<br>   hlev   = (/ 500, 1500, 9000 /)   ; desired interpolation levels<br>   lev_dim = 1                           ; (time,level,lat,lon) => dim number (0,1,2,3)<br>   thgt      = int2p_n_Wrap(hgt,tmp, hlev, linlog, lev_dim)<br>   printVarSummary(thgt)<br><br>=================<br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/int2p_n.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/int2p_n.shtml</a><br><br>================<br></div><div>There are many examples for plotting this type of array.<br><br></div><div>Good Luck<br></div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 16, 2018 at 10:43 PM, dave montgomery <span dir="ltr"><<a href="mailto:montgomery_dave@yahoo.com" target="_blank">montgomery_dave@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:13px"><div id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_3067">Dear Ncl-talk,</div><div id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_4611" dir="ltr">Is there a way to interpolate or convert as well as to plot a 4D (time,level,lat,lon) variable that is already on pressure levels (hPa) to height levels (m)? In other words, my goal is to plot the data in height levels (m) rather than in hPa. I am aware of the "gsn_csm_pres_hgt_vector" function for plotting, however, it only plots pressure levels on the Y-axis left side which it uses as the main axis for its plot and not the Y-axis right side. I even came upon an example from NCL archive (<a href="http://www.hao.ucar.edu/modeling/tgcm/NCL/Graphics/Scripts/zplat_6.ncl" id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_9393" target="_blank">http://www.hao.ucar.edu/<wbr>modeling/tgcm/NCL/Graphics/<wbr>Scripts/zplat_6.ncl</a>) which I used but to no avail. I got error messages: fatal:Undefined identifier: (halfzp) is undefined, can't continue; fatal:["Execute.c":8640]:<wbr>Execute: Error occurred at or near line 16 in file zplat.ncl. I know I don't have the "halfzp" identifier. I thought it's one of NCL resources. I have attached my script and the ".nc" file for your perusal.<br></div><div dir="ltr" id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_9173">Any assistance on the best way to go would be greatly appreciated. <br></div><div dir="ltr" id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_9651">Many thanks.</div><span class="HOEnZb"><font color="#888888"><div dir="ltr" id="m_6015588489607750481yui_3_16_0_ym19_1_1518843635550_9174">Dave  <br></div></font></span></div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>