c--------------------------------- subroutine csa_read implicit real*8 (a-h,o-z) include 'DIMENSIONS' include 'COMMON.CSA' include 'COMMON.BANK' include 'COMMON.IOUNITS' open(icsa_in,file=csa_in,status="old",err=100) read(icsa_in,*) nconf read(icsa_in,*) jstart,jend read(icsa_in,*) nstmax read(icsa_in,*) n1,n2,n3,n4,n5,n6,n7,n8,is1,is2 read(icsa_in,*) nran0,nran1,irr read(icsa_in,*) nseed read(icsa_in,*) ntotal,cut1,cut2 read(icsa_in,*) estop read(icsa_in,*) icmax,irestart read(icsa_in,*) ntbankm,dele,difcut read(icsa_in,*) iref,rmscut,pnccut read(icsa_in,*) ndiff close(icsa_in) return 100 continue return end c--------------------------------- subroutine initial_write implicit real*8 (a-h,o-z) include 'DIMENSIONS' include 'COMMON.CSA' include 'COMMON.BANK' include 'COMMON.IOUNITS' open(icsa_seed,file=csa_seed,status="unknown") write(icsa_seed,*) "seed" close(31) #if defined(AIX) || defined(PGI) open(icsa_history,file=csa_history,status="unknown", & position="append") #else open(icsa_history,file=csa_history,status="unknown", & access="append") #endif write(icsa_history,*) nconf write(icsa_history,*) jstart,jend write(icsa_history,*) nstmax write(icsa_history,*) n1,n2,n3,n4,n5,n6,n7,n8,is1,is2 write(icsa_history,*) nran0,nran1,irr write(icsa_history,*) nseed write(icsa_history,*) ntotal,cut1,cut2 write(icsa_history,*) estop write(icsa_history,*) icmax,irestart write(icsa_history,*) ntbankm,dele,difcut write(icsa_history,*) iref,rmscut,pnccut write(icsa_history,*) ndiff write(icsa_history,*) close(icsa_history) open(icsa_bank1,file=csa_bank1,status="unknown") write(icsa_bank1,*) 0 close(icsa_bank1) return end c--------------------------------- subroutine restart_write implicit real*8 (a-h,o-z) include 'DIMENSIONS' include 'COMMON.IOUNITS' include 'COMMON.CSA' include 'COMMON.BANK' #if defined(AIX) || defined(PGI) open(icsa_history,file=csa_history,position="append") #else open(icsa_history,file=csa_history,access="append") #endif write(icsa_history,*) write(icsa_history,*) "This is restart" write(icsa_history,*) write(icsa_history,*) nconf write(icsa_history,*) jstart,jend write(icsa_history,*) nstmax write(icsa_history,*) n1,n2,n3,n4,n5,n6,n7,n8,is1,is2 write(icsa_history,*) nran0,nran1,irr write(icsa_history,*) nseed write(icsa_history,*) ntotal,cut1,cut2 write(icsa_history,*) estop write(icsa_history,*) icmax,irestart write(icsa_history,*) ntbankm,dele,difcut write(icsa_history,*) iref,rmscut,pnccut write(icsa_history,*) ndiff write(icsa_history,*) write(icsa_history,*) "irestart is: ", irestart write(icsa_history,*) close(icsa_history) return end c---------------------------------