intro story DELWAQ


DELWAQ is the engine of the D-Water Quality and D-Ecology programmes of the Delft3D suite. It is based on a rich library from which relevant substances and processes can be selected to quickly put water and sediment quality models together.

The processes library covers many aspects of water quality and ecology, from basic tracers, dissolved oxygen, nutrients, organic matter, inorganic suspended matter, heavy metals, bacteria and organic micro-pollutants, to complex algae and macrophyte dynamics. High performance solvers enable the simulation of long periods, often required to capture the full cycles of the processes being modelled.

The finite volume approach underlying DELWAQ allows it to be coupled to both the structured grid hydrodynamics of the current Delft3D-FLOW engine and the upcoming D-Flow Flexible Mesh engine (1D-2D-3D) of the Delft3D Flexible Mesh Suite (or even other models such as TELEMAC).

'DELWAQ in open source' is our invitation to all leading experts to collaborate in further development and research in the field of water quality, ecology and morphology using Delft3D. Feel free to post your DELWAQ related questions or comments in this dedicated forum space. If you are new to DELWAQ, the tutorial (in the user manual) is a good place to start. A list of DELWAQ related publications is available here.




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: https://svn.oss.deltares.nl/repos/delft3d/trunk/src/tools_lgpl/matlab/quickplot/progsrc/qpread.m $
% $Id: qpread.m 2498 2013-05-16 12:38:27Z jagers $