intro story D-Flow FM

 

D-Flow Flexible Mesh

D-Flow Flexible Mesh (D-Flow FM) is the new software engine for hydrodynamical simulations on unstructured grids in 1D-2D-3D. Together with the familiar curvilinear meshes from Delft3D 4, the unstructured grid can consist of triangles, pentagons (etc.) and 1D channel networks, all in one single mesh. It combines proven technology from the hydrodynamic engines of Delft3D 4 and SOBEK 2 and adds flexible administration, resulting in:

  • Easier 1D-2D-3D model coupling, intuitive setup of boundary conditions and meteorological forcings (amongst others).
  • More flexible 2D gridding in delta regions, river junctions, harbours, intertidal flats and more.
  • High performance by smart use of multicore architectures, and grid computing clusters.
An overview of the current developments can be found here.
 
The D-Flow FM - team would be delighted if you would participate in discussions on the generation of meshes, the specification of boundary conditions, the running of computations, and all kinds of other relevant topics. Feel free to share your smart questions and/or brilliant solutions! 

 

=======================================================
We have launched a new website (still under construction so expect continuous improvements) and a new forum dedicated to Delft3D Flexible Mesh.

Please follow this link to the new forum: 
/web/delft3dfm/forum

Post your questions, issues, suggestions, difficulties related to our Delft3D Flexible Mesh Suite on the new forum.

=======================================================

** PLEASE TAG YOUR POST! **

 

 

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

 


Message Boards

Boundary condition problem

VY
veli yumruktepe, modified 8 Years ago.

Boundary condition problem

Padawan Posts: 42 Join Date: 2/7/12 Recent Posts
Hi all,

For a long while, I was running my model with "per layer specified" boundary conditions. The domain has 30 layers.

Now I wanted to extend my domain to 60 layers. However as you know, for each boundary, you need to specify values for two ends. So you end up having 121 values for each line.
1(time)+120(condition).

When I start the run, I get several of these errors:
---------------------------------------------
warning: line number 130 of file
60l.bct


possibly too long.
Flow exited abnormally
Check diagnosis file
---------------------------------------------
these line numbers coincide with the initial time step of each boundary

Here is the output of diag file

*** MESSAGE Ocean heat model: Stanton number specified to be 0.13000E-02
*** MESSAGE Ocean heat model: Dalton number specified to be 0.13000E-02
*** MESSAGE Using UNESCO density formulation by default
*** MESSAGE Evaporation taken into account in continuity equation
*** MESSAGE Number of pivot points to convert wind speed to wind drag coef.: 2
*** MESSAGE Momentum solver cyclic method is specified
*** MESSAGE DRYFLP and DPSOPT both specified in MD-file. Using DPSOPT: MAX
*** MESSAGE Transport solver cyclic-method method is specified
*** MESSAGE Upwind advection scheme only near momentum discharges
*** ERROR Reading 60l.bct
Line 1591: table-name 'Boundary Section : 2'
Number of values (146100) in table is not multiple of number of parameters (121)


IS IT SOMETHING TO DO WITH LINE RESTRICTIONS? I REMEMBER SOMEWHERE I SAW MAXIMUN ALLOWABLE UNITS ARE "132" ?
OR AM I DOING SOMETHING WRONG?

BY THE WAY THESE BOUNDARY CONDITION FILES ARE CREATED WITH THE COMMERCIAL VERSION OF FLOW IN WINDOWS, IT DOES NOT WORK THERE EITHER
Adri Mourits, modified 8 Years ago.

RE: Boundary condition problem

Yoda Posts: 1221 Join Date: 1/3/11 Recent Posts
Hi Veli,

This file is read by the so called "table-module" which allows lines to be 1600 characters long.
May be the FLOW-GUI is not able to produce a file with 121 reals on a line (with the pre set format). Can you check that by opening file "60l.bct" in a text editor and counting the number of reals on a data line? Or divide the number of characters on a line by the space occupied by one real.

If this is the problem, a workaround will be to edit the bct-file manually, ensuring that each data line contains the correct number of reals.

Regards,

Adri
VY
veli yumruktepe, modified 8 Years ago.

RE: Boundary condition problem

Padawan Posts: 42 Join Date: 2/7/12 Recent Posts
Hi Adri,

