swash-users — for questions and discussion related to SWASH

[swash-users] Error during compilation of latest version

[Nicolas G. <nic...@ou...>]

Dear devs,

I tried to compile latest version using main branch from gitlab repo but I get this error:

$ cmake --build .
[3/28] Building Fortran object lib/CMakeFiles/swash11.01.dir/SwashBCspectrum.f90.o
FAILED: lib/CMakeFiles/swash11.01.dir/SwashBCspectrum.f90.o
/usr/bin/f95  -I../src -O -w -fno-second-underscore -ffree-line-length-none -Jmod   -fpreprocessed -c lib/CMakeFiles/swash11.01.dir/SwashBCspectrum.f90-pp.f90 -o lib/CMakeFiles/swash11.01.dir/SwashBCspectrum.f90.o
../src/SwashBCspectrum.f90:337:28:

  337 |        call random_seed(put=tseed)
      |                            1
Error: Size of ‘put’ argument of ‘random_seed’ intrinsic at (1) too small (1/33)
[24/28] Building Fortran object lib/CMakeFiles/swash11.01.dir/SwashQuanOutp.f90.o
ninja: build stopped: subcommand failed.

It seems related to this change: https://gitlab.tudelft.nl/citg/wavemodels/swash/-/commit/c4257fc693139cae42d7d95c94a03e633ac32837#2ce7c825b5b808abffd4ec32380f3ce8b24735be

Kind regards,
Nicolas

[swash-users] SWASH usage

[Martin P. <mar...@as...>]

Hello,

I am using SWASH for the first time trying to model propagating of waves in an open water body with some land masses around. I am using a bottom grid file with depths in each point which also includes land points given with nan. In the end I would like to be able to plot a three-hour developed wave state on a map for this particular area in a .vtk/.vts file. I have been struggling to set up the input file with the different commands, especially selecting the timesteps and in the output I do not get any values for Hsig for example over the whole domain. I only get a proper plot for the sea depth and land points. Could someone look into the input file in the attachment and give some directions as to what should I do further?

Best regards,

Martin

[swash-users] Merging output files

[Andrea B. <and...@un...>]

Hi everybody,
I have an issue related to SWASH run in parallel mode not producing the unique MatLab output files after completing the calculations, leaving me with the processor-related output files.
>From the forum of this mailing list, it is clear to me that a solution to this issue is producing the output as netcdf or vtk files. Nonetheless, the last simulation that I launched failed a first time due to a power down and, secondly, resulted in not producing the output files, as described. The problem is that the simulation is a long one, taking some two weeks to terminate.
I then ask if someone is informed about the possibility of merging the processor-related output files with some post-processing procedure/script or whatever, so that I don't have to wait two more weeks to have my results as netcdf or vtk files.
Thank you very much in advance

Andrea Balzano
________________________________________________________________________________
Prof. Ing. Andrea Balzano, PhD
Università di Cagliari
Facoltà di Ingegneria e Architettura
Dipartimento di Ingegneria Civile, Ingegneria Ambientale e Architettura - DICAAR
Via Marengo, 2
09123 Cagliari
Phone.: +39 070 6755304
Fax:    +39 070 6755304
e-mail:  and...@un...<mailto:and...@un...>

Re: [swash-users] SWASH error "the area of the centroid dual is negative..."

[Clayton H. <cl...@ca...>]

It sounds like you have either a zero area element (three colocated
points), or an element in which the points are oriented in the clockwise
direction (yielding a negative area).  Some griding software like Blue
Kenue have tools to find and correct orientation problems.

Cheers,

Clayton

On Wed, Apr 16, 2025, 9:55 a.m. Ana Catarina Zózimo <acz...@ln...>
wrote:

> Dear all,
> A student of mine is currently running the unstructured grid SWASH model
> and encountered an error message: *"The area of the centroid dual is
> negative or zero at the vertex" *after enabling the BREAKING option. This
> issue causes the program to terminate. However, the model operates
> correctly when BREAKING is not activated.
>
> Could you please assist us with this issue?
>
> Thanks and best regards,
>
> --
> Ana Catarina Zózimo
> Doutorada
> Núcleo de Portos e Estruturas Marítimas/Departamento de Hidráulica e Ambiente
> LNEC - Laboratório Nacional de Engenharia CivilAv. do Brasil 101, 1700-066 Lisboa, Portugal <https://www.google.com/maps/search/Av.+do+Brasil+101,+1700-066+Lisboa,+Portugal?entry=gmail&source=g>
> Tel. (+351) 21 844 3756
> ---
> Ana Catarina Zózimo
> Postdoctorate
> Ports and Maritime Structures Unit/Hydraulics and Environment Department
> LNEC - National Laboratory for Civil EngineeringAv. do Brasil 101, 1700-066 Lisboa, Portugal <https://www.google.com/maps/search/Av.+do+Brasil+101,+1700-066+Lisboa,+Portugal?entry=gmail&source=g>
> Tel. (+351) 21 844 3756
>
> _______________________________________________
> swash-users mailing list
> swa...@li...
> https://lists.sourceforge.net/lists/listinfo/swash-users
>

[swash-users] [PLEASE HELP!] Salinity transport in unstructured grids

[LEONORA A. P. Z. (Alumno) <LEO...@us...>]

Hello everyone,
I'm working with the new version 11.01, which includes constituent transport for unstructured grids.
I started with a very simple simulation. First, I ran it without any boundary conditions—everything worked fine. Then I added a small discharge on the left side of the domain—still fine. But when I added salinity to this same setup (with the discharge), the model crashes. According to the .prt file, the time step keeps halving until the water depths at some nodes become negative, which makes no sense since my model is just a box with a flat bathymetry of 10 meters across the entire domain.
Honestly, this is driving me crazy. Could someone please help?
Best regards,
Leonora Pereira Z.
Estudiante de Ingeniería Civil
leo...@us...
Universidad Técnica Federico Santa María

[swash-users] Request for Support Regarding NaN Results in Floating Object Analysis with SWASH

[Semih C. <sem...@do...>]

Dear all,

 

Recently, I have been working on analyzing floating object behavior using
SWASH. I aim to compare the model results with physical experiments. To this
end, I set up a wave flume domain and successfully simulated both regular
and irregular waves without any issues. However, after introducing a fixed
floating object into the model, the "WLEV" outputs across the entire
computational domain turned into "NaN" values. To resolve this issue, I have
tried the following approaches:

 

