intro story Coast / Estuary

Coast / Estuary

Coastal systems are among the most dynamic physical systems on earth and are subject to a large variety of forces. The morphodynamic changes occurring to coastlines worldwide are of great interest and importance. These changes occur as a result of the erosion of sediments, its subsequent transport as bed load or suspended load, and eventual deposition. 
 
Estuaries are partly enclosed water bodies that have an open connection to the coast. Estuaries generally have one or more branching channels, intertidal mudflats and/or salt marshes. Intertidal areas are of high ecological importance and trap sediments (sands, silts, clays and organic matter).
Within the Delft3D modelling package a large variation of coastal and estuarine physical and chemical processes can be simulated. These include waves, tidal propagation, wind- or wave-induced water level setup, flow induced by salinity or temperature gradients, sand and mud transport, water quality and changing bathymetry (morphology). Delft3D can also be used operationally e.g. storm, surge and algal bloom forecasting. 
 
On this discussion page you can post questions, research discussions or just share your experience about modelling coastal and/or estuarine systems with Delft3D FM. 
 

** PLEASE TAG YOUR POST! **

 

 

Sub groups
D-Flow Flexible Mesh
DELWAQ
Cohesive sediments & muddy systems

 

 

Message Boards

Trouble Opening Delft3D-FLOW trih file in Matlab

James Knighton, modified 5 Years ago.

Trouble Opening Delft3D-FLOW trih file in Matlab

Youngling Posts: 12 Join Date: 8/9/13 Recent Posts
Hi,

I am having difficulty opening the Delft3D-FLOW trih output files using the Matlab functions provided with Delft. I can open the files in Quickplot and can verify that the data is present within the file. The Delft manual only really gives examples for trim files. Based on the brief discussion in the Delft Matlab manual, the following is the best I can come up with for the correct syntax:

infile = 'E:\TsunamiResearch\Delft_PacificOcean\GridFiles\MonteCarloDirectory\trih-Flow_Grid_1.dat';
trih_file = qpfopen(infile);
trih_read = qpread(trih_file);
Stations = qpread(trih_file, trih_read(3), 'stations');
AA_WSEL = qpread(trih_file, trih_read(3), 'data', 0, Stations(5));

The following error is returned:
"Caught in qp_getdata/data: Index exceeds matrix dimensions"


Stations(5) = 'AA', if I substitute this string directly into the function as follows:

AA_WSEL = qpread(trih_file, trih_read(3), 'data', 0, 'AA');

The following error is returned:
"Caught in qp_getdata/data: Error using d3d_trihfil (line 355) Element index 1 of 'ZWL' should contain integers in the range of 1:5."

It's not clear to me what the proper syntax should be, or why the two inputs to the function will return different errors. Any help would be much appreciated.

Thank you,
James Knighton
James Knighton, modified 5 Years ago.

RE: Trouble Opening Delft3D-FLOW trih file in Matlab (Answer)

Youngling Posts: 12 Join Date: 8/9/13 Recent Posts
Problem resolved. The 'station' field is the index of the station within the stations array, not the station name.