!---------------------------------------------------------------------
   Subroutine Write_GrADS_Stn_3D (Varnm,X,Nst,Nvar,Ntime,xlon,xlat,Stname,Grads_File) ! (Nst,Nvar,Ntime)
     Character (Len=*) :: Grads_File
     !Real :: X(Nst,Nidx,Nlen,ndr_idx) !Nst,Nmon,NLen,ndr_idx) ! (Nst,2,Nvar)
     Real,    intent (in) :: X(:,:,:)           ! Nst,Nvar,Ntime)
     Real,    intent (in) :: xlon(:), xlat(:)   ! Nst
     Character (Len=*),intent (in) :: Stname(:) ! Nst
     Character (Len=*),intent (in) :: Varnm(:)  ! Nvar
     Integer, intent (in) :: Nst,Nvar,Ntime
     Integer :: ivar, ist, it, ic
     Character (Len=8) :: STID
     Real :: RLAT, RLON, TIM, tmp, Var
     Integer :: NLEV,NFLAG
     !----------------------------------------------------------------------------------
     Open (10,file=trim(Grads_File)//'.dat',FORM='UNFORMATTED',RECORDTYPE='STREAM')
     TIM = 0.0   ! Co the bien doi tu -0.5 .. +0.5; thuong dat bang 0.
     NFLAG = 1   ! Khong doi
        Do it = 1,Ntime   ! tdef 
           NLEV = 1       ! =0 o ban ghi cuoi cung cua 1 timeslice; =1 cho tat ca cac timeslice khac
           ic = 0
           Do ist=1, Nst
              ic = ic+1
              STID = Stname(ist)(1:8)   ! Chi lay 8 ky tu dau tien
              RLAT = xlat(ist); RLON = xlon(ist)
              WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG
              Do ivar=1,Nvar
                 Var = X(ist,ivar,it)
                 WRITE (10) Var   ! Gia tri cua bien
              Enddo
           Enddo
           NLEV = 0
           WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG
        Enddo   ! ivar
     Close(10)
     !---------------------------------------------------------------------!
     !     Write control file Station data                                 !  
     !---------------------------------------------------------------------!
      open  (2001,file=trim(Grads_File)//'.ctl')
      write (2001,'(2A)') 'dset  ','^'//trim(Grads_File)//'.dat'
      write (2001,'(A)')  'dtype station'
      write (2001,'(2A)') 'stnmap ','^'//trim(Grads_File)//'.dat'//'.map'
      write (2001,'(A,F8.1)')  'undef   ',Miss
      write (2001,'(A)')  'title  undefine'
      write (2001,'(A,I3,1x,A)') 'tdef ',Ntime,' linear jan1980 1mo'
      write (2001,'(A,I3)')  'vars  ',Nvar
      Do ivar=1,Nvar
         write (2001,'(A,A)')  trim(Varnm(ivar)),'  0 99 Drough Chars'
      Enddo
      write (2001,'(A)')  'endvars'
      close (2001)
   End Subroutine Write_GrADS_Stn_3D
   !---------------------------------------------------------------------