@@ -65,42 +65,46 @@ void Ions::calc_electron_temperature(Neutrals neutrals, Grid grid, Times time) {
65
65
epsilon = calc_epsilon (neutrals);
66
66
}
67
67
68
- report.print (4 , " Calculating photoelectron heating" );
69
68
70
69
// Photoelectron heating
71
70
if (input.get_do_photoelectron_heating ()) {
71
+ report.print (4 , " Calculating photoelectron heating" );
72
72
Qphe = calc_photoelectron_heating (epsilon);
73
73
}
74
74
75
- report.print (4 , " Calculating ionization heating" );
76
-
77
75
// Ionization heating (includes all ionization sources)
78
76
if (input.get_do_ionization_heating ()) {
77
+ report.print (4 , " Calculating ionization heating" );
79
78
QIonization = calc_ionization_heating (epsilon);
80
79
}
81
80
82
- report.print (4 , " Calculating electron-ion collisions" );
83
81
84
82
// electron-ion collisions
85
83
if (input.get_do_electron_ion_collisional_heating ()) {
84
+ report.print (4 , " Calculating electron-ion collisions" );
85
+
86
86
Qeic = calc_electron_ion_collisions ();
87
87
Qeicp = Qeic[0 ];
88
88
Qeicm = Qeic[1 ];
89
89
Qeic_v = Qeic[2 ]; // Friction
90
90
}
91
91
92
- report.print (4 , " Calculating electron-neutral elastic collisions" );
93
92
94
93
// electron-neutral Elastic collisions
95
94
if (input.get_do_electron_neutral_elastic_collisional_heating ()) {
95
+ report.print (4 , " Calculating electron-neutral elastic collisions" );
96
+
96
97
Qenc = calc_electron_neutral_elastic_collisions (neutrals);
97
98
Qencp = Qenc[0 ];
98
99
Qencm = Qenc[1 ];
99
100
Qenc_v = Qenc[2 ]; // Friction
100
101
}
101
102
103
+
102
104
// electron-neutral inelastic collisions
103
105
if (input.get_do_electron_neutral_inelastic_collisional_heating ()) {
106
+ report.print (4 , " Calculating electron-neutral inelastic collisions" );
107
+
104
108
Qenc_inelastic = calc_electron_neutral_inelastic_collisions (neutrals);
105
109
Qrotm = Qenc_inelastic[0 ];
106
110
Qrotp = Qenc_inelastic[1 ];
0 commit comments