subroutine MREMD
- implicit real*8 (a-h,o-z)
+ implicit none
include 'DIMENSIONS'
include 'mpif.h'
include 'COMMON.CONTROL'
include 'COMMON.VAR'
include 'COMMON.MD'
+#ifdef FIVEDIAG
+ include 'COMMON.LAGRANGE.5diag'
+#else
+ include 'COMMON.LAGRANGE'
+#endif
+ include 'COMMON.QRESTR'
#ifndef LANG0
include 'COMMON.LANGEVIN'
#else
+#ifdef FIVEDIAG
+ include 'COMMON.LANGEVIN.lang0.5diag'
+#else
include 'COMMON.LANGEVIN.lang0'
#endif
+#endif
include 'COMMON.CHAIN'
include 'COMMON.DERIV'
include 'COMMON.GEO'
include 'COMMON.SETUP'
include 'COMMON.MUCA'
include 'COMMON.HAIRPIN'
- integer ERRCODE
+ double precision time00,time01,time02,time03,time04,time05,
+ & time06,time07,time08,time001,tt0
+ double precision scalfac
+ integer i,j,k,il,il1,ii,iex,itmp,i_temp,i_mult,i_iset,i_mset,
+ & i_dir,i_temp1,i_mult1,i_mset1
+ integer ERRCODE,ierr,ierror
double precision cm(3),L(3),vcm(3)
double precision energia(0:n_ene)
double precision remd_t_bath(maxprocs)
external ilen
character*50 tytul
common /gucio/ cm
- integer itime
+ integer itime,i_set_temp,itt,itime_master,irr,i_iset1
+ integer nharp,iharp(4,maxres/3)
cold integer nup(0:maxprocs),ndown(0:maxprocs)
integer rep2i(0:maxprocs),ireqi(maxprocs)
integer icache_all(maxprocs)
integer status(MPI_STATUS_SIZE),statusi(MPI_STATUS_SIZE,maxprocs)
logical synflag,end_of_run,file_exist /.false./,ovrtim
+ double precision t_bath_temp,delta,ene_iex_iex,ene_i_i,ene_iex_i,
+ & ene_i_iex,xxx,tmp,econstr_temp_iex,econstr_temp_i
+ integer iran_num
+ double precision ran_number
cdeb imin_itime_old=0
ntwx_cache=0
c-----------------------------------------------------------------------
subroutine write1rst(i_index)
- implicit real*8 (a-h,o-z)
+ implicit none
include 'DIMENSIONS'
include 'mpif.h'
+ include 'COMMON.CONTROL'
include 'COMMON.MD'
+#ifdef FIVEDIAG
+ include 'COMMON.LAGRANGE.5diag'
+#else
+ include 'COMMON.LAGRANGE'
+#endif
+ include 'COMMON.QRESTR'
include 'COMMON.IOUNITS'
include 'COMMON.REMD'
include 'COMMON.SETUP'
integer*2 i_index
& (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
common /przechowalnia/ d_restart1,d_restart2
+ integer i,j,il1,il,ixdrf
+ integer ierr
t5_restart1(1)=totT
t5_restart1(2)=EK
subroutine write1traj
- implicit real*8 (a-h,o-z)
+ implicit none
include 'DIMENSIONS'
include 'mpif.h'
include 'COMMON.MD'
+ include 'COMMON.QRESTR'
include 'COMMON.IOUNITS'
include 'COMMON.REMD'
include 'COMMON.SETUP'
& p_uscdiff(100*maxprocs)
real p_c(3,(maxres2+2)*maxprocs),r_c(3,maxres2+2)
common /przechowalnia/ p_c
+ integer ii,i,il,j,ixdrf
+ integer ierr
call mpi_bcast(ii_write,1,mpi_integer,
& king,CG_COMM,ierr)
subroutine read1restart(i_index)
- implicit real*8 (a-h,o-z)
+ implicit none
include 'DIMENSIONS'
include 'mpif.h'
+ include 'COMMON.CONTROL'
include 'COMMON.MD'
+#ifdef FIVEDIAG
+ include 'COMMON.LAGRANGE.5diag'
+#else
+ include 'COMMON.LAGRANGE'
+#endif
+ include 'COMMON.QRESTR'
include 'COMMON.IOUNITS'
include 'COMMON.REMD'
include 'COMMON.SETUP'
integer*2 i_index
& (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
common /przechowalnia/ d_restart1
+ integer i,j,il,il1,ixdrf,iret,itmp
+ integer ierr
write (*,*) "Processor",me," called read1restart"
if(me.eq.king)then
end
subroutine read1restart_old
- implicit real*8 (a-h,o-z)
+ implicit none
include 'DIMENSIONS'
include 'mpif.h'
include 'COMMON.MD'
+#ifdef FIVEDIAG
+ include 'COMMON.LAGRANGE.5diag'
+#else
+ include 'COMMON.LAGRANGE'
+#endif
include 'COMMON.IOUNITS'
include 'COMMON.REMD'
include 'COMMON.SETUP'
real d_restart1(3,2*maxres*maxprocs),r_d(3,2*maxres),
& t5_restart1(5)
common /przechowalnia/ d_restart1
+ integer i,j,il,itmp
+ integer ierr
if(me.eq.king)then
open(irest2,file=mremd_rst_name,status='unknown')
read (irest2,*) (i2rep(i),i=0,nodes-1)