Type material (material%eg ...)

Type(NODEMATERIAL) :: material is for saving each node's material variables. There are many material parameters saved in the Type "material"

 material%dope(i)
 material%ea(i)
 material%Eg(i)
 material%impurity(i)
 material%charges(i)
 material%ep(i)
 material%polx(i)
 material%poly(i)
 material%polz(i)
 material%efmass(1:4,i)
 material%trapNT(i)
 material%trapEt(i)
 material%trapdegenercy(i)
 material%traptaun(i)
 material%traptaup(i)
 material%taun(i)
 material%taup(i)

   TYPE(meshformat)  :: mshformat
   TYPE(node)  :: mshnd
   TYPE(elements)  :: mshel
   TYPE(surface)  :: mshsurface
   TYPE(volume)  :: mshvolume
   TYPE(Boundary)  :: mshboun
   TYPE(ELBoundary)  :: mshelboun
   TYPE(sprs2_DP) :: sprsA,sprstem
   TYPE(Thermalcoe) :: mshthermal
   TYPE(Electriccoe) :: mshelectric
   TYPE(NODEMATERIAL) :: material

 do i=1,mshnd%n
     nodetype = mshnd%element(i) ! Get the node number
     x=mshnd%x(i)
     y=mshnd%y(i)
     z=mshnd%z(i)
     material%dope(i)=mshelectric%dope(nodetype)
     material%Eg(i)=mshelectric%Eg(nodetype)
     material%Ea(i)=mshelectric%Ea(nodetype)
     material%impurity(i)=mshelectric%impurity(nodetype)
     material%charges(i)=mshelectric%charges(nodetype)
     material%ep(i)=mshelectric%ep(nodetype)
   !  print*, material%dope(i), material%impurity(i), material%ea(i), material%charges(i), material%ep(i), material%eg(i)
     material%polx(i)=mshelectric%polx(nodetype)
     material%poly(i)=mshelectric%poly(nodetype)
     material%polz(i)=mshelectric%polz(nodetype)

     material%efmass(1:4,i)=mshelectric%efmass(1:4,nodetype) ! OK
     material%trapNT(i)=mshelectric%trapNT(nodetype)  ! OK
     material%trapEt(i)=mshelectric%trapEt(nodetype)   !OK
     material%trapdegenercy(i)=mshelectric%trapdegenercy(nodetype) !OK
     material%traptaun(i)=mshelectric%traptaun(nodetype)  !OK
     material%traptaup(i)=mshelectric%traptaup(nodetype)  !OK
     material%taun(i)=mshelectric%taun(nodetype)   !OK
     material%taup(i)=mshelectric%taup(nodetype)    !OK
     if( ifusegenfunc) then
        generations(i) = mshelectric%taun(nodetype)
     end if
   end do

   do j=mshel%tetraelenumstart,mshel%n
     material%EgEl(j) = mshelectric%Eg(mshel%realvolnum(j)) ! Setup the tetrahydron volnum
     material%polxEl(j)=mshelectric%polx(mshel%realvolnum(j))
     material%polyEl(j)=mshelectric%poly(mshel%realvolnum(j))
     material%polzEl(j)=mshelectric%polz(mshel%realvolnum(j))

!    print*,material%EgEl(j),mshel%realvolnum(j)
   end do