From b5d47585a8d5dac3bbc29bc1dcad161764a09504 Mon Sep 17 00:00:00 2001
From: che <che-email@yandex.ru>
Date: Wed, 15 Jul 2026 09:42:27 +0500
Subject: [PATCH] fix

---
 config.toml      | 6 +++---
 src/equations.rs | 2 +-
 src/main.rs      | 8 ++++----
 3 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/config.toml b/config.toml
index ce304b4..125b3bf 100644
--- a/config.toml
+++ b/config.toml
@@ -1,10 +1,10 @@
 ID = 0
 CreatedAt = 2026-07-11T12:52:12.752206984+05:00
 UpdatedAt = 2026-07-11T12:52:12.752206984+05:00
-rel = 1700.0
-rel_c = 1.0
+rel = 2000.0
+rel_c = 0.0
 le = 1.0
-pr = 1.0
+sc = 1.0
 pe = 1.0
 ma = 0.0
 time = 10.0
diff --git a/src/equations.rs b/src/equations.rs
index b0b3457..8495f4c 100644
--- a/src/equations.rs
+++ b/src/equations.rs
@@ -73,7 +73,7 @@ impl PhiEquation {
         let inv_h = 1.0 / H;
         let inv_h2 = inv_h * inv_h;
         for (i, j) in self.grid.inner_nodes() {
-            let archim = self.le * self.rel * dx!(temp, i, j) + self.rel_c * dx!(conc, i, j);
+            let archim = self.le * self.rel * dx!(temp, i, j) - self.rel_c * dx!(conc, i, j);
             let tmp = self.sc * (lap!(phi, i, j) * inv_h2 + archim * inv_h)
                 + jacobian!(psi, phi, i, j) * inv_h2;
             self.delta[[i, j]] = tmp * dt;
diff --git a/src/main.rs b/src/main.rs
index a613a53..8b47bff 100644
--- a/src/main.rs
+++ b/src/main.rs
@@ -120,12 +120,12 @@ impl System for S {
                 let x = (i as f64) * H;
                 let z = (j as f64) * H;
                 let pi = std::f64::consts::PI;
-                self.temp.f[[i, j]] = 1.0 - z + amp * (2.0 * pi * x / L).sin() * (pi * z).sin();
-                self.phi.f[[i, j]] = 0.0;
+                self.temp.f[[i, j]] = 1.0 - z; //+ amp * (2.0 * pi * x / L).sin() * (pi * z).sin();
+                self.phi.f[[i, j]] = 2.0;
                 self.psi.f[[i, j]] = 0.0;
                 if pe.abs() > 1e-5 {
-                    //self.conc.f[[i, j]] = -f64::exp(z * pe) / (1.0 - f64::exp(pe));
-                    self.conc.f[[i, j]] = 1.0 / pe.abs()
+                    self.conc.f[[i, j]] = -pe * f64::exp(z * pe) / (1.0 - f64::exp(pe));
+                    //self.conc.f[[i, j]] = 1.0
                 } else {
                     self.conc.f[[i, j]] = 0.0;
                 }