檢視 Function getmunbyET(ind,Ex,Ey,T,xp,yp) 的原始碼

function getmunbyET(ind,Ex,Ey,T,xp,yp) is to return the program its mobility according to electric field, temperature, x, y. If some more information is needed, such as carrier  density, it can call the global variable 

    function getmunbyET(ind,Ex,Ey,T,xp,yp)
    use global2d
    implicit none
    integer*4 :: ind   ! element's global index number
    real*8  ::  getmunbyET ! return value of mobility
    integer*4  :: i,j,k,regionindex
    real*8 ::  Ex,Ey , T ,xp,yp, vsat, Esaturation ,Eslope ,mun_low, Esat
    real*8  :: v , x
    real*8  :: A, Asat, Amin, B, Bmax, Bstart,Bslope , nn(3),pp(3) , nave,pave,n0,p0 , orginalmun
    do i=1,3
        nn(i)= me%n( me%element(i,ind) )   ! carrier density is saved in node, not element
        pp(i)= me%p( me%element(i,ind) )   ! mobility is for element. So we need to get the average value
    end do
    regionindex=me%parEL(ind) ! To obtained the region information for this node
    nave=sum(nn)/3.0 ! average carrier density
    pave=sum(pp)/3.0 ! average carrier density
    if(  xp < 1.0d-4 .and. yp < 1.0d-4) then
        x=Ex
        v= B+Bmax*exp((Bstart-Ex)/Bslope)
        getmunbyET=v/Ex  *(n0/nave)
    else
        getmunbyET=mun_low *(n0/nave)
    end if
    print*, " "
    print*, "************************************************************************"
    print*, " Error! This version does not support using external libemufunc.so module"
    print*, " You might accidentally set $useemunfunc  to activate using external module"
    print*, " Please execute 2D-ddcc-dyna.exe instead."
    Print*, " Program Stop! "
    stop
    end function getmunbyET