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. 




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



Message Boards

underlayers and qpread()

Christopher Esposito, modified 3 Years ago.

underlayers and qpread()

Padawan Posts: 35 Join Date: 10/16/12 Recent Posts
Hello. I am trying to use the matlab function qpread() to extract the underlayer grid from my model, but I cannot figure out the syntax.

In the past I had been using vs_let to extract the THLYR parameter, and construct my underlayers by sequentially subtracting thickness from the bed, but I've been noticing some differences between my plots and those made by quickplot. So I think it's best if I switch to using qpread.

Below is the code that I am trying. Either of the qpread calls below will give me an error that says "plot type is not supported for underalyers". I've tried several different similar calls, but haven't gotten it to work. But since I can make underlayer plots in quickplot, I know that this data is accessible somehow...

>> I=qpread(trim);

>> I(62)

ans =

Name: 'sediment volume fraction'
Units: '-'
DimFlag: [1 0 1 1 1]
DataInCell: 1
NVal: 1
VecType: ''
Loc: 'z'
ReqLoc: 'z'
Loc3D: 'b'
Group: 'map-sed-series'
Val2: ''
SubFld: 'sb1'
MNK: 0
UseGrid: 2

>> qpread(trim,I(62),'griddata',22)

>> qpread(trim,1,I(62),'griddata','Sediment_Sand175',22,0,0,0)

And here is the version information on qpread.

% $HeadURL: $
% $Id: qpread.m 2498 2013-05-16 12:38:27Z jagers $