Use NetCDF
      Implicit None

      Real, parameter :: Miss=-99.0
      Character (Len=255) :: Fname
      Integer :: Nst, Nidx, ist, imon, i

      Real, Allocatable :: Lon(:), Lat(:), HWday(:,:)
      Real, Allocatable :: StnID(:), Xmon(:)
      Character (Len=8), Allocatable :: Station(:)

      Nidx=1; Nst=58
      Allocate (Lon(Nst), Lat(Nst), HWday(Nst,Nidx), Station(Nst))
      Fname="cluster_4groups.txt" !"cluster_3groups.txt"
      Open (1, file=Fname,Status="old")
      Read (1,*)
      Do ist=1,Nst
         Read (1,*) i,Lon(ist),Lat(ist),HWday(ist,Nidx)
         write (Station(ist),"(A,I2.2)") "Stn_",ist
      Enddo
      Close (1)
      Call Write_GrADS_Stn (HWday, Lon, Lat, Station, Nst, Nidx, &
      !     "3_clus", "./", "clus", Miss)
           "4_clus", "./", "clus", Miss)

      CONTAINS

      Include "/home/tanpv/Mylibs/Sub_Write_GrADS_Stn.f90"

      END