I have learn ROMs some time. I calculate 15 years,and get some results,
when I check the 15th year result which are shown in the below figure,
I am confused, I want to know Whether my result is reasonable,if not,how to improve my calculation?
I use force,ini,bry and grid files,
my cpp option:
/*
**
** Options for SHIYF.
**
** Application flag: SHIYF_BV
** Input script: ocean_shiyf.in
*/
#define NO_HIS
#undef NETCDF4
#undef PARALLEL_IO
#undef OFFLINE_FLOATS
#define CURVGRID
#define MASKING
#define NONLIN_EOS
#define SOLVE3D
#define SALINITY
#ifdef SOLVE3D
#define SPLINES
#endif
#undef WET_DRY
#undef T_PASSIVE
#ifdef T_PASSIVE
#define ANA_PASSIVE
#define TRC_PSOURCE
#define ANA_TRC_PSOURCE
#define AGE_PASSIVE
#endif
#define NO_WRITE_GRID
#undef OUT_DOUBLE
#define RST_SINGLE
#define AVERAGES
#define AVERAGES2
#ifdef SOLVE3D
#undef AVERAGES_DETIDE
#undef DIAGNOSTICS_TS
#endif
#undef DIAGNOSTICS_UV
#ifdef SOLVE3D
#define DJ_GRADPS
#endif
#define UV_ADV
#define UV_COR
#define UV_LDRAG
#undef UV_SADVECTION
#ifdef SOLVE3D
#define TS_U3HADVECTION
#define TS_C4VADVECTION
#undef TS_MPDATA
#endif
#define UV_VIS2
#undef UV_SMAGORINSKY
#undef VISC_3DCOEF
#define MIX_S_UV
#define VISC_GRID
#undef SPONGE
#ifdef SOLVE3D
#define TS_DIF2
#define MIX_GEO_TS
#define DIFF_GRID
#endif
#ifdef SOLVE3D
#define SOLAR_SOURCE
#define WTYPE_GRID
#define LMD_MIXING
#define LMD_RIMIX
#define LMD_CONVEC
#define LMD_SKPP
#undef LMD_BKPP
#define LMD_NONLOCAL
#define LMD_SHAPIRO
#undef LMD_DDMIX
#endif
#ifdef SOLVE3D
#define NCEP_FLUXES
#endif
#ifdef SOLVE3D
#define SCORRECTION
#undef QCORRECTION
#undef TCLIMATOLOGY
#undef TCLM_NUDGING
#endif
#define RADIATION_2D
#ifdef SOLVE3D
#define ANA_BSFLUX
#define ANA_BTFLUX
#else
#define ANA_SMFLUX
#endif