<div dir="ltr"><div><div><div><div><div><div><div><div><div>Hello,<br><br>There is nothing in the internal code that precludes interpolation to values > 1000hPa.<br><br></div>Please make the file available.<br><br></div>ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><br></div>anonymous<br></div>email<br></div>cd incoming<br></div><div>prompt<br></div>put <span style="color:rgb(0,0,0);background-color:transparent"><a href="http://bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc" target="_blank">bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc</a><br></span></div><span style="color:rgb(0,0,0);background-color:transparent">quit<br><br></span></div><span style="color:rgb(0,0,0);background-color:transparent">Send an email to ncl-talk after the file has been successfully uploaded,<br><br></span></div><span style="color:rgb(0,0,0);background-color:transparent">Cheers<br></span></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, May 17, 2015 at 10:23 PM, Wade Qiang Wei <span dir="ltr"><<a href="mailto:wadew@pku.edu.cn" target="_blank">wadew@pku.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<div><span></span>Hi all,</div><div><br></div><div>I am trying to analyze cam model results in a CCSM3.0 simulation of an exoplanet with P0=1379hPa. It kept giving me weird data after converting from hybrid coordinates to pressure coordinates with vinth2p. So I turned off extrapolation and noticed every output with pnew>=1000hPa is filled with missing values. The original data has levels across 3hpa to 1379hpa and therefore extrapolation should not be necessary from my point of view.</div><div><br></div><div>Could vinth2p deal with <span style="font-size:10.5pt;line-height:1.5;background-color:window">P0 greater than 1000hPa? Has anyone encountered similar problems?</span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><br></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window">Thanks a lot.</span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><br></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window"><br></span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window">Best Rards,</span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window">Wade</span></div>
<div><br></div><div>------------------------------------------------------------------------------</div><span style="color:rgb(0,0,0);background-color:rgba(0,0,0,0)">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br><br>begin<br><br> cam4="~/GI581g/<a href="http://bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc" target="_blank">bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc</a>"<br>;---Read data<br> a = addfile(cam4,"r")<br> T = a->T<br><br></span><span style="color:rgb(0,0,0);background-color:rgba(0,0,0,0)"><span style="white-space:pre-wrap">;original T&lev=xxxxxxxxxxxx</span></span><span style="font-size:10.5pt;line-height:1.5;background-color:window">, </span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">960, 1</span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">086., </span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">1195, </span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">1281<span style="white-space:pre-wrap">, </span></span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">1338, </span><span style="color:rgb(0,0,0);font-size:10.5pt;line-height:1.5;background-color:rgba(0,0,0,0)">1368</span><div><span style="background-color:rgba(0,0,0,0);font-size:10.5pt;line-height:1.5"><br></span></div><div><span style="color:rgb(0,0,0);background-color:rgba(0,0,0,0)"> pnew = (/3,10,20,30,50,70,100,150,200,250,300,400,500,600,700,850,1000,1150,1250,1330,1379/)<br> pnew@units = "hPa" <br> P0mb = 1379 ;mb<br> K = 0.286<br> g=13.5<br><br><br>;===============================================<br> hyam = a->hyam ; get a coefficiants<br> hybm = a->hybm ; get b coefficiants<br> PS = a->PS ; get pressure in Pa<br><br>;************************************************<br> interp = 2 <br>; is extrapolation desired if data is outside the range of PS<br> extrap = False<br>;************************************************<br> T2 = vinth2p(T,hyam,hybm,pnew,PS,interp,P0mb,1,extrap)<br><br> Tout=dim_avg_n_Wrap(T2 (0,:, :,:), 2)<br><br> wks_type="x11"<br> wks = gsn_open_wks(wks_type,"Zonal_Mean_T_Slice"+"_top")<br><br>; -- set resources<br> res = True<br> res@cnFillOn = True ; -- turn on color fill<br> res@cnLineLabelsOn = False ; -- turns off contour line labels<br> res@cnInfoLabelOn = False ; -- turns off contour info label<br> res@cnLinesOn = False<br> res@tiYAxisString = "Pressure"+" [hPa]" ; -- append units to y-axis label<br> res@trYReverse = True ; -- reverses y-axis<br> res@gsnMaximize=True<br><br> plot = gsn_csm_contour(wks,Tout,res)<br><br>end</span>
</div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>