W2 Version 3.12 Release Notes W2 V3.1 Release
Contents
1. 2 Algae mortality by Algae mortality is a constant rate only affected by zooplankton grazing temperature and is not a function of season e g when zooplankton may be grazing Water levels in a If water levels decrease in a waterbody shaped like a bowl bowl the removal of model layers as the water level decreases will cause the model to bomb if an upstream segment dries up Pipes under high The pipes algorithm does not handle high head dynamic head flow conditions in a pipe 5 Time step limitation The time step for stability in a system model is governed in a complex system by the lowest time step for numerical stability If you model have a very dynamic river with several reservoirs the time step for the river will control This can result in very long run times One can still break apart the model and run the pieces separately using the WTHOUT files to provide boundary conditions for downstream waterbodies Partitioning The partitioning coefficient is currently constant for all anic and inorganic compartments W2 Version 3 1 Release Notes 7 O2LIMIT The concentration at which processes are defined as anaerobic is a step function ON OFF switch A Monod kinetic approach will be used in the next release for smoothing the transition from aerobic to anaerobic conditions Internal weir at a Putting an internal weir at a Dam segment does not affect Dam segment the outflow from the selective withdrawal structure One m2. 8 W m C is now correctly close to 0 3 W m C 17 W2 MISCL Miscellaneous minor bug fixes that prevented the 12 27 02 model from running unless the Epiphyton was greater than 1 18 W2 Reading Solar Minor fixes in reading in SRO 1 21 03 Radiation W2 Version 3 1 Release Notes 6 Fix or Description of Bug Enhancement Date Enhance ment Added 19 W2 River A new feature of the W2 code for a river branch 1 21 03 modeling defined as a branch with Slope not equal to zero is for it not to subtract layers even when KB becomes less than KT A fictitious layer of width 0 01 m is added to keep the model running even when KT is more than KB a reflection twice in the older code Epiphyton The GUI did not properly read in the value of this 2 25 03 Burial rate ii N A few further updates to bug fix 19 were made that 2 28 03 Modeling were overlooked earlier redefined KBMIN BAREA and reordered computation of segment widths GUI Saved File Please read the new GUI document in the full model 3 3 03 Names download file for details Now the GUI makes backup copies of the control file bathymetry files and graph file each time the GUI is used to open these files This allows the user to recover from an error and go back to previous versions of these files Derived The derived variables TP and TN did not account for 3 4 03 Variables TP algae and double counted CBODP and CBODN The and TN derived variables TOP and TON did not a
3. affect the magnitude of H the kinetic constants used in the computation of Epiphyton growth 11 W2 Sediment Conversion of the sediment initial concentration from 10 22 02 Initial g m to g m Since the User s Manual has the input Concentration dimensions as g n consistent with the epiphyton input gy units we needed to convert to the internal units of g m W for the first order sediment model 12 PRE EPI 0 The preprocessor bombed when EPI was set to 0 in the 10 28 02 DE ES meaane ee ee that included a leap year were off by 1 day 14 GUI Epiphyton The Epiphyton burial rate EB did not correctly write 11 6 02 Burial rate out the value into the new control file Also a header label was mislabeled 15 GUI Output format The GUI often wrote blanks to the w2_con npt file 12 27 02 and default when adding new Branches or waterbodies from within output the GUI And often zeroes were written rather than variables default values Many of these were corrected and updated Now if one has an existing w2_con npt file when adjusting the waterbodies or branches the new w2_con npt file has many of the new values added automatically This prevents the time of adding or filling in blanks 16 W2 TSED There was a unit error in the code for the sediment 12 27 02 temperature heat exchange The value of CBHE was not divided by RHOWCP This will mean that much higher default values of CBHE will be necessary The older default value of 7E
4. in the bottom layer it had the impact of a TSEDF 0 only for the bottom layer in a multiple layer system FLUX output Flux output formats were corrected and updated 11 25 03 W2 Version 3 1 Release Notes 10 Fix or Description of Bug Enhancement Enhance ment Added Type 62 W2 Epiphyton Added a line of code to prevent negative concentrations 12 3 03 of epiphyton in the case of no epihpyton but continued mortality New line EPI K LJ MAX EPI K I J 0 0 m 63 W2 DZ in KT The value of DZ was overcorrected for stratification 1 1 04 effects if they were evident between the 1 and 2 layers only Old code DZ KT l MAX DZMIN FRAZDZ AZ KT I EXPAZ DZMIN 1 0 EXPAZ Fixed code DZ KT MAX DZMIN FRAZDZ AZO EXPAZ DZMIN 1 0 EXPAZ W2 Reading in A problem was fixed reading in dynamic extinction 1 5 04 fee cticens for me fest ime mera S l for TSR file were being modeled Examples DTRIC In the control files in the example problems and in the 1 28 04 and variable manual the variable DTRIC header in w2_con npt was Manual used twice This was corrected in the headers and in the User Manual The GUI W2 and the preprocessor read these files correctly Manual explanatory figures were added to the User Manual W2 Suspended Earlier revisions introduced a bug in the suspended 2 5 04 _ sediment solids routine only allowing settling from the bottom aa layer W2 Revised IMP A code revision in the Thomas algorithm f
5. was fixed using a 7 21 03 work around from MS Visual Basic The installation error stated An Error Occurred While Registering the file C WINDOWS S YSTEM32 msado25 tbl If the user chose Ignore the installation was fine This error message has been corrected so that the user no longer needs to choose ignore _ Floating Algae Model Enhancement The model correctly handles algae with negative settling velocities i e floating algae This enhancement marks Version 3 12 50 W2 Derived In the derived variable output for the WDOUT file the 7 24 03 variable output derived variable computation for TDG was incorrectly summed 52 W2 TSR output New output variables have been added to the TSR file 8 1 03 O eo These include SRON ET SHADE DEPTH SURFACE WIDTH 53 W2 Auto stepping The way the time step increases after a time step 8 1 03 E aigon vielan has been reviga stigmy nO Interpolation The interpolation variables QINIC HDIC DTRIC and 8 12 03 WDIC were not editable from the GUI These variables were added finally to the GUI W2 Version 3 1 Release Notes 9 i om a meen PRE or Enhance ment Added GUI Type 55 GUI Calendar The calendar function under the Time Control returned 8 12 03 a Julian date that was 1 day off The older Version 2 convention was used where January 1 at 12 noon was 0 5 This was fixed to the Version 3 convention which sets Jan 1 at 12 noon to Julian day 1 5 56 W2 pH Ioni
6. W2 Version 3 12 Release Notes Please address questions about the code to either Tom Cole Scott Wells US Corps of Engineers Department of Civil and Environmental Engineering Waterways Experiments Station Portland State University 601 634 3283 P O Box 751 FAX 601 634 3129 Portland OR 97207 0751 e mail Thomas M Cole erdc usace army mil 503 725 4276 FAX 503 725 5950 e mail scott cecs pdx edu W2 V3 1 Release The current release of the model includes the following 1 The w2 model and preprocessor executables and source code W2V3 uses dynamic array allocation so the user does not need a FORTRAN compiler and one executable works for all problems and an example application The model source codes and executables are in the file w2v3 zip The executables prew2ve exe and w2v3 exe need to be run in the same directory as the input files npt There is no need to set up additional software to run the code except perhaps the Array Viewer software from COMPAQ see next section below If you do not set up the COMPAQ Array Viewer you will get an error when running the w2 executable w2v3 exe The executables should work with the enclosed example problem Also a generic FORTRAN source code is included that does not have the COMPAQ Array Viewer information 2 The GUI preprocessor A postprocessor is under development but is not ready for release The GUI preprocessor is found in the file gui zip There is a s
7. ater level to agree to water level data for a reservoir There is a short write up Word document in the zip file Registration information will be required for model download in order to alert users to bug fixes and enhancements W2 Version 3 1 Release Notes 2 W2 Model Assistance for Array Viewer V2 to V3 Converter and GUI Interface The following list provides information to assist the model user to set up the program and files bal ed a 1 Array Viewer Setup In order to set up the Array Viewer you must install the Compaq Array Viewer Demo Version 1 6 or higher This is available free from the following web site http www digital com fortran aview getviewer html This web address may change This allows the use of the Dynamic Array Visualizer during runs 2 V2 to V3 converter Currently this version works only with the Workshop release version of Version 2 that included a water balance routine 3 GUI Interface Install this using the setup exe file Also if there is a newer w2control exe file in the zip file on the FTP site than installed during setup use this latest w2control exe file by overwriting the one installed during setup The latest file includes the latest bug fixes W2 Known Issues The following list shows known bugs and issues with the current release of the code these are being addressed in the next release Description Restart Capability This has not been tested in Version 3 1 and may not work
8. c strength equation was 11 25 03 IF SALT_WATER JW S2 1 47E 3 1 9885E 2 TDS K 1 3 8E 4 TDS K D TDS K I But was changed to IF SALT_WATER JW S2 1 47E 3 1 9885E 2 TDS K 1 3 8E 5 TDS K D TDS K I Because of an earlier typo Stewart Rounds 57 W2 pH Equilibrium constant for KW as a function of 11 25 03 Temperature was not in exact agreement with tabulated values The old equation was KW 10 0 5242 39 T1K 35 3944 8 350E 3 T1K 11 8261 LOG10 T1K OH The replacement from Stumm and Morgan 3 edition kw 10 0 283 9710 0 05069842 t1k 13323 00 t1k 102 24447 log10 t1k 1119669 t1k tlk OH Stumm and Morgan 3rd edition SR 11 25 03 This only affected results in a series of tests to at most 3 decimal point accuracy Stewart Rounds Initialized the array GRAV 1125703 59 W2 Sun rise sun Changed to the earlier version of the variable HH This 11 25 03 set may affect sunrise and sunset hour depending on your 1 1 04 longitude Old code HH JW 0 261799 HOUR LOCAL STANDARD 0 0666667 EQTNEW 12 0 Fixed code HH JW 0 261799 HOUR LOCAL STANDARD 0 0666667 EQTNEW 12 0 Occurs in 2 places in the code W2 TSEDF Short wave solar radiation penetration to the bed was 11 25 03 not computed correctly if there were multiple vertical layers The one layer system results were still valid This error did not allow the short wave solar that was computed at the bottom of the vertical layer to be available to the water
9. ccount for algae P and N respectively QIN displayed on the run time display was shown as 3 5 03 ask if the value was not used in the code This was fixed by initializing the QIN array to 0 0 This did not affect any model computations Linkage of In rare cases for estuaries either UHS 1 or DHS 1 3 5 03 model the model did not correctly set up the linkage of the branches proper head conditions The model did not run in this situation Code was reorganized in the section setting up the UP_HEAD and DN_HEAD conditions Derived This is a follow up on Bug 24 There were several 3 11 03 Variables TP small additional errors that were fixed for TP TN and TN and TKN _ TKN River A further bug fix for KBMIN for Bug fix 22 3 12 03 modeling kbmin I 1 recomputed 3 13 03 AKBR An incorrect computation of the surface area for 3 31 03 distributed distributing the DTRs resulted in potentially uneven flows distribution of the distributed inflows for a river segment W2 Version 3 1 Release Notes 7 Enhance ment Added Type 30 W2 Spreadsheet The depths in the spreadsheet output were potentially 3 31 03 IO output incorrect for a river system W 31 W KT The model did not constrain KT to always be more 4 1 03 than 1 This bug caused the model to crash when water na levels were above the grid This implies that when the a water level is above the top of the grid the width of layer 2 is used W 3 Epiphyton The snap
10. etup exe routine that installs the Visual Basic W2 V3 1 Model Preprocessor called W2CONTROL Once installed the GUI preprocessor is able to aid the model user in setting up the Control File and in evaluating and changing the bathymetry of the system This preprocessor does not automatically set up the bathymetry of the system nor does it provide post processing support A lot of effort is required to properly set up the model bathymetry prior to using the Bathymetry editor within W2Control Also note that on the W2 web page there may be updates to the file w2control exe If there are copy this new file over the one installed during the setup program There is no need to run the setup program again Note that there is now a separate PDF user manual for the GUI interface 3 User s Manual On the web page the User s Manual is provided in a zipped PDF file 4 A Version 2 to Version 3 1 Converter W2 Version 3 1 Release Notes 1 This utility in the file v2tov31 zip converts CE QUAL W2 V2 files based on the latest Workshop release Version 2 to V3 1 file formats There is currently not a V3 0 to V3 1 converter 5 Waterbalance Utility The code and executable for this are in the file waterbalance zip This utility is released without support at the present time The purpose of this code is to approximate the waterbalance for a reservoir by computing flows positive and negative that will allow the model predicted w
11. heck was added to the preprocessor to check for this condition See Bug Fix 10 13 Multiple dams into Currently the code will allow one dam inflow to a one downstream downstream branch by a user specified outflow file The reach code though does allow multiple dams inflowing to a common downstream branch if the outflow is specified as a hydraulic structure W2 Version 3 1 Release Notes Description 14 Problems reading file Sometimes the control file or bathymetry file cannot be read properly by the GUI interface This can be a result of the text editor used to produce the file You will find that the problem file s look all messed up in NOTEPAD but look OK in the PFE Editor or in WORD and W2 usually can read them OK Sometimes the following will fix the formatting 1 Copy the file to a UNIX workstation and copy it back 2 Load the file in WORD as a Text file add a space somewhere in the file but don t mess up the file formatting then save it as a Text file W2 V3 1 Bug Fixes Enhancements and User Manual Changes Code W2 Fix or Description of Bug Enhancement Date PRE or Enhancement Added 1 Assorted bugs Problem saving MANN or CHEZY 9 19 02 2 WIND Changed dimension from IMX to NWB 9 30 02 structures on the desktop pointed to w2control exe rather than installation w2control31 exe Also a missing installation file that required the user to type Ignore was fixed For rare cases whe
12. ing or splitting segments segments the segment s in the I O statements SNP PRF SPR TSR and WDO were not updated 2 In some cases the output format for new segment widths defaulted to XXEXX after a merge split This has been fixed to allow proper of digits in the F8 field for layer widths 3 When merging segments water levels friction and orientation angles were only averaged globally for the new segments Now there is interpolation of these onto the new segments W2 Version 3 1 Release Notes 8 1 2 3 6 7 8 9 Description of Bug Enhancement Date Added A minor fix in QSS K DHST prior to this DHS was 6 12 03 used instead of DHST This would cause a problem in waterbody to waterbody connections where the downstream waterbody started to loose segments An internal variable CDSUM referring to the order of 7 8 03 the derived variable TDG total dissolved gas was iia incorrect Some internal code was re written since some variables 7 8 03 did not need to have the code allocate memory Hence memory allocation was reduced W2 Shading For full topographic shading the maximum shading 7 16 03 was set to 0 9 from 1 0 allowing for diffuse radiation 47 W2 Hydraulic If a downstream branch received 2 or more hydraulic 7 16 03 elements element inflows and these were specified as DOWN IO in the upstream branch the flow accounting was not al done correctly GUI GUI An installation error on the GUI
13. n inflows were zero the code used 10 7 02 an incorrect value of TIN that may have affected the computation of temperatures at the boundary segments using a higher order transport scheme Problem There was a problem with saving updates to pipes 10 7 02 saving tribs gates pumps spillways internal weirs and tributaries hydraulic from the GUI Also the installation routine short cut epiphyton groups Width of If the width of the structure was greater than the 10 11 02 structure segment width there was a format error in processing Input output errors were fixed when using the time 10 8 02 series I O Also a bug was corrected when of EPI 0 the error Sloping and When a waterbody included a branch with zero slope 10 17 02 flat slope and others with a fixed slope the grid was incorrectly within 1 WB set up within W2 PRE Setting up grid Bug fix 7 was implemented in the Preprocessor so 10 17 02 that the 2 codes evaluated the grid similarly W2 Version 3 1 Release Notes 5 Fix or Description of Bug Enhancement Enhance ment Added Type PRE Exit Box amp The CVF exit box has been removed so that the 10 22 02 Misc program finishes without a cryptic message that it finished with error code 0 This caused confusion a among users and has been eliminated Also some A false errors and I O format errors were corrected 10 W2 Epiphyton The calculation of the sources and sinks had a 10 22 02 oe dimensional error This may
14. only to be net incident short wave solar W2 Planned Enhancements for 2003 The following list shows planned enhancements for the year 2003 lem O Deepon OS Sediment Diagenesis Complex sediment diagenesis model 2 Simultaneous water Currently water surface is solved branch by branch The level solution new technique will involve solving all water surfaces for Yt the system or waterbody simultaneously Water quality and hydrodynamic effect of macrophytes eriphyton and periphyton viscosity formulations into one Ea k e turbulence model k e turbulence model that will collapse all turbulent eddy W2 Version 3 1 Release Notes 12
15. or the 2 14 04 code implicit vertical shear stress computation No effect on aa model results Minor adjustment of DLT computation to allow 2 14 04 smoother transitions of time steps when a violation 71 W2 RESTART Improved Restart option capability added more 2 14 04 option variables to the restart output W2 Version 3 1 Release Notes 11 Fix or Description of Bug Enhancement Date Enhance ment Added The surface width to compute water loss during 2 23 04 evaporation was based on the width of the surface layer KT It is possible under some conditions for the correct width to be B at KTI which is the width of the water surface Old code EVA VPTG T2 KT D TDEW JW FW B KT 1T DLX D 2 45E9 Fixed code EV D VPTG T2 KT D TDEW JW FW B KTI D DLX D 2 45E9 SW 1 23 04 Waterbody For some cases of waterbody waterbody linkages 6 25 04 waterbody based on grid meshing there was a leak of flow linkage unaccounted for Short wave Computation of theoretical clear sky radiation in the 6 25 04 solar radiation W2 code when using term by term approach resulted in reflection 6 being subtracted A search of the actual reference showed that this theoretical term already included reflection and was the net at the water surface Hence only for theoretical computations of SRO the REFL variable was removed For SRO read in from the met file reflection is still accounted for in the code since this is expected
16. shot output for epiphyton was in units of g m 4 1 03 this was fixed to the correct output units of g m 34 GUI ELTRIB and The GUI interpreted elevations for TRIB placement as 4 15 03 TRIC integers rather than real numbers and did not have the TRIC variable in the GUI even though it read and wrote that variable to and from the control file Bug Finder Larry Brower 2 2 2 35 GUI Top View of The connection of the segments was not correct for 4 15 03 Grid cases of branches in series This error was corrected making the top view of the grid more realistic 2 2 2 GUI A Fix or Description of Bug Enhancement Date 36 W Surface and The surface layer fluxes of gases DO and CO2 and 4 28 03 sedimentation the surface layer sedimentation rates for LPOM terms RPOM algae ISS were computed incorrectly for the special case where a surface layer extended above layer KT 37 W Branch inflows W2 allows users to change the inputs water quality to 4 28 03 each branch This was not implemented properly so if a user had different water quality variables in different inflows the code did not handle it properly 38 W Flux Output The flux output variables in the User Manual states 5 20 03 they are in kg day but they were accidentally output in units of g day This has been fixed so that all flux variables are in kg d 3 Merging 1 In the GUI bathymetry editor when the number of 5 28 03 Splitting segments was changed by merg
17. ust limit selective withdrawal rather than use an internal weir at the dam segment Remember the internal weir works for the right hand face of a model layer W2 multiple file error If the model user accidentally enters duplicate file names check for an input file the w2 executable will bomb because it will try to read the file in more than once The first use of the file will lock its availability for the second instance The W2 error message that comes on the screen traceback error should mention the file name that has problems The W2 preprocessor is supposed to now catch this potential error see bug fix 34 10 Raising level of The preprocessor will say there is an error if the user spillway weir above raises the wer spillway gate water level control or any grid other hydraulic element above the current top of the grid The w2 code will still run properly though But more correctly the model user should increase the DZ of the upper most layer to a value that would eliminate this problem But keep in mind that the segment widths from the top layer then extend upward at that same width 11 Conservation of N P Unless the stoichiometry of the organic matter and algae and C are the same the code will conserve organic matter but not N P or C exactly This will be fixed in the next release 12 Internal weirs The internal weir algorithm does not work when all vertical layers of a segment are blocked by the weir A c
Download Pdf Manuals
Related Search