.              Tested various combinations of computational and input domain
sizes (10x40m, 20x40m, 10x100m, 20x100m, 40x100m, 100x100m, etc.).

.              Used different spatial resolutions (1x1m, 0.1x1m, 0.01x1m,
0.1x0.1m, 0.01x0.01m, etc.).

.              Conducted analyses with various time steps (0.1s, 0.01s,
0.005s, 0.001s, 0.0001s).

.              Defined regular and irregular waves approaching perpendicular
to the channel with wave heights between 0.1-1m and periods between 3-6s.

.              Configured the bottom file with both three sides closed and
all sides open boundary conditions.

.              Applied porosity to closed boundaries and compared results
with and without porosity.

.              Added sponge layers at boundaries without incoming waves and
observed the effects with and without sponge layers.

.              Tested both INIT WATLEV and INIT ZERO as initial conditions.

.              Defined water level across the entire domain.

.              Used the SET command to define SWL combinations between 0-5m.

.              Adjusted bathymetry between 0.5m-10m to investigate
bathymetry effects.

.              Defined the "draft" file both across the entire domain and
only at the structure location.

.              Ran analyses using BODY DIM and BODY SOLV commands.

 

Throughout these tests, although I occasionally encountered stability
issues, I was able to resolve them with appropriate adjustments. However,
the issue with NaN "WLEV" output persists specifically when using the DRAFT
and FLOAT commands. Upon reviewing the Draft output file, I can confirm that
the draft file is correctly defined and recognized by the model but despite
all efforts, I continue to receive NaN results.

 

I would greatly appreciate any suggestions or insights you might have
regarding this issue.

 

Thank you for your time and support.

 

Best regards,




Semih CAN
Civil/Coastal Engineer, MSc.

 

Dolfen Engineering&Consultancy

Üsküp Caddesi (Çevre Sokak) 20/3,

Kavaklıdere, 06680 Ankara

T: (312) 468 5930 pbx

F: (312) 468 5941

E:  <mailto:sem...@do...> sem...@do...

W: <http://www.dolfen.com/> www.dolfen.com

 

[swash-users] SWASH error "the area of the centroid dual is negative..."

[Ana C. Z. <acz...@ln...>]

Dear all,

A student of mine is currently running the unstructured grid SWASH model 
and encountered an error message: *"The area of the centroid dual is 
negative or zero at the vertex" *after enabling the BREAKING option. 
This issue causes the program to terminate. However, the model operates 
correctly when BREAKING is not activated.

Could you please assist us with this issue?

Thanks and best regards,

-- 
Ana Catarina Zózimo
Doutorada
Núcleo de Portos e Estruturas Marítimas/Departamento de Hidráulica e Ambiente
LNEC - Laboratório Nacional de Engenharia Civil
Av. do Brasil 101, 1700-066 Lisboa, Portugal
Tel. (+351) 21 844 3756
---
Ana Catarina Zózimo
Postdoctorate
Ports and Maritime Structures Unit/Hydraulics and Environment Department
LNEC - National Laboratory for Civil Engineering
Av. do Brasil 101, 1700-066 Lisboa, Portugal
Tel. (+351) 21 844 3756

[swash-users] Submerged breakwater representation

[DEKNEUDT Q. <que...@ar...>]

Hello everyone,
I would like to represent a submerged breakwater in SWASH 1D to estimate its effectiveness and wave attenuation behind it. What do you think is the best way to represent this structure? Is it directly in the bathymetry .bot or via a porous structure with porosity, height, and grain size?

Thank you for your responses.

I have already conducted a comparison test and the results are presented below.
[cid:image003.png@01DBA8BA.CBB19D40]

Best regards,
Quentin Dekneudt

[swash-users] [HELP!] Issue with Freshwater Appearing in Deactivated Cells in River-Ocean Interaction Simulation

[LEONORA A. P. Z. (Alumno) <LEO...@us...>]

Hello everyone,
I am trying to simulate the interaction between a river and the ocean, specifically to observe the behavior of the plume. For now, I only have a discharge (DISCH) as a boundary condition. My initial conditions are, of course, the bathymetry and a salinity file. My bathymetry is quite simple: the ocean has a steeper slope than the beach slope, and the river is channeled by a 3.5-meter wall. In my initial salinity condition, I set a salinity of 35 throughout the entire ocean domain and the river, and for everything that is "land" (negative value), I set a value of -999 because in my .sws file, I specify that I want the cells containing the value -999 to be deactivated (to avoid unnecessary salinity calculations and reduce simulation time). What happens when I do this is that, at 1 second, freshwater starts appearing at the mouth of my river, specifically around the cells I deactivated, which channel the river (gradual expansion). In my opinion, this is incorrect because freshwater should only appear where the river starts, and therefore where the discharge is. I want to clarify that freshwater does appear at the start of the river, but it also appears in my deactivated cells.
I will attach a photo to make it clearer. Does anyone know why this is happening? Is it inappropriate for me to deactivate the cells that represent my wall? Is it better to leave the entire domain with an initial salinity of 35? Is there any way to speed up the simulation without deactivating the cells?
I am also attaching my .sws, .bot, and .sal files.
The simulation did not crash; it runs fine, but I am disturbed by seeing this issue that shouldn’t be there.
Thanks.

[cid:9bb2b137-fa4e-47b5-8c81-3102a4f6d93f] [cid:157ec24b-b2db-4010-b854-c5d0e43b6530] [cid:36880485-ccec-43ed-a795-c6eda7ffa96f]
I would like to clarify that the blue color represents freshwater and the white color represents saltwater or 35

Leonora Pereira Z.
Estudiante de Ingeniería Civil y MSc. en Ingeniería Civil
Universidad Técnica Federico Santa María
Chile

[swash-users] Problems to register in swash users mailing list

[Ana C. Z. <acz...@ln...>]

Dear Sir/Madam,

A colleague of mine, in CC to this e-mail, is having problems in 
registering in the swash users mailing list. He already filled the form 
to register and updated his profile, but he still can't manage to access 
the group.

Can you please help us on that?

Thank you and best regards,

