fortran-时间的格式化输出

1 声明部分

    implicit none    real :: read_null               # 用于读取时跳过输入文件中不需要的列    integer :: leap_year            # 用于判断是否为闰年。 0是平年；1是闰年    integer,allocatable :: year(:), month(:), day(:), hour(:)       #  时间间隔为1小时debug

2 输出时间格式生成

 # 输出的时间格式生成，起始时段为净雨量文件中的起始时段    allocate(year(size(Q)),month(size(Q)),day(size(Q)),hour(size(Q)))  ! 数组长度等于径流的长度    ! 读取初始时间    filename = trim(Address)//'\'//'data_rs.txt'    open(10,file=filename,access='sequential',status='old')        read(10,*) year(1),month(1),day(1),hour(1)    close(10)    do t=2,size(Q)        leap_year = 0        ! 一般情况下，“小时”+1，而“日，月，年”不变，把它们放在后边处理        hour(t) = hour(t-1)+1        day(t) = day(t-1)        month(t) = month(t-1)        year(t) = year(t-1)        !判断是否为闰年，确定二月的天数        if(mod(year(t),100)==0) then            if(mod(year(t),400)==0) then                leap_year = 1            end if         else            if(mod(year(t),4)==0) then                leap_year = 1            end if        end if        ! 24时统一记为第二天0时        if(hour(t) == 24) then            hour(t)=0            day(t) = day(t)+1        end if        !超出一个月的最大天数，则计入下一个月        select case(month(t))        case(1,3,5,7,8,10,12)            if (day(t)==32)then                day(t) = 1                month(t)= month(t)+1            end if         case(4,6,9,11)            if(day(i)==31)then                day(t) = 1                month(t)= month(t)+1             end if         case(2)            if(leap_year==1 .AND. day(t)==30)then                day(t) = 1                month(t)= month(t)+1             elseif(leap_year==0 .AND. day(t)==29)then                day(t) = 1                month(t)= month(t)+1             end if        end select    end do

3 输出

     filename=trim(Address)//'\'//'output_path'     open(10,file=filename)     !write(10,"(5A6)") "year" ,"month", "day", "hour" ,"Q"     do i=1,size(Q)      write(10,"(4I5,F10.3)")year(i),month(i),day(i),hour(i),Q(i)     end do     close(10)