High density at deeper layers

[jayasankar]

Hello All,

My simulation yields very high densities (around 1045 kg/m3) at deeper layers.
There is no model instability though. Is this acceptable? If not what could be the reasons and how it be taken care of?
I have used CFSR datasets for surface fluxes and SODA data for boundary fluxes.
Thanks in advance.

My model configuration is

/* Model dynamics */
# define SOLVE3D
# define UV_COR
# define UV_ADV
/* Equation of State */
# define SALINITY
# define NONLIN_EOS
# define SPLIT_EOS
/* Lateral Explicit Momentum Mixing */
# undef UV_VIS2
# ifdef UV_VIS2
# define UV_MIX_S
# define UV_VIS_SMAGO
# endif
/* Lateral Tracer Advection (default UP3) */
# define TS_HADV_RSUP3
# undef TS_HADV_UP5
# undef TS_HADV_C4
# undef TS_HADV_WENO5
/* Lateral Explicit Tracer Mixing */
# undef TS_MIX_S
# ifdef TS_HADV_C4
# define TS_DIF2
# undef TS_DIF4
# define TS_DIF_SMAGO
# define TS_MIX_ISO
# endif
/* Sponge layers for UV and TS */
# define SPONGE
/* Vertical Mixing */
# undef BODYFORCE
# undef BVF_MIXING
# define LMD_MIXING
# undef GLS_MIXING
# ifdef LMD_MIXING
# define LMD_SKPP
# define LMD_BKPP
# define LMD_RIMIX
# define LMD_CONVEC
# undef LMD_DDMIX
# define LMD_NONLOCAL
# endif
# ifdef GLS_MIXING
# define GLS_KKL
# undef GLS_KOMEGA
# undef GLS_KEPSILON
# undef GLS_GEN
# undef KANTHA_CLAYSON
# undef CRAIG_BANNER
# undef CANUTO_A
# undef ZOS_HSIG
# endif

# undef ANA_SST
# undef ANA_SSS

[kate]

The model variables are potential temperature and in situ density. If you want potential density, that would be pden, though the model does not (yet) write that out. Anyway, your density does not sound unreasonable for deep ocean basins. How deep are you talking about?

[jayasankar]

Thanks a lot Kate,

The potential density is within range. However the model computed density based on equation of state for seawater(Jackett and McDougall,1995) is around 1045kg/m3 at depths greater than 4000 meters.

[shifei]

Thank you, Kate