-- 
Ana Catarina Zózimo
Doutorada
Núcleo de Portos e Estruturas Marítimas/Departamento de Hidráulica e Ambiente
LNEC - Laboratório Nacional de Engenharia Civil
Av. do Brasil 101, 1700-066 Lisboa, Portugal
Tel. (+351) 21 844 3756
---
Ana Catarina Zózimo
Postdoctorate
Ports and Maritime Structures Unit/Hydraulics and Environment Department
LNEC - National Laboratory for Civil Engineering
Av. do Brasil 101, 1700-066 Lisboa, Portugal
Tel. (+351) 21 844 3756

Re: [swash-users] [PLEASE HELP!] Issues Running SWASH in Parallel with OpenMPI on Ubuntu 24

[Bunnik, T. <t.b...@ma...>]

Hi Leonora,

Could it be that you accidentally build a serial version of swash instead of the mpi version?

Best regards,

Tim


dr. ir. Tim Bunnik
 | Senior Researcher | Research & Development|Monday, 9 am – 5.30 pm​
Tuesday, 9 am – 5.30 pm​
Wednesday, 9 am – 5.30 pm​
Thursday, 9 am – 5.30 pm​
Friday, 9 am – 5.30 am​
T +31 317 49 34 29 | M +31 654 23 42 81 | t.b...@ma... | www.marin.nl
From: LEONORA ANTONIA PEREIRA ZAPATA (Alumno) <LEO...@us...>
Sent: Tuesday, March 18, 2025 3:02 AM
To: swa...@li...
Subject: [swash-users] [PLEASE HELP!] Issues Running SWASH in Parallel with OpenMPI on Ubuntu 24

Algunos contactos que recibieron este mensaje no suelen recibir correos electrónicos de leo...@us...<mailto:leo...@us...>. Por qué es esto importante<https://aka.ms/LearnAboutSenderIdentification>

Hello everyone,



I have been trying to run SWASH in parallel, so I downloaded OpenMPI on Ubuntu version 24.



I ran a simulation with 15 cores, however, I get the impression that it is running the same simulation 15 times, meaning not in parts. In the terminal, the typical percentage that appears when running a simulation shows up 15 times. I also notice that the simulation takes MUCH longer than if I were running it serially.



Has anyone figured out how to resolve this?

[swash-users] [PLEASE HELP!] Issues Running SWASH in Parallel with OpenMPI on Ubuntu 24

[LEONORA A. P. Z. (Alumno) <LEO...@us...>]

Hello everyone,


I have been trying to run SWASH in parallel, so I downloaded OpenMPI on Ubuntu version 24.


I ran a simulation with 15 cores, however, I get the impression that it is running the same simulation 15 times, meaning not in parts. In the terminal, the typical percentage that appears when running a simulation shows up 15 times. I also notice that the simulation takes MUCH longer than if I were running it serially.


Has anyone figured out how to resolve this?

[swash-users] Simulating storm surge with SWASH - issue with spherical coordinates

[Andrea R. <mat...@gm...>]

Dear community,
I would like to use SWASh to simulate storm surge in coastal areas. The
plan is to cover an area of about 40 km wide taking the boundary conditions
from a larger model. The idea is to start with the case a34wndco provided
in the testcases folder.
Step by step:
I run the case as it is. Works fine
I move from 1D to 2D. Works fine
I use Nautical convention. Works fine.
The problem arises when I move from cartesian to spherical coordinates.
With spherical coordinates the model gives a negligible storm surge at the
eastern boundary (wind is blowing from the West at 30 m/s), that should be
on the order of 0.7/0.8 m.

This is the input with cartesian coordinates:
$*************HEADING***********************************
$
PROJ 'A34wnd01' 'A34'
$
$  Wind setup in semi-closed basin, test 1: correction to wind velocity at
surface
$
$ --|--------------------------------------------------------------|--
$   | This SWASH input file is part of the bench mark tests for    |
$   | SWASH. More information about this test can be found in      |
$   | an accompanied document.                                     |
$ --|--------------------------------------------------------------|--
$
$***********MODEL INPUT*********************************
$
SET temp = 4. salinity = 0.1 rhoair = 1.225 NAUTICAL
$
MODE NONST TWOD
$
CGRID 0. 0. 0. 10000. 200. 100 2
$
VERT 1
$
INPGRID BOTTOM 0. 0. 0. 1 1 10000. 200.
READINP BOTTOM 5. 'a34wind2D.bot' 3 0 FREE
$
INIT zero
$
BOU SIDE W CCW BTYPE WLEV CON 0.
$
WIND 30. 270. CONSTANT cd = 0.0043 RELATIVE alpha = 1.
FRIC CHEZY 60.
VISC 0.
$
DISCRET UPW FIRST
DISCRET CORRDEP FIRST
$
$************ OUTPUT REQUESTS *************************
$
POINTS 'GAUGE' 10000. 34.
TABLE  'GAUGE' NOHEAD 'a34wnd012D.tbl' TSEC WATL OUTPUT 000000.000 1 MIN
$
TEST 1,0
COMPUTE 000000.000 10. SEC 080000.000
STOP



This is the input with spherical coordinates:
$*************HEADING***********************************
$
PROJ 'A34wnd01' 'A34'
$
$  Wind setup in semi-closed basin, test 1: correction to wind velocity at
surface
$
$ --|--------------------------------------------------------------|--
$   | This SWASH input file is part of the bench mark tests for    |
$   | SWASH. More information about this test can be found in      |
$   | an accompanied document.                                     |
$ --|--------------------------------------------------------------|--
$
$***********MODEL INPUT*********************************
$
SET temp = 4. salinity = 0.1 rhoair = 1.225 CORIOLIS NAUTICAL
$
MODE NONST TWOD
COORDINATES SPHE
$
CGRID 0. 40. 0. 0.1 0.002 100 2
$
VERT 1
$
INPGRID BOTTOM 0. 40. 0. 1 1 0.1 0.002
READINP BOTTOM 5. 'a34wind2D.bot' 3 0 FREE
$
INIT zero
$
BOU SIDE W CCW BTYPE WLEV CON 0.
$
WIND 30. 270. CONSTANT cd = 0.0043 RELATIVE alpha = 1.
FRIC CHEZY 60.
VISC 0.
$
DISCRET UPW FIRST
DISCRET CORRDEP FIRST
$
$************ OUTPUT REQUESTS *************************
$
POINTS 'GAUGE' 0.1 40.001
TABLE  'GAUGE' NOHEAD 'a34wnd012Dsphe.tbl' TSEC WATL OUTPUT 000000.000 1 MIN
$
TEST 1,0
COMPUTE 000000.000 10. SEC 080000.000
STOP


