A simple question about the t3dmix2_s

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Posts: 39
Joined: Wed Jun 25, 2008 2:49 am
Location: Georgia Institute of Technology

A simple question about the t3dmix2_s

#1 Post by seashellingolds »

Hi All

I happened to look at the code in the t3dmix2_s.h and met something I couldnt understand.

Code: Select all


similar for FE

What I wanted to ask is why cff3 in the last line doesnt need to be divided by Hz? (FX and FE have been multiplied by Hz). Do I miss anything here?

Thanks a lot

User avatar
Posts: 185
Joined: Fri Nov 14, 2003 4:57 pm

Re: A simple question about the t3dmix2_s

#2 Post by shchepet »

Because at this moment t(i,j,k,nnew,itrc) has the meaning of tracer concentration
multiplied by Hz(i,j,k). Since Hz changes because of free surface evolution, tracers
are multiplied by Hz at time step "n" before the barotropic time stepping begins for
the current 3D time step, then all fluxes are added in to the product, and finally the
outcome is divided by Hz, which happens implicitly during the tri-diagonal solver
associated with vertical diffusion.

This algorithm allows to avoid extra storage associated with Hz because only one
time frame is needed for it.

User avatar
Posts: 572
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University

Re: A simple question about the t3dmix2_s

#3 Post by wilkin »

Immediately above the code segment you quote there is the comment:

Code: Select all

! Time-step harmonic, S-surfaces diffusion term (m Tunits).
At this place in the code the tracer array is in flux units of layer thickness (Hz, in m) times Tracer units (Tunits). The return to simple tracer units occurs later.

See also https://www.myroms.org/wiki/index.php/Horizontal_Mixing
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply