[MITgcm-support] Memory issues with regional simulation

[Ruth Musgrave]

Hi everyone,

I’m trying to run a regional model with 1e8 grid points on 64 processors. I have 1Gb of RAM per processor. Sadly with the number of full 3D diagnostics I would like to output, 1Gb of RAM isn’t enough. I think the diagnostics package is allocating double precision arrays for each of my diagnostics, which is what is causing my problems.

I wondered if there is a simple way to reduce the precision of the qdiag array in the diagnostics package. I had a look at it and I can only see very complicated ways of doing this (i.e. I tried and was unsuccessful). Is there something easy that I’m missing?

As a last resort I wondered if I could just run the entire model in single precision by making some changes to CPP_EEMACROS.h. It seems a bit desperate, and I don’t really have a sense of how errors would accumulate over a 1e5 timestep simulation. I’d be very grateful to hear other people’s experience and opinions on doing this.

If anyone has any other suggestions on how to reduce memory usage by the model (other than reduce grid points!) I’d be very happy to hear them.

Thanks in advance,

Ruth