Subroutine Create_file_for_Fcs_DLap (ilt,Out_Dir) !(Lag, Kmon) !------------------------------------------------- !...Structure of the input data file as following (ASCII file): ! Row 1: Title of the file (CHARACTER) ! Row 2: Two positive integer number: N (Number of samples), Mbien (Number of variables) ! Row 3: Mbien values of IDX arrays (IDX(Mbien)) ! Row 4: Name of variables (CHARACTER*15) ! From row 5 until the end of file (N rows), each row contains Mbien values (REAL) Integer, Allocatable :: IDX(:) !=(/3,0,0,0,0,0,0,0,0,0/) ! 1 bien phu thuoc, 9 bien doc lap (3 months x 3 EOF) Integer :: ieigen, imon, iyr, isrc, i, Mvar, ilt Character (Len=7) :: EOFName(1:NMon,1:Nmode) Character (Len=255) :: Fname, Fname1, Fname2 Character (Len=10) :: vnametc Character (Len=*) :: Out_Dir vnametc=trim(adjustl(vartc)) print*,"In Create Files: Nmode,NMon,Nt= ",Nmode,NMon,Nt Do ieigen=1, Nmode Do imon=1, NMon write(EOFName(imon,ieigen),"(A3,I1,A1,I2.2)") "EOF",ieigen,"_",imon Enddo Enddo Do imon=1, NMon !Mvar=Mvar1(imon) write (Fname,"(A,A,I2.2,A)") trim(Out_Dir)//trim(vnametc),"_",imon,".txt" write (Fname1,"(A,A,I2.2,A)") trim(Out_Dir)//"Ann_Out_"//trim(vnametc),"_",imon,".txt" write (Fname2,"(A,A,I2.2,A)") trim(Out_Dir)//"Ann_Save_"//trim(vnametc),"_",imon,".txt" Write (33,"(A)") trim(Fname); Write (44,"(A)") trim(Fname1); Write (55,"(A)") trim(Fname2) Open (3,file=trim(Fname), Status="Unknown") write (3,"(A,I2.2,A)") "Total TC count of month ", imon, " and Predictors" write (3,"(3I6)") Yr1+1, Nt, Nmode+1 write (3,"(20I3)") 3, (0,ivar=1,Nmode) !(IDX(i), i=1,Mvar) !write (3,"(20A15)") "TC_"//MonName(imon),(EOFName(mon_x(k),:,isrc),k=1,Kmon) !write (3,"(20A10)") trim(vnametc),(trim(vname(imon,ivar)),ivar=1,Mvar) write (3,"(20A10)") trim(vnametc),(trim(EOFName(imon,i)),i=1,Nmode) Do iyr=1,Nt ! EOFall(LT,Nmon,Nmode,Nt) write (3,"(20F10.4)") TC(iyr+1,imon), EOFall(ilt,imon,iyr,1:Nmode) Enddo Close (3) Enddo Close (33); Close (44); Close (55) End Subroutine Create_file_for_Fcs_DLap