bad Vertical mixing, are my forcing data wrong?

[shiyf]

Dear ROMS users-

my forcing data is NCEP long term mean data which were averaged from 1981-2010;

Uwind: surface u wind
Vwind: surface v wind
Tair : air temperature at sigma level 0.995, (degC)
Qair: Relative Humidity at sigma level 0.995,(%)
Pair :Surface Pressure,(millibars)
Swrad: net solor shortwave radiation flux,(W/m^2)
lwrad_down:Downward Longwave Radiation Flux ,(W/m^2)
rain : rain fall rate get from pres.surface.mon.ltm.nc file,(kilogram meter-2 second-1),

my KPP option is:

#define UV_ADV
#define UV_COR

#define UV_VIS2
#define UV_QDRAG

#define TS_U3HADVECTION
#define TS_C4VADVECTION

#define TS_DIF4
#define SALINITY

#define NONLIN_EOS


#define SCORRECTION
#define SRELAXATION

#define DJ_GRADPS

#define BULK_FLUXES
#define SOLAR_SOURCE
#define COOL_SKIN
#define LONGWAVE_OUT
#define EMINUSP
#define DIURNAL_SRFLUX

#define SOLVE3D
!#define CURVGRID
#define MASKING
#define PROFILE

#define SPLINES
#define ANA_BSFLUX
#define ANA_BTFLUX

#define VISC_GRID
#define MIX_S_UV
#define DIFF_GRID
#define MIX_GEO_TS

#define LMD_MIXING
#ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LMD_BKPP
# define LMD_NONLOCAL
# define LMD_DDMIX
#endif

#define SPONGE
/*
#define M2CLIMATOLOGY
#define M3CLIMATOLOGY
#define TCLIMATOLOGY
#define ZCLIMATOLOGY

#define M2CLM_NUDGING
#define M3CLM_NUDGING
#define TCLM_NUDGING
*/
#define DEFLATE
#define INLINE_2DIO
#define PERFECT_RESTART

#define AVERAGES_FLUXES
#define AVERAGES

so, can anybody tell why i am getting bad Vertical mixing. Does radiations fields I am using are correct? is anything missing/wrong in cppdefs.h ?