Question about read_md_hist.F90  [SOLVED]

Documentation, Web site and code modifications

Moderators: baguetl, routerov

Locked
denis
Posts: 3
Joined: Sun Oct 06, 2013 2:04 pm

Question about read_md_hist.F90

Post by denis » Sun Oct 06, 2013 2:36 pm

It looks like entropy is not read correctly (ekin is read instead):

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.

mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

Re: Question about read_md_hist.F90  [SOLVED]

Post by mverstra » Sun Oct 06, 2013 4:41 pm

someone fixed this in 7.5 - check release of 7.6 in a few weeks.

Matthieu
Matthieu Verstraete
University of Liege, Belgium

Locked