- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
所以我一直在检查状态方程模拟器值的一些计算问题。经过一些检查后,我发现模块 constantes 没有将值拖到它开头声明的子例程中。我不想每次使用它们时都声明这些变量,我该怎么办?声明一个带有值的子例程,并在每次需要参数时“使用”它?
MODULE Constantes
Implicit None
!Lista de constantes para funciones a y b usadas para calcular Z
!Real,dimension(4)::y
!Real::T,P
Real,Parameter,Dimension(4)::Tc=(/ 190.56, 132.85, 304.12, 126.2 /)
Real,Parameter,Dimension(4)::Pc=(/ 45.992, 34.94, 73.74, 33.98 /)
Real,Parameter,Dimension(4)::w=(/ 0.011, 0.045, 0.225, 0.037 /)
Real,Parameter,Dimension(4)::na=(/ 0.072, 0.608, 0.974, 25765.2 /)
Real,Parameter,Dimension(4)::ma=(/15.68, 9.62, 34.04, 0.00 /)
Real,Parameter,Dimension(4)::ba=(/ 6.48, 66.2, 14.04, 0.000928699 /)
Real,Parameter,Dimension(4)::da=(/ 0.025, 0.03, 0.061, 0.00 /)
Real,Parameter,Dimension(4)::ha=(/30.56, 48.16, 34.67, 30.0 /)
Real,Parameter,Dimension(4)::ua=(/ 13.45, 14.84, 20.79, 0.0 /)
Real,Parameter,Dimension(4)::ka=(/ 1.0, 1.0, 1.0, 1.13 /)
Real,Parameter::c0=0.48,c1=1.574,c2=-0.176,eps=1,sig=0
Double Precision,Parameter::omega=0.08664035,psi=0.42748023
!Variables input para constante A
!Variables input para A y B
Contains
Subroutine CA(y,T,At,ap,aij,D,n)
Implicit None
Double precision::psi
Double precision,intent(out)::At,D
Double Precision,Dimension(4),intent(out)::ap
Double Precision,Dimension(4,4),intent(out)::aij
Real,Dimension(4)::Tc,Pc,w
Real::c0,c1,c2
Double precision,intent(in)::T,n
Double precision,dimension(4),intent(in)::y
Integer::i,j
Do i=1, 4 !Array de constantes a individuales
ap(i)=Tc(i)*Tc(i)
!ap(i)=(psi*(83.14*83.14*Tc(i)*Tc(i))/Pc(i))!*((1+((c0+c1*w(i)+c2*w(i)*w(i))*(1-sqrt(T/Tc(i)))))**(2))
End do
Do i=1,4
Do j=1,4
aij(i,j)=y(i)*y(j)*sqrt(ap(i)*ap(j))!Constantes duales de a
End do
End Do
At=sum(sum(aij,dim=1)) !A total independiente de n
D=At*n*n !A total dependiente de n
End Subroutine CA
Subroutine CB(y,Bt,bi,B,bij,n)
Implicit none
Real,Dimension(4)::Tc,Pc
Double precision, intent(in)::n
Double Precision,Dimension(4),intent(out)::bi
Double precision,dimension(4),intent(in)::y
Double precision::omega,bm
Double precision,intent(out)::Bt,B
Double precision,Dimension(4,4),intent(out)::bij
Integer::i,j
Do i=1,4 !Array de constantes b individuales
bi(i)=omega*83.14*Tc(i)/Pc(i)
End Do
Do i=1,4 !Array constantes duales de b
Do j=1,4
bij(i,j)=0.5*(bi(i)+bi(j))
End Do
End Do
BT=sum(bi*y) !B total independiente de n
B=n*BT !b total dependiente de n
End Subroutine CB
END MODULE Constantes
最佳答案
您遇到的一个问题是您在子例程中重新声明(至少部分)参数数组。例如,行
Real,Dimension(4)::Tc,Pc,w
屏蔽模块变量的存在,例如
Real,Parameter,Dimension(4)::Tc=(/ 190.56, 132.85, 304.12, 126.2 /)
子例程已经可以通过主机关联访问这些参数数组,在这种情况下,在子例程中再次声明它们是错误的。
有些人会争辩说,即使不是绝对必要,您也应该通过参数列表将所有参数传递给过程。
关于arrays - FORTRAN 数组计算错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/42722243/
是的,我知道..,这不是想象的...这是一个真正的 Fortran 问题。 以前的版本是指 Fortran 2003、95、90,甚至 77。 我所说的“向后兼容”是指可以轻松运行为 2008 年以前
我有一个程序,它的变量中有一个值。一旦确定了该值,我想调用另一个程序并使用该变量的值来确定在新程序中的位置。有人知道该怎么做吗? 最佳答案 如果您有 Fortran 2008 编译器,您将拥有标准子例
namelist 是一种有用的 fortran 结构,可以从文件中快速初始化变量。 namelist 有一个名称并包含一组具有已知类型的变量。这使得它类似于 type 结构。 通常情况下,给程序或子例
我正在遍历索引,我正在检查我是否不在第一个循环交互和另一个条件中。如果第一个条件是 .False.,我不想评估第二个条件。 do i = 1, n if ( i /= 1 .and. var(
Fortran 2003 具有用于数组连接的方括号语法,Intel fortran 编译器也支持它。我在这里为矩阵连接写了一个简单的代码: program matrix implicit none r
我正在尝试通过重载类型名称来制作自定义数据类型构造函数。但是,在进行调用时,将调用默认构造函数。我不明白我做错了什么。 这是有问题的代码片段。 module test type, pu
我的最终目标是在 Fortran 中有一个通用的映射函数,即一个接受任意类型 A 的数组和一个 A->B 类型的函数的函数,将此函数应用于给定数组的所有元素并返回一个B 类型的数组。我无法用数组实现它
我正在学习 Fortran,在使用格式编写时发现了一些奇怪的东西(我使用的是 Fortran onlinegdb) Program Hello real, dimension(3,2):: array
Fortran 中的INTERFACE 语句是否使其成为正式实现multiple dispatch 的编程语言? ? (我问是因为所链接的维基百科文章在其看似全面的支持相关范式的示例编程语言列表中并未
我可以使用 Fortran 95 编译器编译 Fortran 90 文件吗? Fortran 95 似乎有很多,但 Fortran 90 没有。 最佳答案 这个可以: NAGWare f95 Comp
嗨,我在 Fortran 中对二维离散化问题强加边界条件时遇到了麻烦。我的离散化网格是一个二维正方形,在 x,y 方向上从 -L 到 L。 我想强加这样的边界条件, 在 x=L 的边界线上,指定了函数
Fortran 是否有与 C assert 等效的标准函数/关键字? ? 我找不到 assert我在Fortran2003标准中提到过。我发现了一些如何使用预处理器的方法,但是在这个 answer建议
我有一系列的作业,使用“;”将它们分配给同一个ike。分开statemnts,但我收到此错误: 1.0;磅(1,9) 1个 错误:(1)处无法分类的陈述 在文件LJ.F90:223中 如果每个语句都在
我正在使用 gfortran -std=f2008。我有一个函数,它返回一个包含可分配数组的派生类型。该函数在返回之前调用allocate()。似乎在分配数组的函数返回之后,数组会自动释放一段时间,并
我制作了这个小型测试程序来“证明”在编译之前(或者如果你让它们可分配),你不能在不指定它们的大小的情况下使用向量。我的观点失败了。我期待本地向量“num”会失败。程序在执行程序之前无法知道它的大小。大
出于优化原因,Fortran 强制子例程或函数的虚拟参数不是别名,即它们不指向相同的内存位置。 我想知道相同的约束是否适用于函数的返回值。 换句话说,对于给定的 myfunc 函数: function
我已经在Fortran 90中编写了一个相当大的程序。它已经运行了一段时间了,但是今天我尝试将其提高一个档次并增加问题的大小(这是研究非标准的有限元求解器,如果那样的话)。可以帮助任何人...)现在,
在 C 和 C++ 中,有许多操作会导致未定义的行为,即允许编译器做任何它想做的事情的情况。 Examples包括在释放变量后使用它,释放变量两次和取消引用空指针。 Fortran 是否也有未定义的行
通常我使用fortran进行数值分析,然后使用matlab、R和python进行后期和前期工作。 我发现 matlab、R 和 python 在终端中提供了命令提示符,以便您可以运行脚本以及从命令行立
在 Fortran 中将变量设置为 +Infinity 的最安全方法是什么?目前我正在使用: program test implicit none print *,infinity() con
我是一名优秀的程序员,十分优秀!