zero.for

来自「有限元计算程序」· FOR 代码 · 共 31 行

C$DEBUG
C$LARGE
	SUBROUTINE ZERO(ELOAD,MELEM,MEVAB,MPOIN,MTOTG,MTOTV,NDOFN,NELEM,
     .                  NEVAB,NGAUS,NSTR1,NTOTG,EPSTN,EFFST,
     .                  NTOTV,NVFIX,STRSG,TDISP,TFACT,
     .                  TLOAD,TREAC,MVFIX)                 !(Page 238)
C*********************************************************************
C
C*** THIS SUBROUTINE INITIALISES VARIOUS ARRAYS TO ZERO
C
C*********************************************************************
	DIMENSION ELOAD(MELEM,MEVAB),STRSG(4,MTOTG),TDISP(MTOTV),
     .            TLOAD(MELEM,MEVAB),TREAC(MVFIX,2),EPSTN(MTOTG),
     .            EFFST(MTOTG)
	TFACT=0.0
	DO 30 IELEM=1,NELEM
	DO 30 IEVAB=1,NEVAB
	ELOAD(IELEM,IEVAB)=0.0
  30  TLOAD(IELEM,IEVAB)=0.0
	DO 40 ITOTV=1,NTOTV
  40  TDISP(ITOTV)=0.0
	DO 50 IVFIX=1,NVFIX
	DO 50 IDOFN=1,NDOFN
  50  TREAC(IVFIX,IDOFN)=0.0
	DO 60 ITOTG=1,NTOTG
	EPSTN(ITOTG)=0.0
	EFFST(ITOTG)=0.0
	DO 60 ISTR1=1,NSTR1
  60  STRSG(ISTR1,ITOTG)=0.0
	RETURN
	END