<html><head><meta http-equiv="content-type" content="text/html; charset=GB2312"><style>body { line-height: 1.5; }body { font-size: 10.5pt; font-family: ΢ÈíÑźÚ; color: rgb(0, 0, 0); line-height: 1.5; }body { font-size: 10.5pt; font-family: ΢ÈíÑźÚ; color: rgb(0, 0, 0); line-height: 1.5; }</style></head><body>
<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/bbS866C200000CPfb8.cam2.h0.1901-1999.avg.nc"<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;">;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;">, </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></body></html>