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 125emin=e(1)do i=1,nif (e(i).lt.emin) emin=e(i)enddodo i=1,ne(i)=(e(i)-emin)*627.5095d0enddoreturnendC Compute ZPEsubroutine zpen(nfreq,n_r,freq,zpe)implicit double precision(a-h,o-z)inclu<strong>de</strong> ’data.h’dimension freq(nconfs,nfreqs)zpe=0.D0do i=1,nfreqif (freq(n_r,i).gt.0.D0) zpe=zpe+0.5D0*freq(n_r,i)enddozpe=zpe*2.859143523D-3C zpe=zpe/349.75d0returnendC Reor<strong>de</strong>r frequenciessubroutine reor<strong>de</strong>rfreq(n,freq,nf)implicit double precision(a-h,o-z)inclu<strong>de</strong> ’data.h’dimension freq(nconfs,nfreqs)common/ska/ scalefC Scale Factordo i=1,nffreq(n,i)=freq(n,i)*scalefenddoC Reor<strong>de</strong>rdo 50 j=1,nfdo 50 k=1,j-1if (freq(n,k) .lt. freq(n,j)) thensalva=freq(n,j)do l=j-1,k,-1freq(n,l+1)=freq(n,l)end dofreq(n,k)=salvaend if50 continuereturnendC Compute the Maxwell-Botzmann weight of each partition functionsubroutine wgt(q,temp,enr)implicit double precision(a-h,o-z)inclu<strong>de</strong> ’data.h’q=q*exp(-enr/(rg*temp))returnend

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

Saved successfully!

Ooh no, something went wrong!