[MITgcm-support] NaNs in the DIC package

[Tudino, Tobia]

Hi Steph,

thanks for the suggestions. I have checked at the other outputs and the salinity seems the only one that reaches unexpected low values in winter, but with a slower time scale in respect to the DIC disappearance.
However, I have checked the input and they seem fine.

Regards,

Tobia

Tobia Tudino
tt282 at exeter.ac.uk<mailto:tt282 at exeter.ac.uk>
PhD Student
Geography
College of Life and Environmental Science
Skype: tobia.tudino



On 25 Aug 2015, at 19:46, Stephanie Dutkiewicz <stephd at ocean.mit.edu<mailto:stephd at ocean.mit.edu>> wrote:


Hi Tobia -
Usually NaN in the DIC fields comes from the pH solver struggling with T or S reaching strange values -- so worth check both those fields to see if there is something wrong there.

You can try defining DIC_BOUNDS in DIC_OPTIONS.h -- this is a hack to try and make sure that the values going into the pH solver are within "reasonable bounds" -- but it is a hack -- it is worth checking why the fields (T,S or any of the DIC, ALK etc) are going out of bound as this might suggest other problems in the simulation.

steph

On Tue, 25 Aug 2015, Tudino, Tobia wrote:

Dear all,
I am currently run the MITgcm at a 2.8 resolution with the GCHEM package included. This run is based on a
previous setting constructed on an idealised forcing field.