* B Vs -VTL

"reinit"
"open /tasdevil/s2/arnottj/research/scripts/dummy_scripts/dummy.ctl"
"set display color white"
"clear"
"set parea 1 10 1 7.5"
*"set parea 1 8.5 1 7.5"

"set string 1"
"set strsiz 0.13 0.13"

"set lat 0 20"
"set lon 260 360"


"set clevs 500"
"set xlopts 1 5 0.2"
"set ylopts 1 5 0.2"
"set grads off"
"d lat"


*Change for desired file type
file ="2002SATL.txt"                                                                                                                                    

colors="1 2 3 4 5 6 7 8 9 10 11 12 13"
ccount=0

rc=1
while (rc != 2)
   data=read(file)
   rc=sublin(data,1)
   if (rc != 2)
      line=sublin(data,2)
      chk=substr(line,1,10)
      if (chk != "----------") 

        lat=subwrd(line,1)
        lon=subwrd(line,2)
       flag=subwrd(line,3)
      color=subwrd(line,4)


        "q w2xy "lon " " lat
         xloc=subwrd(result,3)
         yloc=subwrd(result,6) 
         name=subwrd(ccount,1)
          
          "set line "color" 1 5"
          "draw mark "flag " "xloc  " "yloc  " 0.1"
*           "draw mark 1 "xloc  " "yloc  " 0.1"
      else
         ccount=ccount+1
*         color=subwrd(colors,ccount)
      endif
   endif
endwhile

'draw title 2002 S. Atlantic Tropical Cyclogenesis Points'
"printim 2002SATL.gif"

*"enable print output.gm"
*"print"
*"disable"
*"!gxeps -c -i output.gm -o GEOCLUS.ps"
*rc=close(file)

