sediment settling flux

Sediment modeling collaborators: issues, applications, information exchange

Moderators: arango, robertson, rsignell

Post Reply
Message
Author
Johnny.Zang
Posts: 4
Joined: Tue Dec 20, 2016 3:12 pm
Location: LSU

sediment settling flux

#1 Unread post by Johnny.Zang »

Dear ROMS users,

I am working on the sediment-biological model coupling using COAWST and I found a problem regarding those sediment and biological tracers who have vertical settling velocity.

In sediment transport model (sed_settling.F), sediment tracer concentration is updated after settling flux calculation:
DO i=Istr,Iend
DO k=1,N(ng)
t(i,j,k,nnew,indx)=t(i,j,k,nnew,indx)+(FC(i,k)-FC(i,k-1))
END DO
settling_flux(i,j,ised)=FC(i,0)
END DO

While in biological model (here I use nemuro.h), PON and Opal (these two bio tracers have settling flux like sediment) concentration is updated using
DO k=1,N(ng)
DO i=Istr,Iend
Bio(i,k,ibio)=qc(i,k)+(FC(i,k)-FC(i,k-1))*Hz_inv(i,k)
END DO
END DO

An additional Hz_inv is included to transfer flux (mmol N/m2) to concentration (mmol N/m3). I am wondering why Hz_inv is not included in sediment model? I think these two models use the same method to estimate vertical settling flux, and settling flux should be divided by grid thickness (*Hz_inv(i,k)) to estimate its contribution to the tracer concentration variation at certain grid.

Does anyone have ideas on that?

Best,

John

jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: sediment settling flux

#2 Unread post by jcwarner »

be careful in here.

a quick look i can see the bottom of the bio routine has:
... DO i=Istr,Iend
cff=Bio(i,k,ibio)-Bio_old(i,k,ibio)
t(i,j,k,nnew,ibio)=t(i,j,k,nnew,ibio)+cff*Hz(i,j,k)
so to put the Bio back into the tracer array, the Bio gets multiplied by Hz. so that is probably where the units get back to what is needed.

Johnny.Zang
Posts: 4
Joined: Tue Dec 20, 2016 3:12 pm
Location: LSU

Re: sediment settling flux

#3 Unread post by Johnny.Zang »

Dr. Warner,

Thank you for your reply. Now my understanding is t(i,j,k,nnew,ibio) is bio tracer concentration times Hz rather than bio tracer's concentration, and the units of t(i,j,k,nnew,ibio) and t(i,j,k,nstp,ibio) are different (mmol/m2 vs mmol/m3). Am I correct?

Best,

John

jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: sediment settling flux

#4 Unread post by jcwarner »

the tracer nnew time level is the "tracer * Hz". the tracer nstp time level is just the tracer value.
so yes, be careful.

Post Reply