[Wien] format of structure file

Jutta Rogal rogal at fhi-berlin.mpg.de
Mon Aug 9 11:06:28 CEST 2004


Dear all,
I found an inconsistency for the format of the structure file, 
especially for the R0 and RMT, during the initialization and the main 
programs.  In nn, e.g. we have the following format for read and write
SRC_nn:
nn.f:         READ(20,1050) NAME(JATOM),JRJ(JATOM),R0(JATOM),RMTT(jatom), &
nn.f: 1050 FORMAT(A10,5X,I5,5X,F10.8,5X,F10.4,5x,f5.1)

nn.f:         write(21,1149) 
NAMEN(JATOM)(1:2),jatom,JRJN(JATOM),R0N(JATOM),RMTTN(jatom),zzn(jatom)
nn.f: 1149 FORMAT(A2,i1,7x,' NPT=',I5,'  R0=',F10.8,' RMT=',F10.4,'   
Z:',f5.1)

But in the main programs like lapw2 it is like that
SRC_lapw2:
lapw2.F:         READ(20,1050) 
ANAME(JATOM),JRI(JATOM),R0(JATOM),RMT(JATOM), &
lapw2.F: 1050 FORMAT(A10,5X,I5,5X,F10.9,5X,F10.5,5X,F10.5)

I recongnized this because also in the mini program, where new structure 
files are written, the old format is used, so that the relaxed structure 
files are different from the initial one:
SRC_mini:
haupt.f: 1049 FORMAT(A10,' NPT=',I5,'  R0=',F10.9,' RMT=',F10.5,'   
Z:',f10.5)
haupt.f:        WRITE(21,1049) 
ANAME(JATOM),JRI(JATOM),RNOT(JATOM),RMT(JATOM), &

As far as I could see, at least in my case, the R0 and RMT were still 
read in correctly, but that might depend on how many digits you specify 
after the point.
Regards
Jutta

-- 
--
---------------------------------------------------------------------
   Jutta Rogal

   Fritz-Haber-Institut
   der Max-Planck-Gesellschaft
   Abteilung Theorie
   Faradayweg 4-6
   14195 Berlin/Germany

   Phone: ++49 30 8413 4854
   Fax:   ++49 30 8413 4701
   rogal at fhi-berlin.mpg.de
---------------------------------------------------------------------
----------------------key:1-0.0735-11600-23.05:fhi-------------------





More information about the Wien mailing list