diff --git a/src/processing_mod.F90 b/src/processing_mod.F90 index 291608c5ee789a4204a74d01326e239824b3a6ba..54f99c63ae489810c974c4bf6a0003ebc99634be 100644 --- a/src/processing_mod.F90 +++ b/src/processing_mod.F90 @@ -277,7 +277,7 @@ SUBROUTINE compute_radial_ion_heatflux +tau_i*sqrt_tau_o_sigma_i*imagu*ky_*CONJG(psi(iky,ikx,izgs:izge))& *kernel_i(in,iky,ikx,izgs:izge,0)*(& 0.5_dp*SQRT2*SQRT3*moments_i(ip3_i,in ,iky,ikx,izgs:izge,updatetlevel)& - +1.5_dp*CONJG(moments_i(ip1_i,in ,iky,ikx,izgs:izge,updatetlevel))& !????? + +1.5_dp*moments_i(ip1_i,in ,iky,ikx,izgs:izge,updatetlevel)& +(2._dp*n_dp+1._dp)*moments_i(ip1_i,in ,iky,ikx,izgs:izge,updatetlevel)& -(n_dp+1._dp)*moments_i(ip1_i,in+1,iky,ikx,izgs:izge,updatetlevel)& -n_dp*moments_i(ip1_i,in-1,iky,ikx,izgs:izge,updatetlevel))