SWAN & ROMS coupling problem

[dhson]

I'm trying to run ROMS with SWAN_COUPLING, WET_DRY, UV_PSOURCE, TS_PSOURCE,SEDIMENT and SSW_BBL. I believed they're all toughest modules. I've spent almost two months on playing with ROMS but not yet been successful:
- Without SWAN_COUPLING, I can get reasonable results
- With SWAN_COUPLING, there were a lot of troubles (I guess one of the main reason was that SWAN_COUPLING and WET_DRY can not be together):
+ For the swan.in, if I comment out the followed lines I can get reasonable wave field by SWAN

&& KEYWORD TO CREATE BOTTOM FRICTION GRID &&
INPGRID FRIC CURVILINEAR 0 0 349 499 EXC 9.999990e+005 &
NONSTAT 20061101.000000 25 DAY 20061126.000000

Otherwise Hwave all = 0 like this:
10 304 00:05:00 4.720785E-06 1.765845E+02 1.765845E+02 9.571772E+10
WAV2OCN - (04) imported and (05) exported fields, t = 20061101.000500
- SWAN coupling exchanges wait clock (s):
(Recv= 4.00000000E-02 Send= 0.00000000E+00)
OCN2WAV - (05) imported and (04) exported fields, t = 304 00:05:00
- ROMS coupling exchanges wait clock (s):
(Recv= 6.04280000E+02 Send= 0.00000000E+00)
- ROMS Import: wave direction
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Import: significant wave height
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Import: average wave length
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Import: surface wave relative peak period
(Min= 0.00000000E+00 Max= 9.99999762E-01)
- ROMS Import: bottom wave period
(Min= 0.00000000E+00 Max= 1.06307325E+01)
- ROMS Export: sea surface temperature
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Export: bathymetry
(Min= -2.00000038E+00 Max= 5.49999995E+01)
- ROMS Export: free-surface
(Min= -2.23851289E-02 Max= 2.10000001E+00)
- ROMS Export: vertically integrated u-momentum component
(Min= -1.35769737E-02 Max= 3.25842054E-01)
- ROMS Export: vertically integrated v-momentum component
(Min= -2.68944779E-01 Max= 9.66446239E-02)

+ ROMS blowed up as soon as the first writing out the history file:
11 304 00:05:30 5.445673E-06 1.765847E+02 1.765847E+02 9.571780E+10
12 304 00:06:00 6.296908E-06 1.765850E+02 1.765850E+02 9.571788E+10
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000002
13 304 00:06:30 NaN NaN NaN NaN

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001

Elapsed CPU time (seconds):

Node # 0 CPU: 837.940

*******************************
Here's my header fife:

#define SWAN_COUPLING

#ifdef SWAN_COUPLING

# define MCT_LIB
# define NEARSHORE_MELLOR
#endif

#define UV_ADV
#define UV_COR
#define UV_VIS2
#define MIX_GEO_UV

#define TS_A4HADVECTION
#define TS_A4VADVECTION
#undef TS_MPDATA
#define TS_DIF2
#define MIX_GEO_TS

#define DJ_GRADPS
#undef SALINITY
#define SOLVE3D
#define NONLIN_EOS

#define CURVGRID

#define SPLINES

#define MASKING

#define STATIONS
#define WET_DRY

#define UV_PSOURCE

#define TS_PSOURCE
#undef Q_PSOURCE

#define AVERAGES
#define RADIATION_2D
#define SPONGE

#define WESTERN_WALL


#define RAMP_TIDES
#define SSH_TIDES
#define UV_TIDES
#ifdef SSH_TIDES
# define FSOBC_REDUCED
# define ADD_FSOBC
# define EAST_FSCHAPMAN
# define NORTH_FSCHAPMAN
# define SOUTH_FSCHAPMAN
#else
# define EAST_FSGRADIENT
# define NORTH_FSGRADIENT
# define SOUTH_FSGRADIENT
#endif
#if defined UV_TIDES || defined SSH_TIDES
# define ADD_M2OBC
# define EAST_M2FLATHER
# define NORTH_M2FLATHER
# define SOUTH_M2FLATHER
#else
# define EAST_M2RADIATION
# define NORTH_M2RADIATION
# define SOUTH_M2RADIATION
#endif

#define EAST_M3RADIATION
#define NORTH_M3RADIATION
#define SOUTH_M3RADIATION

#define EAST_TRADIATION
#define NORTH_TRADIATION
#define SOUTH_TRADIATION
#define EAST_TNUDGING
#define NORTH_TNUDGING
#define SOUTH_TNUDGING

/* define only one of the following 5 */
#undef UV_LOGDRAG
#undef UV_QDRAG
#undef MB_BBL
#undef SG_BBL
#define SSW_BBL

#ifdef SG_BBL
# define SG_CALC_ZNOT
# undef SG_LOGINT
#endif
#ifdef MB_BBL
# define MB_CALC_ZNOT
# undef MB_Z0BIO
# undef MB_Z0BL
# undef MB_Z0RIP
#endif
#ifdef SSW_BBL
# define SSW_CALC_ZNOT
# undef SSW_LOGINT
#endif

#if defined MB_BBL || defined SG_BBL || defined SSW_BBL
# ifndef SWAN_COUPLING
# define ANA_WWAVE
# endif
#endif

#ifdef SOLVE3D
# define SEDIMENT
# ifdef SEDIMENT
# define SUSPLOAD
# define BEDLOAD_MPM
# define SED_MORPH
# define RIVER_SEDIMENT
# undef BEDLOAD_SOULSBY
# endif
# if defined SEDIMENT || defined SG_BBL || defined MB_BBL || defined SSW_BBL
# define ANA_SEDIMENT
# endif

# define BULK_FLUXES
# ifdef BULK_FLUXES
# define LONGWAVE
# define ANA_CLOUD
# define ANA_RAIN
# endif


# define ANA_BSFLUX

# define ANA_SSFLUX

# define ANA_BTFLUX

# define ANA_BPFLUX
# define ANA_SPFLUX


# define ATM_PRESS

# define QCORRECTION

# define SOLAR_SOURCE

# define MY25_MIXING
# undef GLS_MIXING
# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# undef CRAIG_BANNER
# undef CHARNOK
# undef ZOS_HSIG
# undef TKE_WAVEDISS
# endif

#endif

[dhson]

Sorry, thanks for any suggestion,