[Wien] A_lm, B_lm and C_lm coefficients
fatemeh.mirjani
f_mirjani at ph.iut.ac.ir
Wed Jun 11 06:52:36 CEST 2008
Dear Sergio;
You should do these changes in l2main.F:
!.......for momentum program
IF(MODUS.EQ.'ALM ') then
write(23,4645)jatom,jri(jatom),r0(jatom), &
dx(jatom),rmt(jatom)
do l=0,lmax2
write(23,*) l
if (l.le.lomax) then
write(23,4646) (RRAD1(jrj,l),RRAD2(jrj,l),RADE1(jrj,l),RADE2(jrj,l), &
a1lo(jrj,1,l),b1lo(jrj,1,l),a1lo(jrj,2,l),b1lo(jrj,2,l),jrj=1,jri(jatom))
else
write(23,4646) (RRAD1(jrj,l),RRAD2(jrj,l),RADE1(jrj,l),RADE2(jrj,l),0.0d0,0.d0,0.0d0,0.d0, &
jrj=1,jri(jatom))
endif
enddo
endif
!................................
!......for momentum program
IF(MODUS.EQ.'ALM ') then
! WRITE(24,*) mu
DO NUM=NEMIN,NEMAX
INDEX=0
DO L=0,LMAX2
MAX=2*L+1
CFAC=IMAG**L
DO M=-l,l
INDEX=INDEX+1
write(24,4893)num,l,m,index,alm(INDEX,NUM) &
,blm(INDEX,NUM),(clm(INDEX,NUM,jlo),jlo=1,2)
! write(24)alm(INDEX,NUM) &
! ,blm(INDEX,NUM),(clm(INDEX,NUM,jlo),jlo=1,2)
ENDDO
ENDDO
ENDDO
ENDIF
!..........................................
--
Best Regards
Fatemeh Mirjani | Email: f_mirjani at ph.iut.ac.ir
Computational Condensed Matter Research Lab. | Tel/Fax Office: +98311-3913746
ICTP Affiliated Center in Isfahan, | Tel Lab.: +98311-3913731
----------Physics Department, Isfahan University of Technology, Iran----------------------
More information about the Wien
mailing list