Merge branch 'master' of https://github.com/loftytopping/PyBox

loftytopping · Oct 23, 2020 · a83e72d · a83e72d
2 parents 5223928 + cb08695
commit a83e72d
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/Parse_eqn_file.py b/Parse_eqn_file.py
@@ -3366,7 +3366,7 @@ def write_partitioning_section_fortran(total_length_y,num_bins,num_species):
     f.write('        Inverse_Kn(:)=1.0D0/Kn(:)\n') 
     f.write('        Correction_part1(:)=(1.33D0+0.71D0*Inverse_Kn(:))/(1.0D0+Inverse_Kn(:))\n') 
     f.write('        Correction_part2(:)=(4.0D0*(1.0D0-alpha_d_org(:)))/(3.0D0*alpha_d_org(:))\n') 
-    f.write('        Correction_part3(:)=1.0E0+(Correction_part2(:)+Correction_part2(:))*Kn(:)\n') 
+    f.write('        Correction_part3(:)=1.0E0+(Correction_part1(:)+Correction_part2(:))*Kn(:)\n') 
     f.write('        Correction(:)=1.0D0/Correction_part3(:)\n') 
     f.write('\n')     	
     f.write('        ! 6) Kelvin factor  \n') 
@@ -3585,7 +3585,7 @@ def write_partitioning_section_fortran_ignore(total_length_y,num_bins,num_specie
     f.write('        Inverse_Kn(:)=1.0D0/Kn(:)\n') 
     f.write('        Correction_part1(:)=(1.33D0+0.71D0*Inverse_Kn(:))/(1.0D0+Inverse_Kn(:))\n') 
     f.write('        Correction_part2(:)=(4.0D0*(1.0D0-alpha_d_org(:)))/(3.0D0*alpha_d_org(:))\n') 
-    f.write('        Correction_part3(:)=1.0E0+(Correction_part2(:)+Correction_part2(:))*Kn(:)\n') 
+    f.write('        Correction_part3(:)=1.0E0+(Correction_part1(:)+Correction_part2(:))*Kn(:)\n') 
     f.write('        Correction(:)=1.0D0/Correction_part3(:)\n') 
     f.write('\n')     	
     f.write('        ! 6) Kelvin factor  \n')