<div dir="ltr"><div><div><div><div>Hi Dennis,<br><br></div>Thanks for your reply. <br></div>I do hope this might be a common issue in the community so that I can get a quick answer. <br><br></div><div>I attached the file you asked (<a href="http://test_hyi2hyo.nc">test_hyi2hyo.nc</a>) and my script (check_ERAint_vert_4.ncl). <br></div><div>Also the files used in my test script. <br><br></div><div>Note that for the test_hyi2hyo I don't know how to write a nc file with different vertical resolution (ERAI vs. CAM).<br></div><div>So for hya_cam and hyb_cam, only (34:59) contains useful information. <br><br></div><div>Thanks for your time again, and see if we can identify the problem.<br> </div></div><div><div><div><br><br></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><span><div><div dir="ltr"><div align="left"><font>
<div align="left"><font><span>
<div><font><font><span><br>Cheers, <br></span></font></font></div></span></font></div>
<div><font><span>Mao-Lin</span></font></div></font></div>=================================<br>Dr. Mao-Lin Shen<br>Researcher<br>Geophysical Institute<br>Faculty of Mathematics and Natural Sciences<br>University of Bergen<br>Allégaten 70, 5007 Bergen, Norway<br>Tel: (Office) +47 <span>55 58 46 54</span><br>Tel: (mobile) +47 94 42 54 55<br><span>E-mail: <a href="mailto:earnestshen@gmail.com" target="_blank">earnestshen@gmail.com</a><br> <a href="mailto:Maolin.Shen@gfi.uib.no" target="_blank">Maolin.Shen@gfi.uib.no</a><br>Homepage</span>: <a href="http://www.uib.no/persons/Maolin.Shen" target="_blank">http://www.uib.no/persons/Maolin.Shen</a><br> <a href="http://folk.uib.no/msh065/" target="_blank">http://folk.uib.no/msh065/</a><br>=================================</div></div></span></div></div></div>
<br><div class="gmail_quote">2015-03-23 23:25 GMT+01:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div>The email does not provide enough information for debugging.<br>I would suggest the following to facilitate debugging.<br><br><br></div>Since "horizontal resolution remains identical in this test", <br>pick one grid point<br><br></div><div>The reason for using {LAT:LAT},{LON:LON} is to force NCL to retain dimensionality.<br></div><div> x(time,lev,lat,lon)<br></div><div>If you use {LAT},{LON}, NCL will eliminate the 'degenerate' dimension<br></div><div> x(time,lev)<br>-----------<br> </div><div><br></div><div><br></div> LAT = 30.0 ; wherever<br></div> LON = 180.0<br><br></div> x = f->ERAI(:,:,{LAT:LAT},{LON:LON}) ; force NCL to retain dimensionality<br></div> printVarSummary(x) ; x(time,level,1,1)<br><br></div> ps = f->PSFC(:,{LAT:LAT},{LON:LON}) ; force NCL to retain dimensionality<div> printVarSummary(ps)<br></div><div><pre> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml" target="_blank"><b>system</b></a>("/bin/rm -f <a href="http://test.nc" target="_blank">test_hyi2hyo.nc</a>") ; remove any pre-existing file
ncdf = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" target="_blank"><b>addfile</b></a>("<a href="http://test.nc" target="_blank">test_hyi2hyo.nc</a>" ,"c") ; open output netCDF file
;===================================================================
; create global attributes of the file (optional)
;===================================================================
fAtt = True ; assign file attributes
fAtt@title = "TEST NCL HYI2HYO: ERA-INTERIM grid point"
fAtt@source_file = "<a href="http://original-file.nc" target="_blank">original-file.nc</a>"
fAtt@Conventions = "None"
fAtt@creation_date = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/systemfunc.shtml" target="_blank"><b>systemfunc</b></a> ("date")
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/fileattdef.shtml" target="_blank"><b>fileattdef</b></a>( ncdf, fAtt ) ; copy file attributes
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml" target="_blank"><b>filedimdef</b></a>( ncdf,"time",-1,True)
;===================================================================
; output variables directly; NCL will call appropriate functions
; to write the meta data associated with each variable
;===================================================================<br> ncdf->hya_erai = hya_erai ; 1D
ncdf->hyb_erai = hyb_erai ; 1D<br> ncdf->hya_cam = hya_cam ; 1D
ncdf->hyb_cam = hyb_cam ; 1D <br><br> ncdf->X_ERAI = x ; 4D <br> ncdf->PS_ERAI = ps ; 3D
<br>===<br></pre><pre>You can send the file and a *clean* simple version of your test script. <br></pre> <br><br> <br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 23, 2015 at 2:55 PM, Mao-Lin Shen <span dir="ltr"><<a href="mailto:maolin.shen@gfi.uib.no" target="_blank">maolin.shen@gfi.uib.no</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<br><div class="gmail_quote"><div dir="ltr"><div class="gmail_quote"><div dir="ltr"><br><div><div><div><br></div>I interpolated data from hybrid ERA-Interim to CAM's hybrid level for nudging.<br></div>But we found a big difference in between when we compared the data on pressure level. <br><br></div>The results show in bellow. <br><div><br>Left
panel is the results applying hyi2hyo to get temperature on CAM's
hybrid level first and then using vinth2p to get temperature on pressure
level. <br><br></div><div>Right panel shows the results interpolated from ERA-Interim by using vinth2p.<br></div><div><br></div><div>Level Temperature (Degree) Temperature (Degree)</div><div>(850hPa) 294.0318218618177 292.5374525210495</div><div>(500hPa) 276.7536469179544 268.8364571841531</div><div>(250hPa) 256.9609800328714 233.3842838572815</div><div>(100hPa) 219.1232631977532 186.2120804720032</div><div><br></div><div>Horizontal resolution remains identical in this test.<br></div><div>The temperature are supposed to be identical, but it is not due to unclear issue. <br></div><div>Does any one have similar problem?<br></div>Or has idea how to tackle this issue?</div></div></div>
</div>
<br>
<br></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>
</blockquote></div><br></div>