这个Fortran变量来自哪里？
有人能告诉我下面的子程序中使用的这个Fortran变量来自哪里？,fortran,fortran-common-block,Fortran,Fortran Common Block,有人能告诉我下面的子程序中使用的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(
costhSUBROUTINE 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
COMMON/XQANG/SUM(300,10),MDIR,COSTH
XQANGCOSTHrealCOMMON/XQANG/SUM(300,10),MDIR,COSTH
XQANGCOSTHreal我猜50岁以下的Fortran程序员中很少有人还在使用它们编写新代码。但是8-80的Fortran程序员仍然在编写包含它们的代码。添加
隐式无costhcommonIMPLICIT NONEcosthcommon