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
   410   411   412   413   414   415   416   417   418   419   420