About imposing Holland Model in ROMS on an idealized domain

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
CeeCee.X
Posts: 20
Joined: Thu Jul 17, 2014 1:31 pm
Location: PNNL, used to work and study at University of Florida

About imposing Holland Model in ROMS on an idealized domain

#1 Post by CeeCee.X » Tue Jan 13, 2015 7:04 pm

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.
Attachments
Capture_blowsup.PNG

User avatar
susonic
Posts: 160
Joined: Tue Aug 21, 2007 5:44 pm
Location: Jeju National University
Contact:

Re: About imposing Holland Model in ROMS on an idealized dom

#2 Post by susonic » Thu Jan 15, 2015 1:52 pm

Did you run the model without tide?

Did you run the model without holland wind input and put constant or 0 wind?

-JH

Post Reply