Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (geom_testpardiso) Browse

fennecs/src/xg_mod.f90
8d7b3488dcccgeom_testpardiso

Last Change
Raw File

xg_mod.f90
View Options

MODULE xg
USE constants
USE basic, ONLY : nr, nz, nrun, zgrid, rgrid, nplasma, itgraph, resfile
USE beam, ONLY: parts
USE fields, ONLY: nrank, vec1, vec2
IMPLICIT NONE
!
DOUBLE PRECISION :: tstep=0.0_db, dnorm=1.0_db
DOUBLE PRECISION :: zb1, zb2, rb1, rb2
DOUBLE PRECISION, ALLOCATABLE :: pzxg(:), prxg(:), pthetxg(:)
INTEGER :: nhist, maxhist
DOUBLE PRECISION, ALLOCATABLE :: thist(:), tepot(:), tekin(:), tetot(:), terror(:)
INTEGER :: x0=250, xsize = 410, y0=0, ysize=335
!
CONTAINS
SUBROUTINE initw
!
! Initialize Xgrafix
!
! Particle phase space
!
ALLOCATE (pzxg(parts%Nptot))
ALLOCATE (prxg(parts%Nptot))
ALLOCATE (pthetxg(parts%Nptot))
!
! Bounds of space domain
!
zb1 = zgrid(0)
zb2 = zgrid(nz)
rb1 = rgrid(0)
rb2 = rgrid(nr)
!
! Time history
!
nhist = 0
maxhist=nrun/itgraph+1
ALLOCATE(thist(0:maxhist))
ALLOCATE(tepot(0:maxhist))
ALLOCATE(tekin(0:maxhist))
ALLOCATE(tetot(0:maxhist))
ALLOCATE(terror(0:maxhist))
!
CALL updt_xg_var
!
CALL xginit(3,'espic2d:'//TRIM(resfile),'espic2d'//TRIM(resfile)//'.o','espic2d'//TRIM(resfile)//'.dump','N4','N5','N6',tstep)
!
! Particle phase space plots
!
pzxg(:) = parts%uz(1:parts%Nptot)
prxg(:) = parts%ur(1:parts%Nptot)
pthetxg(:) = parts%uthet(1:parts%Nptot)
! Es potential
!
CALL xgset2d('linlin','Z','Z-R phase space','open',x0, y0, dnorm, dnorm, &
& .FALSE., .FALSE., ZB1, ZB2, RB1, RB2)
CALL xgscat2d(parts%Z,parts%R,parts%Nptot,1)
!
CALL xgset2d('linlin','PZ','PZ-PR phase space','open',x0+2*xsize, y0, 1._db, 1._db, &
& .TRUE., .TRUE., 0._db, 1._db, 0._db, 1._db)
CALL xgscat2d(pzxg, prxg, parts%Nptot,1)
!
CALL xgset2d('linlin','PZ','PZ-PTHET phase space','open',x0+2*xsize, y0+ysize, 1._db, 1._db, &
& .TRUE., .TRUE., 0._db, 1._db, 0._db, 1._db)
CALL xgscat2d(pzxg, pthetxg, parts%Nptot,1)
!
CALL xgset2d('linlin','Z','Z-PZ phase space','open',x0+xsize, y0, dnorm, 1.0_db, &
& .FALSE., .TRUE., ZB1, ZB2, 0._db, 1._db)
CALL xgscat2d(parts%Z, pzxg, parts%Nptot,1)
!
! Time history
!
CALL xgset2d('linlin','T','Energies', 'open', x0, y0+ysize, 1._db, 1._db, &
& .TRUE., .TRUE., 0._db, 1._db, 0._db, 1._db)
CALL xgcurve(thist, tepot, nhist, 1)
CALL xgcurve(thist, tekin, nhist, 2)
CALL xgcurve(thist, tetot, nhist, 3)
!
CALL xgset2d('linlog','T','Relative error in Etot', 'open', x0+xsize, y0+ysize, 1._db, 1._db, &
& .TRUE., .TRUE., 0._db, 1._db, 0._db, 1._db)
CALL xgcurve(thist, terror, nhist, 1)
!
! Xgrafix control
!
CALL xgevent
CALL xgupdate
END SUBROUTINE initw
!--------------------------------------------------------------------------------
SUBROUTINE updt_xg_var
!
! Update XG variables
!
USE basic, ONLY : cstep, dt
USE basic, ONLY : omegap, omegac, pot
USE beam, ONLY : ekin, epot, etot, etot0
!
! Phase space
tstep = REAL(cstep,db)
pzxg(:) = parts%uz(1:parts%Nptot)
prxg(:) = parts%ur(1:parts%Nptot)
pthetxg(:) = parts%uthet(1:parts%Nptot)
!
! Time history
!
IF(omegac.GT.omegap) THEN
thist(nhist) = nhist*itgraph*dt*omegac
ELSE
thist(nhist) = nhist*itgraph*dt*omegap
END IF
tekin(nhist) = ekin
tepot(nhist) = epot
tetot(nhist) = etot
terror(nhist) = ABS((etot-etot0)/etot) ! To be dispalyed on log scale!
nhist = nhist+1
!
END SUBROUTINE updt_xg_var
END MODULE xg
c4science · Help