@@ -858,7 +858,8 @@ function cea_mixture_per_weight_to_per_mole(mptr, len, per_weight, per_mole) res
858858 per_mole(:ns) = mix% per_mole_from_per_weight(per_weight(:ns))
859859 end function
860860
861- function cea_mixture_chem_eq_ratio_to_of_ratio (mptr , len , oxidant_weights , fuel_weights , chem_eq_ratio , of_ratio ) result(ierr) bind(c)
861+ function cea_mixture_chem_eq_ratio_to_of_ratio (mptr , len , oxidant_weights , fuel_weights , &
862+ chem_eq_ratio , of_ratio ) result(ierr) bind(c)
862863 integer (c_int) :: ierr
863864 type (c_ptr), intent (in ), value :: mptr
864865 integer (c_int), intent (in ), value :: len
@@ -878,7 +879,8 @@ function cea_mixture_chem_eq_ratio_to_of_ratio(mptr, len, oxidant_weights, fuel_
878879 of_ratio = mix% of_from_equivalence(oxidant_weights(:ns), fuel_weights(:ns), chem_eq_ratio)
879880 end function
880881
881- function cea_mixture_weight_eq_ratio_to_of_ratio (mptr , len , oxidant_weights , fuel_weights , weight_eq_ratio , of_ratio ) result(ierr) bind(c)
882+ function cea_mixture_weight_eq_ratio_to_of_ratio (mptr , len , oxidant_weights , fuel_weights , &
883+ weight_eq_ratio , of_ratio ) result(ierr) bind(c)
882884 integer (c_int) :: ierr
883885 type (c_ptr), intent (in ), value :: mptr
884886 integer (c_int), intent (in ), value :: len
@@ -898,7 +900,8 @@ function cea_mixture_weight_eq_ratio_to_of_ratio(mptr, len, oxidant_weights, fue
898900 of_ratio = mix% of_from_phi(oxidant_weights(:ns), fuel_weights(:ns), weight_eq_ratio)
899901 end function
900902
901- function cea_mixture_of_ratio_to_weights (mptr , len , oxidant_weights , fuel_weights , of_ratio , weights ) result(ierr) bind(c)
903+ function cea_mixture_of_ratio_to_weights (mptr , len , oxidant_weights , fuel_weights , &
904+ of_ratio , weights ) result(ierr) bind(c)
902905 integer (c_int) :: ierr
903906 type (c_ptr), intent (in ), value :: mptr
904907 integer (c_int), intent (in ), value :: len
@@ -921,7 +924,8 @@ function cea_mixture_of_ratio_to_weights(mptr, len, oxidant_weights, fuel_weight
921924 of_ratio)
922925 end function
923926
924- function cea_mixture_calc_property (mptr , prop_type , len_weights , weights , temperature , prop_value ) result(ierr) bind(c)
927+ function cea_mixture_calc_property (mptr , prop_type , len_weights , weights , &
928+ temperature , prop_value ) result(ierr) bind(c)
925929 integer (c_int) :: ierr
926930 type (c_ptr), intent (in ), value :: mptr
927931 integer (c_int), intent (in ), value :: prop_type
@@ -953,7 +957,8 @@ function cea_mixture_calc_property(mptr, prop_type, len_weights, weights, temper
953957 end select
954958 end function
955959
956- function cea_mixture_calc_property_multitemp (mptr , prop_type , len_weights , weights , len_temperatures , temperatures , prop_value ) result(ierr) bind(c)
960+ function cea_mixture_calc_property_multitemp (mptr , prop_type , len_weights , &
961+ weights , len_temperatures , temperatures , prop_value ) result(ierr) bind(c)
957962 integer (c_int) :: ierr
958963 type (c_ptr), intent (in ), value :: mptr
959964 integer (c_int), intent (in ), value :: prop_type
@@ -986,7 +991,8 @@ function cea_mixture_calc_property_multitemp(mptr, prop_type, len_weights, weigh
986991 end select
987992 end function
988993
989- function cea_mixture_calc_property_tp (mptr , prop_type , len_weights , weights , temperature , pressure , prop_value ) result(ierr) bind(c)
994+ function cea_mixture_calc_property_tp (mptr , prop_type , len_weights , weights , &
995+ temperature , pressure , prop_value ) result(ierr) bind(c)
990996 integer (c_int) :: ierr
991997 type (c_ptr), intent (in ), value :: mptr
992998 integer (c_int), intent (in ), value :: prop_type
@@ -1019,7 +1025,8 @@ function cea_mixture_calc_property_tp(mptr, prop_type, len_weights, weights, tem
10191025 end select
10201026 end function
10211027
1022- function cea_mixture_calc_property_tp_multitemp (mptr , prop_type , len_weights , weights , len_temperatures , temperatures , pressure , prop_value ) result(ierr) bind(c)
1028+ function cea_mixture_calc_property_tp_multitemp (mptr , prop_type , len_weights , weights , &
1029+ len_temperatures , temperatures , pressure , prop_value ) result(ierr) bind(c)
10231030 integer (c_int) :: ierr
10241031 type (c_ptr), intent (in ), value :: mptr
10251032 integer (c_int), intent (in ), value :: prop_type
@@ -1550,7 +1557,8 @@ function cea_rocket_solver_get_size(sptr, eq_variable, eq_value) result(ierr) bi
15501557 end select
15511558 end function
15521559
1553- function cea_rocket_solver_solve_iac (sptr , slptr , weights , pc , pi_p , n_pi_p , subar , nsubar , supar , nsupar , n_frz , hc_or_tc , use_hc , tc_est , use_tc_est ) result(ierr) bind(c)
1560+ function cea_rocket_solver_solve_iac (sptr , slptr , weights , pc , pi_p , n_pi_p , subar , nsubar , &
1561+ supar , nsupar , n_frz , hc_or_tc , use_hc , tc_est , use_tc_est ) result(ierr) bind(c)
15541562 integer (c_int) :: ierr
15551563 type (c_ptr), intent (in ), value :: sptr
15561564 type (c_ptr), intent (in ), value :: slptr
@@ -1660,7 +1668,8 @@ function cea_rocket_solver_solve_iac(sptr, slptr, weights, pc, pi_p, n_pi_p, sub
16601668 if (.not. solution% converged) ierr = CEA_NOT_CONVERGED
16611669 end function
16621670
1663- function cea_rocket_solver_solve_fac (sptr , slptr , weights , pc , pi_p , n_pi_p , subar , nsubar , supar , nsupar , n_frz , hc_or_tc , use_hc , mdot_or_acat , use_mdot , tc_est , use_tc_est ) result(ierr) bind(c)
1671+ function cea_rocket_solver_solve_fac (sptr , slptr , weights , pc , pi_p , n_pi_p , subar , nsubar , &
1672+ supar , nsupar , n_frz , hc_or_tc , use_hc , mdot_or_acat , use_mdot , tc_est , use_tc_est ) result(ierr) bind(c)
16641673 integer (c_int) :: ierr
16651674 type (c_ptr), intent (in ), value :: sptr
16661675 type (c_ptr), intent (in ), value :: slptr
@@ -2096,9 +2105,11 @@ function cea_shock_solver_solve(sptr, slptr, weights, T0, p0, mach1_or_u1, use_m
20962105 if (refl_froz .eqv. .true. ) reflected_frozen = .true.
20972106
20982107 if (use_mach .eqv. .true. ) then
2099- solution = solver% solve(weights(:nr), T0, p0, mach1= mach1_or_u1, reflected= reflected, incident_frozen= incident_frozen, reflected_frozen= reflected_frozen)
2108+ solution = solver% solve(weights(:nr), T0, p0, mach1= mach1_or_u1, &
2109+ reflected= reflected, incident_frozen= incident_frozen, reflected_frozen= reflected_frozen)
21002110 else
2101- solution = solver% solve(weights(:nr), T0, p0, u1= mach1_or_u1, reflected= reflected, incident_frozen= incident_frozen, reflected_frozen= reflected_frozen)
2111+ solution = solver% solve(weights(:nr), T0, p0, u1= mach1_or_u1, &
2112+ reflected= reflected, incident_frozen= incident_frozen, reflected_frozen= reflected_frozen)
21022113 end if
21032114 if (.not. solution% converged) ierr = CEA_NOT_CONVERGED
21042115 end function
0 commit comments