Function usefielddependent p(indEl, matind, Ex, Ey, T, xp, yp)

Function usefielddependent_n(indEl, matind ,Ex,Ey,T,xp,yp) is activated by the command $usefieldmufunc. This function is similar to command $useemunfunc, $useemupfunc. But this function provide more parameters for use to assign in the input file. Therefore, this function can be continuous expanded to many different mobility models.

   function usefielddependent_p(indEl, matind, Ex, Ey, T, xp, yp)
   use global2Dextra3
   use nrtype
   use global2D, only : KB
   implicit none
   real(dp)  ::  usefielddependent_p
   integer(i4b),intent(in)  :: indEl, matind
   integer(i4b)  :: i,j,k
   real(dp),intent(in) :: Ex, Ey ,T ,xp,yp
   real(dp) :: mu, E, muconst
   if(munpfunctype(matind)==1) then
       if(Ey<0) then
           E=0
       else
           E=Ey
       end if
       !E=abs(Ey)
       if(munpfuncpara(4,matind)*sqrt(E) < 50.0_dp) then
           mu= munpfuncpara(3,matind)*exp(munpfuncpara(4,matind)*sqrt(E))
       else
           mu= munpfuncpara(3,matind)*exp(50.0_dp)
       end if
       usefielddependent_p=mu
   elseif(munpfunctype(matind)==11) then
       if(Ey>0) then
           E=0
       else
           E=abs(Ey)
       end if
       !E=abs(Ey)
       if(munpfuncpara(4,matind)*sqrt(E) < 50.0_dp) then
           mu= munpfuncpara(3,matind)*exp(munpfuncpara(4,matind)*sqrt(E))
       else
           mu= munpfuncpara(3,matind)*exp(50.0_dp)
       end if
       usefielddependent_p=mu
   elseif(munpfunctype(matind)==2) then
       if(Ex<0) then
           E=0
       else
           E=Ex
       end if
       !E=abs(Ey)
       if(munpfuncpara(4,matind)*sqrt(E) < 50.0_dp) then
           mu= munpfuncpara(3,matind)*exp(munpfuncpara(4,matind)*sqrt(E))
       else
           mu= munpfuncpara(3,matind)*exp(50.0_dp)
       end if
       usefielddependent_p=mu
   elseif(munpfunctype(matind)==21) then
       if(Ex>0) then
           E=0
       else
           E=abs(Ex)
       end if
       !E=abs(Ey)
       if(munpfuncpara(4,matind)*sqrt(E) < 50.0_dp) then
           mu= munpfuncpara(3,matind)*exp(munpfuncpara(4,matind)*sqrt(E))
       else
           mu= munpfuncpara(3,matind)*exp(50.0_dp)
       end if
       usefielddependent_p=mu
   end if
   !stop
   end function