通过命令行参数设置数组大小

Question

I know that in Fortran I can declare an array of a number N of elements if N is declared as a parameter before (and so defined to some value).

On the other hand I know I can get a program to accept arguments from the command line with the use of the subroutine get_command_arg .

My question is:

Can I somehow declare an array of a number of elements given by a command when calling the program from the command line?

I'm looking for something like the command line:

./main -30 

where main.f03 would begin with something like:

integer, parameter :: N = get_command_arg(1)
real :: x(N) ...

I'm trying not to define the arrays as allocatable.

Answer 1

There are certain languages where you can initialize variables and named constants from command line (like chapel), but Fortran is not one of them.

You say "I am trying not to define the arrays as allocatable." but that is the problem. You simply have to. There is no other way.

Your code

integer, parameter :: N = get_command_arg(1)

real :: x(N) ...

is illegal for several reasons.

You cannot put GET_COMMAND_ARGUMENT() into a constant expression because it is not among allowed functions. It does not return compile-time constant values. And parameter initializers must be set at the compile time.

GET_COMMAND_ARGUMENT() is a subroutine, not a function. It can return more stuff, not just the value, but also the length and status. It is not pure and the Fortran standard is trying to use only pure functions. Other things, like RANDOM_NUMBER() , are subroutines. It is a good style to follow in your own programs too.

The only way in Fortran to create arrays that change from run to run is to make the array allocatable or pointer . There are also automatic arrays for local arrays.

Answer 2

You can do it this way without allocatables. You just have to pass the size (after converting it to an integer) into a subroutine. But really, I see no reason not to use allocatables for something like this.

program main

implicit none

integer :: n,arg_len,istat
character(len=100) :: arg

call get_command_argument(1,value=arg,status=istat)
if (istat/=0) error stop 'error: cannot read first arg'

read(arg,'(I100)',iostat=istat) n
if (istat/=0) error stop 'error: first arg not an integer'

call real_main(n)

contains

    subroutine real_main(n)
    integer,intent(in) :: n
    integer,dimension(n) :: ival
    ival = 1
    write(*,*) ival
    end subroutine real_main

end program main

Example use:

> ./main 1
       1
> ./main 2
       1           1
> ./main 3
       1           1           1