subroutine bassel(rlat,rlon,ibasin)
c     This routine determines which basin a given 
c     lat/lon point is in
c
c     Input: 
c       rlat - latitude deg N
c       rlon - longitude deg E (0-360 convention)
c 
c     Output
c       ibasin = 1 for Atlantic
c              = 2 for east Pacific
c              = 3 for west Pacific
c              = 4 for southern hemisphere
c              = 5 for N. Indian Ocean 
c
c     Modified June 2013 by M. DeMaria to change definitions of ibasin=4 and 5
c
      ibasin = 1
      if     (rlat                      .ge.  0.0) then
         if     (rlat                      .le.  8.0) then
            if (rlon .lt.  285.0) ibasin=2
         elseif (rlat .gt.  8.0 .and. rlat .le. 15.0) then
            if (rlon .lt.  275.0) ibasin=2
         elseif (rlat .gt. 15.0 .and. rlat .le. 17.5) then
            if (rlon .lt.  270.0) ibasin=2
         elseif (rlat .gt. 17.5                     ) then
            if (rlon .lt.  260.0) ibasin=2
         endif
      endif
c
      if ((rlon .lt. 180.0) .and. (rlat .ge. 0.0))  ibasin=3
c
      if (                        (rlat .lt. 0.0))  ibasin=4
c
      if ((rlon .lt. 100.0) .and. (rlat .ge. 0.0))  ibasin=5
      
      return
      end