Use NetCDF
      Implicit None

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

      Real, Allocatable :: Lon(:), Lat(:), HWday(:,:)
      Real, Allocatable :: StnID(:), Xmon(:)

      Fname="HW.dat"
      Open (1, file=Fname, Status="Old")
      Read (1,*)
      Read (1,*) Nst,Nmon
      Allocate (Lon(Nst),Lat(Nst),HWday(Nst,Nmon))
      Do ist=1,Nst
         Read (1,*) Lat(ist), Lon(ist), (HWday(ist,imon),imon=1,Nmon)
      Enddo

      Allocate (StnID(Nst), Xmon(Nmon))
      Do ist=1,Nst; StnID(ist)=ist; Enddo
      Do imon=1,Nmon; Xmon(imon)=imon; Enddo
      Fname="HWday.nc"
      Call wr_ncdf_2d (Fname,StnID,Xmon,Nst,Nmon,HWday,"hwday", &
           "HW days","day/mon","stn","mon",Miss)
      Call wr_ncdf_1d_apd (Fname,Lon,"lon",Nst,"Longitude","deg-E","stn",Miss)
      Call wr_ncdf_1d_apd (Fname,Lat,"lat",Nst,"Latitude","deg-N","stn",Miss)
 
      CONTAINS
      Include "/home/tanpv/Mylibs/sub_wr_ncdf_Tan.F90"
      END