Fortran命令行输入规范部分

Question

我是fortran的新手，在尝试通过命令行传递参数时遇到问题。例如，我的工作代码有以下几行：

!experimental parameters
real (kind=8), parameter :: rhot                =1.2456
!density of top fluid
real (kind=8), parameter :: rhob                = 1.3432
!density of bottom fluid
real (kind=8), parameter :: rhof                = (rhot-rhob)
!generic function of rhot rhob

我想通过上交

./fulltime_2009_05_15_Fortran9 [Value rhot] [Value rhob]

使用类似的东西

call get_command_argument(1,input1)
call get_command_argument(2,input2)
read(input1,*) rhot
read(input2,*) rhob

问题是，我声明了像rhof这样的参数，这些参数取决于输入的值。因此，我希望立即应用用户输入的值，以便所有相关参数都可以使用这些值。但是，如果将我的代码修改为：

real (kind=8) :: rhot,  rhoB
call get_command_argument(1,input1)
call get_command_argument(2,input2)
read(input1,*) rhot
read(input2,*) rhob
real (kind=8), parameter :: rhof    = (rhot-rhob)

我得到错误：规范语句不能出现在可执行部分中。

我对如何解决这个问题有什么想法或建议吗？

Answer 1

命令行参数无法更改编译时常数(parameter(。它在编译时是固定的。

此代码：

real (kind=8) :: rhot,  rhoB
real (kind=8) :: rhof
call get_command_argument(1,input1)
call get_command_argument(2,input2)
read(input1,*) rhot
read(input2,*) rhob
rhof    = (rhot-rhob)

会编译得很好。但是不能在一个普通语句之后有一个变量声明。

你可以有

real (kind=8) :: rhot,  rhoB
call get_command_argument(1,input1)
call get_command_argument(2,input2)
read(input1,*) rhot
read(input2,*) rhob
block
real (kind=8) :: rhof    = (rhot-rhob)

在Fortran 2008中，但是不能使用非常量表达式定义编译时常量。

某些编程语言确实允许在运行配置期间设置一类常量，然后再进行固定。Fortran不是其中之一。