The file with the bathymetry is in both cases:
 1.0 1.0
 1.0 1.0


Does anybody have a clue about what is going wrong?
Cheers,
Andrea

Re: [swash-users] Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

[<mer...@do...>]

Leonora,

I have checked the sws file and there may be a problem in your output
commands. In the "Compute" part you have;

'COMPUTE 000000.0 0.09 SEC 043200.0'

This means that start the computation at time 0 and compute with intervals
of 0.09 seconds, end it at time 4 hours 32 minutes 00 seconds (If you do not
change the default lock-up settings) . You can check the notation for this
part. 

 





 

 

Saygılarımla,

Best regards,

Mert Yaman
İnşaat Yüksek Mühendisi
Civil/Coastal Engineer, MSc.

 

 

 

Dolfen Danışmanlık Mühendislik A.Ş.

Üsküp Caddesi (Çevre Sokak) 20/3,

Kavaklıdere, 06680 Ankara

T: (312) 468 5930 pbx

F: (312) 468 5941

E:  <mailto:mer...@do...> mer...@do...

W:  <http://www.dolfen.com/> www.dolfen.com

 

From: LEONORA ANTONIA PEREIRA ZAPATA (Alumno) <LEO...@us...> 
Sent: Monday, March 10, 2025 6:03 AM
To: swa...@li...
Subject: [swash-users] Issues with Simulating Salinity Plume in SWASH: Model
Crashes When Approaching the Coast

 

Hello SWASH community,

I am trying to simulate the interaction between the ocean and a river to
analyze the salinity plume. Currently, I am working only with a river
discharge flowing into the sea (no waves or tides involved). However, when
the plume starts approaching the coast, the model crashes. Honestly, I don't
fully understand why this is happening.

My domain is 1200x999 meters with cells of 3x3 meters each. The simulation
time is 12 hours, but it crashes at 6 hours. According to the SWASH manual,
when dealing with salinity and sediment analysis, it is suggested to use TVD
as the numerical scheme. I applied it, but the model still crashes.

I will attach the .sws, .bot, and .sal files in case anyone can help me.

Could it be the bathymetry? In my opinion, it shouldn't cause any issues.

The .sws file contains the following:

$*************HEADING***********************************

$

PROJ 'tester' '01'

$

$***********MODEL INPUT*********************************

$

MODE NONSTationary TWOD

COORDinates CARTesian

$

CGRID 0. 0. 0. 999 1200 333 400

VERT 2

$

INPGRID BOTTOM  0. 0. 0. 333 400 3.0 3.0

READINP BOTTOM 1. 'bathymetry.bot' 1 0 FREE

$

INITial ZERO

$

INPTRANS SALINITY 0. 0. 0. 333 400 3.0 3.0 EXCEPTION -999

READTRANS SALINITY 1. 'salinity.sal' 1 0 FREE

$

BOU SEGM XY 999 321 999 303 BTYPE DISCH SMOO 10 MIN CON FOUR -1

BOU SIDE W CCW BTYPE RADIATION

BOU SIDE N CCW BTYPE RADIATION

BOU SIDE S CCW BTYPE RADIATION

$

VISC H CON 0.001

FRIC MANN 0.019

NONHYDROSTATIC STAndard 1

$

DISCRET UPW MOM

DISCRET UPW WMOM

DISCRET CORRDEP

DISCRET MINMOD

$

TIMEI METH EXPL 0.3 0.6

$

$************ OUTPUT REQUESTS *************************

$

QUANTITY VELK

QUANTITY SALK

$

GROUP 'OUTPUT' SUBGRID 1 334 1 401

BLOCK 'OUTPUT' NOHEAD 'salinidad_limiter.mat' LAY-OUT 4 SALK OUTPUT 000000.0
60 SEC

BLOCK 'OUTPUT' NOHEAD 'velocidad_limiter.mat' LAY-OUT 4 VELK OUTPUT 000000.0
60 SEC

$

TEST 1,0

COMPUTE 000000.0 0.09 SEC 043200.0

STOP

I hope someone can help me.

Best regards,
Leonora Pereira
+56976494317
leo...@us... <mailto:leo...@us...> 
Civil Engineering Student, UTFSM

 

 

Re: [swash-users] Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

[<mer...@do...>]

Hello Leonora,

 

If you can share .prt file too, I think it can be helpfull. Also, if swash
generates an error file at the end, it may contain the reason of the crash.
Only thing comes to my mind is that you can try to increase the number of
vertical layers. There is a section discusses vertical layers in the SWASH
manual, you can find some information there too. Hope you can solve it in no
time. 

 

Saygılarımla,

Best regards,

Mert Yaman
İnşaat Yüksek Mühendisi
Civil/Coastal Engineer, MSc.

 

 

 

Dolfen Danışmanlık Mühendislik A.Ş.

Üsküp Caddesi (Çevre Sokak) 20/3,

Kavaklıdere, 06680 Ankara

T: (312) 468 5930 pbx

F: (312) 468 5941

E:  <mailto:mer...@do...> mer...@do...

W:  <http://www.dolfen.com/> www.dolfen.com

 

From: LEONORA ANTONIA PEREIRA ZAPATA (Alumno) <LEO...@us...> 
Sent: Monday, March 10, 2025 6:03 AM
To: swa...@li...
Subject: [swash-users] Issues with Simulating Salinity Plume in SWASH: Model
Crashes When Approaching the Coast

 

Hello SWASH community,

I am trying to simulate the interaction between the ocean and a river to
analyze the salinity plume. Currently, I am working only with a river
discharge flowing into the sea (no waves or tides involved). However, when
the plume starts approaching the coast, the model crashes. Honestly, I don't
fully understand why this is happening.

My domain is 1200x999 meters with cells of 3x3 meters each. The simulation
time is 12 hours, but it crashes at 6 hours. According to the SWASH manual,
when dealing with salinity and sediment analysis, it is suggested to use TVD
as the numerical scheme. I applied it, but the model still crashes.

I will attach the .sws, .bot, and .sal files in case anyone can help me.

