File: streamlines.pro

package info (click to toggle)
cloc 2.04-1
  • links: PTS, VCS
  • area: main
  • in suites: trixie
  • size: 7,776 kB
  • sloc: perl: 29,368; cpp: 1,219; ansic: 334; asm: 267; makefile: 240; sh: 186; sql: 144; java: 136; ruby: 111; cs: 104; python: 84; pascal: 52; lisp: 50; cobol: 35; f90: 35; haskell: 35; objc: 25; php: 22; javascript: 15; fortran: 9; ml: 8; xml: 7; tcl: 2
file content (43 lines) | stat: -rw-r--r-- 700 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
 
; http://www.harrisgeospatial.com/docs/streamlines.html

; Read the data.

RESTORE, FILEPATH('globalwinds.dat', SUBDIR=['examples','data'])

 

; Set up the map projection, grid, and continents.

map = MAP('Equirectangular', POSITION=[0.1,0.1,0.9,0.9], $

    LIMIT=[0,-160,80,-50], TITLE='Wind Streamlines')

 

; Change some map grid properties.

grid = map.MAPGRID

grid.LINESTYLE = "dotted"

grid.ANTIALIAS = 0

grid.LABEL_POSITION = 0

grid.LABEL_ANGLE = 0

grid.FONT_SIZE=11

 

cont = MAPCONTINENTS(FILL_COLOR="light gray")

 

; Display the streamlines on top of the map.

stream = STREAMLINE(u, v, x, y, /OVERPLOT, $

    STREAMLINE_STEPSIZE=0.05, $

   RGB_TABLE=33, AUTO_COLOR=1, THICK=5)