I have a 0.6x0.6 degrees simple grid setup with closed boundaries and constant initial fields for tracers (temp=17C, salt=20ppt).
I also use constant values for the BULK_FLUXES (Tair=23.567 C).
There is no wind or surface momentum flux and no salt or freshwater input from the surface.
When I check the his file I see high temperature values at the surface.
In the following figure red color is for the initial value and the purple is from the latest timestep (after 7 days).
- surf-temp-problem.gif (4.9 KiB) Viewed 1704 times
What might be the cause for this ?
I have the version 594 and my application.h file is as follows.
#undef UV_ADV /* no advection */
#undef UV_COR /* no coriolis */
#define MIX_GEO_TS /* mixing on geopotential (constant Z) surfaces */
#define MIX_GEO_UV /* mixing on geopotential (constant Z) surfaces */
#define UV_QDRAG /* quadratic bottom friction */
#undef UV_VIS2 /* turn ON harmonic horizontal mixing */
#undef TS_DIF2 /* turn ON harmonic horizontal mixing */
#define SALINITY /* salinity */
#define ANA_SSFLUX /* Analytical surface salinity flux */
#define SOLVE3D
#define SPHERICAL
#define MASKING
#define ANA_INITIAL /* analytical initial conditions */
/*#undef ANA_VMIX*/
#define BULK_FLUXES /* Calculate surface wind stress and surface net heat flux */
#ifdef BULK_FLUXES /* @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ */
#define ANA_PAIR
#define ANA_TAIR
#define ANA_HUMIDITY
#define ANA_CLOUD
#define ANA_RAIN
#define ANA_WINDS /* analytical surface winds */
#define ALBEDO /* Calculate shortwave radiation (define it before ana_srflux) */
#define ANA_SRFLUX /* Analytical surface shortwave radiation flux */
#define LONGWAVE /* use Berliand relation */
/*#undef ANA_SMFLUX /* Analytical surface momentum stress */
#else /* @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ */
#define ANA_SMFLUX /* Analytical surface momentum stress */
#define ANA_STFLUX /* Analytical surface temperature flux */
#endif /* @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ */
#define ANA_BSFLUX /* analytical bottom salinity flux */
#define ANA_BTFLUX /* analytical bottom temperature flux */
Regards,
ozgur