Could it be the bathymetry? In my opinion, it shouldn't cause any issues.

The .sws file contains the following:

$*************HEADING***********************************

$

PROJ 'tester' '01'

$

$***********MODEL INPUT*********************************

$

MODE NONSTationary TWOD

COORDinates CARTesian

$

CGRID 0. 0. 0. 999 1200 333 400

VERT 2

$

INPGRID BOTTOM  0. 0. 0. 333 400 3.0 3.0

READINP BOTTOM 1. 'bathymetry.bot' 1 0 FREE

$

INITial ZERO

$

INPTRANS SALINITY 0. 0. 0. 333 400 3.0 3.0 EXCEPTION -999

READTRANS SALINITY 1. 'salinity.sal' 1 0 FREE

$

BOU SEGM XY 999 321 999 303 BTYPE DISCH SMOO 10 MIN CON FOUR -1

BOU SIDE W CCW BTYPE RADIATION

BOU SIDE N CCW BTYPE RADIATION

BOU SIDE S CCW BTYPE RADIATION

$

VISC H CON 0.001

FRIC MANN 0.019

NONHYDROSTATIC STAndard 1

$

DISCRET UPW MOM

DISCRET UPW WMOM

DISCRET CORRDEP

DISCRET MINMOD

$

TIMEI METH EXPL 0.3 0.6

$

$************ OUTPUT REQUESTS *************************

$

QUANTITY VELK

QUANTITY SALK

$

GROUP 'OUTPUT' SUBGRID 1 334 1 401

BLOCK 'OUTPUT' NOHEAD 'salinidad_limiter.mat' LAY-OUT 4 SALK OUTPUT 000000.0
60 SEC

BLOCK 'OUTPUT' NOHEAD 'velocidad_limiter.mat' LAY-OUT 4 VELK OUTPUT 000000.0
60 SEC

$

TEST 1,0

COMPUTE 000000.0 0.09 SEC 043200.0

STOP

I hope someone can help me.

Best regards,
Leonora Pereira
+56976494317
leo...@us... <mailto:leo...@us...> 
Civil Engineering Student, UTFSM

 

 

[swash-users] Issues with Simulating Salinity Plume in SWASH: Model Crashes When Approaching the Coast

[LEONORA A. P. Z. (Alumno) <LEO...@us...>]

Hello SWASH community,
I am trying to simulate the interaction between the ocean and a river to analyze the salinity plume. Currently, I am working only with a river discharge flowing into the sea (no waves or tides involved). However, when the plume starts approaching the coast, the model crashes. Honestly, I don't fully understand why this is happening.
My domain is 1200x999 meters with cells of 3x3 meters each. The simulation time is 12 hours, but it crashes at 6 hours. According to the SWASH manual, when dealing with salinity and sediment analysis, it is suggested to use TVD as the numerical scheme. I applied it, but the model still crashes.
I will attach the .sws, .bot, and .sal files in case anyone can help me.
Could it be the bathymetry? In my opinion, it shouldn't cause any issues.
The .sws file contains the following:

$*************HEADING***********************************
$
PROJ 'tester' '01'
$
$***********MODEL INPUT*********************************
$
MODE NONSTationary TWOD
COORDinates CARTesian
$
CGRID 0. 0. 0. 999 1200 333 400
VERT 2
$
INPGRID BOTTOM  0. 0. 0. 333 400 3.0 3.0
READINP BOTTOM 1. 'bathymetry.bot' 1 0 FREE
$
INITial ZERO
$
INPTRANS SALINITY 0. 0. 0. 333 400 3.0 3.0 EXCEPTION -999
READTRANS SALINITY 1. 'salinity.sal' 1 0 FREE
$
BOU SEGM XY 999 321 999 303 BTYPE DISCH SMOO 10 MIN CON FOUR -1
BOU SIDE W CCW BTYPE RADIATION
BOU SIDE N CCW BTYPE RADIATION
BOU SIDE S CCW BTYPE RADIATION
$
VISC H CON 0.001
FRIC MANN 0.019
NONHYDROSTATIC STAndard 1
$
DISCRET UPW MOM
DISCRET UPW WMOM
DISCRET CORRDEP
DISCRET MINMOD
$
TIMEI METH EXPL 0.3 0.6
$
$************ OUTPUT REQUESTS *************************
$
QUANTITY VELK
QUANTITY SALK
$
GROUP 'OUTPUT' SUBGRID 1 334 1 401
BLOCK 'OUTPUT' NOHEAD 'salinidad_limiter.mat' LAY-OUT 4 SALK OUTPUT 000000.0 60 SEC
BLOCK 'OUTPUT' NOHEAD 'velocidad_limiter.mat' LAY-OUT 4 VELK OUTPUT 000000.0 60 SEC
$
TEST 1,0
COMPUTE 000000.0 0.09 SEC 043200.0
STOP

I hope someone can help me.
Best regards,
Leonora Pereira
+56976494317
leo...@us...
Civil Engineering Student, UTFSM

[swash-users] SWAN MPI Setup

[<mer...@do...>]

Dear Swan Community, 

 

Is there any step by step guide to install SWAN with MPI support? Could you
please send a link for it if it exists? Thank you in advance. 

 

Saygılarımla,

Best regards,

Mert Yaman
İnşaat Yüksek Mühendisi
Civil/Coastal Engineer, MSc.

 

 

 

Dolfen Danışmanlık Mühendislik A.Ş.

Üsküp Caddesi (Çevre Sokak) 20/3,

Kavaklıdere, 06680 Ankara

T: (312) 468 5930 pbx

F: (312) 468 5941

E:  <mailto:mer...@do...> mer...@do...

W:  <http://www.dolfen.com/> www.dolfen.com

 

[swash-users] Wave Impact Forces

[Laura S. <l.s...@qu...>]

Hello!

I see that SWASH can output wave impact forces on floating bodies (e.g., moored ships).

Is there a way to output wave impact forces on a static body (e.g., wall)? Or would this have to be calculated from other standard outputs?

Thank you,

Laura Szczyrba, Ph.D. (she/her)

Earth Scientist

BGC Engineering



M: 703-853-2616   E: lsz...@bg...<mailto:lsz...@bg...>

Suite 408 – 27 Princess Street, Kingston, ON K7L 1A3

www.bgcengineering.ca<https://www.bgcengineering.ca/>

