When I constructed the initial and lateral boundary conditions for the model, HYCOM data were interpolated (hot start).
However, the surface temperature of the model was more than 4 degrees Celsius higher than HYCOM after only 5 days of simulation.
I tried turning off all the forcing files (Qair, Pair, Tair, winds, etc.) to find out why the surface temperature was increasing (I actually turned off bulk_fluxes, solar_source, etc. in the .h file), but the surface temperature still increased regardless of the forcing files.
I also simulated setting all initial temperatures to 16 degrees Celsius or 8 degrees Celsius, but this also generated high temperatures that did not come out of the boundary.
All lateral boundary conditions were well interpolated and applied as seen in the vertical profiles (boundary condition type: Clamped).
I have uploaded the input file (roms_ccs.in) andt the header file (ccs.h) I used.
Can someone please tell me why I am getting high temperature even though it is not applied to the boundary and what I should do to fix it?
Best regards,
Bo-Kyung Kim.
Code: Select all
/*
** svn $Id: scb.h 1099 2022-01-06 21:01:01Z arango $
*******************************************************************************
** Copyright (c) 2002-2022 The ROMS/TOMS Group **
** Licensed under a MIT/X style license **
** See License_ROMS.txt **
*******************************************************************************
**
** Options for Southern California Bight.
**
** Application flag: CCS
** Input script: roms_ccs.in
*/
#undef AFT_EIGENMODES /* Adjoint Finite Time Eigenmodes */
#undef CORRELATION /* Background-error Correlation Check */
#undef FORCING_SV /* Forcing Singular Vectors */
#undef FT_EIGENMODES /* Finite Time Eigenmodes */
#undef I4DVAR /* Incremental, strong constraint 4D-Var */
#define NLM_DRIVER /* Nonlinear Basic State trajectory */
#undef OPT_PERTURBATION /* Optimal perturbations */
#undef PICARD_TEST /* Picard Iterations Test */
#undef RBL4DVAR /* Strong/Weak constraint RBL4D-Var */
#undef R4DVAR /* Strong/Weak constraint R4D-Var */
#undef R_SYMMETRY /* Representer Matrix Symmetry Test */
#undef SANITY_CHECK /* Sanity Check */
#undef SO_SEMI /* Stochastic Optimals: Semi-norm */
#undef TLM_CHECK /* Tangent Linear Model Check */
/*
**-----------------------------------------------------------------------------
** Nonlinear basic state tracjectory.
**-----------------------------------------------------------------------------
*/
#if defined NLM_DRIVER
# define UV_U3ADV_SPLIT
# define UV_COR
# define UV_LDRAG
# define UV_VIS2
# define MIX_S_UV
# undef SPLINES_VDIFF
# undef SPLINES_VVISC
# define DJ_GRADPS
# define TS_U3HADVECTION
# define TS_C4VADVECTION
# undef TS_U3ADV_SPLIT
# define TS_DIF2
# define MIX_S_TS
# undef MIX_GEO_TS
# define SALINITY
# define NONLIN_EOS
# define CURVGRID
# define MASKING
# define SOLVE3D
# define AVERAGES
# define AVERAGES_QUADRATIC
# define DIAGNOSTICS_UV
# define DIAGNOSTICS_TS
# define LMD_MIXING
# undef MY25_MIXING
# ifdef LMD_MIXING
# undef SOLAR_SOURCE
# define LMD_RIMIX
# define LMD_CONVEC
# undef LMD_NONLOCAL
# undef LMD_SKPP
# define RI_SPLINES
# endif
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# undef ANA_RAIN
# undef FORWARD_MIXING
# undef FORWARD_READ
# undef FORWARD_WRITE
# define OUT_DOUBLE
# undef BULK_FLUXES
# undef LONGWAVE_OUT
# include "globaldefs.h"
#else
/*
**-----------------------------------------------------------------------------
** Adjoint-based drivers.
**-----------------------------------------------------------------------------
*/
# define UV_ADV
# define UV_COR
# define UV_LDRAG
# define UV_VIS2
# define MIX_S_UV
# define SPLINES_VDIFF
# define SPLINES_VVISC
# define DJ_GRADPS
# define TS_U3HADVECTION
# define TS_C4VADVECTION
# define TS_DIF2
# define MIX_S_TS
# undef MIX_GEO_TS
# define SALINITY
# define NONLIN_EOS
# define CURVGRID
# define MASKING
# define SOLVE3D
# define LMD_MIXING
# ifdef LMD_MIXING
# define SOLAR_SOURCE
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_NONLOCAL
# define LMD_SKPP
# define RI_SPLINES
# endif
# define ANA_BSFLUX
# define ANA_BTFLUX
# if defined RBL4DVAR || defined R4DVAR
# define CONVOLVE
# endif
# define VCONVOLUTION
# define IMPLICIT_VCONV
# undef FULL_GRID
# undef FORWARD_MIXING
# undef FORWARD_READ
# undef FORWARD_WRITE
# define OUT_DOUBLE
#endif
