Demo entry 6314808

fortrn

   

Submitted by anonymous on Oct 30, 2016 at 04:17
Language: Fortran. Code size: 970 Bytes.

!Normalization for gyrokinetic case
  if(iload/=9)then
! 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).