[swash-users] Vegetation vertical distribution

[Ignacio <bar...@gm...>]

Dear SWASH community,

I am setting up some simulations with a number of vertical layers. In
the simulations I will have vegetation in a given area. The vegetation
properties vary vertically (I am setting up two vertical layers, and I
would like to match them to the vertical variation of the vegetation).

At the moment I am using the command VEGEtation like this to create
two vertical layers for the vegetation:
VEGE 1. 0.4 4 1.2 0.5 0.4 5 1.2

My question is, what is the order and reference point for the
vegetation heights? Is it relative to the bottom like the Porosity? Is
their order from the top to the bottom like the Vertical command? If I
understand correctly, given my VEGE command above, I will have a 0.5 m
height layer at the bottom with 5 stems/m2 and a 1.0 m height layer on
top of it with 4 stems/m2 (total 1.5m), right?

It seems to me a bit counterintuitive, I would start writing from the
top and then each consecutive layer on top of it, but based on how
SWASH works with similar commands, I would say it works like this for
consistency. I did not find an answer in the manual.

I will appreciate any insights on this topic! Thanks in advance and
Happy New Year!

Ignacio

Re: [swash-users] SWASH curvilinear grid and bathymetry reading.

[Hafizul I. <ton...@ya...>]

 Hello Zhiling,
Thank you very much for the response. Unfortunately, in my case, even after changing the grid format, I keep on receiving the same error message. Can you please mention the SWASH version you are using. 
And please share any other thoughts or suggestions that you think might help resolve the issue.

Thanks,Hafizul  


    On Friday, December 6, 2024 at 03:30:44 PM GMT+1, <tx....@gm...> wrote:  
 
 <!--#yiv7341332093 filtered {}#yiv7341332093 filtered {}#yiv7341332093 filtered {}#yiv7341332093 filtered {}#yiv7341332093 filtered {}#yiv7341332093 filtered {}#yiv7341332093 p.yiv7341332093MsoNormal, #yiv7341332093 li.yiv7341332093MsoNormal, #yiv7341332093 div.yiv7341332093MsoNormal {margin:0cm;font-size:12.0pt;font-family:"Aptos", sans-serif;}#yiv7341332093 span.yiv7341332093EmailStyle18 {font-family:"Aptos", sans-serif;color:windowtext;}#yiv7341332093 .yiv7341332093MsoChpDefault {font-size:10.0pt;}#yiv7341332093 filtered {}#yiv7341332093 div.yiv7341332093WordSection1 {}-->
Hi Hafizul,

  

I can successfully run the test case after changing IDLA to 5 for ‘READ COOR’ line

And change the xy_grid.grd to sth like this

  



Basically two 10 by 4 blocks: upper one for x coord and lower for y

  

Cheers,

Zhiling

  

  

From: Hafizul Islam via swash-users <swa...@li...> 
Sent: 06 December 2024 13:48
To: swa...@li...
Subject: [swash-users] SWASH curvilinear grid and bathymetry reading.

  

Hello,

  

I am a new user of SWASH and have been struggling quite a bit to understand SWASH input file formats (for grid and bathymetry). I have been trying to provide a curvilinear grid for xy dimension and a bathymetry for z in SWASH using a .grd and a .dep file. I keep getting errors while SWASH tries to read the bathymetry file. To understand the problem better, I tried to create a simple 10x4 grid and a sloped bathymetry. I get the following error. Can someone please help me by explaining the issue?   

  

My input files are in the attachment. 

  

Thanks in advance,

Hafizul

  

  

CGRID CURVI 9 3 EXCEPTION -999.000 -999.000

    

 READ  COOR  1. 'xy_grid.grd' IDLA=6 0  FREE

    

 $

    

 VERT 3

    

 $$

    

 INPGRID BOTTOM CURV STAGgered

    

 READ BOT 1. 'bathymetry.dep' 0  FREE

 ** Error            : 3 comp. grid points on one line

   0  0         0.0000         5.0000   1  0         0.0000         5.0000   1  1         2.0000         2.0000

 ** Error            : grid angle <0 or >180 degrees

   1  0         0.0000         5.0000   1  1         2.0000         2.0000   0  1         0.0000         0.0000

 ** Error            : grid angle <0 or >180 degrees

   1  1         2.0000         2.0000   0  1         0.0000         0.0000   0  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  1         0.0000         0.0000   0  0         0.0000         5.0000   1  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  1         2.0000         2.0000   1  2         0.0000         5.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  2         0.0000         5.0000   0  2         0.0000         5.0000   0  1         0.0000         0.0000

 ** Error            : 3 comp. grid points on one line

   0  2         0.0000         5.0000   1  2         0.0000         5.0000   1  3         3.0000         3.0000

 ** Error            : grid angle <0 or >180 degrees

   1  2         0.0000         5.0000   1  3         3.0000         3.0000   0  3         1.0000         1.0000

 ** Error            : grid angle <0 or >180 degrees

   1  3         3.0000         3.0000   0  3         1.0000         1.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  3         1.0000         1.0000   0  2         0.0000         5.0000   1  2         0.0000         5.0000

  
  

Re: [swash-users] SWASH curvilinear grid and bathymetry reading.

[<tx....@gm...>]

Hi Hafizul,

 

I used version 9.01.  Unfortunately, I cannot think of any measures other than the two changes I mentioned.

 

Good luck : )

 

Best,

Zhiling

 

 

From: Hafizul Islam <ton...@ya...> 
Sent: 06 December 2024 15:18
To: swa...@li...; tx....@gm...
Subject: Re: [swash-users] SWASH curvilinear grid and bathymetry reading.

 

Hello Zhiling,

 

Thank you very much for the response. Unfortunately, in my case, even after changing the grid format, I keep on receiving the same error message. Can you please mention the SWASH version you are using. 

 

And please share any other thoughts or suggestions that you think might help resolve the issue.

 

 

Thanks,

Hafizul  

 

 

 

On Friday, December 6, 2024 at 03:30:44 PM GMT+1, <tx....@gm... <mailto:tx....@gm...> > wrote: 

 

 

Hi Hafizul,

 

I can successfully run the test case after changing IDLA to 5 for ‘READ COOR’ line

And change the xy_grid.grd to sth like this

 



Basically two 10 by 4 blocks: upper one for x coord and lower for y

 

Cheers,

