<div dir="ltr">Hi Eric,<div><br></div><div>Your inputs are valid. Thanks for your help.</div><div><br></div><div>On another case, with the same data set I am using the below command to display the u_velocity and v_velocity vars as vector this time: </div>
<div><br></div><div><div><b>reinit</b></div><div><b>sdfopen <a href="http://rtofs_glo_2ds_f024_daily_prog.nc">rtofs_glo_2ds_f024_daily_prog.nc</a></b></div><div><b>set vpage 0 11 0 8.5</b></div><div><b>set parea 0 11 0 8.5</b></div>
<div><b>set lon 199.0503 291.3354</b></div><div><b>set lat -24.15176601231298 25.77516058680343</b></div><div><b>set mproj scaled</b></div><div><b>set mproj off</b></div><div><b>set display color white</b></div><div><b>set datawarn off</b></div>
<div><b>set grads off</b></div><div><b>set timelab off</b></div><div><b>set xlab off</b></div><div><b>set ylab off</b></div><div><b>set frame off</b></div><div><b>set t 1</b></div><div><b>set cint 40</b></div><div><b>set display color white</b></div>
<div><b>set background 0 </b></div><div><b>set gxout vector</b></div><div><b>set cthick 6</b></div><div><b>set ccolor 1</b></div><div><b>set arrscl 0.11 </b></div><div><b><br></b></div><div><b><font color="#ff0000" style="background-color:rgb(255,255,255)">d skip(u_velocity*1.95,3,6);skip(v_velocity*1.95,3,6)</font></b></div>
</div><div><b><br></b></div><div>And the output is as per expectation and is attached.</div><div><br></div><div>In display command, Am using skip factor to clearly visualize the arrow display otherwise it goes sluggish. What I need is to display the clear output but the value must be average of the other skipped value. So that user can choose the area, where intensity is high from full plot (withgout going down to greater zoom on plots). </div>
<div><br></div><div>I point am explaining is, I guess, areal average but how can i achieve in my requirement.</div><div><br></div><div><br></div><div>--Anil</div><div><br></div><div><br></div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Fri, Apr 4, 2014 at 12:38 AM, Eric Altshuler <span dir="ltr"><<a href="mailto:ela@cola.iges.org" target="_blank">ela@cola.iges.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div style="font-size:12pt;font-family:arial,helvetica,sans-serif">Hello Anil,<div><br></div><div>You need to do 'set gxout stream' before displaying the variables. By default, the gxout type for displaying two variables ('display U;V') is 'vector', so you'll get a vector plot unless you tell grads you want something different. To control streamline appearance, use 'set strmden …' The plot you attached was probably made by manually drawing the currents, and grads can't do this automatically.</div>
<div><br></div><div>Best regards,<br><br><div><span name="x"></span>Eric L. Altshuler<br>Research Scientist<br>Center for Ocean-Land-Atmosphere Studies<br>253 Research Hall, Mail Stop 6C5<div>George Mason University</div>
<div>4400 University Drive<br>Fairfax, VA 22030 USA<br><br>E-mail: <a href="mailto:ela@cola.iges.org" target="_blank">ela@cola.iges.org</a><br>Phone: (703) 993-5725<br>Fax: (703) 993-5770</div><span name="x"></span><br></div>
<br><hr><b>From: </b>"Anil Kumar Sharma" <<a href="mailto:aanil649@gmail.com" target="_blank">aanil649@gmail.com</a>><br><b>To: </b>"GrADS Users Forum" <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
<b>Sent: </b>Thursday, April 3, 2014 3:13:38 AM<br><b>Subject: </b>[gradsusr] To Plot Ocean Current streams<div><div class="h5"><br><br><div dir="ltr">Hi Everyone,<div><br></div><div>I have NetCDF data downloaded for Ocean Current from NOAA website and is attached for reference. Below is the NCDump</div>
<div><br></div><div><div><b>netcdf rtofs_glo_2ds_f024_daily_prog {</b></div>
<div><b>dimensions:</b></div><div><b> lat = 360 ;</b></div><div><b> lev = 1 ;</b></div><div><b> lon = 720 ;</b></div><div><b> time = 1 ;</b></div><div><b>variables:</b></div><div><b> double lat(lat) ;</b></div>
<div><b> lat:grads_dim = "y" ;</b></div><div><b> lat:grads_mapping = "linear" ;</b></div><div><b> lat:grads_size = "2160" ;</b></div><div><b> lat:units = "degrees_north" ;</b></div>
<div><b> lat:long_name = "latitude" ;</b></div><div><b> lat:minimum = -90. ;</b></div><div><b> lat:maximum = 89.90947 ;</b></div><div><b> lat:resolution = 0.08333f ;</b></div>
<div><b> double lev(lev) ;</b></div><div><b> lev:grads_dim = "z" ;</b></div><div><b> lev:grads_mapping = "levels" ;</b></div><div><b> lev:units = "millibar" ;</b></div>
<div><b> lev:long_name = "altitude" ;</b></div><div><b> lev:minimum = 1. ;</b></div><div><b> lev:maximum = 1. ;</b></div><div><b> double lon(lon) ;</b></div><div>
<b> lon:grads_dim = "x" ;</b></div><div><b> lon:grads_mapping = "linear" ;</b></div><div><b> lon:grads_size = "4320" ;</b></div><div><b> lon:units = "degrees_east" ;</b></div>
<div><b> lon:long_name = "longitude" ;</b></div><div><b> lon:minimum = 74.16 ;</b></div><div><b> lon:maximum = 434.06227 ;</b></div><div><b> lon:resolution = 0.08333f ;</b></div>
<div><b> double time(time) ;</b></div><div><b> time:grads_dim = "t" ;</b></div><div><b> time:grads_mapping = "linear" ;</b></div><div><b> time:grads_size = "9" ;</b></div>
<div><b> time:grads_min = "00z13mar2014" ;</b></div><div><b> time:grads_step = "1dy" ;</b></div><div><b> time:units = "days since 1-1-1 00:00:0.0" ;</b></div>
<div><b> time:long_name = "time" ;</b></div><div><b> time:minimum = "00z13mar2014" ;</b></div><div><b> time:maximum = "00z21mar2014" ;</b></div><div>
<b> time:resolution = 1.f ;</b></div><div><b> float u_velocity(time, lev, lat, lon) ;</b></div><div><b> u_velocity:_FillValue = 1.26765e+030f ;</b></div><div><b> u_velocity:missing_value = 1.26765e+030f ;</b></div>
<div><b> u_velocity:long_name = "eastward_sea_water_velocity (m/s) " ;</b></div><div><b> float v_velocity(time, lev, lat, lon) ;</b></div><div><b> v_velocity:_FillValue = 1.26765e+030f ;</b></div>
<div><b> v_velocity:missing_value = 1.26765e+030f ;</b></div><div><b> v_velocity:long_name = "northward_sea_water_velocity (m/s) " ;</b></div><div><b><br></b></div><div><b>// global attributes:</b></div>
<div><b> :title = "2D Sfc Daily Prognostic 00Z13mar2014: Forecast, downlo</b></div><div><b>aded Mar 13 15:18 UTC" ;</b></div><div><b> :Conventions = "COARDS\n",</b></div>
<div>
<b> "GrADS" ;</b></div><div><b> :dataType = "Grid" ;</b></div><div><b> :history = "Thu Mar 13 11:21:41 EDT 2014 : imported by GrADS Data Server 2.0" ;</b></div>
<div><b> :NCO = "20140314" ;</b></div><div><b>}</b></div></div><div><br></div><div><br></div><div>I need to Plot the current Streams as output, Attached image can be treated as output sample.<br>
currently, I am using below command for displaying u_velocity and v_velocity variables </div>
<div><br></div><div><div><b>d skip(u_velocity*1.95,14,14);skip(v_velocity*1.95,14,14)</b></div></div><div><b><br></b></div><div>with skip function I am only able to plot arrow but no stream (we can say axis of current)</div>
<div><br></div><div>I really don't know how to get this working. Any help on it please</div><div><br></div><div>-- Anil</div><div><br> </div>
</div>
<br></div></div>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>