[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