include line break

gafusion · Dec 4, 2023 · 8ea308a · 8ea308a
1 parent d4d1635
commit 8ea308a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git 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)