读无格式文件’1948.dat’, ‘1949.dat’, …, ‘2013.dat’,把每个月的数据单独写成有格式顺序文件,以月份的前三个字母命名,分别为’Jan.dat’, ’Feb.dat’, …, ’Dec.dat’, 每个文件2列66行,第一列为年份,第二列为温度数据。(要求:写12个文件时,要使用循环结构,不要复制12段重复性代码)。
program main
implicit none
integer::i,s=1948
character(len=5)::dirname
character(len=8)::filename
real,dimension(12)::a
open(10,file='TI.dat',form='unformatted',access='direct',recl=48)
do i=1,66
read(10,rec=i)a
print*,a
write(dirname,'(i4,"s")')s
write(filename,'(i4,".dat")')s
call system("mkdir "//dirname)
open(s,file=trim(dirname)//'/'//filename,form='unformatted',access='direct',recl=48)
write(s,rec=1)a
close(s)
s=s+1
enddo
endprogram
character(len=7),dimension(12)::monname
do i=1,66
write(dirname,'(i4,"s")')dir
write(filename,'(i4,".dat")')dir
open(i,file=trim(dirname)//'/'//filename,form='unformatted',access='direct',recl=4)
do j=1,12
read(i,rec=j)a(i,j)
enddo
close(i)
dir=dir+1
enddo
do j=1,12
monname=['Jan.dat','Feb.dat','Mar,dat','Apr.dat','May.dat','Jun.dat','Jul.dat','Aug.dat','Sep.dat','Oct.dat','Nov.dat','Dec.dat']
open(k,file='/'//monname(j),status='new')
do i=1,66
write(k,'(i4,f4.2)')s,a(i,j)
s=s+1
enddo
close(k)
k=k+1
s=1948
enddo
endprogram
1.PNG (73.41 KB, 下载次数: 326)
2.PNG (23.79 KB, 下载次数: 316)
fcode 发表于 2020-6-12 08:05
文件权限有问题,检查文件夹及文件的归属权。不行的话可以设置成 chmod 666
fcode 发表于 2020-6-13 08:00
我不清楚在你的计算机上发生了什么。但是现象来看,是这个原因。
论坛里没有神仙,都是凡人,我们不会掐指 ...
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |