From b70158fa97d6854c19c1190354f92e683206913d Mon Sep 17 00:00:00 2001
From: George Faykin <pofnutsy@gmail.com>
Date: Fri, 6 Jun 2025 19:06:53 +0300
Subject: [PATCH] =?UTF-8?q?=D0=98=D1=81=D0=BF=D1=80=D0=B0=D0=B2=D0=BB?=
=?UTF-8?q?=D0=B5=D0=BD=D0=B0=20=D0=BE=D1=88=D0=B8=D0=B1=D0=BA=D0=B0=20?=
=?UTF-8?q?=D0=B2=20=D1=80=D0=B0=D1=81=D1=87=D1=91=D1=82=D0=B5=20=D0=B0?=
=?UTF-8?q?=D0=B4=D0=B2=D0=B5=D0=BA=D1=86=D0=B8=D0=B8:=20=D0=BE=D1=88?=
=?UTF-8?q?=D0=B8=D0=B1=D0=BA=D0=B0=20=D0=B3=D1=80=D0=B0=D0=BD=D0=B8=D1=87?=
=?UTF-8?q?=D0=BD=D0=BE=D0=B3=D0=BE=20=D0=B7=D0=BD=D0=B0=D1=87=D0=B5=D0=BD?=
=?UTF-8?q?=D0=B8=D1=8F=20=D1=8D=D0=BB=D0=B5=D0=BC=D0=B5=D0=BD=D1=82=D0=B0?=
=?UTF-8?q?=20=D1=81=D0=B5=D1=82=D0=BA=D0=B8=20=D0=BF=D0=B5=D1=80=D0=B2?=
=?UTF-8?q?=D0=BE=D0=B3=D0=BE=20=D1=8D=D0=BB=D0=B5=D0=BC=D0=B5=D0=BD=D1=82?=
=?UTF-8?q?=D0=B0?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
advection.f90 | 8 +++++---
neodn2.f90 | 12 +++++++-----
2 files changed, 12 insertions(+), 8 deletions(-)
diff --git a/advection.f90 b/advection.f90
index 2c502b2..018d5b6 100644
--- a/advection.f90
+++ b/advection.f90
@@ -8,7 +8,7 @@ program advection_solver
! Параметры сетки
integer, parameter :: N = 100 ! Число узлов по пространству
- integer, parameter :: M = 1000 ! Число шагов по времени
+ integer, parameter :: M = 100 ! Число шагов по времени
real, parameter :: dt = 0.1 ! Шаг по времени
real :: Hsoil = 1.E+3 ! Размер области по пространству
real :: dz ! Шаг по пространству
@@ -58,7 +58,7 @@ program advection_solver
! Вычисление аргумента сетки
dz = Hsoil/real(N)
! Постоянная сетка по пространству
- !forall (i=1:N) z_cen(i) = i*dz
+ forall (i=1:N) z_cen(i) = i*dz
! Неоднородная сетка по пространству
!forall (i=1:N) z_cen(i) = (i/5.)**2.
@@ -67,7 +67,9 @@ program advection_solver
! Вычисление граничного значения ячейки
forall (i=2:N) z_edg(i) = (z_cen(i-1) + z_cen(i))/2.
- z_edg(1) = -z_edg(2) - z_cen(1)
+ if (z_cen(1) == 0.) z_edg(1) = z_cen(1) - z_edg(2)
+
+ z_edg(1) = 2.*z_cen(1) - z_edg(2)
z_edg(N+1) = 2.*z_cen(N) - z_edg(N)
! Вычисление шага по пространственной сетке
diff --git a/neodn2.f90 b/neodn2.f90
index 21a6326..46b0881 100644
--- a/neodn2.f90
+++ b/neodn2.f90
@@ -106,17 +106,19 @@ program neodn2
! ----------------------------------------------- расчёт сетки --------------------------------------------------------
! Вычисление аргумента сетки
dz = Hsoil/real(N)
- ! Постоянная сетка по пространству
+ ! Постоянная сетка по пространству
forall (i=1:N) z_cen(i) = i*dz
! Неоднородная сетка по пространству
- !forall (i=0:N) z_edg(i) = (i/5.)**2.
- forall (i=0:N) z_edg(i) = Hsoil * ((1.0 + 12.0)**((real(i) - 1.0)/(real(N) - 1.0)) - 1.0) / 12.0
+ !forall (i=1:N) z_cen(i) = (i/5.)**2.
+ forall (i=1:N) z_cen(i) = Hsoil * ((1.0 + 12.0)**((real(i) - 1.0)/(real(N) - 1.0)) - 1.0) / 12.0
- ! Вычисление центрального значения ячейки
+ ! Вычисление граничного значения ячейки
forall (i=2:N) z_edg(i) = (z_cen(i-1) + z_cen(i))/2.
- z_edg(1) = z_edg(2) - z_cen(1)
+ if (z_cen(1) == 0.) z_edg(1) = z_cen(1) - z_edg(2)
+
+ z_edg(1) = 2.*z_cen(1) - z_edg(2)
z_edg(N+1) = 2.*z_cen(N) - z_edg(N)
! Вычисление шага по пространственной сетке
--
GitLab