Skip to content
Snippets Groups Projects
Commit f7125b3d authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

comment the correction operator for heat flux

parent c35f0204
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 2 deletions
src/processing_mod.F90
+ 2
− 2
View file @ f7125b3d
...@@ -370,8 +370,8 @@ SUBROUTINE compute_radial_electron_heatflux ...@@ -370,8 +370,8 @@ SUBROUTINE compute_radial_electron_heatflux
ENDDO ENDDO
ENDIF ENDIF
! Add polarisation contribution ! Add polarisation contribution
integrant(izs:ize) = integrant(izs:ize) + tau_e*imagu*ky_& ! integrant(izs:ize) = integrant(izs:ize) + tau_e*imagu*ky_&
*CONJG(phi(iky,ikx,izs:ize))*phi(iky,ikx,izs:ize) * HF_phi_correction_operator(iky,ikx,izs:ize) ! *CONJG(phi(iky,ikx,izs:ize))*phi(iky,ikx,izs:ize) * HF_phi_correction_operator(iky,ikx,izs:ize)
! Integrate over z ! Integrate over z
integrant(izgs:izge) = Jacobian(izgs:izge,0)*integrant(izgs:izge) integrant(izgs:izge) = Jacobian(izgs:izge,0)*integrant(izgs:izge)
call simpson_rule_z(integrant,integral) call simpson_rule_z(integrant,integral)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment