I tried two cases, one with 'define bulk_flux', the other one without bulk_flux. All the atmospheric forcing are from ECMWF-ERA-interim. Simulated salinity of two cases are similar.
The simulated surface temperature without 'bulk_flux' is similar with SST of JPL data. But the simulated temperature with 'bulk_flux' is overestimated (see figures). The only difference of the two headfiles is that 'bulk_flux' was defined in 'with bulk_fulx' case. So the problem should be related to the 'bulk_flux' formula. What should I do about it?
I saw a post viewtopic.php?f=14&t=3066 about overestimating the temperature. Should I add 'define ALBEDO' in the 'with bulk_flux' case?
Thanks.
headfile of 'with bulk_flux' case:
Code: Select all
/* %%BIOLOGICAL OPTIONS%% */
#define BIO_FENNEL
#ifdef BIO_FENNEL
# define SOLAR_SOURCE
# define DIURNAL_SRFLUX
# undef ANA_CLOUD
# define OXYGEN
# undef PO4
# undef BIO_SOD
# define BIO_SEDIMENT
# define DENITRIFICATION
# define RIVER_BIOLOGY
# define RIVER_PON
# define DIAGNOSTICS_BIO
# define ANA_SPFLUX
# define ANA_BPFLUX
#endif
#define UV_ADV
#define UV_COR
#define UV_VIS2
#define MIX_S_UV
#define TS_MPDATA
#define DJ_GRADPS
#define TS_DIF2
#define MIX_S_TS
#define SALINITY
#define MASKING
#define SOLVE3D
#define AVERAGES
#define CURVGRID
#define NONLIN_EOS
#define UV_QDRAG
/* %%%%%%%%%%% added; use bulk_fluxes%%%%%*/
#define VISC_GRID
#define DIFF_GRID
#define SPLINES
#define BULK_FLUXES
#define LONGWAVE
#define EMINUSP
#define ANA_BTFLUX
#define ANA_BSFLUX
#define ANA_NUDGCOEF /* %% nudging coefficients, ana_nudgcoef.h %% */
#undef ANA_SPONGE /* %% sponge coefficients, ana_sponge.h, ini_hmixcoef.F %% */
#define QCORRECTION
#undef SCORRECTION
/* %%%%%%%%%%%%% vertical turbulent mixing %%%%%%%%%%%%*/
#define MY25_MIXING
# define N2S2_HORAVG
/* %%%%%%%%%%%%%%%%%% tide %%%%%%%%%%%%%%%%%*/
# define SSH_TIDES
# define UV_TIDES
# ifdef SSH_TIDES /* if turn on SSH_TIDES */
# define ADD_FSOBC
# else
# define ANA_FSOBC
# endif
# ifdef UV_TIDES /* if turn on UV_TIDES */
# define ADD_M2OBC
# else
# define ANA_M2OBC
# endif
/* boundary condition */
#define RADIATION_2D
#define PERFECT_RESTART
#undef DIAGNOSTICS_UV
#undef DIAGNOSTICS_TS