Index

Rundeck: P1SoM40

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

de_solver_main

Module:

File: ATURB_E1.f

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

---

Declaration:

subroutine de_solver_main(x, x0, p1, p4, rhoebydz, bydzerho, flux_bot, flux_top, dtime, n, qlimit)

Variables:

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

the unknown to be solved (at main 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.

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

coeff. of the d.e.

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

(j) rhoe(j+1)/dz(j)

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

(j) 1/(dze(j)*rho(j))

flux_bot : real*8, intent(in)

flux from the bottom

flux_top : real*8, intent(in)

flux at the top

dtime : real*8, intent(in)

time step

n : integer, intent(in)

number of main grid

qlimit : logical, intent(in)

true if tracer must be positive definite

tridiag : used from tridiag_mod

sub : real*8, dimension(n)

dia : real*8, dimension(n)

sup : real*8, dimension(n)

rhs : real*8, dimension(n)

j : integer

loop variable

alpha : real*8