请教高手这个子函数是什么意思？

laodao325

下面是子函数newxyz中的例子
请高手帮忙解释一下最后部分中的
VCORN(3,I)=VCORN(3,I)-DELALF*(VCORN(3,I)-ZREF)/(ZMAX-ZREF)是什么意思啊？
C     Sample coding: Vertex movement calculation for a reciprocating
C                    piston in a fixed cylinder, where the piston is
C                    activated by a rotating crank mechanism, whose para-
C                    meters are:
C                              C = connecting rod length
C                              A = crank radius
C                              OMEGA = rotational speed
C                              DELTHE = crank angle covered per time step
C                              ZREF = z coordinate below which the mesh
C                                        remains unchanged
C
C      PI=3.14159265
C      C=40.0
C      A=9.0
C      OMEGA=100.
C      ZREF=1.
C      DELTHE=PI*OMEGA/30.*DT
C      THETA=REAL(ITER)*DELTHE
C      TN=SQRT(C**2-A**2*(SIN(THETA-DELTHE)**2))
C      ALFN=A*(1.-COS(THETA-DELTHE))-C+TN
C      TNP1=SQRT(C**2-A**2*(SIN(THETA)**2))
C      ALFNP1=A*(1.-COS(THETA))-C+TNP1
C      DELALF=ALFNP1-ALFN
CC---- LOCATE MAX Z
C      ZMAX=VCORN(3,1)
C      DO 20 I=2,NDMAXU
C        IF(VCORN(3,I).GT.ZMAX) ZMAX=VCORN(3,I)
C   20 CONTINUE
CC---- CALCATULE NEW Z valueS
C      DO 10 I=1,NDMAXU
C        IF(VCORN(3,I).LE.ZREF) GO TO 10
C        VCORN(3,I)=VCORN(3,I)-DELALF*(VCORN(3,I)-ZREF)/(ZMAX-ZREF)
C   10 CONTINUE
C-------------------------------------------------------------------------C
     RETURN
     END
C