!## Reformat data to apply in R-CDF
!## 15-07-2016
PROGRAM Conv_2_column
   implicit none
   integer, parameter   :: NST = 151, SYR = 2012, EYR = 2014, NYR = EYR - SYR + 1, &
                           NMN = 12
   integer              :: ist, iyr, imn, itmp, NDY, idy
   character (4)        :: cyr, cmn, cdy
   character (15)       :: stname(NST)
   character (100)      :: filenm, ofile
   character (500)      :: cmd
   real                 :: dat(NST,NYR,NMN,31), rtmp

   !## Read station name
   open (1, file = "lonlat.info",status = "old")
      do ist = 1, NST
         read (1,*) stname(ist), rtmp, rtmp
      enddo
   close (1)

   !## Read data
   dat = -99.
   NDY = 31
   do ist = 1, NST
      filenm = "Daily_Converted/output/R_"//trim(stname(ist))//".txt"
      open (2,file = filenm)
         do iyr = 1, NYR
            read (2,*)
            do idy = 1, NDY
               read (2,*) itmp, (dat(ist,iyr,imn,idy), imn = 1, NMN)
            enddo
         enddo
      close (2)
   enddo

   !## Rewrite data
   do iyr = 1, NYR
      itmp = SYR + iyr - 1
      write (cyr,'(i4)') itmp
      do imn = 1, NMN
         write (cmn,'(i2)') imn
         if (imn < 10) write (cmn,'A,(i1)') "0",imn
         ofile = "Daily_Converted_Column/"//trim(cyr)//trim(cmn)//"/"
         cmd = "if [ -d "//trim(ofile)//" ]; then echo ...; else mkdir -p "//trim(ofile)//"; fi" 
         Call system(cmd)
         do ist = 1, NST
            filenm = trim(ofile)//"R_"//trim(stname(ist))//".txt"
            open (3, file = filenm)
               NDY = cNDAY(itmp,imn)
               do idy = 1, NDY
                  write (3,'(F8.3)') dat(ist,iyr,imn,idy)
               enddo
            close (3)
         enddo
      enddo
   enddo
CONTAINS
include "/work/users/thanhnx/Func_Sub/func_cNDAY.F90" 
END