'reinit' cyear = '2010' data = '_ERA' fname = 'TC_freq_' cyear data 'sdfopen ' fname '.nc' 'set parea 1 10.5 0.4 8.5' 'set display color white' 'set grid off' 'c' 'set mpdset hires' 'set xlint 5' 'set xlopts 1 5 0.15' 'set ylopts 1 5 0.15' 'set ylint 5' 'set gxout shaded' 'color 0 6 1 -kind white->blue->cyan->yellow->red' 'd tc_freq' *'set gxout grid' *'set ccolor 0' *'set dignum 0' *'d tc_freq' 'draw title ERA_TC_' cyear 'cbar' 'printim ' fname '.png'