Indeed my file contains
for each real 1-space and 15 number characters by multiplying 121 with 16 you exceed 1600 character restriction. I wil try to format the file myself. Does it have to be in E^(+)(-1) format
for example 1.0000E+001
VY
veli yumruktepe, modified 8 Years ago.

RE: Boundary condition problem

Padawan Posts: 42 Join Date: 2/7/12 Recent Posts
Hi Adri,

When I try to format the file myself, I get this error in diag file. Maybe there is a certain format option that I am missing.

*** MESSAGE Ocean heat model: Stanton number specified to be 0.13000E-02
*** MESSAGE Ocean heat model: Dalton number specified to be 0.13000E-02
*** MESSAGE Using UNESCO density formulation by default
*** MESSAGE Evaporation taken into account in continuity equation
*** MESSAGE Number of pivot points to convert wind speed to wind drag coef.: 2
*** MESSAGE Momentum solver cyclic method is specified
*** MESSAGE DRYFLP and DPSOPT both specified in MD-file. Using DPSOPT: MAX
*** MESSAGE Transport solver cyclic-method method is specified
*** MESSAGE Upwind advection scheme only near momentum discharges
*** ERROR Reading bound.bct
Line 130: 0.0000E+00 -3.0558E-02 -3.0489E-02 -3.0419E-02 -3.0350E-02 -3.0246E-02 -3.0108E-02 -2.9900E-02 -2.9623E-02 -2.9346E-02 -2.9166E-02 -2.9031E-02 -2.8896E-02 -2.8749E-02 -2.8652E-02 -2.8502E
Line parse error: too many sub-fields
*** ERROR Flow exited abnormally
Adri Mourits, modified 8 Years ago.

RE: Boundary condition problem

Yoda Posts: 1221 Join Date: 1/3/11 Recent Posts
Hi Veli,

The format (with or without exponent) does not matter, as long as they are all reals, separated with one or more spaces (not tabs).

The error message suggests that there are too many reals on line 130 in your bct file. Your first error, two posts above, contains an error in line 1591. May be you accidentally changed line 130 wrongly.
Can you try starting with the original file and start changing the input on line 1591?

Regards,

Adri
VY
veli yumruktepe, modified 8 Years ago.

RE: Boundary condition problem

Padawan Posts: 42 Join Date: 2/7/12 Recent Posts
Hi Adri,

I am not sure I understand your suggestion. At line 1591, a new table starts for the boundary no:2
so basically my tables are like 129 lines of header
1461 lines of data
Adri Mourits, modified 8 Years ago.

RE: Boundary condition problem

Yoda Posts: 1221 Join Date: 1/3/11 Recent Posts
Hi Veli,

It looks like the data related to the first boundary (line 129 to 1590) was OK and the data related to the second boundary (line 1719 to ...) contains errors. You can check this by using your original input and remove all boundaries except the first one. That should work.

My suggestion was to leave the data of boundary 1 unchanged and start editing at the data of the second boundary.

In case boundary 1 and boundary 2 need exactly the same data format, you have another option: copy the data(-format) from boundary 1 (line 129 to 1590) to boundary 2 (replacing line 1719 to 3180).

Regards, Adri
LS
LUIS HENRIQUE SILVA, modified 7 Years ago.

RE: Boundary condition problem

Youngling Posts: 12 Join Date: 9/30/11 Recent Posts
Adri Mourits:
Hi Veli,

It looks like the data related to the first boundary (line 129 to 1590) was OK and the data related to the second boundary (line 1719 to ...) contains errors. You can check this by using your original input and remove all boundaries except the first one. That should work.

My suggestion was to leave the data of boundary 1 unchanged and start editing at the data of the second boundary.

In case boundary 1 and boundary 2 need exactly the same data format, you have another option: copy the data(-format) from boundary 1 (line 129 to 1590) to boundary 2 (replacing line 1719 to 3180).

Regards, Adri


Hi Adri,

After to have the error message

too many sub-fields

I carried out several tests increasing the number of levels. I got this error message when the
number of levels was greater than 50 levels. This means that a line of my bcc file has 101
numbers (1 + 50 + 50).

This message is generated by the subroutine scannr, where the number of sub-strings in each line of the file bcc can not be greater than maxfld (each sub-string is a group of digits separated by space, or ', or ").

In the original code maxfld = 100, that is the cause of this error message. I changed
to maxfld = 120 in the module tables.f90, and the Delft3D is running very well.

Regards

Luis