“…Where x, y, and z are the coordinates of the Cartesian system; 𝑢𝑢, 𝑣𝑣, and 𝑤𝑤 are the velocities in the directions of the Cartesian system, p is the water pressure, t is time, ε xx , ε xy , ε xz , ε yx , ε yy , ε yz , ε zx , ε zy , ε zz are the eddy turbulence coefficients, g is the acceleration due to gravity, 𝐷𝐷 𝑥𝑥 , 𝐷𝐷 𝑌𝑌 , and 𝐷𝐷 𝑧𝑧 are the eddy diffusion coefficients, 𝜌𝜌 is the density of water, Γ x , Γ y , and Γ z are the external forces; 𝜃𝜃 is the salinity, and 𝜃𝜃𝜃𝜃 is the salinity source/sink. The RMA11 model, which uses the finite element method, has been used successfully to model water quality in estuaries, bays, lakes, and rivers [15][16][17][18]. It receives the velocity, temperature, and salinity fields from the RMA10 model and uses them to solve the advection-diffusion constituent transport equations [15].…”