Ocean Modeling Discussion

Both is the case as NDIA=NRST and NDIA is defined anyway as an integer factor applied to DT.

This reply only partly applies to my setup. Indeed, my NAVG=0 (=> no output) differs from my NRST=1440 (=> daily restart writing). However, AVERAGES is switched off and DIAGNOSTICS_TS and DIAGNOSTICS_BIO are switched on, though, with NDIA=NRST=1440 . Hence, I would expect that the diagnostics output ...

It seems the DEF_HIS logic for determining the output record (see below) also has a flaw. ! ! Set unlimited time record dimension to the appropriate value. ! IF (ndefHIS(ng).gt.0) THEN HIS(ng)%Rindex=((ntstart(ng)-1)- & & ndefHIS(ng)*((ntstart(ng)-1)/ndefHIS(ng)))/ & & nHIS(ng) ELSE HIS(ng)%Rindex ...

Virivarola wrote: Specifically, in fennel_mod.h the text headings tell us they are [mg_C/(mg_Chl W m-2 day)]

Which version of ROMS do you use? In my fennel_mod.h (downloaded in August this year) the unit is [1/(W m-2 day)] consistently, incl. the heading.

Best,
Fabian

Hello, There seems to be a bug in the DEF_DIAGS at the point where the output time step for the diagnostics file is determined after a restart of ROMS. In my code (which is from August this year) lines 721-734 read as: ! ! Set unlimited time record dimension to the appropriate value. ! IF (nRST(ng ...

Hello, I am not exactly sure whether this is a real bug, but I would consider the following problem at least as an avoidable inconsistency. I use ROMS version 854 with the Fennel biology on the northern Gulf of Mexico shelf. The Mississippi and Atchafalaya rivers are implemented as horizontal point ...