do kblock = 1, nblock c data output and output frequence unit=1000+kblock outdir='' c 输出时间间隔 outputstep=0.02 noutputstep=period/outputstep c 工作目录调用 CALL VGETOUTDIR( OUTDIR, LENOUTDIR ) do i =1,noutputstep time_tol=abs(i*outputstep-time(1)) if ((time_tol).gt.(0.0001)) then continue else c 按单元号生成输出文本 write(cfile,'(i8.5,a4)') jelem(kblock),'.txt' cfile=TRIM(ADJUSTL(outdir))//'\'//TRIM(ADJUSTL(cfile)) open(unit, file=cfile, action="WRITE") write(unit,*) time(1) endif enddo end do
来源:https://www.cnblogs.com/structurer/p/11793129.html