Integer, Parameter :: Nmon=12, Yr1=1982, Yr2=2009, Nyrs=(Yr2-Yr1)+1
      Integer, Parameter :: Ntime=Nyrs*Nmon, LT=6
      Real, Parameter :: Miss=-99.0, missing_value = 9.999e+20
      Real, Parameter :: Dlon=1.0, Dlat=1.0, Lon1=60.0, Lon2=290.0, Lat1=-20.0, lat2=60.0

      Integer :: Nst,Nlon,Nlat, ilon,ilat,ist, ilt
      Integer :: ix1, ix2, iy1, iy2
      Real :: x1,x2,y1,y2

      ilon1=1; ilat1=1

      x1=130; x2=150; y1=20; y2=30

      ix1=Int((x1-lon1)/dlon)+1
      ix2=Int((x2-lon1)/dlon)+1
      iy1=Int((y1-lat1)/dlon)+1
      iy2=Int((x2-lat1)/dlon)+1

      print*,x1,x2,y1,y2
      print*,ix1, ix2, iy1, iy2

      End