Abnormal low temperature in bottom layer

[zhiguoMEI]

Hi everyone ! Recently I got a problem with my ROMS-only model.The temperature in bottom layer was below -10 deg C and this abnormal temp happened only in the slope region. If I change my CPP define TS_C4HADVCTION into TS_U3HADVTION this problem just disappear but other problems appear.Could you please give me some help about how to solve this abnormal temp problem while using TS_C$HADVECTION CPP define ? Thank you so much!

Below is my CPP define
/*Physics options*/
#define UV_ADV
#define UV_COR
#define UV_VIS2

#define VISC_GRID
#define DIFF_GRID
#define TS_DIF2
#define MIX_S_UV
#define MIX_GEO_TS
#define SALINITY

/*Sediment*/

/*Numeric options*/
#define CURVGRID
#define SPHERICAL
#define MASKING
#define SOLVE3D
#define SPLINES_VVISC
#define SPLINES_VDIFF
#define NONLIN_EOS
#define DJ_GRADPS

#define TS_U3HADVECTION
#undef TS_C4HADVECTION
#undef TS_MPDATA
#undef WET_DRY
#undef TS_U3ADV_SPLIT
#define UV_U3ADV_SPLIT
/*bottom*/
#define UV_QDRAG

/*Output*/
#define AVERAGES
#define PERFECT_RESTART
#define RST_SINGLE
/*Surface and bottom boundary options*/
#define BULK_FLUXES
#ifdef BULK_FLUXES
# define SOLAR_SOURCE
# define LONGWAVE_OUT
# define EMINUSP
# define WIND_MINUS_CURRENT
#else
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SRFLUX
#endif

#define ANA_BTFLUX
#undef ANA_SSFLUX
#define ANA_BSFLUX
#define ANA_SPFLUX
#define ANA_BPFLUX
# undef LIMIT_VDIFF
# undef LIMIT_VVISC

/* Vertical subgridscale turbulence closure */
# define GLS_MIXING

# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# define RI_SPLINES
# undef CRAIG_BANNER
# undef CHARNOK
# endif

# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_BKPP
# define LMD_SKPP
# undef LMD_NONLOCAL
# endif

Here is model result :

[gresendiz]

Hello! If the area is very shallow the problem can be the heat flux, you can try use the forcing and Qcorrection instead of Bulk Fluxes. I had a very similar problem in a shallow domain, the temperature was very high or very low.

[zhiguoMEI]

Thank you for your help ! I check the bathymetry where this problem happen. Below is a temperature slice along 25.7N near Taiwan island.. The depth of low temperature region is below 1600 meters , I don't think this is too shallow.

[zhiguoMEI]

Hello! If the area is very shallow the problem can be the heat flux, you can try use the forcing and Qcorrection instead of Bulk Fluxes. I had a very similar problem in a shallow domain, the temperature was very high or very low.

Thank you for your help ! I check the bathymetry where this problem happen. Below is a temperature slice along 25.7N near Taiwan island.. The depth of low temperature region is below 1600 meters , I don't think this is too shallow.