Commit c3cc1bc6 authored by Juan Caneses Marin (nfc)'s avatar Juan Caneses Marin (nfc)
Browse files

Completed the removal of padding suceesfully, results are correct based on postprocessing

parent e7893b3c
Loading
Loading
Loading
Loading
+2 −2
Original line number Diff line number Diff line
&in_nml
! Simulation name:
! ===============
in%fileDescriptor = '2021_02_03a',
in%fileDescriptor = '2021_02_04a',

! Magnetic field input data:
! =========================
@@ -12,7 +12,7 @@ in%nz = 501, ! N

! Simulation conditions:
! =====================
in%Nparts = 500000,                                                              ! Total number of particles
in%Nparts = 50000,                                                              ! Total number of particles
in%Nsteps = 5000,                                                               ! Total number of time steps
in%dt     = 0.5E-7,                                                             ! Time step in [s]
in%zmax = +5.0,
+1 −2
Original line number Diff line number Diff line
@@ -44,7 +44,7 @@ REAL(r8) :: df
! simulation time:
REAl(r8) :: tp
! Particle position (zp), kinetic energy (kep), pitch angle (xip):
REAL(r8), DIMENSION(:,:)  , ALLOCATABLE :: xip, zp, kep
REAL(r8), DIMENSION(:)  , ALLOCATABLE :: xip, zp, kep
! subset of zp, kep and xip to save:
REAL(r8), DIMENSION(:,:), ALLOCATABLE :: zp_hist, kep_hist, xip_hist
! Subset of tp:
@@ -280,7 +280,6 @@ TimeStepping: do j = 1,in%Nsteps
            		  WRITE(*,*) ''
            		  WRITE(*,*) '*********************************************************************'
            		  WRITE(*,*) "Number of threads given: ", in%threads_given
            		  WRITE(*,*) "Padding: ", in%padding, " Real(r8)"
            		  WRITE(*,*) '*********************************************************************'
            		  WRITE(*,*) ''
            		end if
+1 −1
Original line number Diff line number Diff line
@@ -7,7 +7,7 @@ INPUT_FILE="xp_IonSlowingDown.in"

# Set number of threads:
# ===================================================
export OMP_NUM_THREADS=1
export OMP_NUM_THREADS=4

# Set processor binding for openMP:
# ===================================================