<div dir="ltr"><div dir="ltr"><div dir="ltr"><a href="http://www.ncl.ucar.edu/Applications/pot_vort.shtml"><b>http://www.ncl.ucar.edu/Applications/pot_vort.shtml</b></a></div><div dir="ltr"><br></div><div>[1]</div><div>"Previously, I computed PV using data from ERA-Interim on a <span style="color:rgb(0,0,255)">gaussian grid </span><b>without any troubles</b>."</div><div><br></div><div>[2] <br></div><div>"Now I try the same, but with ERA-Interim data on a<b> regular 1x1 grid.</b></div><div><b><br></b></div><div><b>---</b></div><div><b><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/pot_vort_hybrid.shtml">pot_vort_hybrid</a></b>: This function works with both gaussian and regular grids. Same code, just a switch for the grid type<br></div><div><br></div><div>---</div><div>I <b>speculate</b> that<b> </b>the ECMWFspectral  interpolation from the source gaussian grid to 1x1 is the issue. <b>Spectral interpolation can result in "ringing" due to Gibbs phenomenon.</b></div><div><b><br></b></div><div><b>D</b><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><b></b> </div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 6, 2019 at 10:31 AM De Vries  Andries <<a href="mailto:andries.devries@env.ethz.ch">andries.devries@env.ethz.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear NCL team & NCL users,<br>
<br>
I have an issue with the computation of potential vorticity (PV) on hybrid levels using the function "pot_vort_hybrid". Previously, I computed PV using data from ERA-Interim on a gaussian grid without any troubles. Now I try the same, but with ERA-Interim data on a regular 1x1 grid. I found that values toward the pole are fluctuating with respect to PV that has been provided directly by the ECMWF. (The reason for computing PV myself is to obtain PV isentropic surfaces with higher intervals than provided by ECMWF). In other words, the PV that I compute show increasing, positive & negative alternating values as compared to the PV that is provided by the ECMWF. <br>
<br>
Any idea if I do something wrong, and/or how to fix this issue?<br>
<br>
Please, find attached the two scripts for the computation of PV on a Gaussian grid and on a regular grid, as well as a short script that prints the numerical differences between the two output files. The needed data files to compute PV on a guassian and regular grid, as well as files with PV from ECMWF are uploaded.<br>
<br>
Looking forward for any feedback!<br>
Many thanks!<br>
<br>
Cheers,<br>
Andries<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>