<div dir="ltr">Thank you Dennis and Adam, finally the code is running. I am grateful for your inputs<div><br></div><div>Best Regards</div><div>Ipshita</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Sep 10, 2016 at 6:09 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</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"><div>diri = "/Users/ipshita/Documents/<wbr>grads_NCL/"<br><br>c=addfile(diri+"<a href="http://vwnd.2008.nc" target="_blank">vwnd.2008.nc</a>",<wbr>"r")<br>e=addfile(diri+"<a href="http://pres.sfc.mon.mean.nc" target="_blank">pres.sfc.mon.<wbr>mean.nc</a>","r")<br><br>v = c->vwnd ; (time,level,lat,lon) <br>ps = e->pres ; (time,lat,lon) hPa<br><br>printVarSummary(ps)<br>print("---")<br>printVarSummary(v)<br>print("---")<br>print(v&level) <wbr> ; hPa (1000-10) => bottom-to-top<br>print("---")<br>; =========================<br></div>; function requires top-to-bottom; reverse order<br>;==========================<br><div>v = v(:,::-1,:,:)<br><br>printVarSummary(v) <br>print("---")<br>print(v&level) <wbr> ; hPa (10-1000) => top-to-bottom<br>print("---")<br><br>; =========================<br></div><div>; function requires Pa ... not hPa/mb<br>;==========================<br>ps =(/ ps*100 /) ; function requires Pa<br>ps@units = "Pa"<br><br>v&level = (/ 100*v&level /) ; function requires Pa<br>v&level@units = "Pa" <br><br>zmpsi = zonal_mpsi_Wrap(v, v&lat, v&level, ps)<br>printVarSummary(zmpsi)<br><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 9, 2016 at 5:10 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</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 Ipshita,<div>The printVarSummary output shows what the problem is:</div><span><div><div style="font-size:12.8px">Type: float</div><div style="font-size:12.8px">Total Size: 714816 bytes</div><div style="font-size:12.8px"> 178704 values</div><div style="font-size:12.8px">Number of Dimensions: 3</div><div style="font-size:12.8px">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[level | 17] x [lat | 73] x [lon | 144]</div><div style="font-size:12.8px">Coordinates: </div><div style="font-size:12.8px"> <span style="background-color:rgb(255,217,102)"> level: [10..1000]</span></div></div><div style="font-size:12.8px"><span style="background-color:rgb(255,217,102)"><br></span></div></span><div style="font-size:12.8px">The error message that you are recieving states that the levels need to be in Pa:</div><span><div style="font-size:12.8px"><span style="font-size:12.8px">fatal:zonal_mpsi: The pressure array must be between the values of 500 and 100500 (exclusive), and monotonically increasing.</span><br></div><div style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div></span><div style="font-size:12.8px"><span style="font-size:12.8px">Thus, you need to convert your pressure levels to Pa. I am unsure whether zonal_mpsi uses the level coordinate variable attached to your v1 array, so it is probably better to change that (as opposed to just altering the v1&level array itself.</span></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">level = v1&level</div><div style="font-size:12.8px">delete(v1&level)</div><div style="font-size:12.8px">level = level*100.</div><div style="font-size:12.8px">level&level = level</div><div style="font-size:12.8px">level@units = "Pa"</div><div style="font-size:12.8px">v1&level = level</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">That should fix the issue..</div><span><font color="#888888"><div style="font-size:12.8px">Adam</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><span style="font-size:12.8px"> </span></div></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 9, 2016 at 12:01 PM, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</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 ,<div><br></div><div>Thank you for your input. I did try the above mentioned correction. I still get an error :</div><div><br></div><div>fatal:zonal_mpsi: The pressure array must be between the values of 500 and 100500 (exclusive), and monotonically increasing<br></div><div><br></div><div>Here is the printVarSummary and print for v1 Variable: v1</div><div><br></div><div>Type: float</div><div>Total Size: 714816 bytes</div><div> 178704 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[level | 17] x [lat | 73] x [lon | 144]</div><div>Coordinates: </div><div> level: [10..1000]</div><div> lat: [90..-90]</div><div> lon: [ 0..357.5]</div><div>Number Of Attributes: 14</div><div> long_name :<span style="white-space:pre-wrap"> </span>Monthly mean v wind</div><div> units :<span style="white-space:pre-wrap"> </span>m/s</div><div> precision :<span style="white-space:pre-wrap"> </span>2</div><div> least_significant_digit :<span style="white-space:pre-wrap"> </span>1</div><div> var_desc :<span style="white-space:pre-wrap"> </span>v-wind</div><div> level_desc :<span style="white-space:pre-wrap"> </span>Multiple levels</div><div> statistic :<span style="white-space:pre-wrap"> </span>Mean</div><div> parent_stat :<span style="white-space:pre-wrap"> </span>Other</div><div> missing_value :<span style="white-space:pre-wrap"> </span>-9.96921e+36</div><div> valid_range :<span style="white-space:pre-wrap"> </span>( -125, 160 )</div><div> dataset :<span style="white-space:pre-wrap"> </span>NCEP Reanalysis Derived Products</div><div> actual_range :<span style="white-space:pre-wrap"> </span>( -68.57001, 68.56999 )</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.96921e+36</div><div> time :<span style="white-space:pre-wrap"> </span>1297320</div><div><br></div><div><br></div><div>Variable: level (coordinate)</div><div>Type: float</div><div>Total Size: 68 bytes</div><div> 17 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[level | 17]</div><div>Coordinates: </div><div>Number Of Attributes: 7</div><div> units :<span style="white-space:pre-wrap"> </span>millibar</div><div> long_name :<span style="white-space:pre-wrap"> </span>Level</div><div> positive :<span style="white-space:pre-wrap"> </span>down</div><div> GRIB_id :<span style="white-space:pre-wrap"> </span>100</div><div> GRIB_name :<span style="white-space:pre-wrap"> </span>hPa</div><div> actual_range :<span style="white-space:pre-wrap"> </span>( 1000, 10 )</div><div> axis :<span style="white-space:pre-wrap"> </span>Z</div><div>(0)<span style="white-space:pre-wrap"> </span>10</div><div>(1)<span style="white-space:pre-wrap"> </span>20</div><div>(2)<span style="white-space:pre-wrap"> </span>30</div><div>(3)<span style="white-space:pre-wrap"> </span>50</div><div>(4)<span style="white-space:pre-wrap"> </span>70</div><div>(5)<span style="white-space:pre-wrap"> </span>100</div><div>(6)<span style="white-space:pre-wrap"> </span>150</div><div>(7)<span style="white-space:pre-wrap"> </span>200</div><div>(8)<span style="white-space:pre-wrap"> </span>250</div><div>(9)<span style="white-space:pre-wrap"> </span>300</div><div>(10)<span style="white-space:pre-wrap"> </span>400</div><div>(11)<span style="white-space:pre-wrap"> </span>500</div><div>(12)<span style="white-space:pre-wrap"> </span>600</div><div>(13)<span style="white-space:pre-wrap"> </span>700</div><div>(14)<span style="white-space:pre-wrap"> </span>850</div><div>(15)<span style="white-space:pre-wrap"> </span>925</div><div>(16)<span style="white-space:pre-wrap"> </span>1000</div><div><br></div><div>Variable: PsPa</div><div>Type: float</div><div>Total Size: 42048 bytes</div><div> 10512 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[73] x [144]</div><div>Coordinates: </div><div>Number Of Attributes: 1</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>-9.96921e+36</div><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 9, 2016 at 8:10 AM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</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 Ipshita,<div>You are altering the order of your level dimension here:</div><div><span style="font-size:12.8px">v1=v_jan_1948(::-1,:,:) </span><br></div><div><span style="font-size:12.8px">and that is completely fine. However, you are </span><span style="font-size:12.8px">then flipping the level dimension again here when you specify v1&level:</span></div><div><span style="font-size:12.8px">zmpsi = zonal_mpsi(v1,v1&lat,v1&level(</span><span style="font-size:12.8px"><wbr>::-1),PsPa)</span><span style="font-size:12.8px"><br></span></div><div><br></div><div>Try changing this:</div><div><span style="font-size:12.8px">zmpsi = zonal_mpsi(v1,v1&lat,v1&level(</span><span style="font-size:12.8px"><wbr>::-1),PsPa)</span><br></div><div><span style="font-size:12.8px">to this:</span></div><div><span style="font-size:12.8px">zmpsi = zonal_mpsi(v1,v1&lat,v1&level</span><span style="font-size:12.8px">,<wbr>PsPa)</span><span style="font-size:12.8px"><br></span></div><div><br></div><div>If that does not fix the issue please add printVarSummary and print statements to your script and examine the output. The end of your script should look like this:</div><div><span style="font-size:12.8px">v1=v_jan_1948(::-1,:,:)</span><br></div><div><span style="font-size:12.8px">PsPa=ps_jan_1948*100</span><span style="font-size:12.8px"><br></span></div><div><div>printVarSummary(v1)</div><div>print(v1&level)</div></div><div><span style="font-size:12.8px">zmpsi = zonal_mpsi(v1,v1&lat,v1&level</span><span style="font-size:12.8px">,<wbr>PsPa)</span><br></div><div><br></div><div>If you still are having trouble please respond to the ncl-talk email list with the printVarSummary and print output.</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Thu, Sep 8, 2016 at 5:02 PM, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><br clear="all"><div>Dear Users,</div><div><br></div><div>I want to plot stream function of wind (v) and I am using zonal_mpsi to calculate the stream function. I keep getting error messages.</div><div><br></div><div>fatal:zonal_mpsi: The pressure array must be between the values of 500 and 100500 (exclusive), and monotonically increasing.<br></div><div><br></div><div>I am not sure how to fix this problem.</div><div><br></div><div>I am attaching my code.</div><div><br></div><div>I read in the V wind and surface pressure</div><div><br></div><div>v-wind has 4 dimensions (time,lev,lat,lon)</div><div><br></div><div>Surface pressure has 3 dimensions (time,lat,lon)</div><div><br></div><div>Here is the code I wrote</div><div><br></div><div>;*****************************<wbr>******************************<wbr>***</div><div><div>;*****************************<wbr>**************</div><div>load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl"</div><div>;*****************************<wbr>**************</div><div><br></div><div>;*****************************<wbr>***************</div><div>;This is to read wind speed and wind direction</div><div>;from 1948-2016</div><div>;*****************************<wbr>***************</div><div><br></div><div>a=addfile("/Users/ipshita/Docu<wbr>ments/grads_NCL/<a href="http://wspd.mon.mean.nc" target="_blank">wspd.mon.mean.<wbr>nc</a>","r")</div><div>b=addfile("/Users/ipshita/Docu<wbr>ments/grads_NCL/<a href="http://uwnd.mon.mean.nc" target="_blank">uwnd.mon.mean.<wbr>nc</a>","r")</div><div>c=addfile("/Users/ipshita/Docu<wbr>ments/grads_NCL/<a href="http://vwnd.mon.mean.nc" target="_blank">vwnd.mon.mean.<wbr>nc</a>","r")</div><div>d=addfile("/Users/ipshita/Docu<wbr>ments/grads_NCL/<a href="http://omega.mon.mean.nc" target="_blank">omega.mon.mean<wbr>.nc</a>","r")</div><div>e=addfile("/Users/ipshita/Docu<wbr>ments/grads_NCL/<a href="http://pres.sfc.mon.mean.nc" target="_blank">pres.sfc.mon.m<wbr>ean.nc</a>","r")</div><div><br></div><div>v=c->vwnd ;Wind in horizontal direction in m/s</div><div><br></div><div>level=c->level ;Pressure level from 1000hpa to 10 hpa</div><div><br></div><div>lat=c->lat</div><div><br></div><div>lon=c->lon</div><div><br></div><div>ps=e->pres</div><div><br></div><div>ps_jan_1948=ps(0,:,:)</div><div><br></div><div>printMinMax(ps_1948_2015,False<wbr>)</div><div><br></div><div>v_jan_1948=v(0,:,:,:)</div><div>;*****************************<wbr>****************</div><div>;Going from low to high in terms of level</div><div>;*****************************<wbr>****************</div><div><br></div><div>v1=v_jan_1948(::-1,:,:)</div><div>;*****************************<wbr>****************</div><div>;Converting Pressure into Pascals</div><div>;*****************************<wbr>****************</div><div><br></div><div>PsPa=ps_jan_1948*100</div><div><br></div><div>printMinMax(PsPa,False)</div><div><br></div><div>;*****************************<wbr>****************</div><div>zmpsi = zonal_mpsi(v1,v1&lat,v1&level(<wbr>::-1),PsPa)</div></div><div>;*****************************<wbr>******************************<wbr>********</div><div><br></div><div>I will be grateful if you can help me resolve this issue.</div><div><br></div><div>Best Regards</div><span><font color="#888888"><div>Ipshita</div><div><div dir="ltr"><br></div></div>
</font></span></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><span><font color="#888888"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a> </font></span><span><font color="#888888"><a href="tel:303-497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></font></span></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><div>Ipshita Majhi<br></div>PhD Candidate<br></div>University of Alaska , Fairbanks<br></div>Atmospheric Science Department<br></div><a href="tel:%28907%29978-4220" value="+19079784220" target="_blank">(907)978-4220</a> <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a> </font></span><span><font color="#888888"><a href="tel:303-497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>
</div></div><br>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><div>Ipshita Majhi<br></div>PhD Candidate<br></div>University of Alaska , Fairbanks<br></div>Atmospheric Science Department<br></div>(907)978-4220 <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</div>