Page 415 - Handbook Of Multiphase Flow Assurance
P. 415
414 10. Research methods in flow assurance
COMMON /NATOMS/ Norgn, Nwatrs, M1, M2, M3, M4
COMMON /SPACE / xms(1000), yms(1000), zms(1000)
COMMON /iori / iori
COMMON /offset/ xo,yo,zo,phi,theta,xi
COMMON /ZSAVE / zrun(100),epolymer,Erun
COMMON /chain / monodata(100,100),monocount(99),nbackbone
COMMON /types / iatomtype(1000,2),nbondfound,natom,ngroups,ncall
COMMON /STEPS / Deltrn, Delrot, Delpiv, Acctrn, Accrot, Accpiv
COMMON /rans / randseed
COMMON /POTPAM/ Casqof,Co,Rsqon,Rsqoff,Rcut,Pi,Tk,Temp
COMMON /RATIOS/ Ratiot, Ratior, Ratiop
open(unit=1,file='SAVEFILE.',status='unknown')
write (1,'(A8)') fname
write (1,'(I9)') ngroups
write (1,*) Temp
write (1,'(I12)') ncalls+1
write (1,*) randseed
write (1,'(I9)') ndead
write (1,'(I9)') ncycle
write (1,'(I9)') kprint
write (1,'(I9)') kdprnt
write (1,'(I2)') a
write (1,'(I2)') ap
write (1,*) b
write (1,*) c
write (1,*) xo,yo,zo
if(Iori.eq.1) then
write (1,*) 1
write (1,*)phi,theta,xi
else
write (1,*) 0
endif
write (1,'(2I9)') KCYC,Neq
write (1,'(4F32.18)') EOLD,Epold,ERUN,epolymer
C E segments
do 50 i=1,nbackbone-2
50 write (1,'(99F23.18)') (Epoly(i,j),j=i+2,nbackbone)
write (1,'(99F23.18)') (Esurf(i),i=1,nbackbone)
write (1,'(3I16)') Kacct, Kaccr, Kaccp
write (1,'(6F24.18)') Deltrn,Delrot,Delpiv, Ratiot, Ratior, Ratiop
write (1,'(99I10)') (nhb(i),i=1,nbackbone)
write (1,'(99F27.18)') (zrun(i),i=1,100)
write (1,'(2999F22.18)') (xms(i),yms(i),zms(i),i=1,natom)
c write blank spaces to fill the write buffer and force a write to disk

