Blow up

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
javadi
Posts: 115
Joined: Mon Jan 27, 2020 5:55 pm
Location: UoT

Blow up

#1 Unread post by javadi »

Dear users

I have applied all files including grid, initial, forcing, rivers... to my model. After that, the model began to run. The model ran for 26 days without any problem but after that, I have faced a blow-up. Based on my previous experiences and this forum, I have done some actions:

1) I decrease the time step from 100 to 5 s
2) I have revised my .h file. for example, using LMD_MIXING

In all cases, the model has faced blow-up in 26 days after beginning to run, exactly at the same time.

I have checked my log file. This blowing up refers to a grid that is land(!) as the model output shows the grid out of the domain.


How to solve this problem?

Here is the log file
Attachments
log.txt
(44.2 KiB) Downloaded 149 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blow up

#2 Unread post by kate »

This is when you need to become a detective. Which field is blowing up. You say it is in the land? Is that the only place going bad?

javadi
Posts: 115
Joined: Mon Jan 27, 2020 5:55 pm
Location: UoT

Re: Blow up

#3 Unread post by javadi »

Thanks(!)

I do not understand what you mean ''which field is blowing up''.

Initially, I guessed that this point is related to the masking process as I have edited grids(masking water or land when making a grid.). For this reason, I have checked again this point. This grid is land and I have not changed it.

About the result of the model, the model output is rational when checking temperature or salinity...

To be honest, I cannot answer this only fields the model is unstable but I think that most of the instability can be as a result of land.
I have monitored all time steps when running the model. The model's maximum speed varies from 0.4 to 1.5 in all-time steps.

I do not understand why models run equations for grids that are land(!).

MY .h is here:


#define UV_ADV
#define UV_COR
#define UV_LDRAG
#define UV_VIS2
#undef MIX_GEO_UV
#define MIX_S_UV
#define SPLINES_VDIFF
#define SPLINES_VVISC
#define DJ_GRADPS
#define TS_DIF2
#undef TS_DIF4
#undef MIX_GEO_TS
#define MIX_S_TS

# define LMD_MIXING
# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LI_FOX_KEMPER
# undef LMD_BKPP
# define LMD_NONLOCAL
# define LMD_SHAPIRO
# define LMD_DDMIX
# define LIMIT_VDIFF
# endif
#define SALINITY
#define SOLVE3D
#define AVERAGES
#define DIAGNOSTICS_TS
#define DIAGNOSTICS_UV
#undef ANA_GRID
#undef ANA_INITIAL
#undef ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_BTFLUX
#define ANA_BSFLUX
#undef ANA_PSOURCE
#undef ANA_SPONGE
#if defined BIO_FENNEL || defined ECOSIM || \
defined NPZD_POWELL || defined NEMURO
# define ANA_BIOLOGY
# define ANA_SPFLUX
# define ANA_BPFLUX
# define ANA_SRFLUX
#endif
#if defined NEMURO
# define HOLLING_GRAZING
# undef IVLEV_EXPLICIT
#endif
#ifdef BIO_FENNEL
# define CARBON
# define DENITRIFICATION
# define BIO_SEDIMENT
# define DIAGNOSTICS_BIO
#endif

#ifdef PERFECT_RESTART
# undef AVERAGES
# undef DIAGNOSTICS_BIO
# undef DIAGNOSTICS_TS
# undef DIAGNOSTICS_UV
# define OUT_DOUBLE

Post Reply