Abnormal warm on seawater temperature.

[cyhsu]

Hi, All,

I tried to use the atmospheric flux results from one of my previous runs to activate a new simulation. I had checked the shflux, swrad, sustr, svstr from either of my new run and my old run, they are fully equal. I can't check the EminusP with swflux, because the new run did not output it.

However, the conversion from EminusP to swflux is

swflux (cm day-1) = EminusP (m s-1) * 100/86400.

Since the all of the atmospheric forcing to the new simulation is from old run, sustr and svstr are in u and v coordinate.

This is my header file

Code: Select all

/* Model Configuration Options */
# define SOLVE3D
# undef CURVGRID
# define MASKING
# define SPLINES            
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_TS

# define AVERAGES
# undef RST_SINGLE
# undef QCORRECTION
# define SCORRECTION
# define MODIFY_SALT_TNUDG
# undef SRELAXATION
# undef SPONGE

/* Pressure Gradient Algorithm */
# define DJ_GRADPS

/* Tracer Equation Options */
# undef T_PASSIVE
# define SALINITY
# define NONLIN_EOS

/* Tracer advection */
# undef TS_MPDATA /* Reduces checkerboard pattern compared to  C4HADVECTION */
# define TS_A4HADVECTION /* Hernan's Recomendation */
# define TS_A4VADVECTION /* Hernan's Recomendation */

/* Tracer horizontal mixing */
# define DIFF_GRID
# define TS_DIF2
# define MIX_GEO_TS

/* Momemtum horizontal mixing */
# define VISC_GRID
# define MIX_S_UV
# define UV_VIS2

/* Momentum Equation Options */
# define UV_ADV
# define UV_COR
# define UV_QDRAG

/* turbulence mixing scheme */
# define LMD_MIXING   /* KPP turned on */
# ifdef LMD_MIXING
#  define LMD_RIMIX
#  define LMD_CONVEC
#  define LMD_SKPP
#  define LMD_BKPP
#  define LMD_NONLOCAL
# endif
# undef MY25_MIXING   /* MY25 turned off */
# ifdef MY25_MIXING
#  define N2S2_HORAVG
#  define KANTHA_CLAYSON
# endif

# undef  GLS_MIXING
# ifdef GLS_MIXING
#  define KANTHA_CLAYSON
#  undef  CANUTO_A
#  define N2S2_HORAVG
# endif

/* surface forcing */
# undef BULK_FLUXES  /*# define BULK_FLUXES*/ /* original set up */
# ifdef BULK_FLUXES
#  define LONGWAVE_OUT
#  define EMINUSP
#  define SOLAR_SOURCE
#  undef  DIURNAL_SRFLUX
# else
#  define ANA_SSFLUX
#  define ANA_SMFLUX
#  define ANA_STFLUX
# endif

/* define River flows */
# undef  UV_PSOURCE
# undef  TS_PSOURCE

/* Select  Biological model option */
# undef RIVER_BIOLOGY

# undef NPZD_POWELL
# ifdef NPZD_POWELL
#  define ANA_SPFLUX
#  define ANA_BPFLUX
#  undef CONST_PAR
#  define SPITZ
# endif

/*define boundary conditon */
# define RADIATION_2D

/* bottom surface ANA Flux */
# define ANA_BSFLUX
# define ANA_BTFLUX

/* define water mass relaxation */
# undef TCLIMATOLOGY
# undef TCLM_NUDGING

/* Jaisons Modifications */
# undef RELATIVE_WIND

/* Memory Management */
# undef INLINE_2DIO

Here I attached one time step of the ROMS simulation output on "sea surface seawater temperature (temp)"
https://www.dropbox.com/s/xddyd8zkt4szj ... M.png?dl=0


It turns out, If I let model has a longer run, I can see the sea surface temperature is reaching to 80 degree Celsius. Could anyone give me suggestions how to solve this abnormal warming?

[hwusklec]

Over-heating or over-cooling might occur occasionally when using prescribed heat flux to drive the model, especially in shallow area, to my own experience. Let's say a positive heat flux is prescribed, i.e. the sea receives the heat from air. In reality the ocean temperature will rise thus reduce the temperature contract with air, as a consequence the sensible heat flux will decrease, hence the over-heating will not happen. Or the warmed water is transported to elsewhere thus over-heating also unlikely happens.
However, if the heat flux is fixed and the flushing time of ocean water is long, the water temperature will increase continuously. Especially if the water depth is shallow the over-heating will be severe.

May be you can use bulk formula to calculate the heat flux internally, this will probably resolve your problem. Sorry I cannot see your dropbox file, blocked by the fire wall ..., so don't know if I have answered your question.