Initdip.f90

SUBROUTINE Initdip(trel,kon,koff,tswitch,dip,N,Nc)
        IMPLICIT NONE
                INTEGER*8 :: i, j, N, Nc, dip(N,Nc)
                REAL*8 :: con, trel(N,Nc), tswitch(N,Nc),kon,koff
                REAL*4 :: rand
! Loop over chains
                DO i = 1,Nc

! Loop over links
                        DO j = 1,N

! Assign initial state randomly according to on/off rates
! Calculate time till next switch (Poisson process)
                                con = kon/(kon+koff)
                                        IF(rand(0) .LT. con)THEN
                                                dip(j,i) = 1
                                                trel(j,i) = 0.d0
                                                tswitch(j,i) = -log(rand(0))/koff
                                        ELSE
                                                dip(j,i) = 0
                                                trel(j,i) = 0.d0
                                                tswitch(j,i) = -log(rand(0))/kon
                                        ENDIF
                        ENDDO
                ENDDO
END SUBROUTINE Initdip