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

Cleaned up code after succesfully adding fueling constrained re-injection

parent b8433828

InputFiles/xp_IonSlowingDown.in

+1 −0
Original line number Diff line number Diff line
@@ -15,6 +15,7 @@ params%nz = 501, 
params%NC = 10000,                                                              ! Total number of particles

params%NS = 50000,                                                              ! Total number of time steps

params%dt = 1.5E-7,                                                             ! Time step in [s]

params%G  = 15.0E+21,                                                           ! Real particle injection rate



! Time steps to record:

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

InputFiles/xp_TestCollisions.in

+1 −0
Original line number Diff line number Diff line
@@ -15,6 +15,7 @@ params%nz = 501, ! 
params%NC = 100000,                                                             ! Total number of particles

params%NS = 5000,                                                               ! Total number of time steps

params%dt = 0.5E-7,                                                             ! Time step in [s]

params%G  = 15.0E+21,                                                           ! Real particle injection rate



! Time steps to record:

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

src/linFP.f90

+3 −7
Original line number Diff line number Diff line
@@ -255,10 +255,8 @@ NS_loop: DO j = 1,params%NS 
    ! ==============================================================================

    dNR  = a_new*(uN1 + uN2)*dt - (N1 + N2)*dt

    dNSP = dNR/alpha

!    dER = (E3 + E5 - E4 - E1 - E2)*dt

    NR  = NR  + dNR

    NSP = NSP + dNSP

 !   ER  = ER + dER



    ! Assign new NR and NSP:

    ! ============================================================================
@@ -266,8 +264,6 @@ NS_loop: DO j = 1,params%NS 
    output%NSP(j) = NSP

    output%ER(j)  = ER

    

    !WRITE(*,*) 'ER: ', output%ER(j)



    ! Assign particle and energy rates in physical units [P/s] and [J/s]

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

    output%Ndot1(j) = N1