<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div>*** The most important rule in data processing is "look at your data". ***<br>------<br><br></div>The reason that ncl-talk recommends using <br><br></div>printVarSummary(...)<br><br></div>and<br><br></div>printMinMax(..., ...)<br>------<br><br></div>is so that *users* can examine a variable's<br><br>(a) data type<br>(b) dimension sizes and dimension names<br></div></div>(c) units<br></div>etc<br><br></div>You have<br><br> <br> T = in->t ; select variable to ave<br> W = in->w<br> W=-W<br> V = in->v<br> U = in->u<br> Q = in->q<br><br> printVarSummary(T)<br> printVarSummary(W)<br> printVarSummary(V)<br> printVarSummary(U)<br> printVarSummary(hyam)<br> printVarSummary(hybm)<br> printVarSummary(PS)<br><br><br></div>then you immediately do the vertical interpolation<br><br>;---Define other arguments required by vinth2p<br> interp = 2<br> pnew = (/1000,975,950,900,850,800,750,700,600,500,400,300,250,200,150,100,50/)<br> pnew@units = "mb"<br><br>;---Interpolate to pressure levels on pressure levels<br> t = vinth2p(T,hyam,hybm,pnew,PS,interp,P0mb,1,False)<br><br>=========<br></div>The initial<br><br></div>printVarSummary(T) <br><div><br><div>Variable: T</div><div><b>Type: short</b></div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 7</div><div> standard_name : air_temperature</div><div> long_name : Temperature</div><div> units : K</div><div><b> add_offset : 249.3816473297584</b></div><div><b> scale_factor : 0.002009905912815991<br><br></b></div><div>Before use, you ***must***<b> UNPACK the variable<br><br>---<br></b></div><div>Please read the documentation of the following:<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml">https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml</a><b><br><br> T = short2flt(in->t) ; select variable: unpack; later ave<br></b></div><div><b> printVarSummary(T)<br></b></div><div><b> printMinMax(T,0)<br><br></b></div><div>Please look at the output.<br></div><div><b></b></div><div><b><br></b>---<br></div><div>So, unpack each variable that is type short on the file. Thendo the vertical interpolation, etc<br><br><br></div><div>Good Luck<br></div><div><b><br><br><br></b></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 27, 2018 at 5:25 AM, Jenny Taylor <span dir="ltr"><<a href="mailto:jenny63taylor91@gmail.com" target="_blank">jenny63taylor91@gmail.com</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">My thanks to Dennis Shea and Adam Philips for valuable advice.<div>I am trying to plot ERA-Interim model data at different vertical pressure levels. I have downloaded the 'hyam' and 'hybm' for model levels (attached below) as specified in the ERA-Interim link; </div><div><br></div><div><a href="https://rda.ucar.edu/datasets/ds627.0/docs/Eta_coordinate/index.html" target="_blank">https://rda.ucar.edu/datasets/<wbr>ds627.0/docs/Eta_coordinate/<wbr>index.html</a><br></div><div><br></div><div>Downloaded the daily data (time step=4) for all 60 levels in the separate file (4 GB approx). Then I merged these two files in CDO and time mean the merged data.</div><div><br></div><div>Summary of the final processed data is given below;</div><div><div><br></div><div><div>Variable: PS1</div><div>Type: short</div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 6</div><div> long_name : Logarithm of surface pressure</div><div> units : ~</div><div> add_offset : 11.19559352539406</div><div> scale_factor : 1.015380172027476e-05</div><div> _FillValue : -32767</div><div> missing_value : -32767</div><div><br></div><div>Variable: T</div><div>Type: short</div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 7</div><div> standard_name : air_temperature</div><div> long_name : Temperature</div><div> units : K</div><div> add_offset : 249.3816473297584</div><div> scale_factor : 0.002009905912815991</div><div> _FillValue : -32767</div><div> missing_value : -32767</div><div><br></div><div>Variable: W</div><div>Type: short</div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 7</div><div> standard_name : lagrangian_tendency_of_air_<wbr>pressure</div><div> long_name : Vertical velocity</div><div> units : Pa s**-1</div><div> add_offset : -1.799464318323461</div><div> scale_factor : 0.0001232090956517211</div><div> _FillValue : -32767</div><div> missing_value : -32767</div><div><br></div><div>Variable: V</div><div>Type: short</div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 7</div><div> standard_name : northward_wind</div><div> long_name : V component of wind</div><div> units : m s**-1</div><div> add_offset : -0.1744413644931722</div><div> scale_factor : 0.001333290509781928</div><div> _FillValue : -32767</div><div> missing_value : -32767</div><div><br></div><div>Variable: U</div><div>Type: short</div><div>Total Size: 2724120 bytes</div><div> 1362060 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [level | 60] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> level: [ 1.. 60]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 7</div><div> standard_name : eastward_wind</div><div> long_name : U component of wind</div><div> units : m s**-1</div><div> add_offset : -14.17403101359716</div><div> scale_factor : 0.002084243991189315</div><div> _FillValue : -32767</div><div> missing_value : -32767</div><div><br></div><div>Variable: hyam</div><div>Type: float</div><div>Total Size: 240 bytes</div><div> 60 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes: [lvl | 60]</div><div>Coordinates:</div><div> lvl: [1..60]</div><div>Number Of Attributes: 2</div><div> long_name : a model ave</div><div> units : Pa</div><div><br></div><div>Variable: hybm</div><div>Type: float</div><div>Total Size: 240 bytes</div><div> 60 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes: [lvl | 60]</div><div>Coordinates:</div><div> lvl: [1..60]</div><div>Number Of Attributes: 2</div><div> long_name : b model ave</div><div> units : Pa Pa**-1</div><div><br></div><div>Variable: PS</div><div>Type: float</div><div>Total Size: 90804 bytes</div><div> 22701 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes: [1] x [161] x [141]</div><div>Coordinates:</div><div>Number Of Attributes: 1</div><div> _FillValue : -32767</div><div><br></div><div>Variable: t</div><div>Type: float</div><div>Total Size: 1543668 bytes</div><div> 385917 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 17] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> lev_p: [1000..50]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 1</div><div> _FillValue : -32767</div><div><br></div><div>Variable: u</div><div>Type: float</div><div>Total Size: 1543668 bytes</div><div> 385917 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 17] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> lev_p: [1000..50]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 1</div><div> _FillValue : -32767</div><div><br></div><div>Variable: v</div><div>Type: float</div><div>Total Size: 1543668 bytes</div><div> 385917 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 17] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> lev_p: [1000..50]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 1</div><div> _FillValue : -32767</div><div><br></div><div>Variable: w</div><div>Type: float</div><div>Total Size: 1543668 bytes</div><div> 385917 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 17] x [latitude | 161] x [longitude | 141]</div><div>Coordinates:</div><div> time: [969705..969705]</div><div> lev_p: [1000..50]</div><div> latitude: [40.. 0]</div><div> longitude: [65..100]</div><div>Number Of Attributes: 1</div><div> _FillValue : -32767</div></div></div><div><br></div><div>My output figure is attached below. Vertical motion is not shaped in proper form. The sample data used in this plot are for August 2010. Kindly inform me what I am doing wrong in this case. <br></div><div><br></div><div>Thanking you,</div><div>Jenny</div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 27, 2018 at 3:46 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><div>Hybrid levels are a form of terrain following vertical coordinate formulation. <br>Unlike isobaric levels where the vertical levels are the same at every grid point, <br>hybrid levels differ at each grid point <b>(j,i) </b>because they are a function of surface pressure<br></div>which varies with time and location. The climate atmospheric model uses the following<br></div>formulation:<br><div><div><br> <b>
<i>p([t,] k,j,i) = a(k)*p0 + b(k)*ps([t,] k,j,i) </i></b>; t=>time, k=>level, j=>latitude; i=longitude<br><br>Each calculated pressure level is called the '<b>m</b>id-level' (hya<b>m</b> and hyb<b>m</b>). Each <br>'<b>m</b>id-level' has an upper and lower boundary called the<b> i</b>nterface levels (hya<b>i</b>, hyb<b>i</b>). <br><br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/pres_hybrid_ccm.shtml" target="_blank">http://www.ncl.ucar.edu/Docume<wbr>nt/Functions/Built-in/pres_<wbr>hybrid_ccm.shtml</a><br><br></div><div>There are other hybrid level formulations so be aware of what is being used.<br><br></div><div>A detailed description with some pictures is here:<br></div><div><br><a href="https://rda.ucar.edu/datasets/ds627.0/docs/Eta_coordinate/index.html" target="_blank">https://rda.ucar.edu/datasets/<wbr>ds627.0/docs/Eta_coordinate/in<wbr>dex.html</a><br><br>====<br></div><div>You should talk with an advisor about these topics. ncl-talk is really for language issues.<br><br></div><div>Good luck<br></div></div></div><div class="m_3418026016588114415HOEnZb"><div class="m_3418026016588114415h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 26, 2018 at 2:23 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 Jenny,<div>With regards to whether the plot looks correct: It looks like you took NCL's vector_5 example script, and modified it slightly. That script uses the <a href="http://atmos.nc" target="_blank">atmos.nc</a> NetCDF file, which according to the global history attribute is an old NCAR-CSM file. The data may be October averages. You're selecting 20E as the longitude to plot, but are basing your weights on 160E. All that being said, i<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">t is honestly up to you to investigate the file you are reading in, and to interpret the results. I am not an expert at looking at this kind of diagnostic.</span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">With regards to your question about hyam and hybm, that information is given within the NetCDF file, which you can see by doing a ncdump <a href="http://atmos.nc" target="_blank">atmos.nc</a>. More information on hyam/hybm can be found by entering either term within the NCL website search engine. </span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">If you have any further questions please respond to the ncl-talk email list.</span></div><span class="m_3418026016588114415m_-5384143621164482691HOEnZb"><font color="#888888"><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Adam </span></div></font></span></div><div class="m_3418026016588114415m_-5384143621164482691HOEnZb"><div class="m_3418026016588114415m_-5384143621164482691h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 26, 2018 at 4:33 AM, Jenny Taylor <span dir="ltr"><<a href="mailto:jenny63taylor91@gmail.com" target="_blank">jenny63taylor91@gmail.com</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>Hello,</div><div><br></div>I am sending you the script based on which it is made. Kindly find the script attached below. I have one question on this matter. In the dataset (<a href="http://atmos.nc" target="_blank">atmos.nc</a>) what are those variables named as 'hyam' and 'hybm'? <div>It is showing as 'hybrid A coefficient at layer midpoints' and 'hybrid B coefficient at layer midpoints'.</div><div><br><div>Kindly also inform me what I am doing wrong as per your observation over NH.</div><div><br></div><div>Thanks,</div><div>Jenny<br><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 26, 2018 at 12:16 PM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</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>What did you plot? Is that omega or vertical velocity? I am wondering why the northern hemisphere winds are not convergent with the positive vertical motion values like they are in the southern hemisphere.</div><div><br></div><div>Barry</div></div><div class="gmail_extra"><div><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756h5"><br><div class="gmail_quote">On Mon, Feb 26, 2018 at 8:32 AM, Jenny Taylor <span dir="ltr"><<a href="mailto:jenny63taylor91@gmail.com" target="_blank">jenny63taylor91@gmail.com</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>Thank you very much, Adam Phillips. I follow your instruction, please look at my attached output figure. Kindly inform me it is in proper shape or not. I think this time it makes some sense.</div><div><br></div><div>Thanks again,</div><div>Jenny</div><div><div><br></div></div></div><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498HOEnZb"><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498h5"><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Feb 24, 2018 at 2:14 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 Jenny,<div>You are reading in OMEGA. Note that OMEGA has opposite signs and different units from the conventional W variable. </div><div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><br>For upward motion: W is positive and OMEGA is negative<br>for downward motion: W is negative and OMEGA is positive<br><br>A basic, simple relationship<br> <span class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992gmail-m_-4504464062761879542gmail-m_-4724114898966346179gmail-il">omega</span><span> </span>= -rho*grav*w ; Pa/s<br>Hence:<br> w = -<span class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992gmail-m_-4504464062761879542gmail-m_-4724114898966346179gmail-il">omega</span>/(rho*grav) ; (say) m/s<br><br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">Note the sign ....<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">-----<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">I believe all you have to do is flip the sign of your W variable:<br> W = in->OMEGA<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"> W = -W <span> </span><br>assuming you want upward arrows for upward motion.</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><br></div>Also note that NCL has a omega_to_w function:</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/omega_to_w.shtml" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Contributed/omeg<wbr>a_to_w.shtml</a><br></div><div><br></div><div>Hope that all makes sense! If not or if you have further questions please respond to the ncl-talk email list.</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><div><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908h5"><br><div class="gmail_quote">On Fri, Feb 23, 2018 at 1:24 AM, Jenny Taylor <span dir="ltr"><<a href="mailto:jenny63taylor91@gmail.com" target="_blank">jenny63taylor91@gmail.com</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">Mary,<div><br></div><div>I was trying with the following attached script. But vertical motion near the equator, as well as the 30S and 30N, shows rising motion in the output which is not correct in real term. Please kindly look into my script and figure and inform me what I am missing or doing any wrong in this.</div><div><br></div><div>Data structure is below;</div><div><br></div><div><div>Variable: t</div><div>Type: float</div><div>Total Size: 655360 bytes</div><div> 163840 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 20] x [lat | 64] x [lon | 128]</div><div>Coordinates:</div><div> time: [ 791.. 791]</div><div> lev_p: [950..50]</div><div> lat: [-87.8638..87.8638]</div><div> lon: [ 0..357.1875]</div><div>Number Of Attributes: 1</div><div> _FillValue : 9.96921e+36</div><div><br></div><div>Variable: u</div><div>Type: float</div><div>Total Size: 655360 bytes</div><div> 163840 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 20] x [lat | 64] x [lon | 128]</div><div>Coordinates:</div><div> time: [ 791.. 791]</div><div> lev_p: [950..50]</div><div> lat: [-87.8638..87.8638]</div><div> lon: [ 0..357.1875]</div><div>Number Of Attributes: 1</div><div> _FillValue : 9.96921e+36</div><div><br></div><div>Variable: v</div><div>Type: float</div><div>Total Size: 655360 bytes</div><div> 163840 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 20] x [lat | 64] x [lon | 128]</div><div>Coordinates:</div><div> time: [ 791.. 791]</div><div> lev_p: [950..50]</div><div> lat: [-87.8638..87.8638]</div><div> lon: [ 0..357.1875]</div><div>Number Of Attributes: 1</div><div> _FillValue : 9.96921e+36</div><div><br></div><div>Variable: w</div><div>Type: float</div><div>Total Size: 655360 bytes</div><div> 163840 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes: [time | 1] x [lev_p | 20] x [lat | 64] x [lon | 128]</div><div>Coordinates:</div><div> time: [ 791.. 791]</div><div> lev_p: [950..50]</div><div> lat: [-87.8638..87.8638]</div><div> lon: [ 0..357.1875]</div><div>Number Of Attributes: 1</div><div> _FillValue : 9.96921e+36</div></div><div><br></div><div><br></div><div>Thanking you,</div><div>Jenny</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 22, 2018 at 1:17 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">Jenny,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In order to do plot overlays of different variables, you need to use the "overlay" procedure in NCL.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">--------<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">------<wbr>--</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">----<wbr>----</span></div><br></div><div class="gmail_default" style="font-size:small">The overlay procedure works one of two ways:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[1] It overlays data from one plot (called the "overlay" plot) into the data space of another plot (called the "base" plot). </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The two plots must have intersecting data spaces for this to work. By "data space", I mean the range of the X and Y axes of the two plots must intersect. Otherwise, you will not see the overlay plot at all.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For an example, see overlay_8.ncl at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/overlay.shtml#ex8" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/overlay.shtml#ex8</a><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[2] It simply "lines up" the two plots in the same rectangular region and draws them. It's not as common to do overlays this way, because it requires that you do all the work of making sure the two plots are in the same data space.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For an example, see overlay_9.ncl at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><a href="http://www.ncl.ucar.edu/Applications/overlay.shtml#ex9" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/overlay.shtml#ex9</a><br></div><div class="gmail_default" style="font-size:small"><br></div></div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">------<wbr>--</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">--------</span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;background-color:rgb(255,255,255);float:none;display:inline">----<wbr>----</span><div class="gmail_default" style="font-size:small"><br class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992m_6400160184359177392m_8469251774671623798gmail-Apple-interchange-newline">It doesn't matter whether you are doing contour plots, vector plots, or XY plots: you can overlay any of these on another one using the "overlay" procedure.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The "overlay" examples page has several other scripts showing how to use the overlay procedure. In particular, you might want to look at overlay_1.ncl, overlay_6.ncl, and/or wrf_gsn_5.ncl:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><a href="http://www.ncl.ucar.edu/Applications/overlay.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/overlay.shtml</a></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">If you continue to have problems, please post back to ncl-talk and include your script and a "printVarSummary" of the data you are trying to plot.</div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">Thanks,</div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">--Mary</div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992m_6400160184359177392h5">On Fri, Feb 16, 2018 at 9:59 AM, Jenny Taylor <span dir="ltr"><<a href="mailto:jenny63taylor91@gmail.com" target="_blank">jenny63taylor91@gmail.com</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 class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992m_6400160184359177392h5"><div dir="auto">Dear NCL Users,<div dir="auto"><br></div><div dir="auto">I am trying to plot vector wind field at different pressure level for latitute range 40N-40S. I am following the below mention script;</div><div dir="auto"><a href="https://www.ncl.ucar.edu/Applications/Scripts/h_lat_7.ncl" target="_blank">https://www.ncl.ucar.edu/Appli<wbr>cations/Scripts/h_lat_7.ncl</a><br></div><div dir="auto"><br></div><div dir="auto">I am trying to overlay vertical motion (rising/sinking) in vector arrows and zonal wind in color scheme.</div><div dir="auto"><br></div><div dir="auto">Can anyone suggest me, in that case which funtion I need to use in it?</div><div dir="auto">Any help will be appriciated.</div><div dir="auto"><br></div><div dir="auto">Thanks,</div><div dir="auto">Jenny</div></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><br></div>
</blockquote></div><br></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><br clear="all"><div><br></div></div></div><span class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908HOEnZb"><font color="#888888">-- <br><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498m_-2543755846130901908m_3198375659086582992gmail_signature" 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)%20497-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></div>
</blockquote></div><br></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><br clear="all"><div><br></div>-- <br></div></div><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384m_7314380712568923756m_-351116948212892498gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_3418026016588114415m_-5384143621164482691m_-3801804134061286384gmail_signature" 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)%20497-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></div>
</div></div></blockquote></div><br></div>