Parameters in my model

[xwbin]

Hi,
I wanna to make my model closer to the hycom in order to run some sensitive experiments.I use the hycom data as the boundary and initial condition.The only force is wind from CFRS.
The outcome is bad,and I suspect that I don't set the parameters well because most of them are difficult to undertand. I also upload the model output as a picture which compares to the hycom.You can see the Kuroshio intrusion is so strong even form an eddy in the western of Luzon strait.
Here are my parameters.I wonder that someone could give me advice and point out some obvious mistakes.

5.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
5.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
-1.000 gls_p GLS stability exponent.
0.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
5.5500E-01 gls_c1 GLS shear production coefficient.
8.3300E-01 gls_c2 GLS dissipation coefficient.
-6.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
2.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
2.0000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
3.0000E-03 Zob Bottom roughness (m).
1.0000E-03 Zos Surface roughness (m).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
300.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1028.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.00 time_ref Reference time for units attribute (yyyymmdd.dd)
2.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
2.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
1.0000E-03 Znudg Nudging/relaxation time scale (days)
for free-surface.
1.0000E-03 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
1.0000E-03 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
5.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
10.000 T0 Background potential temperature (C) constant.
30.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).

Here is my .h file.
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#undef UV_LOGDRAG
#undef UV_SMAGORINSKY
#define UV_VIS2
#ifdef UV_VIS2
# define VISC_GRID
# define MIX_S_UV
#endif

#define TS_MPDATA
#undef TS_DIF2
#ifdef TS_DIF2
# define DIFF_GRID
# define MIX_S_TS
#endif

#define DJ_GRADPS
#define NONLIN_EOS
#define SALINITY
#define SOLVE3D
#define SPLINES
#define MASKING
#define CURVGRID
#define SPHERICAL

#define AVERAGES
#define AVERAGES_AKV
#define AVERAGES_AKT
#define AVERAGES_AKS
#define AVERAGES_FLUXES

#undef DIAGNOSTICS_UV
#undef DIAGNOSTICS_TS

#undef FSOBC_REDUCED
#undef EAST_FS_REDUCED
#undef WEST_FS_REDUCED
#undef NORTH_FS_REDUCED

#define PERFECT_RESTART

#undef TCLIMATOLOGY
#undef TCLM_NUDGING

#define RADIATION_2D

#define GLS_MIXING
#undef MY25_MIXING
#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# undef CANUTO_A
# define N2S2_HORAVG
#endif

#undef SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
# undef SED_DENS
# undef wsed_power
# undef no_deposition
#endif

#undef ANA_PSOURCE
#undef ANA_GRID
#undef ANA_INITIAL
#undef ANA_SEDIMENT
#undef ANA_SMFLUX
#define ANA_STFLUX
#define ANA_BTFLUX
#define ANA_SSFLUX
#define ANA_BSFLUX
#define ANA_SPFLUX
#define ANA_BPFLUX
#undef ANA_FSOBC
#undef ANA_M2OBC
#undef ANA_M3OBC
#undef ANA_TOBC
#undef ANA_TCLIMA
#undef ANA_VMIX

Any suggestions are welcome!
Thanks