Problem in small scale and large scale

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
liangliang
Posts: 61
Joined: Mon Jul 20, 2009 2:41 pm
Location: Port And Costal Engineering Laboratory

Problem in small scale and large scale

#1 Post by liangliang » Fri Nov 22, 2013 10:04 am

Hi all :
I have met a strange problem. My first case is based on shoreface(you can find some information
https://www.myroms.org/wiki/index.php/SHOREFACE_CASE). The domain size is 181*8(Lm*Mm),bathymetric slope is 1/50. I generate a forcing file from swan ,which the value in Y direction are all same(I get some help from J.C Warner , thanks very much !).The case is run well , and I get some result.
Then I extended the domain size to 181*599(Lm *Mm), also extended the forcing file, and not change other things. But it blow up at the 1300st time step. My boundary in south and north is NS_PERIODIC. I also try to change it to NORTHERN_WALL SOUTHERN_WALL ,compile ROMS again and run two case. The problem is the same, run well in small scale case , blow up in the large scale case . Can anyone lay some hints on the reason for that? Thank you very much!
--chenzhen

Time-Stepping parameters.
NTIMES == 7200
DT == 1d0
NDTFAST == 30

Vertical parameters
N == 20
Vtransform == 1 Vstretching == 1
THETA_S == 0.0d0
THETA_B == 0.0d0
TCLINE == 0.0d0

cpp.h
#define ROMS_MODEL
#undef WEC_MELLOR
#define WEC_VF
#define WDISS_CHURTHOR
#define BOTTOM_STREAMING
#define SURFACE_STREAMING

#define UV_VIS2
#define MIX_S_UV
#define DIAGNOSTICS_UV
#define AVERAGES
#define AVERAGES_WEC
#define WET_DRY
#define OUT_DOUBLE
#define UV_ADV
#define UV_C2ADVECTION
#undef TS_MPDATA
#define TS_U3HADVECTION
#define DJ_GRADPS
#undef SALINITY
#define SOLVE3D
#define SPLINES
#define NS_PERIODIC
#define EASTERN_WALL
#define WEST_FSCHAPMAN
#define WEST_M2FLATHER
#define WEST_M3GRADIENT


#define MASKING
#ifdef MASKING
# define ANA_MASK
#endif
#undef ANA_GRID
#define ANA_INITIAL
#define ANA_FSOBC
#define ANA_M2OBC
#define ANA_SMFLUX
#define UV_QDRAG

#ifdef SOLVE3D
# undef SSW_BBL
# ifdef SSW_BBL
# define SSW_CALC_ZNOT
# undef SSW_LOGINT
# endif

# undef SEDIMENT
# ifdef SEDIMENT
# undef SED_BIOMASS
# undef SED_MORPH
# define SUSPLOAD
# define BEDLOAD_MPM
# undef BEDLOAD_SOULSBY
# define AVERAGES_BEDLOAD
# 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
# undef ANA_VMIX

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

#endif
bath.png
Hwave.png
ubar.png

Post Reply