<div dir="ltr"><div><div><div><div><div><div><div>Hi Kacie,<br><br></div>I see the issue.<br><br></div>As noted in the documentation:<br> <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/pot_vort_isobaric.shtml">http://www.ncl.ucar.edu/Document/Functions/Contributed/pot_vort_isobaric.shtml</a><br><br> "Compute potential vorticity on isobaric levels and a global rectilinear grid."<br>----<br></div>After you regrid, you have a rectilinear grid. However, it is not "global." EG:<br><br>Variable: tkReGr
<p class="MsoNormal">Number of Dimensions: 3</p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]</p>
<p class="MsoNormal">Coordinates:</p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]</p>
<p class="MsoNormal"> lat: [ 1.. 85] <==== Regional grid<br></p>
lon: [ 150..358.5] <====<br><br></div>The 'pot_vort_isobaric' function uses spherical harmonic functions to calculate the required meridional and zonal gradients. The spherical harmonic functions are "highly accurate" however they require global grids.<br><br>===<br></div>As I recall, in the summer NCL Workshop, we discussed this issue. :-)<br><br></div>As an exercise, I suggested that you copy and modify 'pot_vort_isobaric' to use the 6.4.0 'grad_latlon_cfd' function to replace the spherical harmonic gradient functions. <br><br> <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/grad_latlon_cfd.shtml">http://www.ncl.ucar.edu/Document/Functions/Contributed/grad_latlon_cfd.shtml</a><br><br></div>Some gradient examples are at: <a href="http://www.ncl.ucar.edu/Applications/gradients.shtml">http://www.ncl.ucar.edu/Applications/gradients.shtml</a><br><div><br></div><div>Actually, in your script, you have a commented line:<br> ;load "/grid-scratch/kshourd/code/potVortCFD.ncl"<br><br></div><div>So .... it looks like you have a created a function but you did not use it??? <br><br>===<br></div><div>I decided to create a function 'pot_vort_isobaric_cfd'. It has slightly different arguments.<br></div><div>See attached.<br><br></div><div>Likely, I will include 'pot_vort_isobaric_cfd' in the 6.4.0 release as part of contributed.ncl<br></div><div><br></div><div>D<br></div><div><div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 7, 2016 at 12:31 PM, Kacie Shourd <span dir="ltr"><<a href="mailto:Kacie.Shourd@dri.edu" target="_blank">Kacie.Shourd@dri.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="#0563C1" vlink="#954F72" lang="EN-US">
<div>
<p class="MsoNormal">Hi all,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am encountering some problems with my potential vorticity (PV) and subsequent isentropic potential vorticity calculations (IPV). The patterns I am getting look right, but the magnitude of the values does not make sense. I assume pot_vort_isobaric()
function outputs in raw units rather than PVUs? <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have checked my units/script over and over and over again, and cannot figure out why the values are a few orders of magnitude off. It seems as though maybe my pressure values are being input into the equations wrong (or some other variable
is). The patterns of PV make sense for the event, but the magnitude is most certainly wrong.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Any insight into what I am doing wrong or any other reasons the calculated PV values are coming out really large. Scripts, output, and necessary files to run script can be found here:
<a href="ftp://pubfiles.dri.edu/pub/kshourd/nclTalkIPVHelp/" target="_blank">ftp://pubfiles.dri.edu/pub/<wbr>kshourd/nclTalkIPVHelp/</a> , and variable summaries are as follows:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: tkReGr<u></u><u></u></p>
<p class="MsoNormal">Type: float<u></u><u></u></p>
<p class="MsoNormal">Total Size: 32485452 bytes<u></u><u></u></p>
<p class="MsoNormal"> 8121363 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 3<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<u></u><u></u></p>
<p class="MsoNormal"> lat: [ 1.. 85]<u></u><u></u></p>
<p class="MsoNormal"> lon: [ 150..358.5]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 15<u></u><u></u></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<u></u><u></u></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<u></u><u></u></p>
<p class="MsoNormal"> long_name : Temperature<u></u><u></u></p>
<p class="MsoNormal"> units : K<u></u><u></u></p>
<p class="MsoNormal"> _FillValue : 1e+20<u></u><u></u></p>
<p class="MsoNormal"> level_indicator : 100<u></u><u></u></p>
<p class="MsoNormal"> grid_number : 221<u></u><u></u></p>
<p class="MsoNormal"> parameter_table_version : 131<u></u><u></u></p>
<p class="MsoNormal"> parameter_number : 11<u></u><u></u></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<u></u><u></u></p>
<p class="MsoNormal"> forecast_time : 0<u></u><u></u></p>
<p class="MsoNormal"> forecast_time_units : hours<u></u><u></u></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<u></u><u></u></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<u></u><u></u></p>
<p class="MsoNormal"> missing_value : 1e+20<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: uN<u></u><u></u></p>
<p class="MsoNormal">Type: float<u></u><u></u></p>
<p class="MsoNormal">Total Size: 32485452 bytes<u></u><u></u></p>
<p class="MsoNormal"> 8121363 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 3<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<u></u><u></u></p>
<p class="MsoNormal"> lat: [ 1.. 85]<u></u><u></u></p>
<p class="MsoNormal"> lon: [ 150..358.5]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 15<u></u><u></u></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<u></u><u></u></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<u></u><u></u></p>
<p class="MsoNormal"> long_name : u-component of wind<u></u><u></u></p>
<p class="MsoNormal"> units : m/s<u></u><u></u></p>
<p class="MsoNormal"> _FillValue : 1e+20<u></u><u></u></p>
<p class="MsoNormal"> level_indicator : 100<u></u><u></u></p>
<p class="MsoNormal"> grid_number : 221<u></u><u></u></p>
<p class="MsoNormal"> parameter_table_version : 131<u></u><u></u></p>
<p class="MsoNormal"> parameter_number : 33<u></u><u></u></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<u></u><u></u></p>
<p class="MsoNormal"> forecast_time : 0<u></u><u></u></p>
<p class="MsoNormal"> forecast_time_units : hours<u></u><u></u></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<u></u><u></u></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<u></u><u></u></p>
<p class="MsoNormal"> missing_value : 1e+20<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: vN<u></u><u></u></p>
<p class="MsoNormal">Type: float<u></u><u></u></p>
<p class="MsoNormal">Total Size: 32485452 bytes<u></u><u></u></p>
<p class="MsoNormal"> 8121363 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 3<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<u></u><u></u></p>
<p class="MsoNormal"> lat: [ 1.. 85]<u></u><u></u></p>
<p class="MsoNormal"> lon: [ 150..358.5]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 15<u></u><u></u></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<u></u><u></u></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<u></u><u></u></p>
<p class="MsoNormal"> long_name : v-component of wind<u></u><u></u></p>
<p class="MsoNormal"> units : m/s<u></u><u></u></p>
<p class="MsoNormal"> _FillValue : 1e+20<u></u><u></u></p>
<p class="MsoNormal"> level_indicator : 100<u></u><u></u></p>
<p class="MsoNormal"> grid_number : 221<u></u><u></u></p>
<p class="MsoNormal"> parameter_table_version : 131<u></u><u></u></p>
<p class="MsoNormal"> parameter_number : 34<u></u><u></u></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<u></u><u></u></p>
<p class="MsoNormal"> forecast_time : 0<u></u><u></u></p>
<p class="MsoNormal"> forecast_time_units : hours<u></u><u></u></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<u></u><u></u></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<u></u><u></u></p>
<p class="MsoNormal"> missing_value : 1e+20<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: p<u></u><u></u></p>
<p class="MsoNormal">Type: integer<u></u><u></u></p>
<p class="MsoNormal">Total Size: 116 bytes<u></u><u></u></p>
<p class="MsoNormal"> 29 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 1<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 2<u></u><u></u></p>
<p class="MsoNormal"> long_name : isobaric level<u></u><u></u></p>
<p class="MsoNormal"> units : Pa<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: thetaN<u></u><u></u></p>
<p class="MsoNormal">Type: float<u></u><u></u></p>
<p class="MsoNormal">Total Size: 32485452 bytes<u></u><u></u></p>
<p class="MsoNormal"> 8121363 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 3<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<u></u><u></u></p>
<p class="MsoNormal"> lat: [ 1.. 85]<u></u><u></u></p>
<p class="MsoNormal"> lon: [ 150..358.5]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 5<u></u><u></u></p>
<p class="MsoNormal"> _FillValue : 1e+20<u></u><u></u></p>
<p class="MsoNormal"> long_name : potential temperature<u></u><u></u></p>
<p class="MsoNormal"> units : K<u></u><u></u></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<u></u><u></u></p>
<p class="MsoNormal"> missing_value : 1e+20<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">--- (FROM tkReGr&lv_ISBL3)<u></u><u></u></p>
<p class="MsoNormal">Variable: p<u></u><u></u></p>
<p class="MsoNormal">Type: integer<u></u><u></u></p>
<p class="MsoNormal">Total Size: 116 bytes<u></u><u></u></p>
<p class="MsoNormal"> 29 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 1<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 2<u></u><u></u></p>
<p class="MsoNormal"> long_name : isobaric level<u></u><u></u></p>
<p class="MsoNormal"> units : Pa<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal">Variable: ipv<u></u><u></u></p>
<p class="MsoNormal">Type: float<u></u><u></u></p>
<p class="MsoNormal">Total Size: 4480752 bytes<u></u><u></u></p>
<p class="MsoNormal"> 1120188 values<u></u><u></u></p>
<p class="MsoNormal">Number of Dimensions: 3<u></u><u></u></p>
<p class="MsoNormal">Dimensions and sizes: [thlev | 4] x [lat | 337] x [lon | 831]<u></u><u></u></p>
<p class="MsoNormal">Coordinates:<u></u><u></u></p>
<p class="MsoNormal"> thlev: [325..370]<u></u><u></u></p>
<p class="MsoNormal"> lat: [ 1.. 85]<u></u><u></u></p>
<p class="MsoNormal"> lon: [ 150..358.5]<u></u><u></u></p>
<p class="MsoNormal">Number Of Attributes: 4<u></u><u></u></p>
<p class="MsoNormal"> _FillValue : 1e+20<u></u><u></u></p>
<p class="MsoNormal"> long_name : Isentropic Potential Vorticity at 325K<u></u><u></u></p>
<p class="MsoNormal"> short_name : IPV<u></u><u></u></p>
<p class="MsoNormal"> units : K m2/kg/s<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks!!<u></u><u></u></p>
<p class="MsoNormal">Best,<u></u><u></u></p>
<p class="MsoNormal">Kacie<u></u><u></u></p>
<p class="MsoNormal">---<u></u><u></u></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:black">Kacie N. Shourd</span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Malgun Gothic",sans-serif;color:black">드림</span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u> <u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:#222222">Graduate Research Assistant</span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:#222222">Desert Research Institute</span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:#222222"><a href="mailto:Kacie.Shourd@dri.edu" target="_blank"><span style="color:blue">Kacie.Shourd@dri.edu</span></a></span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:#222222">Work: <a href="tel:%28775%29%20674-7049" value="+17756747049" target="_blank">(775) 674-7049</a></span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<br>
<font face="Arial" size="2">PUBLIC RECORDS NOTICE: In accordance with NRS Chapter 239, this email and responses, unless otherwise made confidential by law, may be subject to the Nevada Public Records laws and may be disclosed to the public upon request.</font>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>