;*******************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;*******************************************

;***********************************************
;read in the data

do j=1980,1981

data=  readAsciiTable("~/Documents/NCL_files/Storm_tracks/data/OND_1979_2014/OND_"+j+".txt", 5, "float", 0)

time= data(:,0)      ; storm time
cnt = data(:,1)     ; storm number
lon = data(:,2)     ; longitude
lat = data(:,3)     ; latitude 

t=dimsizes(cnt)

sn=cnt(0)-cnt(t-1)

count_all=new(35,"float")
count=0
do i=0,sn

ix=ind(time.eq.274.)   
x=cnt(ix(0))+i

y=ind(cnt.eq.x)

lat_st=lat(y(0))
lat_end=lat(dimsizes(y))
lon_st=lon(y(0))
lon_end=lon(dimsizes(y))

if (lat_end.gt.lat_st.and.lat_st.ge.30.and.lat_end.ge.60) then
       count=count+1
end if

count_all=count

print(y)

delete(y)
delete(lat_st)
delete(lat_end)
delete(lon_st)
delete(lon_end)

end do
  
 delete(data) 
 delete(time)
 delete(cnt)
 delete(lon)
 delete(lat)

end do