[Wien] Some interesting observation with "runsp_lapw -so -orb"on Mac OSX
Gavin Abo
gsabo at crimson.ua.edu
Mon Dec 24 20:37:26 CET 2012
Your statement in lapwso.f never runs:
! find out whether orb potential should be added (iorbpot=1)
if(iunit.eq.12)then
iorbpot=1
endif
It gets skipped over by the goto 11 in:
DO 10 I = LEN(fname), 1, -1
IF (fname(I:I) .EQ. '.') THEN
IF (fname(i+1:i+4) .eq. 'in1c') flcom=.true.
goto 11
ENDIF
10 CONTINUE
You could put the statement before, but not after the DO 10 section.
On 12/23/2012 7:38 PM, Zhu, Jianxin wrote:
> Dear Peter,
>
>
> In my previous message, I sent you the two files with the changes based on
> what you suggested.
> For a given case with given values of U and J, "runsp_lapw -so -orb"
> converges nicely.
> Just now, when I changed to use other values of U and J, the above run
> seems to not see the use of new U and J values and just stops at the first
> iteration.
>
>
>
> Cheers,
>
> Jianxin
More information about the Wien
mailing list