Re: Streamline Density

Thanks Tom,

I was hoping that the McIDAS streamlines and the VisAD ones were the same, and the latter just required a bit of tweaking.

The McIDAS streamlines certainly look better as the density improves (and if one forces the density of the VisAD streamlines to be greater than 2.0, then strange things happen such as streamlines crossing each other... probably why the density is limited to 2.0 in the code!).

So in the longer term, we'd be interested discussing getting the your (McIDAS) routine into VisAD.



Tom Whittaker wrote:

I cannot speak for VisAD, however in McIDAS the algorithm for drawing
streamlines is the one described in:  Monthly Weather Review: Vol.
105, No. 6, pp. 786–788.

The "density" value is used to determine in what grid "boxes"
streamlines can start.  Highest density says that a streamline will
start in every grid box and be drawn to its termination (rules for
stopping drawing are described in the paper).  Streamlines are always
started in the center of a grid box, thus no more than one will ever
be started in a box -- that does not, of course, prevent more than one
streamline from passing through a grid box.

Hope that helps.


On 3/29/06, James Kelly <J.Kelly@xxxxxxxxxx> wrote:

I've been looking at the quality of streamline analyses (of atmospheric
wind fields) lately.

I was experimenting with streamline displays in the McIDAS system, and
thought (IMHO) that the streamline analyses looked better as the density
of the streamlines increased.

I then tried to increase the density of the streamlines in a visad
display, and but wasn't able to generate the same density as the McIDAS
Should I be able to do this? If so how?

Tried changing "streamlineDensity" from the default of 1.0 to higher
values of 10.0, 100.0, 1000.0 and 1000.0, but didn't notice much change
beyond 100.0.

Further details:
1) In McIDAS I did this:
-100 CINT=0 COLOR=15

Image is attached as "mcidas_streamline.png"

2) In Visad/jython I did this:

maps = (xMap, yMap, rngUMap, rngVMap)

disp = subs.makeDisplay(maps)
flowcontrol = rngUMap.getControl()
strmDensity = flowcontrol.getStreamlineDensity()
stepFactor = flowcontrol.getStepFactor()
print " streamline density " , strmDensity, " stepFactor ", stepFactor

drLAPSuv=subs.addData("windbarbs", uv, disp)

Image is attached as "visad_streamline.png"



James Kelly
Australian Bureau of Meteorology

Tom Whittaker
University of Wisconsin-Madison
1225 W. Dayton Street
Madison, WI  53706  USA
ph: 608-262-2759

To unsubscribe visad, visit:

James Kelly Meteorological Systems Section Bureau of Meteorology
PO Box 1289K                                  Melbourne 3001, Australia
Phone: 61-3-9669-4724 Fax: 61-3-9669-4128     Email: J.Kelly@xxxxxxxxxx

To unsubscribe visad, visit:

  • 2006 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the visad archives: