12.07.2015 Views

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

C. Programas 111read(92,30) (grad(i),i=1,3*iat)do i=1,3*iatgrad(i)=grad(i)*8.432971931d-4enddoendifendifC Gradient type 2if (iop.eq.2) thenif (0.ne.in<strong>de</strong>x(string,’CARTESIAN COORDINATE DERIVATIVES’)) thenread(92,*)read(92,*)read(92,*)read(92,50) (grad(i),i=1,3*iat)do i=1,3*iatgrad(i)=grad(i)*8.432971931d-4enddoendifendifC Hessianif (iop.eq.2) thenif (0.ne.in<strong>de</strong>x(string,’FULL FORCE MATRIX, INVOKED BY "DFORCE"’))& thenread(92,*)read(92,*)read(92,*)read(92,*)call hessendifendif99 continueclose(92)10 format(a80)19 format(36x,g16.5)20 format(29x,g12.6)30 format(50x,g12.6)40 format(16x,d9.6)50 format(15x,f12.6,1x,f12.6,1x,f12.6)returnendCsubroutine hessimplicit double precision(a-h,o-z)parameter (natoms=100)common/gg/ energy,grad(3*natoms),gh(3*natoms,3*natoms)common/geom/x(natoms),y(natoms),z(natoms),is(natoms),iat,ich,imultopen(unit=92,file=’input.out’,status=’old’)k=1imax=7c=1j=1l=1icurt=43do while (j.lt.3*iat)do i=c,3*iat

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!