**Forcing Flux**input files to simulate my model. Comparing with SODA Climatology, The temperature results are not so good(shown below).

Then be advised to using

**Bulk Flux**input files, I open #define BULK_FLUXES in .h files. The result in the first month is close to SODA, but it gets worse in following months gradually(shown below).

Expect for #define BULK_FLUXES, what else cpps should I define?

Following is Activated C-preprocessing Options in my .h file.

oceanM Northwest Pacific Model

ANA_BSFLUX Analytical kinematic bottom salinity flux.

ANA_BTFLUX Analytical kinematic bottom temperature flux.

ASSUMED_SHAPE Using assumed-shape arrays.

AVERAGES Writing out time-averaged nonlinear model fields.

BULK_FLUXES Surface bulk fluxes parameterization.

CURVGRID Orthogonal curvilinear grid.

DIFF_GRID Horizontal diffusion coefficient scaled by grid size.

DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.

DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).

DOUBLE_PRECISION Double precision arithmetic.

EMINUSP Compute Salt Flux using E-P.

INLINE_2DIO Processing 3D IO level by level to reduce memory needs.

KANTHA_CLAYSON Kantha and Clayson stability function formulation.

MASKING Land/Sea masking.

MIX_GEO_TS Mixing of tracers along geopotential surfaces.

MIX_S_UV Mixing of momentum along constant S-surfaces.

MPI MPI distributed-memory configuration.

MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.

NONLINEAR Nonlinear Model.

NONLIN_EOS Nonlinear Equation of State for seawater.

NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files.

N2S2_HORAVG Horizontal smoothing of buoyancy and shear.

POWER_LAW Power-law shape time-averaging barotropic filter.

QCORRECTION Surface net heat flux correction.

K_GSCHEME Third-order upstream advection of TKE fields.

!RST_SINGLE Double precision fields in restart NetCDF file.

SALINITY Using salinity.

SOLVE3D Solving 3D Primitive Equations.

SPLINES Conservative parabolic spline reconstruction.

TS_U3HADVECTION Third-order upstream horizontal advection of tracers.

TS_C4VADVECTION Fourth-order centered vertical advection of tracers.

TS_DIF2 Harmonic mixing of tracers.

UV_ADV Advection of momentum.

UV_COR Coriolis term.

UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.

UV_C4VADVECTION Fourth-order centered vertical advection of momentum.

UV_QDRAG Quadratic bottom stress.

UV_VIS2 Harmonic mixing of momentum.

VAR_RHO_2D Variable density barotropic mode.

VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Any suggestion are welcomed, thank you in advance.