Subroutine that normalizes the simulation variables with the previously computed characteristic scales.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(KORC_PARAMS), | intent(inout) | :: | params | Core KORC simulation parameters. |
||
| type(SPECIES), | intent(inout), | DIMENSION(:), ALLOCATABLE | :: | spp | An instance of KORC's derived type SPECIES containing all the information of different electron species. See korc_types. |
|
| type(FIELDS), | intent(inout) | :: | F | @param[in,out] F An instance of KORC's derived type FIELDS containing all the information about the fields used in the simulation. See korc_types and korc_fields. |
||
| type(PROFILES), | intent(inout) | :: | P | @param[in,out] P An instance of KORC's derived type PROFILES containing all the information about the plasma profiles used in the simulation. See korc_types.f90 and korc_profiles.f90. |
subroutine normalize_variables(params,spp,F,P)
!! @note Subroutine that normalizes the simulation variables with
!! the previously computed characteristic scales. @endnote
TYPE(KORC_PARAMS), INTENT(INOUT) :: params
!! Core KORC simulation parameters.
TYPE(SPECIES), DIMENSION(:), ALLOCATABLE, INTENT(INOUT) :: spp
!! An instance of KORC's derived type SPECIES containing all
!! the information of different electron species. See [[korc_types]].
TYPE(FIELDS), INTENT(INOUT) :: F
!! @param[in,out] F An instance of KORC's derived type FIELDS
!! containing all the information about the fields used in the simulation.
!! See [[korc_types]] and [[korc_fields]].
TYPE(PROFILES), INTENT(INOUT) :: P
!! @param[in,out] P An instance of KORC's derived type PROFILES containing all the information about the plasma profiles used in
!! the simulation. See korc_types.f90 and korc_profiles.f90.
INTEGER :: ii
!! @param ii Interator of spp array.
! Normalize params variables
params%dt = params%dt/params%cpp%time
params%simulation_time = params%simulation_time/params%cpp%time
params%snapshot_frequency = params%snapshot_frequency/params%cpp%time
params%minimum_particle_energy = params%minimum_particle_energy/params%cpp%energy
params%init_time=params%init_time/params%cpp%time
! Normalize particle variables
do ii=1_idef,size(spp)
spp(ii)%q = spp(ii)%q/params%cpp%charge
spp(ii)%m = spp(ii)%m/params%cpp%mass
spp(ii)%Eo = spp(ii)%Eo/params%cpp%energy
spp(ii)%Eo_lims = spp(ii)%Eo_lims/params%cpp%energy
spp(ii)%wc = spp(ii)%wc*params%cpp%time
spp(ii)%wc_r = spp(ii)%wc_r*params%cpp%time
spp(ii)%vars%X = spp(ii)%vars%X/params%cpp%length
spp(ii)%vars%V = spp(ii)%vars%V/params%cpp%velocity
spp(ii)%vars%Rgc = spp(ii)%vars%Rgc/params%cpp%length
spp(ii)%Ro = spp(ii)%Ro/params%cpp%length
spp(ii)%Zo = spp(ii)%Zo/params%cpp%length
spp(ii)%r_inner = spp(ii)%r_inner/params%cpp%length
spp(ii)%r_outter = spp(ii)%r_outter/params%cpp%length
spp(ii)%sigmaR = spp(ii)%sigmaR/params%cpp%length
spp(ii)%sigmaZ = spp(ii)%sigmaZ/params%cpp%length
spp(ii)%falloff_rate = spp(ii)%falloff_rate*params%cpp%length
spp(ii)%Xtrace = spp(ii)%Xtrace/params%cpp%length
spp(ii)%Spong_b = spp(ii)%Spong_b/params%cpp%length
spp(ii)%Spong_w = spp(ii)%Spong_w/params%cpp%length
spp(ii)%dR = spp(ii)%dR/params%cpp%length
spp(ii)%dZ = spp(ii)%dZ/params%cpp%length
end do
! Normalize electromagnetic fields and profiles
F%Bo = F%Bo/params%cpp%Bo
F%Eo = F%Eo/params%cpp%Eo
F%AB%Ero = F%AB%Ero/params%cpp%Eo
F%AB%rmn=F%AB%rmn/params%cpp%length
F%AB%sigmamn=F%AB%sigmamn/params%cpp%length
F%Ro = F%Ro/params%cpp%length
F%Zo = F%Zo/params%cpp%length
F%E_dyn = F%E_dyn/params%cpp%Eo
F%E_pulse=F%E_pulse/params%cpp%time
F%E_width=F%E_width/params%cpp%time
F%t0_2x1t=F%t0_2x1t/params%cpp%time
F%circumradius=F%circumradius/params%cpp%length
if (F%E_profile.eq.'MST_FSA') then
F%AB%a = F%AB%a/params%cpp%length
F%AB%Ro = F%AB%Ro/params%cpp%length
end if
P%a = P%a/params%cpp%length
P%R0 = P%R0/params%cpp%length
P%Z0 = P%Z0/params%cpp%length
P%R0_RE = P%R0_RE/params%cpp%length
P%Z0_RE = P%Z0_RE/params%cpp%length
P%neo = P%neo/params%cpp%density
P%n_ne = P%n_ne/params%cpp%density
P%n_shelf = P%n_shelf/params%cpp%density
P%Teo = P%Teo/params%cpp%temperature
P%n_REr0=P%n_REr0/params%cpp%length
P%n_tauion=P%n_tauion/params%cpp%time
P%n_tauin=P%n_tauin/params%cpp%time
P%n_tauout=P%n_tauout/params%cpp%time
P%n_shelfdelay=P%n_shelfdelay/params%cpp%time
P%n_lamfront=P%n_lamfront/params%cpp%length
P%n_lamback=P%n_lamback/params%cpp%length
P%n_lamshelf=P%n_lamshelf/params%cpp%length
if (params%profile_model .EQ. 'ANALYTICAL') then
if (params%field_eval.eq.'interp') then
if (ALLOCATED(P%X%R)) P%X%R = P%X%R/params%cpp%length
if (ALLOCATED(P%X%Z)) P%X%Z = P%X%Z/params%cpp%length
if (ALLOCATED(P%ne_2D)) P%ne_2D = P%ne_2D/params%cpp%density
if (ALLOCATED(P%Te_2D)) P%Te_2D = P%Te_2D/params%cpp%temperature
end if
else if (params%profile_model(1:8) .EQ. 'EXTERNAL') then
if (ALLOCATED(P%X%R)) P%X%R = P%X%R/params%cpp%length
if (ALLOCATED(P%X%Z)) P%X%Z = P%X%Z/params%cpp%length
if (ALLOCATED(P%ne_2D)) P%ne_2D = P%ne_2D/params%cpp%density
if (ALLOCATED(P%Te_2D)) P%Te_2D = P%Te_2D/params%cpp%temperature
if (ALLOCATED(P%ne_3D)) P%ne_3D = P%ne_3D/params%cpp%density
if (ALLOCATED(P%Te_3D)) P%Te_3D = P%Te_3D/params%cpp%temperature
if (params%profile_model(10:10).eq.'H') then
if (ALLOCATED(P%nRE_2D)) P%nRE_2D = P%nRE_2D/params%cpp%density
if (ALLOCATED(P%nAr0_2D)) P%nAr0_2D = P%nAr0_2D/params%cpp%density
if (ALLOCATED(P%nAr1_2D)) P%nAr1_2D = P%nAr1_2D/params%cpp%density
if (ALLOCATED(P%nAr2_2D)) P%nAr2_2D = P%nAr2_2D/params%cpp%density
if (ALLOCATED(P%nAr3_2D)) P%nAr3_2D = P%nAr3_2D/params%cpp%density
if (ALLOCATED(P%nD_2D)) P%nD_2D = P%nD_2D/params%cpp%density
if (ALLOCATED(P%nD1_2D)) P%nD1_2D = P%nD1_2D/params%cpp%density
endif
end if
if (params%field_model(1:10) .EQ. 'ANALYTICAL') then
F%AB%Bo = F%AB%Bo/params%cpp%Bo
F%AB%lambda = F%AB%lambda/params%cpp%length
F%AB%Bpo = F%AB%Bpo/params%cpp%Bo
if (.not.(F%E_profile.eq.'MST_FSA')) then
F%AB%a = F%AB%a/params%cpp%length
F%AB%Ro = F%AB%Ro/params%cpp%length
end if
if (params%field_eval.eq.'interp') then
if (ALLOCATED(F%B_2D%R)) F%B_2D%R = F%B_2D%R/params%cpp%Bo
if (ALLOCATED(F%B_2D%PHI)) F%B_2D%PHI = F%B_2D%PHI/params%cpp%Bo
if (ALLOCATED(F%B_2D%Z)) F%B_2D%Z = F%B_2D%Z/params%cpp%Bo
if (ALLOCATED(F%E_2D%R)) F%E_2D%R = F%E_2D%R/params%cpp%Eo
if (ALLOCATED(F%E_2D%PHI)) F%E_2D%PHI = F%E_2D%PHI/params%cpp%Eo
if (ALLOCATED(F%E_2D%Z)) F%E_2D%Z = F%E_2D%Z/params%cpp%Eo
if (ALLOCATED(F%B_3D%R)) F%B_3D%R = F%B_3D%R/params%cpp%Bo
if (ALLOCATED(F%B_3D%PHI)) F%B_3D%PHI = F%B_3D%PHI/params%cpp%Bo
if (ALLOCATED(F%B_3D%Z)) F%B_3D%Z = F%B_3D%Z/params%cpp%Bo
if (ALLOCATED(F%E_3D%R)) F%E_3D%R = F%E_3D%R/params%cpp%Eo
if (ALLOCATED(F%E_3D%PHI)) F%E_3D%PHI = F%E_3D%PHI/params%cpp%Eo
if (ALLOCATED(F%E_3D%Z)) F%E_3D%Z = F%E_3D%Z/params%cpp%Eo
if (ALLOCATED(F%PSIp)) F%PSIp = F%PSIp/ &
(params%cpp%Bo*params%cpp%length**2)
F%PSIP_min = F%PSIP_min/ &
(params%cpp%Bo*params%cpp%length**2)
if (.not.params%field_model.eq.'M3D_C1') then
F%PSIp_lim = F%PSIp_lim/ &
(params%cpp%Bo*params%cpp%length**2)
end if
F%X%R = F%X%R/params%cpp%length
! Nothing to do for the PHI component
F%X%Z = F%X%Z/params%cpp%length
if (params%orbit_model(3:5).eq.'pre') then
if (ALLOCATED(F%gradB_2D%R)) F%gradB_2D%R = F%gradB_2D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_2D%PHI)) F%gradB_2D%PHI = F%gradB_2D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_2D%Z)) F%gradB_2D%Z = F%gradB_2D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%curlb_2D%R)) F%curlb_2D%R = F%curlb_2D%R/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_2D%PHI)) F%curlb_2D%PHI = F%curlb_2D%PHI/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_2D%Z)) F%curlb_2D%Z = F%curlb_2D%Z/ &
(1./params%cpp%length)
if (ALLOCATED(F%gradB_3D%R)) F%gradB_3D%R = F%gradB_3D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_3D%PHI)) F%gradB_3D%PHI = F%gradB_3D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_3D%Z)) F%gradB_3D%Z = F%gradB_3D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%curlb_3D%R)) F%curlb_3D%R = F%curlb_3D%R/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_3D%PHI)) F%curlb_3D%PHI = F%curlb_3D%PHI/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_3D%Z)) F%curlb_3D%Z = F%curlb_3D%Z/ &
(1./params%cpp%length)
end if
end if
else if (params%field_model(1:8) .EQ. 'EXTERNAL') then
if (ALLOCATED(F%B_3D%R)) F%B_3D%R = F%B_3D%R/params%cpp%Bo
if (ALLOCATED(F%B_3D%PHI)) F%B_3D%PHI = F%B_3D%PHI/params%cpp%Bo
if (ALLOCATED(F%B_3D%Z)) F%B_3D%Z = F%B_3D%Z/params%cpp%Bo
if (ALLOCATED(F%E_3D%R)) F%E_3D%R = F%E_3D%R/params%cpp%Eo
if (ALLOCATED(F%E_3D%PHI)) F%E_3D%PHI = F%E_3D%PHI/params%cpp%Eo
if (ALLOCATED(F%E_3D%Z)) F%E_3D%Z = F%E_3D%Z/params%cpp%Eo
if (ALLOCATED(F%PSIp)) F%PSIp = F%PSIp/(params%cpp%Bo*params%cpp%length**2)
if (ALLOCATED(F%PSIp3D)) F%PSIp3D = F%PSIp3D/ &
(params%cpp%Bo*params%cpp%length**2)
F%PSIP_min = F%PSIP_min/(params%cpp%Bo*params%cpp%length**2)
F%PSIp_lim = F%PSIp_lim/(params%cpp%Bo*params%cpp%length**2)
if (ALLOCATED(F%B_2D%R)) F%B_2D%R = F%B_2D%R/params%cpp%Bo
if (ALLOCATED(F%B_2D%PHI)) F%B_2D%PHI = F%B_2D%PHI/params%cpp%Bo
if (ALLOCATED(F%B_2D%Z)) F%B_2D%Z = F%B_2D%Z/params%cpp%Bo
if (ALLOCATED(F%E_2D%R)) F%E_2D%R = F%E_2D%R/params%cpp%Eo
if (ALLOCATED(F%E_2D%PHI)) F%E_2D%PHI = F%E_2D%PHI/params%cpp%Eo
if (ALLOCATED(F%E_2D%Z)) F%E_2D%Z = F%E_2D%Z/params%cpp%Eo
if (params%orbit_model(3:5).eq.'pre') then
if (ALLOCATED(F%gradB_2D%R)) F%gradB_2D%R = F%gradB_2D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_2D%PHI)) F%gradB_2D%PHI = F%gradB_2D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_2D%Z)) F%gradB_2D%Z = F%gradB_2D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%curlb_2D%R)) F%curlb_2D%R = F%curlb_2D%R/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_2D%PHI)) F%curlb_2D%PHI = F%curlb_2D%PHI/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_2D%Z)) F%curlb_2D%Z = F%curlb_2D%Z/ &
(1./params%cpp%length)
if (ALLOCATED(F%gradB_3D%R)) F%gradB_3D%R = F%gradB_3D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_3D%PHI)) F%gradB_3D%PHI = F%gradB_3D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%gradB_3D%Z)) F%gradB_3D%Z = F%gradB_3D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%curlb_3D%R)) F%curlb_3D%R = F%curlb_3D%R/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_3D%PHI)) F%curlb_3D%PHI = F%curlb_3D%PHI/ &
(1./params%cpp%length)
if (ALLOCATED(F%curlb_3D%Z)) F%curlb_3D%Z = F%curlb_3D%Z/ &
(1./params%cpp%length)
end if
if (F%dBfield) then
if (ALLOCATED(F%dBdR_2D%R)) F%dBdR_2D%R = F%dBdR_2D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdR_2D%PHI)) F%dBdR_2D%PHI = F%dBdR_2D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdR_2D%Z)) F%dBdR_2D%Z = F%dBdR_2D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_2D%R)) F%dBdPHI_2D%R = F%dBdPHI_2D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_2D%PHI)) F%dBdPHI_2D%PHI = F%dBdPHI_2D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_2D%Z)) F%dBdPHI_2D%Z = F%dBdPHI_2D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_2D%R)) F%dBdZ_2D%R = F%dBdZ_2D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_2D%PHI)) F%dBdZ_2D%PHI = F%dBdZ_2D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_2D%Z)) F%dBdZ_2D%Z = F%dBdZ_2D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdR_3D%R)) F%dBdR_3D%R = F%dBdR_3D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdR_3D%PHI)) F%dBdR_3D%PHI = F%dBdR_3D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdR_3D%Z)) F%dBdR_3D%Z = F%dBdR_3D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_3D%R)) F%dBdPHI_3D%R = F%dBdPHI_3D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_3D%PHI)) F%dBdPHI_3D%PHI = F%dBdPHI_3D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdPHI_3D%Z)) F%dBdPHI_3D%Z = F%dBdPHI_3D%Z/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_3D%R)) F%dBdZ_3D%R = F%dBdZ_3D%R/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_3D%PHI)) F%dBdZ_3D%PHI = F%dBdZ_3D%PHI/ &
(params%cpp%Bo/params%cpp%length)
if (ALLOCATED(F%dBdZ_3D%Z)) F%dBdZ_3D%Z = F%dBdZ_3D%Z/ &
(params%cpp%Bo/params%cpp%length)
end if
if (F%B1field) then
if (params%field_model(10:13).eq.'MARS') then
if (ALLOCATED(F%B1Re_2D%R)) F%B1Re_2D%R = F%B1Re_2D%R/ &
params%cpp%Bo
if (ALLOCATED(F%B1Re_2D%PHI)) F%B1Re_2D%PHI = F%B1Re_2D%PHI/ &
params%cpp%Bo
if (ALLOCATED(F%B1Re_2D%Z)) F%B1Re_2D%Z = F%B1Re_2D%Z/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2D%R)) F%B1Im_2D%R = F%B1Im_2D%R/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2D%PHI)) F%B1Im_2D%PHI = F%B1Im_2D%PHI/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2D%Z)) F%B1Im_2D%Z = F%B1Im_2D%Z/ &
params%cpp%Bo
else if (params%field_model(10:14).eq.'AORSA') then
if (ALLOCATED(F%B1Re_2DX%X)) F%B1Re_2DX%X = F%B1Re_2DX%X/ &
params%cpp%Bo
if (ALLOCATED(F%B1Re_2DX%Y)) F%B1Re_2DX%Y = F%B1Re_2DX%Y/ &
params%cpp%Bo
if (ALLOCATED(F%B1Re_2DX%Z)) F%B1Re_2DX%Z = F%B1Re_2DX%Z/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2DX%X)) F%B1Im_2DX%X = F%B1Im_2DX%X/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2DX%Y)) F%B1Im_2DX%Y = F%B1Im_2DX%Y/ &
params%cpp%Bo
if (ALLOCATED(F%B1Im_2DX%Z)) F%B1Im_2DX%Z = F%B1Im_2DX%Z/ &
params%cpp%Bo
endif
end if
if (F%E1field) then
if (ALLOCATED(F%E1Re_2DX%X)) F%E1Re_2DX%X = F%E1Re_2DX%X/ &
params%cpp%Eo
if (ALLOCATED(F%E1Re_2DX%Y)) F%E1Re_2DX%Y = F%E1Re_2DX%Y/ &
params%cpp%Eo
if (ALLOCATED(F%E1Re_2DX%Z)) F%E1Re_2DX%Z = F%E1Re_2DX%Z/ &
params%cpp%Eo
if (ALLOCATED(F%E1Im_2DX%X)) F%E1Im_2DX%X = F%E1Im_2DX%X/ &
params%cpp%Eo
if (ALLOCATED(F%E1Im_2DX%Y)) F%E1Im_2DX%Y = F%E1Im_2DX%Y/ &
params%cpp%Eo
if (ALLOCATED(F%E1Im_2DX%Z)) F%E1Im_2DX%Z = F%E1Im_2DX%Z/ &
params%cpp%Eo
end if
F%X%R = F%X%R/params%cpp%length
! Nothing to do for the PHI component
F%X%Z = F%X%Z/params%cpp%length
if (F%Dim2x1t) then
F%X%PHI=F%X%PHI/params%cpp%time
end if
! end if
end if
end subroutine normalize_variables