Fortran 2003 中的非标准类型声明答案

【问题标题】：Nonstandard type declaration in Fortran 2003Fortran 2003 中的非标准类型声明
【发布时间】：2015-03-04 12:39:38
【问题描述】：

我编写了一个科学的 fortran 代码，没有使用任何特定的 fortran 标准。但我现在必须声明我使用的是哪个 fortran 标准。

我说我使用的是 fortran 2003，因为我需要 get_command_argument 和 command_argument_count 内部函数。但是，当使用标志 -std=f2003 检查代码标准时，编译会失败。

我在某些部分收到关于实数类型声明的错误。例如当我在模块中声明变量时：

module innout
implicit none
real*8,parameter               :: nan=-1.
real*8,allocatable,save        :: windU(:),windV(:)
real*8,allocatable,save        :: input_param(:,:),input_rad(:,:)
real*8,allocatable,save        :: prein(:),input(:),ref_lev(:)
character(30),allocatable,save :: sceneclass(:)
end module innout

我在所有实变量中都收到消息“非标准类型声明 REAL*8”。

有人知道发生了什么吗？

【问题讨论】：

尝试改用real64。
谢谢，但我认为这是 Fortran 2008
始终使用标签fortran 并且仅在必要时添加版本以区分您的问题是具体的。例如，您不能使用 Fortran 2008，而只能使用 Fortran 90。

标签： fortran standards variable-declaration fortran2003

【解决方案1】：

real*8 不是，也从来不是 Fortran 标准类型声明。如今，声明 64 位实数的最简单方法可能是从内部模块 iso_fortan_env 导入 named-constant real64，如下所示：

use, intrinsic :: iso_fortran_env
...
real(real64) :: my_var

还有其他方法，涉及selected_real_kind和其他机制，但是如果您想使用IEEE浮点类型进行编程，那么real64和real32是一个不错的选择。

正如@AlexanderVogt 在评论中指出的那样，这些标准命名常量已添加到 2008 年标准中的语言中。我使用过的最新编译器版本已经实现了它们。

【讨论】：

REAL64 在 Fortran 2003 中不适用于 ISO_Fortran_env！命名常量是在 2008 年添加的。
好吧，好吧，好吧，我以为它们是在 2003 年添加的。
见"The new features of Fortran 2008," Ch. 10.12, "Named Constants" ;-)

【解决方案2】：

kind 说明符是要走的路...

如果你限制自己使用 Fortran 2003 Standard，那么你需要先使用kind() 或selected_real_kind() 函数来确定对应的种类：

module innout
  implicit none
  integer,parameter                  :: REAL64 = kind(1.d0)

  real(kind=REAL64),parameter        :: nan=-1._REAL64
  real(kind=REAL64),allocatable,save :: windU(:),windV(:)
  real(kind=REAL64),allocatable,save :: input_param(:,:),input_rad(:,:)
  real(kind=REAL64),allocatable,save :: prein(:),input(:),ref_lev(:)
  character(30),allocatable,save     :: sceneclass(:)
end module innout

如果您被允许/您的编译器支持 Fortran 2008，我会推荐模块 ISO_Fortran_env 和预定义常量 REAL64：

module innout
  use,intrinsic :: ISO_Fortran_env, only: REAL64
  implicit none
  real(kind=REAL64),parameter         :: nan=-1._REAL64
  real(kind=REAL64),allocatable,save  :: windU(:),windV(:)
  real(kind=REAL64),allocatable,save  :: input_param(:,:),input_rad(:,:)
  real(kind=REAL64),allocatable,save  :: prein(:),input(:),ref_lev(:)
  character(30),allocatable,save      :: sceneclass(:)
end module innout

【讨论】：

或者，为了尴尬，integer, parameter :: real64 = kind(1.d0) 和迁移到 2008 标准会容易得多:-)
我已经检查过了，我也可以使用 2008 标准，所以我将使用该选项。谢谢！！

【解决方案3】：

我已经找到答案了。

在变量声明中使用以下内容，似乎可以正常工作：

integer, parameter :: dp = selected_real_kind(15, 307)

来自http://fortranwiki.org/fortran/show/Real+precision

【讨论】：