from rsf.proj import *

# this exercise illustrates fetching data,
# calling simple functions, and various 2d plots. 

# Download data
Fetch('austin-w.HH','bay')

# People without internet, turn off fetch, and instead link data
#datapath = '../DATA/bay/austin-w.HH'
# Flow('austin-w.HH',datapath,'cp')

# convert to rsf format, and declare what dimensions
Flow('austin','austin-w.HH',
   '''
   dd form=native | 
   put d1=0.01 unit1=km d2=0.01 unit2=km|
   window min1=0 max1=12 min2=0 max2=12
   ''')

# plotting function, play around with parameters!
plotting = '''
grey title="Digital Elevation Map" mean=n color=i
scalebar=n barlabel=Elevation barunit=m
'''

# display digital elevation map
Result('austin',plotting)

# apply smoothing and display map
# uncomment and change values below!!!
#Flow('austin-s','austin','smooth rect1=30 rect2=30')
#Result('austin-s',plotting)

# for loop creating contour map
# add 'austin-s' to the loop list!!!
for dem in ('austin',):
    # plot the map as a .vpl
#    Plot(dem,plotting)
    # generate contour file
    contour = dem+'-cont'
#    Plot(contour,dem,'contour nc=6 dc=50 c0=100 wanttitle=n')
    #overlay the contours on the map, display
#    Result(contour,[dem,contour],'Overlay')

End()