From 8ea308a5e4200868f457dee72bc328d66f0baada Mon Sep 17 00:00:00 2001 From: Harry Dudding Date: Mon, 4 Dec 2023 08:31:55 +0000 Subject: [PATCH] include line break --- tglf/src/tglf_multiscale_spectrum.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/tglf/src/tglf_multiscale_spectrum.f90 b/tglf/src/tglf_multiscale_spectrum.f90 index e08953f77..8df8ebaeb 100644 --- a/tglf/src/tglf_multiscale_spectrum.f90 +++ b/tglf/src/tglf_multiscale_spectrum.f90 @@ -443,10 +443,11 @@ SUBROUTINE get_multiscale_spectrum if(YTs(i)==0.0 .OR. kP>=kT)then field_spectrum_out(2,j,i) = 0.0 else if (ky0 <= kP) then ! initial quadratic - sig_ratio = (aoverb * (ky0 ** 2) + ky0 + coverb) / (aoverb * (k0 ** 2) + k0 + coverb) + sig_ratio = (aoverb * (ky0 ** 2) + ky0 + coverb) / (aoverb * (k0 ** 2) + k0 + coverb) field_spectrum_out(2,j,i) = Ys(i) * (sig_ratio ** c_1) * Fky *(eigenvalue_spectrum_out(1,j,i)/gamma0)**(2 * expsub) else ! connecting quadratic - doversig0 = ((Ys(i) / YTs(i))**(1.0/abs(c_1)))-((aoverb*(kP**2)+kP+coverb-((kP-kT)*(2*aoverb*kP+1)))/(aoverb*(k0**2)+k0+coverb)) + doversig0 = ((Ys(i) / YTs(i))**(1.0/abs(c_1)))- & + ((aoverb*(kP**2)+kP+coverb-((kP-kT)*(2*aoverb*kP+1)))/(aoverb*(k0**2)+k0+coverb)) doversig0 = doversig0 * (1.0/((kP-kT)**2)) eoversig0 = - 2 * doversig0 * kP + ((2 * aoverb * kP + 1)/(aoverb * (k0 ** 2) + k0 + coverb)) foversig0 = ((Ys(i) / YTs(i))**(1.0/abs(c_1))) - eoversig0 * kT - doversig0 * (kT ** 2)