I working on NEMURO model. In the last loop of the code, in the "Update global tracer variables" part, we have:

DO itrc=1,NBT

ibio=idbio(itrc)

DO k=1,N

DO i=Istr,Iend

t(i,j,k,nnew,ibio)=MAX(t(i,j,k,nnew,ibio)+(Bio(i,k,ibio)-Bio_bak(i,k,ibio))* Hz(i,j,k), 0.d0)

END DO

END DO

END DO

DO itrc=1,NBT

ibio=idbio(itrc)

DO k=1,N

DO i=Istr,Iend

t(i,j,k,nnew,ibio)=MAX(t(i,j,k,nnew,ibio)+(Bio(i,k,ibio)-Bio_bak(i,k,ibio))* Hz(i,j,k), 0.d0)

END DO

END DO

END DO

With:

t is the matrix of our variables,

Bio and Bio_bak temporary matrix of our variables.

Unities are concentration.

So I don't understand why we times Bio-Bio_bak by Hz, unit of Hz in meter. All ecosystem model coupled to ROMS Rutgers do the same (and no those coupled to ROMS Agrif).

Thank you in advance for your answer.

Cheers

Fanny