Commit eef47507 authored by Caneses Marin, Juan Francisco's avatar Caneses Marin, Juan Francisco
Browse files

Tested the RF heating operator, saw x24 speed up with openMP and all physics... 

Tested the RF heating operator, saw x24 speed up with openMP and all physics modeuls: particle push, reinjection, collisions and RF heating
parent 0c8356db

InputFiles/xp_TestDrag.in

+9 −9
Original line number Diff line number Diff line 
&in_nml

! Simulation name:

! ===============

in%fileDescriptor = '2021_02_11b',

in%fileDescriptor = '2021_02_11d',



! Magnetic field input data:

! =========================
@@ -21,7 +21,7 @@ in%zmin = -5.0, 
in%iSave = .true.,

in%iPush = .true.,                                                             ! Enable RK4 particle pusher

in%iColl = .true.,                                                             ! Collisions or no collisions

in%iHeat = .false.,                                                            ! Enable RF heating operator

in%iHeat = .true.,                                                            ! Enable RF heating operator

in%iPotential = .false.,                                                       ! Enable electric potential

in%iDrag = .true.,                                                            ! Enable plasma drag
@@ -62,13 +62,13 @@ in%IC_xip = 0.707, ! 


! RF heating operator conditions:

! ===============================

in%f_RF = 28.0E+9                                                              ! RF frequency

in%zRes1 = 2.4                                                                 ! Defines interval where RF is present

in%zRes2 = 3.0                                                                 ! Defines interval where RF is present

in%kper = 29920.                                                               ! Perpendicular wave number of EBW

in%kpar = 5864.                                                                ! Parallel wave number of EBW

in%Ew   = 10000                                                                ! EBW E-minus wave electric field magnitude

in%n_harmonic = 4                                                              ! Cyclotron harmonic

in%f_RF  = 12.195E+6                                                              ! RF frequency

in%zRes1 = 0.                                                                 ! Defines interval where RF is present

in%zRes2 = +1.5                                                                 ! Defines interval where RF is present

in%kper  = 125.                                                               ! Perpendicular wave number of EBW

in%kpar  = 30.                                                                ! Parallel wave number of EBW

in%Ew    = 4000                                                                ! EBW E-minus wave electric field magnitude

in%n_harmonic = 2                                                              ! Cyclotron harmonic



! Electric potential conditions:

! =============================

src/linFP.f90

+1 −0
Original line number Diff line number Diff line
@@ -260,6 +260,7 @@ AllTime: do j = 1,in%Nsteps 
           CALL CyclotronResonanceNumber(zp(i),kep(i),xip(i),f1,in,spline_B)

           df = dsign(1.d0,f0*f1)

           if (df .LT. 0 .AND. zp(i) .GT. in%zRes1 .AND. zp(i) .LT. in%zRes2)  then

              !WRITE(*,*) 'Resonance at zp: ', zp(i)

              CALL RFHeatingOperator(zp(i),kep(i),xip(i),ecnt3,pcnt3,in,spline_B,spline_dB,spline_ddB,spline_dV)

           end if

        end if