[Wien] dosplot2_lapw and DOS in interstitial space
Ulrich Wedig
U.Wedig at fkf.mpg.de
Fri Feb 20 16:56:30 CET 2009
With the release of WIEN2k 9.1 I took the opportunity to update my old
WIEN installation. I managed to install the new version on my system
(openSUSE 10.2, Intel-Compiler 11.0, MKL 10.1) without any problems,
considering the recommendations given in the emails of Peter Blaha (5
Feb 2009) and Jean-Claude le Bossé (10 Feb 2009) in this mailing list.
Thanks a lot to all the WIEN2k developers for their work.
Concerning the dosplot2_lapw script I have a remark: The partial DOS of
the interstitial space is not notified and labeled correctly. I have
modified $WIENROOT/dosplot2_lapw and also
$WIENROOT/SRC_w2web/htdocs/exec/extract-dos-cases_lapw (in order to be
effective also in w2web) as follows (lines with <===== added):
.
.
printlineint:
if ( $numdos > 0)then
@ grepline = $grepline + 3
set atom = `head -$grepline <$file.int|tail -1|awk '{print $1}'`
set qtl = `head -$grepline <$file.int|tail -1|awk '{print $2}'`
set atomname = `grep Z: $file.struct|head -$atom|tail -1|cut -f1 -d" " `
set nat = `grep Z: $file.struct| wc -l ` <======
if ($atom == 0) then
set label = "total DOS"
set natom = $file
else if ($atom > $nat) then <======
set label = "interstitial" <======
set natom = $file <======
set atomname = "" <======
else
.
.
With these modifications the scripts work fine for me.
Best regards
Ulrich Wedig
--
-----------------------------------------------------------------
Dr. Ulrich Wedig Tel. 0711/6891535
Max-Planck-Institut fuer Festkoerperforschung FAX 0711/6891502
Heisenbergstr. 1
70569 Stuttgart U.Wedig at fkf.mpg.de
-----------------------------------------------------------------
More information about the Wien
mailing list