Subroutine parameter mod

subroutine parameter_mod is a powerful module. This subroutine is called after the program reads in all parameters and all node's material properties. After all parameters are assigned and users used this command $assignparbyfunc, The program will call this subroutine if 2D-ddcc-dyna.exe in linux version is used. You don't need to assigned or modify all parameters listed below. You only need to modify what you want to change.

This subroutine is only for experts. Any wrong assignment will lead the program to crash

subroutine parameter_mod
use nrtype
use global2d
use global2dextra
use utility, only : findinarray,lintp ,ifintriangleperiod !,callgeneration
implicit none
integer(I4B)  :: i,j,k,i1,i2 ,ksk
real(dp) :: const,temp1,temp2,temp3,KBT
real(dp)  :: xl, xr, yb ,yt , me_1,me_2,mhh,mlh, Incom, AA,BB,CC,DD, tempEg,TempM(4),tempPE(2),tempdi(2)
logical :: notfind , errorstop

 KBT=KB*par%t
 const=2.0_dp*((M0*KB*par%t)/(TWOPI_D*HBAR*HBAR))**1.5
 print *, "loading parameters to nodes"   
 do i = 1,me%nnode
     if( (me%x(i)>xl) .and. (me%x(i)<xr) .and. (me%y(i)>yb) .and. (me%y(i)<yt) ) then
               me%eg(i) = tempEg    ! bandgap  (eV)
               me%dope(i) = 1.0e18   ! doping  (1/cm^3)
               me%ea(i) = 0.025         ! activation energy (eV)
               me%ea(i) = 0.025         ! activation energy (eV)
               me%efmass(1,i) = tempM(1)   ! effective mass of electron
               me%efmass(2,i) = tempM(2)   ! effective mass of electron perp
               me%efmass(3,i) = tempM(3)   ! effective mass of HHnano
               me%efmass(4,i) = tempM(4)   ! effective mass of LF
               me%pe(1,i) =tempPE(1)           ! polarization of z
               me%pe(2,i) =0.0           ! polarization of x
               me%xnd(i)= 1.0e0         ! impurity doping of material
               me%taun(i) = 1.0e-8                     ! nonradiative lifetime electron
               me%taup(i) = 1.0e-8                     ! nonradiative lifetime hole
               me%mun(i) = 300.0  ! electron mobility
               me%mup(i) = 10.0    ! hole mobility
               me%Ecoff(i) = 0.5  ! band offset
               me%par(i) = 1  ! assign this region is equal to original input file's region number (If you don't know what it means, marked it  !me%par(i)
               me%xnc(i)= const*(me%efmass(1,i)*me%efmass(1,i)*me%efmass(2,i))**0.5  ! don't modify this part
               me%generation(i)=real(i) ! generation rate here is simply a demo
               affinity(i)=4.0_dp  ! material affinity
               affinitypar(me%par(i))=29.0_dp
            end if
    end do
   ! Here is to assign element's parameter
   do i = 1,me%nelement
      if( (me%xp(i) >= xl) .and. (me%xp(i) < xr) .and. &
          (me%yp(i)<= yt) .and. (me%yp(i)> yb) ) then
               me%diEL(1:2,i) = 12.90*8.85d-14  ! dielectri constant along x(1)  y(2)
               me%egEL(i) = 3.3    ! assign the element's eg (eV)
               me%peEl(1,i) = 1.0e13   ! assign the element's polarization field in z
               me%peEl(2,i) = 0.0   ! assign the element's polarization field in x
               me%munEl(i) =300.0    ! electron mobility
               me%mupEl(i) =10.0    ! hole mobility
               me%parEl(i) = 1  ! assign this region is equal to original input file's region number (If you don't know what it means, marked it  !me%parEl(i)
               affinityEL(i) = 4.0_dp ! assign the affinity of element
       end if
     end do
     print*,"This is a demo of libmodpar.f90. You should not use this one without proper changes. The program will be stopped unless you unmark stop command"
     stop

end subroutine