Delft3D-NESTWQ
Offline nesting of water quality models
User Manual
April, 2004
User Manual Delft3D-NESTWQ
Published and printed by:
WL | Delft Hydraulics
main office Rotterdamseweg 185 p.o. box 177 2600 MH Delft The Netherlands telephone fax
internet www
For support contact: e-mail: www: telephone: fax:
Copyright © 2004 WL | Delft Hydraulics
All rights reserved. No part of this document may be reproduced in any form by print, photo print, photo copy, microfilm or any other means, without written permission from the publisher: WL | Delft Hydraulics.
delft3d.support@wldelft.nl http://www.wldelft.nl/d3d/ +31 15 285 85 55 +31 15 285 85 82 +31 15 285 85 85 +31 15 285 85 82 info@wldelft.nl http://www.wldelft.nl
Delft3D-NESTWQ Contents User Manual April, 2004 Version 1.31
Contents
1
2
3 4 5 6 WL | Delft Hydraulics
A guide to this manual..........................................................................................1–1 1.1 Introduction................................................................................................1–1 1.2 Manual version and revisions....................................................................1–1 1.3 Typographical conventions........................................................................1–2 1.4 Changes with respect to previous versions................................................1–2
Introduction to Delft3D-NESTWQ.....................................................................2–1 2.1 System name and subsystems....................................................................2–1 2.2 Functionality..............................................................................................2–1 2.3
Overview subsystems................................................................................2–2
Getting started.......................................................................................................3–1 3.1 Starting Delft3D.........................................................................................3–1 3.2 Getting into Delft3D-NESTWQ................................................................3–3 3.3
Exploring some menu options...................................................................3–6
Input and output...................................................................................................4–1 4.1 Input NESTWQ1.......................................................................................4–1 4.2 Output NESTWQ1.....................................................................................4–2 4.3 Input NESTWQ2.......................................................................................4–2 4.4
Output NESTWQ2.....................................................................................4–3
Tutorial..................................................................................................................5–1 5.2 Introduction................................................................................................5–3 5.3 Input and output NESTWQ1.....................................................................5–5 5.4
Input and output NESTWQ2.....................................................................5–7
Limitations, warnings and errors........................................................................6–1
i
April, 2004 Version 1.31 Contents Delft3D-NESTWQ
User Manual
6.1 6.2 6.3 7 8
Limitations................................................................................................6–1 Warnings....................................................................................................6–2 Errors.........................................................................................................6–2
Description files....................................................................................................7–1 Theoretical background.......................................................................................8–1 8.1 8.2
Boundary segments...................................................................................8–1 Nest segments and weights........................................................................8–1
9
References.............................................................................................................9–1
ii
WL | Delft Hydraulics
Delft3D-NESTWQ User Manual A guide to this manual
April, 2004 Version 1.31
1
1.1
A guide to this manual
Introduction
This User Manual concerns the offline water quality nesting module, Delft3D-NESTWQ, of the Delft3D software suite. To make this manual more accessible we will briefly describe the contents of each chapter and appendix.
If this is your first time to start working with Delft3D-NESTWQ we suggest you to read and practice the getting started of Chapter 3 and the tutorial of Chapter 5. These chapters explain the user input and guide you through the definition of your first nested simulation.
Chapter 2: Introduction to Delft3D-NESTWQ, provides specifications of Delft3D-NESTWQ, such as the areas of applications, the functionality provided and an overview of the system.
Chapter 3: Getting started, explains the use of the overall menu program, which gives access to all Delft3D modules (for which you have authorisation) and to the pre- and post-processing tools. Delft3D-NESTWQ being one of them as part of the WAQ module.
Chapter 4: Input and output, provides detailed information on the parameters in the input file and the output files from Delft3D-NESTWQ.
Chapter 5: Tutorial, emphasis at giving you some first hands-on experience in using the Delft3D-NESTWQ tool.
Chapter 6: Limitations, warnings and errors, discusses the underlying assumptions and warnings and errors that might occur.
Chapter 7: Description files, explains in detail the contents of input and output files.
Chapter 8: Theoretical background, provides information about boundary segments, nest segments and weights.
Chapter 9: References, provides a list of related Delft3D modules and background information.
1.2 Manual version and revisions
A manual applies to a certain release of the related numerical program. This manual, version 1.30, applies to Delft3D-NESTWQ version 2.00.03.
WL | Delft Hydraulics
1–1
April, 2004 Version 1.31 A guide to this manual
Delft3D-NESTWQ
User Manual
1.3 Typographical conventions
Throughout this manual, the following conventions help you to distinguish between different elements of text to help you learn about Delft3D and NESTWQ.
Example Description Title of a window or sub-window. Water quality
Sub-windows are displayed in the WAQ Module window
and can not be moved. Tools
Item from a menu, title of a push button or the name of a user interface input field.
Upon selecting this item (click or in some case double click with the left mouse button on it) a related action will be executed; in most cases it will result in displaying some other (sub-)window.
In case of an input field you are supposed to enter input data of the required format and in the required domain.
Directory names, filenames, and path names are expressed between smaller than, <, and larger than, >, signs. For the Linux and UNIX environment a forward slash (/) is used instead of the backward slash (\\) for PCs.
Data to be typed by you into the input fields are displayed between double quotes.
Selections of menu items, option boxes etc. are described as such: for instance ‘select Save and go to next window’. Commands to be typed by you are given in the font Courier New, 10 points.
Units are given between square brackets when used next to the formulae. Leaving them out might result in misinterpretation.
<\utorial\\waq\\nestwq> “27 08 1999” delft3d-menu [m/s] [-] 1.4 Changes with respect to previous versions Version Description 1.31 Operation of Change working directory updated in Chapter 3. Extension of Delft3D-MENU with Delft3D-QUICKPLOT. 1.30 Reference version for these change notes. 1–2 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Introduction to Delft3D-NESTWQ April, 2004 Version 1.31 2 Introduction to Delft3D-NESTWQ The transfer of data from an encompassing or 'overall' numerical model to an embedded or 'nested' numerical model is called nesting. In general the overall model has a coarse resolution of grid cells, whereas the nested model has a higher resolution. At the boundary locations of the nested model the results from the overall model are required as boundary conditions for the nested model. The boundary conditions can be water levels, currents, fluxes or discharges in case of hydrodynamic models, and water quality parameters in case of water quality models. 2.1 System name and subsystems The procedure of nesting through concentrations between Delft3D-WAQ (or PART) models is performed by the system NESTWQ. In this procedure two steps can be distinguished which are handled by separate subsystems: • NESTWQ1, for the determination of nest segments and nest weights in the overall model. The concentrations at these segments are used by the next subsystem. • NESTWQ2, for the generation of boundary conditions for the boundary segments in the nested model from the results at the nest segments in the overall model. 2.2 Functionality The functionality of the subsystem NESTWQ1 is: • determine for each boundary segment in the nested model which nest segments in the overall model are required for the nesting of results to these boundary segments • determine for each nest segment the weight factors for the interpolation from the surrounding nest segments to the boundary segment in the nested model The functionality of the subsystem NESTWQ2 is: • interpolate the concentration time-series from the nest segments in the overall model to boundary conditions for the boundary segments in the nested model • average these boundary conditions in case of aggregated boundary segments in the nested model • generate a Delft3D-WAQ boundary conditions file for the nested model WL | Delft Hydraulics 2–1 April, 2004 Version 1.31 Introduction to Delft3D-NESTWQ Delft3D-NESTWQ User Manual 2.3 Overview subsystems Data flow NestWQnestwq1.inp *.lga overall *.lga detail *.cco overall *.cco detail NESTWQ1Administration *.admDiagnostic *.dianestwq2.inp *.map overall nr of substances NESTWQ2Diagnostic *.diaBoundaries for detail model *.bcc run detail model Figure 2-1 Data flow diagram of Delft3D-NESTWQ 2–2 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Getting started April, 2004 Version 1.31 3 3.1 Getting started Starting Delft3D To start Delft3D: • On an MS Windows platform: select Delft3D in the Applications menu or click on the Delft3D icon on the desktop. • On Linux and UNIX machines: type delft3d-menu on the command line. Next the title window of Delft3D is displayed, Figure 3-1: Figure 3-1 Title window of Delft3D After a short while the main window of the Delft3D-MENU appears, Figure 3-2. WL | Delft Hydraulics 3–1 April, 2004 Version 1.31 Getting started Delft3D-NESTWQ User Manual Figure 3-2 Main window Delft3D-MENU Several menu options are shown. In Figure 3-2 all options are sensitive. If options are insensitive this indicates that you are not authorised to use this option; check your license file ( • Click on the Exit push button. The window will be closed and you are back in the Windows Desk Top screen for PCs or on the command line for Linux and UNIX workstations. You can exit Delft3D-MENU in two other ways: • Select File - Exit in the toolbar menu or • Click the ‘close’ button (x) in the top-right corner of the menu bar. The first method is to be preferred. Remark: • In this and the following chapters several windows are shown to illustrate the presentation of Delft3D-MENU and Delft3D-NESTWQ. These windows are grabbed from the PC-platform. For Linux and UNIX workstations the content of the windows is the same, but the colours may be different. On the PC-platform you can set your preferred colours by using the Display Properties. 3–2 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Getting started April, 2004 Version 1.31 3.2 Getting into Delft3D-NESTWQ To continue restart the menu program as indicated above. To access the nesting tool between an overall water quality model an a detailed water quality model, select the far-field water quality module. • Click the Water Quality button. Next the selection window for Far-field water quality (all modules) is displayed, see Figure 3-3. Figure 3-3 Selection window for Far-field water quality (all modules) The nesting functionality is exactly the same for all four water quality modules, therefore: • Click the General button. The selection window for Water quality (general) is displayed (Figure 3-4), in which you can convert and/or aggregate hydrodynamic results, prepare a substance file using the Process Library Configuration Tool, prepare a water quality input file, run the pre-processor, execute a computation in foreground or in background, inspect the monitoring files with information on the execution, visualise the results, and access additional tools. WL | Delft Hydraulics 3–3 April, 2004 Version 1.31 Getting started Delft3D-NESTWQ User Manual Figure 3-4 Selection window for Water Quality (general) Before continuing with any of the selections of this Water Quality window, you must select the directory in which you are going to nest water quality results: • Click the Change working directory button. Next the Browse for Folder window, Figure 3-5, is displayed (your current directory may differ, depending on the location of your Delft3D installation). Figure 3-5 Browse for Folder window • Navigate to and open the 3–4 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Getting started April, 2004 Version 1.31 • Open the • Select the Figure 3-6 Browse for Folder window to set the working directory to Next the Water Quality window is re-displayed, but now the changed current working directory is displayed in the title bar, see Figure 3-7. Figure 3-7 Current working directory The nesting tools are part of the Additional tools, hence: • Click on Tools. The Nesting programs for Water Quality window, Figure 3-8, contains the supporting programs for the various steps in the nesting of water quality models. WL | Delft Hydraulics 3–5 April, 2004 Version 1.31 Getting started Delft3D-NESTWQ User Manual Figure 3-8 Selection window for Nesting Water Quality • Click on Nest Input (1). The input file for NESTWQ1 is opened, see Figure 3-9. Figure 3-9 Input file for NESTWQ1 • Select File – Exit, to close the input file. • Select Return, to close the Nesting programs for Water Quality window. • Select Return, and Return once again, to enter the main window of Delft3D-MENU. • Click Exit. The window is closed and the control is returned to the desk top or the command line. 3.3 Exploring some menu options To guide you through some menu options please follow the example in Chapter 4, Tutorial. 3–6 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Input and output April, 2004 Version 1.31 4 Input and output NESTWQ1 and NESTWQ2 have each their own input file. The input files 4.1 Input NESTWQ1 The input file locations overall model locations nested model segments in Delft3D-WAQ monitor format and weight factors required for the nesting procedure in NESTWQ2 The filenames between quotes are free to choose with a maximal length of 128 characters. Remarks: • The required information about the vertical structure of the models is contained in the grid files. If the overall model is 2DH and the nested model is 3D, NESTWQ2 will generate a uniform 3D profile for the nested model. If the overall model is 3D and the nested model is 2DH, NESTWQ2 will generate depth-averaged boundary conditions for the nested model. • If both models are 3D then the number of layers in each model must be equal. Example input file '..\\hkcrs.3d\\com-o3d.cco' ; X,Y co-ordinates overall model '..\\hkcrs.3d\\com-o3d.lga' ; grid table with active cells overall model '..\\siulam.3d\\com-d3d.cco' ; X,Y co-ordinates nested model '..\\siulam.3d\\com-d3d.lga' ; grid table with active cells nested model '3d-3d.dia' ; diagnostics from NESTWQ1 '3d-3d.adm' ; nest segments and weight factors overall model WL | Delft Hydraulics 4–1 April, 2004 Version 1.31 Input and output Delft3D-NESTWQ User Manual 4.2 Output NESTWQ1 The output files of NESTWQ1 are: Delft3D-WAQ monitor format required for the nesting procedure in NESTWQ2 Remarks: • It is not necessary to include the required nest segments in the simulation with the overall model as the nesting is performed with the results in the map file. The map file contains the concentrations at every grid point, thus including the required nest segments. • The administration file contains the required information about the vertical structure of each model. See Chapter 6 for an example of these files. 4.3 Input NESTWQ2 The input file segment the overall nest segments and weight factors required for the nesting procedure in NESTWQ2 boundary segments of the nested model The filenames between quotes are free to choose with a maximal length of 128 characters. Remark: • The first ‘nosys’ active substances will be transferred as boundary conditions for the nested model. Example input file '..\\nestwq1.v20\\3d-3d.adm' ; nest segments and weight factors overall model 'o3d-waq.map' ; map file overall model 4 ; number of active substances to nest '3d-3d.dia' ; diagnostics from NESTWQ2 '3d3dwaq.bcc' ; boundary conditions for nested model 4–2 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Input and output April, 2004 Version 1.31 4.4 Output NESTWQ2 The output files of NESTWQ2 are: the nested model See Chapter 6 for an example of the diagnostics file. WL | Delft Hydraulics 4–3 Delft3D-NESTWQ Tutorial User Manual April, 2004 Version 1.31 5 Tutorial At this moment no graphical user interface is available to perform the nesting of water quality models. In order to perform the nesting, correctly, the following steps are required: Edit/Create You have to specify the names of the <*.lga> and <*.cco> files of the overall and detailed model. Specify the administration file. On forehand the hydrodynamics of both the overall and nested model must be converted to formats required for Delft3D-WAQ (Couple menu, see Delft3D-WAQ User Manusl). Create administration file (execute NESTWQ1) Within this step a diagnostic file and an administration file is created which is required input for the second step of the of the nesting (NESTWQ2). Edit/Create You must specify the name of the administration file (created by step 2), the name of the water quality map file of the overall model , the number of active substances to nest and the name of the binary boundary output file for the detailed model. Create boundary conditions (execute NESTWQ2) Within this step a diagnostic file and a boundary file is created which is required input for running the water quality detailed model. Examine diagnostic files of NESTWQ1 and NESTWQ2 After executing NESTWQ1 and/or NESTWQ2 you can inspect the if the modules have run correctly. Edit nested WAQ input file In this step you must add the created boundary filename (step 4) by editing block 4 of the water quality input file: • goto block #4 • change ‘1 : information in this file’ to ‘-2 : information in binary file’ • write on new record the name of boundary file between quotes • delete all other information in the following records in block #4 (until block #5) In order to demonstrate the use of Delft3D-NESTWQ the following tutorial can be applied: WL | Delft Hydraulics 5–1 April, 2004 Version 1.31 Tutorial Delft3D-NESTWQ User Manual Available data files: 1 Start Delft3D-MENU 2 Select Water Quality from the main menu 3 Select General from water quality model 4 Activate in WAQ-GUI data group Hydrodynamics and select 7 Activate in WAQ-GUI data group Hydrodynamics and select 9 Activate data group ‘Substance’ and select file 12 Add two discharges at (32,28) and (44,28) and assign a load of 1000 g/s of CTR1 to both discharges (data group ‘Discharges’) by button ‘edit data’ 13 Activate ‘Output Option’ data group and set output timers for the map file to 1-1-2000 (start) ; 3-1-2000 (stop) ; 1 hour time step 14 Save scenario as <*.scn> file (e.g. 16 Select 3 ‘Run the pre-processors’ (file 19 Visualise the results with Delft3D-GPP Repeat the above steps for the detailed model ( 20 Select Additional tools (NESTWQ) in the Delft3D Far-field modelling menu. The following menu will appear: 1. Edit/Create You have to specify the names of the <*.lga> and <*.cco> files of the overall and detailed model conform the documentation (e.g. 2. Create administration file (run NESTWQ1) Within this step a diagnostic file and an administration file is created which is required input for the second step of the of the nesting (NESTWQ2) 3. Edit/Create You must specify the name of the administration file (created by step 2), the name of the water quality map file of the overall model (e.g 5–2 WL | Delft Hydraulics Delft3D-NESTWQ Tutorial User Manual April, 2004 Version 1.31 number of active substances to nest and the name of the binary boundary output file for the detailed model (e.g. <2d-2d.bcc>). 4. Create boundary conditions (run NESTWQ2) Within this step a diagnostic file and an boundary file is created which is required input for running the water quality detailed model. 5. Examine diagnostic files of NESTWQ1 and NESTWQ2 6. Edit nested WAQ input file In this step you must add the created boundary file name (step 4) by editing block 4 of the water quality input file (e.g. - change ‘1 : information in this file’ to ‘-2 : information in binary file’ - write on new record the name of boundary file between quotes e.g. ‘ 2d-2d.bcc ’ - delete all other information in the following records in block #4 (until block #5) 21 Return to Delft3D-far field modelling menu and select 3 ‘Run the preprocessors’ and select the adapted 22 Run the detailed model and evaluate the results. 5.2 Introduction The Upgraded WAHMO Model is a fine-grid numerical model for the Hong Kong waters. It is based on the shallow water equations in curvilinear co-ordinates. The spatial resolution is typical 300-1500 m. The model was developed in view of the need for a detailed and accurate tidal and water quality model that could be used for accurate computations. A further interest was provided by transport and water quality simulations on seasonal and annual scales. For this purpose the fine-grid water quality model was aggregated. The SIU LAM model is a fine-grid numerical model for the waters north of Lantau Island. It uses a curvilinear fine grid. The hydrodynamic model was developed to investigate the effects of various lay-outs for a proposed typhoon shelter at Siu Lam. As an example to demonstrate the nesting of water quality models a water quality model based upon the hydrodynamic SIU LAM model was set up. For the generation of water quality boundary conditions for the WQ SIU LAM model, the model is nested in a coarse version of the WQ UPGRADE WAHMO model. Figure 5-1 shows the coarse grid UPGRADE WAHMO together with the fine grid SIU LAM. The locations of the boundary segments of the SIU LAM model are presented in Figure 5-2. WL | Delft Hydraulics 5–3 April, 2004 Version 1.31 Tutorial Delft3D-NESTWQ User Manual Figure 5-1 Upgrade and Siu Lam models Figure 5-2 Siu Lam model with 2 open boundaries (east and west) 5–4 WL | Delft Hydraulics Delft3D-NESTWQ Tutorial User Manual April, 2004 Version 1.31 In the following application grid cells have not been aggregated in both models. Both 3D models have 5 layers. 5.3 Input and output NESTWQ1 Input file '..\\hkcrs.3d\\com-o3d.cco' ; X,Y co-ordinates overall model '..\\hkcrs.3d\\com-o3d.lga' ; grid table with active cells overall model '..\\siulam.3d\\com-d3d.cco' ; X,Y co-ordinates nested model '..\\siulam.3d\\com-d3d.lga' ; grid table with active cells nested model '3d-3d.dia' ; diagnostics NESTWQ1 '3d-3d.adm' ; nest segments and weight factors overall model Diagnostics file <3d-3d.dia>: *** Delft3D-WAQ utility NESTWQ1 *** v2.0, Mar 1998 *** >> Input from nestwq1.inp X,Y co-ordinates overall model : ..\\hkcrs.3d\\com-o3d.cco Grid table with active cells overall model : ..\\hkcrs.3d\\com-o3d.lga X,Y co-ordinates nested model ; ..\\siulam.3d\\com-d3d.cco Grid table with active cells nested model : ..\\siulam.3d\\com-d3d.lga Diagnostics and required nest segments : 3d-3d.dia Nest segments and weight factors overall model : 3d-3d.adm *** WARNING *** negative 3rd weight for -9 1 *** WARNING *** negative 4th weight for -9 1 *** WARNING *** negative 3rd weight for -10 1 *** WARNING *** negative 4th weight for -10 1 *** WARNING *** negative 3rd weight for -53 1 *** WARNING *** negative 4th weight for -53 1 >> List of required nest segments in DELWAQ monitor format Number of required nest segments : 51 1719 'Nest segment 1719 ' 1720 'Nest segment 1720 ' 1782 'Nest segment 1782 ' 1781 'Nest segment 1781 ' 1780 'Nest segment 1780 ' 1718 'Nest segment 1718 ' 2897 'Nest segment 2897 ' 2959 'Nest segment 2959 ' 1779 'Nest segment 1779 ' 1717 'Nest segment 1717 ' 2958 'Nest segment 2958 ' 2896 'Nest segment 2896 ' WL | Delft Hydraulics 5–5 April, 2004 Version 1.31 Tutorial Delft3D-NESTWQ User Manual 1778 'Nest segment 1778 ' 1716 'Nest segment 1716 ' 2957 'Nest segment 2957 ' 2895 'Nest segment 2895 ' 1654 'Nest segment 1654 ' 1655 'Nest segment 1655 ' 1715 'Nest segment 1715 ' 1653 'Nest segment 1653 ' 1714 'Nest segment 1714 ' 1652 'Nest segment 1652 ' 1713 'Nest segment 1713 ' 1651 'Nest segment 1651 ' 2956 'Nest segment 2956 ' 2894 'Nest segment 2894 ' 1712 'Nest segment 1712 ' 1650 'Nest segment 1650 ' 1711 'Nest segment 1711 ' 1649 'Nest segment 1649 ' ....... ....... 1766 'Nest segment 1766 ' 1890 'Nest segment 1890 ' 1889 'Nest segment 1889 ' 1827 'Nest segment 1827 ' 1888 'Nest segment 1888 ' 1826 'Nest segment 1826 ' >> Number of warnings : 6 Nest administration data file <3d-3d.adm> (only part is given): *** Delft3D-WAQ utility NESTWQ1 *** v2.0, Mar 1998 *** >> Files used to determine nest characteristics X,Y co-ordinates overall model : ..\\hkcrs.3d\\com-o3d.cco Grid table with active cells overall model : ..\\hkcrs.3d\\com-o3d.lga Number of segments overall model per layer : 4658 Number of layers overall model : 5 X,Y co-ordinates nested model : ..\\siulam.3d\\com-d3d.cco Grid table with active cells nested model : ..\\siulam.3d\\com-d3d.lga Number of segments nested model per layer : 7257 Number of layers nested model : 5 >> List of boundary segments, required nest segments and weight factors Nest model Overall model ---------- -------------- SEGNR #NST SEGNR WEIGHT -1 4 1719 0.47121 1720 0.28688 1782 0.09155 1781 0.15037 SEGNR #NST SEGNR WEIGHT -2 4 1719 0.68195 1720 0.09224 1782 0.02690 1781 0.19890 ........... SEGNR #NST SEGNR WEIGHT -6 2 2897 0.83103 2959 0.16897 SEGNR #NST SEGNR WEIGHT -7 4 1718 0.62568 5–6 WL | Delft Hydraulics Delft3D-NESTWQ Tutorial User Manual April, 2004 Version 1.31 1719 0.20044 1781 0.04219 1780 0.13169 ....... ...... -68 4 1827 0.43794 1889 0.25714 1888 0.11280 1826 0.19211 5.4 Input and output NESTWQ2 Input file '..\\nestwq1.v20\\3d-3d.adm' ; nest segments and weights overall model 'o3d-waq.map' ; map file overall model 4 ; number of active substances to nest '3d-3d.dia' ; diagnostics NESTWQ2 '3d3dwaq.bcc' ; boundary conditions for nested model The diagnostics file <3d-3d.dia>: *** Delft3D-WAQ utility NESTWQ2 *** v2.0, Mar 1998 *** >> Input from nestwq2.inp Nest segments and weight factors overall model : ..\\nestwq1.v20\\3d-3d.adm Map file overall model : o3d-waq.map Number of active substances to nest : 4 Diagnostics NESTWQ2 : 3d-3d.dia Boundary conditions for nested model : 3d3dwaq.bcc *** Delft3D-WAQ utility NESTWQ1 *** v2.0, Mar 1998 *** >> Files used to determine nest characteristics X,Y co-ordinates overall model : ..\\hkcrs.3d\\com-o3d.cco Grid table with active cells overall model : ..\\hkcrs.3d\\com-o3d.lga Number of segments overall model per layer : 4658 Number of layers overall model : 5 X,Y co-ordinates nested model : ..\\siulam.3d\\com-d3d.cco Grid table with active cells nested model : ..\\siulam.3d\\com-d3d.lga Number of segments nested model per layer : 7257 Number of layers nested model : 5 >> Delft3D-WAQ map file characteristics WL | Delft Hydraulics 5–7 April, 2004 Version 1.31 Tutorial Delft3D-NESTWQ User Manual Model and run identification : 1- Delft3D-WAQ ----------------------- 1 2- Upgrade WAHMO coarse model -------- 2 3- test ------------------------------ 3 4- map results to nest SIU LAM ------- 4 Total number of variables : 7 Variable names on map file : Substance 1 Substance 2 Substance 3 Substance 4 Inactive subst 5 Inactive subst 6 Extra variable 7 Selected variables to nest : Substance 1 Substance 2 Substance 3 Substance 4 Total number of segments : 23290 >> Number of warnings : 0 >> Number of error messages : 0 5–8 WL | Delft Hydraulics Delft3D-NESTWQ User Manual 6 6.1 WL | Delft Hydraulics Limitations, warnings and errors April, 2004 Version 1.31 Limitations, warnings and errors Limitations • The X,Y co-ordinates of the overall grid and the nested grid are stored in so-called (binary) cco-files. These files are platform dependent. • The overall and nested grid have the same co-ordinate system. • The time-series boundary conditions for the (aggregated) nested model are derived from the binary Delft3D-WAQ map file from the (aggregated) overall model. • The first NOSYS active substances on the map file of the overall model are transferred to the nested model; selection of substances is not possible. • The location of a boundary segment in the nested model is obtained through mirroring from the interior domain. • The centre of a boundary segment in the nested model before aggregation, and the centres of the surrounding overall nest segments determine the weight factors. • The nest procedure transfers 2DH and 3D overall model results to 2DH and 3D boundary conditions for the nested model. In case both models are 3D then the number of layers in both models must be equal. 2DH overall model results will be transformed to uniform 3D boundary conditions for the nested model. 3D overall model results will be transformed to depth-averaged boundary conditions if the nested model is 2DH. • The time-series boundary conditions file for the nested model is binary and only suitable for Delft3D-WAQ. • The time span and time interval for the time-series boundary conditions of the nested model are determined by the map time span and interval of the map file from the overall model. • The nesting procedure Delft3D-NESTWQ does not account for drying in the overall model. • If there are missing boundary segment numbers (gaps in the numbering) they will get the value 0.0 as concentration in the boundary conditions file. 6–1 April, 2004 Version 1.31 Limitations, warnings and errors Delft3D-NESTWQ User Manual 6.2 Warnings Warnings from subsystem NESTWQ1 can be: *** WARNING *** centre of boundary segment -xxxxx, (M,N) = (xxx,xxx) not inside a nest segment the nearest nest segment will be taken When the centre of a boundary segment is not inside the polygon spanned by the centres of the nest segments, extrapolation takes place. This may cause one or more negative weight factors. Consequently, this may lead to negative concentrations. To avoid this, NESTWQ2 will reset negative concentrations to 0.0. The following four warnings may occur: *** WARNING *** negative 1st weight for -xxx xxx *** WARNING *** negative 2nd weight for -xxx xxx *** WARNING *** negative 3rd weight for -xxx xxx *** WARNING *** negative 4th weight for -xxx xxx When the sequence of boundary segment numbers is interrupted, the following warning is issued. In the time-series file for the nested model a missing boundary segment will have zero concentrations. *** WARNING *** boundary segment -xxx is missing NESTWQ2 gives no warnings. 6.3 Errors Errors from subsystem NESTWQ1 can be: *** ERROR *** The file nestwq1.inp does not exist The following error is related to missing grid or segment numbers files. *** ERROR *** non-existing file If the format or structure of the input, <*.cco> or *** ERROR *** while reading nestwq1.inp *** ERROR *** premature EOF nestwq1.inp *** ERROR *** while reading 6–2 WL | Delft Hydraulics Delft3D-NESTWQ User Manual Limitations, warnings and errors April, 2004 Version 1.31 *** ERROR *** premature EOF The following dimensions may be too small; they can be increased in the main program. *** ERROR *** increase MNMAX to xxxxx *** ERROR *** increase MAXBND to xxxxx *** ERROR *** increase MAXAGR to xxxxx *** ERROR *** increase MAXNST to xxxxx When the <*.cco> or <*.lga> file mismatch, the errors are: *** ERROR *** inconsistent MMAX MMAX-cco = xxxxx MMAX-lga = xxxxx *** ERROR *** inconsistent NMAX NMAX-cco = xxxxx NMAX-lga = xxxxx *** ERROR *** inconsistent LAYT LAYT-cco = xxxxx LAYT-lga = xxxxx The next two errors may not occur (when they occur your data files were corrupted). *** ERROR *** increase NMAX in pinpol to xxxx *** ERROR *** number of nest segments : xxxxx at boundary segment -xxx xxx Failures in the determination of the weight factors may cause the errors: *** ERROR *** singular points for ichoic = x *** ERROR *** no convergence newton-raphson iteration ichoic = xxx Errors from subsystem NESTWQ2 can be: *** ERROR *** The file nestwq2.inp does not exist WL | Delft Hydraulics 6–3 April, 2004 Version 1.31 Limitations, warnings and errors Delft3D-NESTWQ User Manual The following error is related to missing nest administration data from NESTWQ1 or Delft3D-WAQ map file. *** ERROR *** non-existing file If the format or structure of the input, nest administration or Delft3D-WAQ map file is incorrect, the following errors will be given. *** ERROR *** while reading nestwq2.inp *** ERROR *** premature EOF nestwq2.inp *** ERROR *** while reading The following dimensions may be too small; they can be increased in he main program. *** ERROR *** increase MAXVAR to xxxxx *** ERROR *** increase MAXSYS to xxxxx *** ERROR *** increase MAXBND to xxxxx *** ERROR *** increase MAXAGR to xxxxx *** ERROR *** increase MAXSEG to xxxxx When there is an inconsistency between the files used for NESTWQ1 and the overall map file used in NESTWQ2 the error will be: *** ERROR *** inconsistent number of segments NOSEG from NESTWQ1 = xxxxx NOSEG from mapfile = xxxxx An inconsistency in the number of layers for 3D models will cause the error: *** ERROR *** 3D-3D nesting with different number of layers is NOT implemented 6–4 WL | Delft Hydraulics Delft3D-NESTWQ Description files User Manual April, 2004 Version 1.31 7 Description files The files used by NESTWQ1 and NESTWQ2 are described as follows. Name : Name Type Access Contents 2 Name Type Format Contents : the - - - Name Type Format Contents : the - - - Name Type Access Contents 2 3 WL | Delft Hydraulics : transparent empty record 2 2 integers; 3 reals, 2 integers; MMAX, NMAX, XCOR1, YCOR1, 0.0, 0, LAYT 9 records with 0.0 NMAX*MMAX reals; (X(k), k=1,NMAX*MMAX) NMAX*MMAX reals; (Y(k), k=1,NMAX*MMAX) : : free-formatted file contains information about: subsystem name and version number the names of the input and output files of NESTWQ1 the required (overall) segments in Delft3D-WAQ monitor format - warnings and errors occuured while executing NESTWQ1 : : free-formatted file contains information about: subsystem name and version number the filenames which have been used in the nesting procedure - the number of segments per layer in both models the number of layers in both models - per boundary segment the required (overall) nest segments and their weight factors : 7–1 April, 2004 Version 1.31 Description files Delft3D-NESTWQ User Manual 4 for each time step, an integer and NOVAR * NOSEG reals; itime, ((CONC(ivar,iseg),ivar=1,NOVAR),iseg=1,NOSEG) Name : Format : free-formatted Contents : the file contains information about: - subsystem name and version number - the names of input and output files used by NESTWQ2 - the number of active substances - - - - - - Name Type Access Contents itime, 7–2 the names of input files used by NESTWQ1 model and run identification of the simulation with the overall model total number of variables on the map file variable names on the map file selected variables to nest warnings and errors occurred while executing NESTWQ2 : : for each time step, an integer and NOSYS * NOBND reals; ((CONC(isys,ibnd),isys=1,NOSYS),ibnd=1,NOBND) WL | Delft Hydraulics Delft3D-NESTWQ User Manual Theoretical background April, 2004 Version 1.31 8 8.1 Theoretical background Boundary segments Boundary segments are segments with a negative segment number. In case of an aggregated grid, there may be more than one boundary segment with the same negative number. Normally, when aggregation is done, interior grid cells will be aggregated; the boundary segments will not be aggregated. In the nesting procedure the time-series from the nest segments are interpolated to the centre of each (individual) boundary segment. If boundary segments are aggregated these interpolated time-series are then averaged over the boundary segments with the same segment number. 8.2 Nest segments and weights The centre of each boundary segment determines in which (overall) nest segment the centre is located. Depending on the quadrant within this nest segment, maximal three other adjacent nest segments are determined. It may be possible that 3, 2, 1 or no other nest segments can be specified. For instance, when the first nest segment is at the edge of the overall grid, only 1 additional nest segment may be identified. When the centre of a boundary segment is not included in any overall segment then the nearest overall segment is taken as the only nest segment. The criterion for the nearest segment is that segment with the shortest distance between its centre and the centre of the boundary segment. The centres of the nest segments and the centre of the boundary segment determine the weights for each nest segment. The options are linear (2 and 3 nest segments) and bi-linear (4 nest segments) interpolation. WL | Delft Hydraulics 8–1 Delft3D-NESTWQ References User Manual April, 2004 Version 1.31 9 References Delft Hydraulics, 2003a. Delft3D-WAQ User Manual, version 4.00. Delft Hydraulics, 2003b. Delft3D-PART User Manual, version 2.00. Delft Hydraulics, 2004c. Delft3D-GPP User Manual, version 2.11. Hulsen, L.J.M., 1994. Definition study Coupling between the CSM-8 and the KUSTSTROOK model, Docid t1222t01.wp4 (in Dutch). Hulsen, L.J.M., 1995. Functional and technical design Coupling between the CSM-8 and the KUSTSTROOK model, Docid t1222t02.wp8 (in Dutch). Mooiman, J., 1994. Linear interpolation in 1, 2 and 3 dimensions (in Dutch). WL | Delft Hydraulics 9–1 因篇幅问题不能全部显示，请点此查看更多更全内容