search

jagurs-admin / jagurs

JAGURS is a tsunami simulation code solving linear/nonlinear long-wave/Boussinesq equations with/without effects of elastic deformation of the Earth due to tsunami load and vertical profile of seawater density stratification. This code was paralleled by using OpenMP and MPI.
Other
50 stars 18 forks source link

Error upon computing with Boussinesq dispersion #13

Closed amg1199 closed 6 months ago

amg1199 commented 1 year ago

Hello,

We are currently comparing non-dispersive and dispersive computation in JAGURS. While we haven't experienced any issues with non-dispersive computing, when using dispersive or Boussinesq computing, the execution was halted due to this error:

=== ERROR!!!
=== Maximum wave height exceeds 1 million
=== in domain ID   1!
=== Probable instability!
=== Execution is stopped at step    43318 and finalized!

We used GMT for resampling grids until the 2nd daughter grid (3 grids in total) with resolutions: 15 arcsec, 5 arcsec, and 1.66666667 arcsec. I have attached the bathymetries and parameters for your reference.

Hoping for your response! Thank you.

General_Santos_City_Philippines_Bathymetries_Parameters.zip

jagurs-admin commented 1 year ago

Hi

Thank you for using JAGURS and question. JAGURS dispersion calculations often cause numerical instability. Unfortunately, we don't know any fundamental solutions to avoid the instability. When numerical instability occurred during dispersion calculations, we checked the bathymetric data at the point of the instability. The instability is more likely to appear in the bathymetric data containing anomalous values, the seafloor slope changes abruptly, or the water depth is too deep for the grid spacing.

Thank you.

From: 'amg1199' via jagurs_admin @.> Sent: Friday, September 29, 2023 3:52 PM To: jagurs-admin/jagurs @.> Cc: Subscribed @.***> Subject: [jagurs-admin/jagurs] Error upon computing with Boussinesq dispersion (Issue #13)

Hello,

We are currently comparing non-dispersive and dispersive computation in JAGURS. While we haven't experienced any issues with non-dispersive computing, when using dispersive or Boussinesq computing, the execution was halted due to this error:

=== ERROR!!!

=== Maximum wave height exceeds 1 million

=== in domain ID 1!

=== Probable instability!

=== Execution is stopped at step 43318 and finalized!

We used GMT for resampling grids until the 2nd daughter grid (3 grids in total) with resolutions: 15 arcsec, 5 arcsec, and 1.66666667 arcsec. I have attached the bathymetries and parameters for your reference.

Hoping for your response! Thank you.

General_Santos_City_Philippines_Bathymetries_Parameters.ziphttps://github.com/jagurs-admin/jagurs/files/12761708/General_Santos_City_Philippines_Bathymetries_Parameters.zip

— Reply to this email directly, view it on GitHubhttps://github.com/jagurs-admin/jagurs/issues/13, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AEKCJZFD2V3LVLSBM6YCMLLX4ZVZHANCNFSM6AAAAAA5MAO7UI. You are receiving this because you are subscribed to this thread.Message ID: @.**@.>>

To unsubscribe from this group and stop receiving emails from it, send an email to @.**@.>.