Index

Rundeck: P1SoM40

Created: Thu May 9 03:30:17 EDT 2024

de_solver_edge

Module:

File: ATURB_E1.f

Summary: differential eqn solver for x using tridiagonal method d/dt x = d/dz (P1 d/dz x) - P3 x + P4
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine de_solver_edge(x, x0, p1, p3, p4, rhobydze, bydzrhoe, x_surf, dtime, n)

Variables:

x : real*8, dimension(n), intent(out)

the unknown to be solved (at edge drid)

x0 : real*8, dimension(n), intent(in)

x at previous time step

p4 : real*8, dimension(n), intent(in)

coeff. of the d.e.

p3 : real*8, dimension(n), intent(in)

coeff. of the d.e.

p1 : real*8, dimension(n), intent(in)

coeff. of the d.e.

rhobydze : real*8, dimension(n), intent(in)

(j) rho(j)/dze(j)

bydzrhoe : real*8, dimension(n), intent(in)

(j) 1/(dz(j-1)*rhoe(j))

x_surf : real*8, intent(in)

surface value of x

dtime : real*8, intent(in)

time step

n : integer, intent(in)

number of vertical edge grid

tridiag : used from tridiag_mod

teeny : used from constant

j : integer

loop variable

sub : real*8, dimension(n)

dia : real*8, dimension(n)

sup : real*8, dimension(n)

rhs : real*8, dimension(n)