Demo entry 6314822

fortrn

   

Submitted by anonymous on Oct 30, 2016 at 05:03
Language: Fortran. Code size: 764 Bytes.

! add diagonal term for poisson operator
if(idiagonal==1)then
  do i=mpsilow,mpsihigh
    if(feload==0)then
      diagonal(i)=(1.0-real(magnetic))*meshne(i)*qelectron*qelectron*rho0*rho0/meshte(i)
    else
      diagonal(i)=(1.0-real(magnetic))*(meshne(i)*qelectron*qelectron+meshnfe(i)*qfaste*qfaste)*rho0*rho0/meshte(i)
    endif
    if(fload==0)then
      ddum(i)=qion*qion*meshni(i)*rho0*rho0/meshti(i)!+diagonal(i)
    else
      ddum(i)=qion*qion*(meshni(i)+meshnf(i)*afast/aion)*rho0*rho0/meshti(i)+diagonal(i)
    endif
    do j=0,mtheta(i)
       ij=igrid(i)+j
       do k=1,nindex0(ij)
          ring0(k,ij)=0.0_wp-ring0(k,ij)*ddum(i)*rhom(ij)*rhom(ij)
       enddo
       ring0(6,ij)=ring0(6,ij)!+diagonal(i)
    enddo
  enddo
endif

This snippet took 0.00 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).