[Wien] A_lm, B_lm and C_lm coefficients
Sergio Yanuen Rodriguez
srodriguez at physics.tamu.edu
Wed Jun 11 19:20:35 CEST 2008
Dear Fatemeh
Thank you very much for your help, I did the changes and it compiles
without error. Just one more question, I have to run
x lapw2 -alm
or how do I activate the MODUS ALM to get the coefficients and in which
file are they stored?
Sergio Y. Rodriguez
Physics Department
Texas A&M University
> 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----------------------
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
>
More information about the Wien
mailing list