Question about read_md_hist.F90
Posted: Sun Oct 06, 2013 2:36 pm
It looks like entropy is not read correctly (ekin is read instead):
Is it expected or a bug? Anyway, it looks like histEnt is not used now, so it doesn't affect ABINIT results.
Code: Select all
!Ekin
ncerr = nf90_get_var(ncid, ekin_id,hist%histEk(:))
if(ncerr /= NF90_NOERR) call handle_ncerr(ncerr," read variable ekin")
!Ekin
ncerr = nf90_get_var(ncid, ekin_id,hist%histEnt(:))
if(ncerr /= NF90_NOERR) call handle_ncerr(ncerr," read variable entropy")Is it expected or a bug? Anyway, it looks like histEnt is not used now, so it doesn't affect ABINIT results.