Page 414 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 413

100 CONTINUE
C
RETURN
END
C
C-------------------------------------------------------------------
C
SUBROUTINE XFORMB(xc,yc,zc,xx,yy,zz)
IMPLICIT REAL*8(A-H,O-Z)
COMMON /EANEW / Cphin,Sphin,Cthen,Sthen,Cxin,Sxin
C
a11 = Cphin*Cthen
a21 = Cthen*Sphin
a31 = Sthen
C
a12 = -Sphin*Cxin - Sthen*Cphin*Sxin
a22 = Cphin*Cxin - Sthen*Sphin*Sxin
a32 = Cthen*Sxin
C
a13 = Sphin*Sxin - Sthen*Cphin*Cxin
a23 = -Cphin*Sxin - Sthen*Sphin*Cxin
a33 = Cthen*Cxin
C
xx = a11*xc + a12*yc + a13*zc
yy = a21*xc + a22*yc + a23*zc
zz = a31*xc + a32*yc + a33*zc
C
RETURN
END
C
C-------------------------------------------------------------------
C
SUBROUTINE CRASHSAVE
IMPLICIT REAL*8(A-H,O-Z)
real*4 randseed
integer a, ap
character*8 fname
COMMON /fname / fname
COMMON /CRASH / ncrash, ncalls
COMMON /ESAVE / Eold,Epold,Epoly(99,99),Ep(99,99),Esurf(99),Es(99)
COMMON /NCC / Kacct, Kaccr, Kaccp
COMMON /NHB / nhb(100),nhbpoly,nhbcurr(100)
COMMON /NRUN / Ncycle,Ndead,Kprint,Kdprnt,Kcyc,Neq
COMMON /SURF1 / a, ap
COMMON /SURF2 / b, c