#------------------------------------------------------------------------------------------------
# Create canvas ---
canvas = vcs.init(geometry=(1200,800))
canvas.open()
template = canvas.createtemplate()
template.blank(['title','mean','min','max','dataname','crdate','crtime','units']) ## Turn off additional information
#template.list() ## This commend could give list of items 

# Set ploting range ---
lat1=-60
lat2=80
lon1=0
lon2=360

# Wind speed ---
data['wnd'] = MV.sqrt(data['ua']**2+data['va']**2)
iso = canvas.createisofill()
iso.datawc_x1 = lon1
iso.datawc_x2 = lon2
iso.datawc_y1 = lat1
iso.datawc_y2 = lat2
canvas.setcolormap('blue_to_orange')
canvas.plot(data['wnd'],iso,template)

# Geopotential height field ---
lines1 = vcs.createisoline()
lines1.datawc_x1 = lon1
lines1.datawc_x2 = lon2
lines1.datawc_y1 = lat1
lines1.datawc_y2 = lat2
lines1.linecolors = ['black']
canvas = vcs.init(geometry=(1200, 800))
canvas.open()
template = canvas.createtemplate()
template.blank(
    ['title', 'mean', 'min', 'max', 'dataname', 'crdate', 'crtime',
     'units'])  ## Turn off additional information
#template.list() ## This commend could give list of items

# Set ploting range ---
lat1 = -60
lat2 = 80
lon1 = 0
lon2 = 360

# Wind speed ---
data['wnd'] = MV.sqrt(data['ua']**2 + data['va']**2)
iso = canvas.createisofill()
iso.datawc_x1 = lon1
iso.datawc_x2 = lon2
iso.datawc_y1 = lat1
iso.datawc_y2 = lat2
canvas.setcolormap('blue_to_orange')
canvas.plot(data['wnd'], iso, template)

# Geopotential height field ---
lines1 = vcs.createisoline()
lines1.datawc_x1 = lon1
lines1.datawc_x2 = lon2
lines1.datawc_y1 = lat1
lines1.datawc_y2 = lat2
lines1.linecolors = ['black']