Browse code

Отступы.

Dmitrii Kashin authored on 04/09/2014 09:00:15
Showing 1 changed files
... ...
@@ -1,7 +1,7 @@
1 1
 /*
2 2
   Здесь содержатся установки начального поля для того, чтобы
3 3
   исследовать устойчивость метода на СУ. 
4
-  Для сравнения с Власенко.
4
+  Для сравнения с Власенко: test-1
5 5
 */
6 6
 
7 7
 #include "grid-init.h"
... ...
@@ -9,26 +9,26 @@
9 9
 void
10 10
 set_initial_field(Grid1D& grid)
11 11
 {
12
-double phi_L = 0.999999;
13
-double rho_L = 1; // [kg/m^3]
14
-double Ug_L = 400; // [m/s]
15
-double Ul_L = 400; // [m/s]
16
-
17
-double phi_R = phi_L;
18
-double Ug_R = sos*sos/(kappa*Ug_L);
19
-double Ul_R = sos*sos/(kappa*Ul_L);
20
-double rho_R = rho_L * Ug_L/Ug_R;
21
-
22
-ConservVars lv = primitive2conserv((PrimitiveVars){phi_L,rho_L,Ug_L,Ul_L});
23
-ConservVars rv = primitive2conserv((PrimitiveVars){phi_R,rho_R,Ug_R,Ul_R});
24
-
25
-grid.lbc.vars = lv;
26
-grid.rbc.vars = rv;
27
-
28
-for (unsigned int i=0; i<grid.N; i++)
29
-  if (i<grid.N/2)
30
-    grid.cells[i].vars = lv;
31
-  else
32
-    grid.cells[i].vars = rv;
12
+  double phi_L = 0.999999;
13
+  double rho_L = 1; // [kg/m^3]
14
+  double Ug_L = 400; // [m/s]
15
+  double Ul_L = 400; // [m/s]
16
+  
17
+  double phi_R = phi_L;
18
+  double Ug_R = sos*sos/(kappa*Ug_L);
19
+  double Ul_R = sos*sos/(kappa*Ul_L);
20
+  double rho_R = rho_L * Ug_L/Ug_R;
21
+  
22
+  ConservVars lv = primitive2conserv((PrimitiveVars){phi_L,rho_L,Ug_L,Ul_L});
23
+  ConservVars rv = primitive2conserv((PrimitiveVars){phi_R,rho_R,Ug_R,Ul_R});
24
+  
25
+  grid.lbc.vars = lv;
26
+  grid.rbc.vars = rv;
27
+  
28
+  for (unsigned int i=0; i<grid.N; i++)
29
+    if (i<grid.N/2)
30
+      grid.cells[i].vars = lv;
31
+    else
32
+      grid.cells[i].vars = rv;
33 33
 }
34 34