Skip to content

Commit d110bf7

Browse files
committed
Tprog convert for sum of tendencies on rhs
1 parent c5c1e79 commit d110bf7

File tree

3 files changed

+5
-5
lines changed

3 files changed

+5
-5
lines changed

src/Continuity.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -62,7 +62,7 @@ function continuity_itself!(Diag::DiagnosticVars{T,Tprog},
6262

6363
@inbounds for j 1:n
6464
for i 1:m
65-
dη[i+1,j+1] = -(Tprog(dUdx[i,j+1]) + Tprog(dVdy[i+1,j]))
65+
dη[i+1,j+1] = -(Tprog(dUdx[i,j+1] + dVdy[i+1,j]))
6666
end
6767
end
6868
end

src/Diffusion.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -234,7 +234,7 @@ function add_drag_diff_tendencies!( u::Array{Tprog,2},
234234

235235
@inbounds for j 1:n
236236
for i 1:m
237-
u[i+2,j+2] += Δt_diff*(Tprog(Bu[i+1-ep,j+1]) + Tprog(LLu1[i,j+1]) + Tprog(LLu2[i+1-ep,j]))
237+
u[i+2,j+2] += Δt_diff*(Tprog(Bu[i+1-ep,j+1] + LLu1[i,j+1] + LLu2[i+1-ep,j]))
238238
end
239239
end
240240

@@ -245,7 +245,7 @@ function add_drag_diff_tendencies!( u::Array{Tprog,2},
245245

246246
@inbounds for j 1:n
247247
for i 1:m
248-
v[i+2,j+2] += Δt_diff*(Tprog(Bv[i+1,j+1]) + Tprog(LLv1[i,j+1]) + Tprog(LLv2[i+1,j]))
248+
v[i+2,j+2] += Δt_diff*(Tprog(Bv[i+1,j+1] + LLv1[i,j+1] + LLv2[i+1,j]))
249249
end
250250
end
251251
end

src/rhs.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -296,7 +296,7 @@ function momentum_u!(Diag::DiagnosticVars{T,Tprog},S::ModelSetup) where {T,Tprog
296296

297297
@inbounds for j 1:n
298298
for i 1:m
299-
du[i+2,j+2] = Tprog(qhv[i,j]) - Tprog(dpdx[i+1-ep,j+1]) + Tprog(Fx[i,j])
299+
du[i+2,j+2] = Tprog(qhv[i,j] - dpdx[i+1-ep,j+1] + Fx[i,j])
300300
end
301301
end
302302
end
@@ -316,7 +316,7 @@ function momentum_v!(Diag::DiagnosticVars{T,Tprog},S::ModelSetup) where {T,Tprog
316316

317317
@inbounds for j 1:n
318318
for i 1:m
319-
dv[i+2,j+2] = -Tprog(qhu[i,j]) - Tprog(dpdy[i+1,j+1]) + Tprog(Fy[i,j])
319+
dv[i+2,j+2] = Tprog(-qhu[i,j] - dpdy[i+1,j+1] + Fy[i,j])
320320
end
321321
end
322322
end

0 commit comments

Comments
 (0)