Ecosystem Model - ROMS

[fannychenillat]

Hello,

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

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

[kate]

The two versions of ROMS have different timestepping schemes. It's just part of the tracer timestep.

[fannychenillat]

The two versions of ROMS have different timestepping schemes. It's just part of the tracer timestep.

Actually, I try to work with NEMURO coupled to ROMS Agrif,

So If I change
(1) the loading of tracers arrays to
Bio(i,k,ibio)= max ( t(i,j,k,nnew,itrc), 0.0)
and
(2)the extracting of Bio to tracers arrays to
t(i,j,k,nnew,itrc)= min(t(i,j,k,nnew,itrc),0.0) + Bio(i,k,ibio)
i.e. as the same way used in ecosystem model of ROMS Agrif,
I should find the same results as NEMURO in ROMS Rutgers?

I mean: do not used t(i,j,k,nstp,ibio) in my ecosystem model is not a big deal for physics?

Thanks

Fanny

[kate]

You should do what is consistent for the model you are using.

The one true test is to try it...

[fannychenillat]

Hey

Thank you Kate. I have tried it and my nitrogen budget is constant.
So I guess I use the right way.

Cheers,

Fanny