Zhiling

 

 

From: Hafizul Islam via swash-users <swa...@li... <mailto:swa...@li...> > 
Sent: 06 December 2024 13:48
To: swa...@li... <mailto:swa...@li...> 
Subject: [swash-users] SWASH curvilinear grid and bathymetry reading.

 

Hello,

 

I am a new user of SWASH and have been struggling quite a bit to understand SWASH input file formats (for grid and bathymetry). I have been trying to provide a curvilinear grid for xy dimension and a bathymetry for z in SWASH using a .grd and a .dep file. I keep getting errors while SWASH tries to read the bathymetry file. To understand the problem better, I tried to create a simple 10x4 grid and a sloped bathymetry. I get the following error. Can someone please help me by explaining the issue?   

 

My input files are in the attachment. 

 

Thanks in advance,

Hafizul

 

 

CGRID CURVI 9 3 EXCEPTION -999.000 -999.000

    

 READ  COOR  1. 'xy_grid.grd' IDLA=6 0  FREE

    

 $

    

 VERT 3

    

 $$

    

 INPGRID BOTTOM CURV STAGgered

    

 READ BOT 1. 'bathymetry.dep' 0  FREE

 ** Error            : 3 comp. grid points on one line

   0  0         0.0000         5.0000   1  0         0.0000         5.0000   1  1         2.0000         2.0000

 ** Error            : grid angle <0 or >180 degrees

   1  0         0.0000         5.0000   1  1         2.0000         2.0000   0  1         0.0000         0.0000

 ** Error            : grid angle <0 or >180 degrees

   1  1         2.0000         2.0000   0  1         0.0000         0.0000   0  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  1         0.0000         0.0000   0  0         0.0000         5.0000   1  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  1         2.0000         2.0000   1  2         0.0000         5.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  2         0.0000         5.0000   0  2         0.0000         5.0000   0  1         0.0000         0.0000

 ** Error            : 3 comp. grid points on one line

   0  2         0.0000         5.0000   1  2         0.0000         5.0000   1  3         3.0000         3.0000

 ** Error            : grid angle <0 or >180 degrees

   1  2         0.0000         5.0000   1  3         3.0000         3.0000   0  3         1.0000         1.0000

 ** Error            : grid angle <0 or >180 degrees

   1  3         3.0000         3.0000   0  3         1.0000         1.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  3         1.0000         1.0000   0  2         0.0000         5.0000   1  2         0.0000         5.0000

 

Re: [swash-users] SWASH curvilinear grid and bathymetry reading.

[<tx....@gm...>]

Hi Hafizul,

 

I can successfully run the test case after changing IDLA to 5 for ‘READ COOR’ line

And change the xy_grid.grd to sth like this

 



Basically two 10 by 4 blocks: upper one for x coord and lower for y

 

Cheers,

Zhiling

 

 

From: Hafizul Islam via swash-users <swa...@li...> 
Sent: 06 December 2024 13:48
To: swa...@li...
Subject: [swash-users] SWASH curvilinear grid and bathymetry reading.

 

Hello,

 

I am a new user of SWASH and have been struggling quite a bit to understand SWASH input file formats (for grid and bathymetry). I have been trying to provide a curvilinear grid for xy dimension and a bathymetry for z in SWASH using a .grd and a .dep file. I keep getting errors while SWASH tries to read the bathymetry file. To understand the problem better, I tried to create a simple 10x4 grid and a sloped bathymetry. I get the following error. Can someone please help me by explaining the issue?   

 

My input files are in the attachment. 

 

Thanks in advance,

Hafizul

 

 

CGRID CURVI 9 3 EXCEPTION -999.000 -999.000

    

 READ  COOR  1. 'xy_grid.grd' IDLA=6 0  FREE

    

 $

    

 VERT 3

    

 $$

    

 INPGRID BOTTOM CURV STAGgered

    

 READ BOT 1. 'bathymetry.dep' 0  FREE

 ** Error            : 3 comp. grid points on one line

   0  0         0.0000         5.0000   1  0         0.0000         5.0000   1  1         2.0000         2.0000

 ** Error            : grid angle <0 or >180 degrees

   1  0         0.0000         5.0000   1  1         2.0000         2.0000   0  1         0.0000         0.0000

 ** Error            : grid angle <0 or >180 degrees

   1  1         2.0000         2.0000   0  1         0.0000         0.0000   0  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  1         0.0000         0.0000   0  0         0.0000         5.0000   1  0         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  1         2.0000         2.0000   1  2         0.0000         5.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   1  2         0.0000         5.0000   0  2         0.0000         5.0000   0  1         0.0000         0.0000

 ** Error            : 3 comp. grid points on one line

   0  2         0.0000         5.0000   1  2         0.0000         5.0000   1  3         3.0000         3.0000

 ** Error            : grid angle <0 or >180 degrees

   1  2         0.0000         5.0000   1  3         3.0000         3.0000   0  3         1.0000         1.0000

 ** Error            : grid angle <0 or >180 degrees

   1  3         3.0000         3.0000   0  3         1.0000         1.0000   0  2         0.0000         5.0000

 ** Error            : 3 comp. grid points on one line

   0  3         1.0000         1.0000   0  2         0.0000         5.0000   1  2         0.0000         5.0000

 

[swash-users] SWASH curvilinear grid and bathymetry reading.

[Hafizul I. <ton...@ya...>]

Hello,
I am a new user of SWASH and have been struggling quite a bit to understand SWASH input file formats (for grid and bathymetry). I have been trying to provide a curvilinear grid for xy dimension and a bathymetry for z in SWASH using a .grd and a .dep file. I keep getting errors while SWASH tries to read the bathymetry file. To understand the problem better, I tried to create a simple 10x4 grid and a sloped bathymetry. I get the following error. Can someone please help me by explaining the issue?   
My input files are in the attachment. 
Thanks in advance,Hafizul

