<html><head><meta http-equiv="content-type" content="text/html; charset=UTF-8"><style>body { line-height: 1.5; }blockquote { margin-top: 0px; margin-bottom: 0px; margin-left: 0.5em; }div.foxdiv20150519002546125606 { }body { font-size: 10.5pt; font-family: 微软雅黑; color: rgb(0, 0, 0); line-height: 1.5; }</style></head><body>
<div><span></span><div><span style="background-color: rgba(0, 0, 0, 0);">Hello,</span></div><div><span style="background-color: rgba(0, 0, 0, 0);"><br></span></div><div><span style="background-color: rgba(0, 0, 0, 0);">The netCDF file （</span><span style="font-size: 10.5pt; line-height: 1.5; background-color: window;">bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc）</span><span style="font-size: 10.5pt; line-height: 1.5; background-color: rgba(0, 0, 0, 0);"> has been uploaded to </span><a href="http://ftp.cgd.ucar.edu" target="_blank" style="font-size: 10.5pt; line-height: 1.5; background-color: window;">ftp.cgd.ucar.edu</a><span style="font-size: 10.5pt; line-height: 1.5; background-color: window;">/incoming </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;">Thank you.</span></div><div><br></div><div><br></div><div>Regards,</div><div>Wade</div><blockquote> </blockquote></div>
<div><br></div><div><br></div>
<blockquote style="margin-top: 0px; margin-bottom: 0px; margin-left: 0.5em;"><div> </div><div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"><div style="PADDING-RIGHT: 8px; PADDING-LEFT: 8px; FONT-SIZE: 12px;FONT-FAMILY:tahoma;COLOR:#000000; BACKGROUND: #efefef; PADDING-BOTTOM: 8px; PADDING-TOP: 8px"><div><b>From:</b> <a href="mailto:shea@ucar.edu">Dennis Shea</a></div><div><b>Date:</b> 2015-05-18 23:36</div><div><b>To:</b> <a href="mailto:wadew@pku.edu.cn">Wade Qiang Wei</a></div><div><b>CC:</b> <a href="mailto:ncl-talk@ucar.edu">ncl-talk</a></div><div><b>Subject:</b> Re: [ncl-talk] Strange output with vinth2p when P0>1000hPa</div></div></div><div><div class="FoxDiv20150519002546125606"><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: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: 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>
</div></div></blockquote>
</body></html>