有人能告诉我以下子程序中使用的costh
变量吗?子程序从哪里获取该变量的值?是一些错误还是内置函数?
完整的子程序如下。
SUBROUTINE MOMENT(I,A,AP,MODE,EP,NCASC,ID)
COMMON/MOM/VY(99996),VZ(99996),VX(99996)
DIMENSION ERES(99996),FCT(99996)
COMMON/XQANG/SUM(300,10),MDIR,COSTH
EQUIVALENCE(ERES(1),VY(1)),(FCT(1),VZ(1))
IF(ID-2)10,30,40
C INITIALIZATION
10 N=1
NN=N+NCASC-1
DO 20 J=N,NN
VY(J)=EP*MDIR
VX(J)=0.
20 VZ(J)=EP*(1-MDIR)
GO TO 60
C
C CALCULATION OF MOMENT
C
CKM 30 RN4=RANF(0.)
30 RN4=RAN(ISEED) !KM
PHI=6.28318*RN4
SOX=2.*AP*EP/(A**2*931.5)
C BUG FIX A*(AP+A)REPLACED BY A**2 12/15/82
C IF(SOX.LT.0.)WRITE(*,943)AP,EP,A
VT=SQRT(SOX)
SOX=1.-COSTH**2
IF(SOX.LT.0.)SOX=0.
SINTH=SQRT(SOX)
VZSE=VT*COSTH
VYSE=VT*SINTH*SIN(PHI)
VXSE=VT*SINTH*COS(PHI)
VZ(I)=VZ(I)+VZSE
VY(I)=VY(I)+VYSE
VX(I)=VX(I)+VXSE
VZP=VZ(I)-VZSE*(A/AP+1.)
VYP=VY(I)-VYSE*(A/AP+1.)
VXP=VX(I)-VXSE*(A/AP+1.)
VPP2=VZP*VZP+VYP*VYP+VXP*VXP
EPART=AP*VPP2*469.
SOX=0.
IF(VPP2.GT.0.0)SOX=VZP/SQRT(VPP2)
C DO NOT USE QUICK FUNCTIONS HERE
ANG=ACOS(SOX)*180./3.1415927
CALL OUTEM(2,MODE,EPART,ANG)
GO TO 60
C
C END CALCULATION
C
40 VFTS=VX(I)**2+VY(I)**2+VZ(I)**2
ERES(I)=0.5*A*VFTS*931.5
IF(VFTS.NE.0.)GO TO 50
FCT(I)=0.0
GO TO 60
50 FCT(I)=ACOS(VZ(I)/SQRT(VFTS))*180./3.1415927
60 RETURN
END
答案 0 :(得分:2)
这一行
COMMON/XQANG/SUM(300,10),MDIR,COSTH
通知子例程一个名为XQANG
的公共块,它有一个名为COSTH
的元素。如果没有其他信息,并且在该年份的代码中,这很可能是real
变量。
公共块是早期的Fortran机制,用于跨程序单元共享变量,例如主程序和子程序之间。在简单的使用中,相同的公共块声明将在多个位置找到,具有相同的变量列表。每个声明都引用相同的变量。
虽然有一个转折,但要小心谨慎。公共块实际上是共享的内存块,并不要求声明的每个实例都标识相同的变量。常见块的一个常见用法是在一个位置声明一个100个实数的数组,但是在其他地方声明两个50个实数的数组 - 相同的内存,不同的变量。
更好的是,它们也可以用来改变变量的类型(排序)。公共块的一种用法可能包含占用4个字节的实变量,而同一块的另一种用法可能在同一位置包含一个4字节的整数变量 - 相同的位,不同的解释。
这些曲折是常见块被广泛弃用的原因之一。他们弃用的另一个原因是它们模糊了变量的共享,现在我们大多数人更喜欢通过参数列表显式地将参数传入和传出子例程。
我猜很少有50岁以下的Fortran程序员仍在使用它们编写新代码。但是,8-80岁的Fortran程序员仍在研究包含它们的代码。