<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=ks_c_5601-1987">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:"\@Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi all,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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? <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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/">ftp://pubfiles.dri.edu/pub/kshourd/nclTalkIPVHelp/</a> , and variable summaries are as follows:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: tkReGr<o:p></o:p></p>
<p class="MsoNormal">Type: float<o:p></o:p></p>
<p class="MsoNormal">Total Size: 32485452 bytes<o:p></o:p></p>
<p class="MsoNormal"> 8121363 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 3<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<o:p></o:p></p>
<p class="MsoNormal"> lat: [ 1.. 85]<o:p></o:p></p>
<p class="MsoNormal"> lon: [ 150..358.5]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 15<o:p></o:p></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<o:p></o:p></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<o:p></o:p></p>
<p class="MsoNormal"> long_name : Temperature<o:p></o:p></p>
<p class="MsoNormal"> units : K<o:p></o:p></p>
<p class="MsoNormal"> _FillValue : 1e+20<o:p></o:p></p>
<p class="MsoNormal"> level_indicator : 100<o:p></o:p></p>
<p class="MsoNormal"> grid_number : 221<o:p></o:p></p>
<p class="MsoNormal"> parameter_table_version : 131<o:p></o:p></p>
<p class="MsoNormal"> parameter_number : 11<o:p></o:p></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<o:p></o:p></p>
<p class="MsoNormal"> forecast_time : 0<o:p></o:p></p>
<p class="MsoNormal"> forecast_time_units : hours<o:p></o:p></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<o:p></o:p></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<o:p></o:p></p>
<p class="MsoNormal"> missing_value : 1e+20<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: uN<o:p></o:p></p>
<p class="MsoNormal">Type: float<o:p></o:p></p>
<p class="MsoNormal">Total Size: 32485452 bytes<o:p></o:p></p>
<p class="MsoNormal"> 8121363 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 3<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<o:p></o:p></p>
<p class="MsoNormal"> lat: [ 1.. 85]<o:p></o:p></p>
<p class="MsoNormal"> lon: [ 150..358.5]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 15<o:p></o:p></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<o:p></o:p></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<o:p></o:p></p>
<p class="MsoNormal"> long_name : u-component of wind<o:p></o:p></p>
<p class="MsoNormal"> units : m/s<o:p></o:p></p>
<p class="MsoNormal"> _FillValue : 1e+20<o:p></o:p></p>
<p class="MsoNormal"> level_indicator : 100<o:p></o:p></p>
<p class="MsoNormal"> grid_number : 221<o:p></o:p></p>
<p class="MsoNormal"> parameter_table_version : 131<o:p></o:p></p>
<p class="MsoNormal"> parameter_number : 33<o:p></o:p></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<o:p></o:p></p>
<p class="MsoNormal"> forecast_time : 0<o:p></o:p></p>
<p class="MsoNormal"> forecast_time_units : hours<o:p></o:p></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<o:p></o:p></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<o:p></o:p></p>
<p class="MsoNormal"> missing_value : 1e+20<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: vN<o:p></o:p></p>
<p class="MsoNormal">Type: float<o:p></o:p></p>
<p class="MsoNormal">Total Size: 32485452 bytes<o:p></o:p></p>
<p class="MsoNormal"> 8121363 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 3<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<o:p></o:p></p>
<p class="MsoNormal"> lat: [ 1.. 85]<o:p></o:p></p>
<p class="MsoNormal"> lon: [ 150..358.5]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 15<o:p></o:p></p>
<p class="MsoNormal"> sub_center : The North American Regional Reanalysis (NARR) Project<o:p></o:p></p>
<p class="MsoNormal"> center : US National Weather Service - NCEP (WMC)<o:p></o:p></p>
<p class="MsoNormal"> long_name : v-component of wind<o:p></o:p></p>
<p class="MsoNormal"> units : m/s<o:p></o:p></p>
<p class="MsoNormal"> _FillValue : 1e+20<o:p></o:p></p>
<p class="MsoNormal"> level_indicator : 100<o:p></o:p></p>
<p class="MsoNormal"> grid_number : 221<o:p></o:p></p>
<p class="MsoNormal"> parameter_table_version : 131<o:p></o:p></p>
<p class="MsoNormal"> parameter_number : 34<o:p></o:p></p>
<p class="MsoNormal"> model : North American Regional Reanalysis (NARR)<o:p></o:p></p>
<p class="MsoNormal"> forecast_time : 0<o:p></o:p></p>
<p class="MsoNormal"> forecast_time_units : hours<o:p></o:p></p>
<p class="MsoNormal"> initial_time : 09/26/2010 (18:00)<o:p></o:p></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<o:p></o:p></p>
<p class="MsoNormal"> missing_value : 1e+20<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: p<o:p></o:p></p>
<p class="MsoNormal">Type: integer<o:p></o:p></p>
<p class="MsoNormal">Total Size: 116 bytes<o:p></o:p></p>
<p class="MsoNormal"> 29 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 1<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 2<o:p></o:p></p>
<p class="MsoNormal"> long_name : isobaric level<o:p></o:p></p>
<p class="MsoNormal"> units : Pa<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: thetaN<o:p></o:p></p>
<p class="MsoNormal">Type: float<o:p></o:p></p>
<p class="MsoNormal">Total Size: 32485452 bytes<o:p></o:p></p>
<p class="MsoNormal"> 8121363 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 3<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29] x [lat | 337] x [lon | 831]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> lv_ISBL3: [100..1000]<o:p></o:p></p>
<p class="MsoNormal"> lat: [ 1.. 85]<o:p></o:p></p>
<p class="MsoNormal"> lon: [ 150..358.5]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 5<o:p></o:p></p>
<p class="MsoNormal"> _FillValue : 1e+20<o:p></o:p></p>
<p class="MsoNormal"> long_name : potential temperature<o:p></o:p></p>
<p class="MsoNormal"> units : K<o:p></o:p></p>
<p class="MsoNormal"> remap : remapped via ESMF_regrid_with_weights: Bilinear<o:p></o:p></p>
<p class="MsoNormal"> missing_value : 1e+20<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">--- (FROM tkReGr&lv_ISBL3)<o:p></o:p></p>
<p class="MsoNormal">Variable: p<o:p></o:p></p>
<p class="MsoNormal">Type: integer<o:p></o:p></p>
<p class="MsoNormal">Total Size: 116 bytes<o:p></o:p></p>
<p class="MsoNormal"> 29 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 1<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [lv_ISBL3 | 29]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 2<o:p></o:p></p>
<p class="MsoNormal"> long_name : isobaric level<o:p></o:p></p>
<p class="MsoNormal"> units : Pa<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">---<o:p></o:p></p>
<p class="MsoNormal">Variable: ipv<o:p></o:p></p>
<p class="MsoNormal">Type: float<o:p></o:p></p>
<p class="MsoNormal">Total Size: 4480752 bytes<o:p></o:p></p>
<p class="MsoNormal"> 1120188 values<o:p></o:p></p>
<p class="MsoNormal">Number of Dimensions: 3<o:p></o:p></p>
<p class="MsoNormal">Dimensions and sizes: [thlev | 4] x [lat | 337] x [lon | 831]<o:p></o:p></p>
<p class="MsoNormal">Coordinates:<o:p></o:p></p>
<p class="MsoNormal"> thlev: [325..370]<o:p></o:p></p>
<p class="MsoNormal"> lat: [ 1.. 85]<o:p></o:p></p>
<p class="MsoNormal"> lon: [ 150..358.5]<o:p></o:p></p>
<p class="MsoNormal">Number Of Attributes: 4<o:p></o:p></p>
<p class="MsoNormal"> _FillValue : 1e+20<o:p></o:p></p>
<p class="MsoNormal"> long_name : Isentropic Potential Vorticity at 325K<o:p></o:p></p>
<p class="MsoNormal"> short_name : IPV<o:p></o:p></p>
<p class="MsoNormal"> units : K m2/kg/s<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks!!<o:p></o:p></p>
<p class="MsoNormal">Best,<o:p></o:p></p>
<p class="MsoNormal">Kacie<o:p></o:p></p>
<p class="MsoNormal">---<o:p></o:p></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"><o:p></o:p></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"><o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><o:p> </o:p></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"><o:p></o:p></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"><o:p></o:p></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"><o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Courier New";color:#222222">Work: (775) 674-7049</span><span style="font-size:9.5pt;font-family:"Arial",sans-serif;color:#222222"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></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>
</body>
</html>