Hi,
I was wondering if anyone has the experience of imposing Holland wind field in ROMS.
I have an idealized rectangular domain (250 km * 350 km, resolution = 2km) with tidal forcing propagating in from the east boundary and west boundary closed. I created a forcing netcdf file using Holland model and tried to add it into my tidal simulation. However, since the netcdf requires realistic domain, I have to assign lat. and long. to my idealized domain. Now the Lat/lon of the four edges of my domain are 28.2-29.80/5.00-7.80. I also set spherical=T.
Yet, from now on, with or without the wind forcing, the simulation blows up at the first step, exactly at 10min40sec and almost at the same cell. I attached to figure to illustrate the blows up pattern from my restart file.
Here are my cpp option:
#define ROMS_MODEL
#define UV_VIS2
#define MIX_S_UV
#define MASKING
#define UV_ADV
#define UV_COR
#define TS_MPDATA
#define DJ_GRADPS
#define WET_DRY
#undef SOLVE3D
#define SPLINES
#undef BODYFORCE
#define DIAGNOSTICS_UV
#define ANA_INITIAL
#undef ANA_WINDS
#undef ANA_SMFLUX
#define ANA_FSOBC
#define ANA_M2OBC
/* define only one of the following 5 */
#undef UV_LOGDRAG
#define UV_QDRAG
#undef MB_BBL
#undef SG_BBL
# undef SSW_BBL
#ifdef SSW_BBL
# define SSW_CALC_ZNOT
#endif
#ifdef SOLVE3D
# define GLS_MIXING
# ifdef GLS_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif
# undef SEDIMENT
# ifdef SEDIMENT
# define SUSPLOAD
# undef BEDLOAD_SOULSBY
# undef BEDLOAD_MPM
# define SED_MORPH
# endif
# if defined SEDIMENT || defined SG_BBL || defined MB_BBL || defined SSW_BBL
# define ANA_SEDIMENT
# endif
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_BPFLUX
# define ANA_BTFLUX
# define ANA_BSFLUX
# define ANA_SPFLUX
# define ANA_SRFLUX
#endif
And I also tried to the combinations of (Chapman + Flather or Claped + Reduced) for my east boundary condition, and for south and north boundary I used gradient option.
However, neither of them works.
I was wondering if anyone can give me some suggestions to how to determine and solve the blows up problem.
About imposing Holland Model in ROMS on an idealized domain
Re: About imposing Holland Model in ROMS on an idealized dom
Did you run the model without tide?
Did you run the model without holland wind input and put constant or 0 wind?
-JH
Did you run the model without holland wind input and put constant or 0 wind?
-JH
Joonho Lee