CGRID CURVI 9 3 EXCEPTION -999.000 -999.000     READ  COOR  1. 'xy_grid.grd' IDLA=6 0  FREE     $     VERT 3     $$     INPGRID BOTTOM CURV STAGgered     READ BOT 1. 'bathymetry.dep' 0  FREE ** Error            : 3 comp. grid points on one line   0  0         0.0000         5.0000   1  0         0.0000         5.0000   1  1         2.0000         2.0000 ** Error            : grid angle <0 or >180 degrees   1  0         0.0000         5.0000   1  1         2.0000         2.0000   0  1         0.0000         0.0000 ** Error            : grid angle <0 or >180 degrees   1  1         2.0000         2.0000   0  1         0.0000         0.0000   0  0         0.0000         5.0000 ** Error            : 3 comp. grid points on one line   0  1         0.0000         0.0000   0  0         0.0000         5.0000   1  0         0.0000         5.0000 ** Error            : 3 comp. grid points on one line   1  1         2.0000         2.0000   1  2         0.0000         5.0000   0  2         0.0000         5.0000 ** Error            : 3 comp. grid points on one line   1  2         0.0000         5.0000   0  2         0.0000         5.0000   0  1         0.0000         0.0000 ** Error            : 3 comp. grid points on one line   0  2         0.0000         5.0000   1  2         0.0000         5.0000   1  3         3.0000         3.0000 ** Error            : grid angle <0 or >180 degrees   1  2         0.0000         5.0000   1  3         3.0000         3.0000   0  3         1.0000         1.0000 ** Error            : grid angle <0 or >180 degrees   1  3         3.0000         3.0000   0  3         1.0000         1.0000   0  2         0.0000         5.0000 ** Error            : 3 comp. grid points on one line   0  3         1.0000         1.0000   0  2         0.0000         5.0000   1  2         0.0000         5.0000

Re: [swash-users] Simulate wave interactions with floating body

[Shai A. <sha...@cv...>]

Hi,

You can't use both 'float' and 'body dim'. One's for 0dof and the other's for 6dof.
Also, you should increase the number of iterations in both the body and time integration solvers.
Cheers,

Get Outlook for Android<https://aka.ms/AAb9ysg>
________________________________
From: theofilos koronakis <tko...@gm...>
Sent: Friday, November 1, 2024 7:24:19 PM
To: swa...@li... <swa...@li...>
Subject: [swash-users] Simulate wave interactions with floating body

Hello everyone,

I have a question. In order to simulate wave interaction with a floating body , I have to input a draft file along with commands FLOAT , BODY DIM and BODY SOLVER ?

If yes , I already tried it and it doesn't seem to be working. The scenario im using contains a harbour area as well
Here's the code i'm using :

PROJ 'test' 'SM'
MODE DYN TWOD
COORD CART
set depmin = 0.001
CGRID REG 0. 0. 0. 1000.0 1000.0 250 200 rep y
VERT 1
INPGRID BOTTOM 0. 0. 0. 250 200 4 5
READINP BOTTOM -1.0 'depth.bot' 1 0 FREE
INPGRID POROSITY 0. 0. 0. 250 200 4 5
READINP POROSITY 1.0 'porosity.n' 1 0 FREE
INPGRID DRAFT 0. 0. 0. 250 200 4 5
READINP DRAFT 1.0 'draft.n' 1 0 FREE
SPONGE EAST 220
FLOAT alpha=0.5 theta=0.5
BODY DIMENSION l=1 mass=285.42 Ix=1000 Iy=1000 Iz=1000 cogx=370 cogy=420 cogz=0
BODY SOLVER CH COUPLING KBC
POROSITY 0.1 height=2 wper=8
INIT ZERO
BOU SIDE W CCW BTYPE WEAK CON REG 2 8 0
BRE alpha=0.6 beta=0.3
NONHYD box PREC ILUD PROJECTION
DISCRET UPW UMOM MOM H
TIMEI METH IMPLICIT
QUANTITY HSIG HRMS dur 1 MIN
GROUP 'test' 1 251 1 101
TABLE 'test' HEAD 'test_SM' TIME XP YP HRMS HSIG WATL PRESS NHPRES
BLOCK 'COMPGRID' NOHEAD 'output.mat' LAY 1 XP YP HRMS HSIG WATL PRESS NHPRES
BLOCK 'COMPGRID' NOHEAD 'BOTL.mat' XP YP BOTL
TEST 1,0
COMPUTE 000000.000 0.02 SEC 000500.00 STOP

Best regards ,
Theofilos

External e-mail, be judicious when opening attachments or links

[swash-users] Simulate wave interactions with floating body

[theofilos k. <tko...@gm...>]

Hello everyone,

I have a question. In order to simulate wave interaction with a floating
body , I have to input a draft file along with commands FLOAT , BODY DIM
and BODY SOLVER ?

If yes , I already tried it and it doesn't seem to be working. The scenario
im using contains a harbour area as well
Here's the code i'm using :

PROJ 'test' 'SM'
MODE DYN TWOD
COORD CART
set depmin = 0.001
CGRID REG 0. 0. 0. 1000.0 1000.0 250 200 rep y
VERT 1
INPGRID BOTTOM 0. 0. 0. 250 200 4 5
READINP BOTTOM -1.0 'depth.bot' 1 0 FREE
INPGRID POROSITY 0. 0. 0. 250 200 4 5
READINP POROSITY 1.0 'porosity.n' 1 0 FREE
INPGRID DRAFT 0. 0. 0. 250 200 4 5
READINP DRAFT 1.0 'draft.n' 1 0 FREE
SPONGE EAST 220
FLOAT alpha=0.5 theta=0.5
BODY DIMENSION l=1 mass=285.42 Ix=1000 Iy=1000 Iz=1000 cogx=370 cogy=420
cogz=0
BODY SOLVER CH COUPLING KBC
POROSITY 0.1 height=2 wper=8
INIT ZERO
BOU SIDE W CCW BTYPE WEAK CON REG 2 8 0
BRE alpha=0.6 beta=0.3
NONHYD box PREC ILUD PROJECTION
DISCRET UPW UMOM MOM H
TIMEI METH IMPLICIT
QUANTITY HSIG HRMS dur 1 MIN
GROUP 'test' 1 251 1 101
TABLE 'test' HEAD 'test_SM' TIME XP YP HRMS HSIG WATL PRESS NHPRES
BLOCK 'COMPGRID' NOHEAD 'output.mat' LAY 1 XP YP HRMS HSIG WATL PRESS NHPRES
BLOCK 'COMPGRID' NOHEAD 'BOTL.mat' XP YP BOTL
TEST 1,0
COMPUTE 000000.000 0.02 SEC 000500.00 STOP

Best regards ,
Theofilos