0

谁能告诉我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
4

1 回答 1

2

这条线

COMMON/XQANG/SUM(300,10),MDIR,COSTH

通知子例程关于一个名为的公共块XQANG，该块具有一个名为 的元素COSTH。在没有其他信息的情况下，并且在那个年份的代码中，这很可能是一个real变量。

公共块是早期的 Fortran 机制，用于在程序单元之间共享变量，例如在主程序和子程序之间。在直接使用中，相同的公共块声明将在多个位置找到，具有相同的变量列表。每个声明都引用相同的变量。

不过有一个转折点，需要小心谨慎。公共块实际上是一个共享的内存块，并不要求声明的每个实例都标识相同的变量。公共块的一个常见用途是在一个位置声明一个包含 100 个实数的数组，但在其他地方声明两个数组，每个数组包含 50 个实数——相同的内存，不同的变量。

更好的是，它们也可以用来改变变量的类型（有点）。公共块的一种用法可能包含一个占用 4 个字节的实变量，而同一块的另一种用法可能在同一位置包含一个 4 字节整数变量 - 相同的位，不同的解释。

这些扭曲是常见块被广泛弃用的原因之一。他们弃用的另一个原因是他们掩盖了变量的共享，现在我们大多数人更喜欢通过他们的参数列表显式地将参数传入和传出子例程。

我猜很少有 50 岁以下的 Fortran 程序员仍在使用它们编写新代码。但是 8 - 80 岁的 Fortran 程序员仍在研究包含它们的代码。

于 2018-02-23T07:14:19.850 回答