技术

一个Fortran 读文件并带参数的模板

Fortran 作为科学计算常用语言，自然免不了各种文件读写的操作，以下是一个简单的模板

参见 gist bdb32fffcd08b9ab21ea:

program gaugeselect
implicit none
character(len=80):: filename
character(len=80):: arg1 ,arg2
character(len=80):: buff
logical:: flag=.false.
!================================================================!
call getarg(1,arg1)
call getarg(2,arg2)
if ( len_trim(arg1)==0 )then
print *,"please specify the argument of running file!!"
print *," usage: ./gaugeselect filename "
stop
else
filename=trim(arg1)
inquire(file=trim(filename),exist=flag)
if(flag)then
open(5,file=trim(filename))
do while(.true.)
read(5,*,end=100)buff
enddo
100 close(5)
print *,'file: ',trim(filename), " processed"
else
print *,"oops........"
print *,'can NOT open file: ',trim(filename)
print *," "
endif
endif
end program gaugeselect
view raw fortranargs.f90 hosted with ❤ by GitHub

或者，直接参见下面：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
program gaugeselect
implicit none
character(len=80):: filename
character(len=80):: arg1 ,arg2
character(len=80):: buff
logical:: flag=.false.
!================================================================!
call getarg(1,arg1)
call getarg(2,arg2)
if (len_trim(arg1)==0 )then
  print *,"please specify the argument of running file!!"
  print *,"usage: ./gaugeselect filename"
  stop
else
  filename=trim(arg1)
  inquire(file=trim(filename),exist=flag)
  if(flag)then
    open(5,file=trim(filename))
    do while(.true.)
      read(5,*,end=100)buff
    enddo
    100 close(5)
    print *,'file:',trim(filename), "processed" 
  else 
    print *,"oops........"
    print *,'can NOT open file:',trim(filename) 
    print *," "
  endif
endif
end program gaugeselect