We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent e8dd429 commit d6a791fCopy full SHA for d6a791f
profile.c
@@ -291,7 +291,7 @@ void Profile(AtomicLine *line)
291
phi_delta * sin2_gamma * atmos.cos_2chi[mu][k] * sv[k];
292
phi_U[k] +=
293
phi_delta * sin2_gamma * atmos.sin_2chi[mu][k] * sv[k];
294
- phi_V[k] += sign *
+ phi_V[k] += sign * line->c_fraction[n] *
295
0.5*(phi_sp - phi_sm) * atmos.cos_gamma[mu][k] * sv[k];
296
297
if (input.magneto_optical) {
@@ -303,7 +303,7 @@ void Profile(AtomicLine *line)
303
psi_delta * sin2_gamma * atmos.cos_2chi[mu][k] * sv[k];
304
psi_U[k] +=
305
psi_delta * sin2_gamma * atmos.sin_2chi[mu][k] * sv[k];
306
- psi_V[k] += sign *
+ psi_V[k] += sign * line->c_fraction[n] *
307
0.5 * (psi_sp - psi_sm) * atmos.cos_gamma[mu][k] * sv[k];
308
}
309
0 commit comments