Home

W2 Version 3 Release Notes - Civil & Environmental Engineering

1. POAER K I PO4ER K D ERR GIHIESEPGUCUJEISEPQET END DO NEW CODE IF EPIPHYTON CALC JW JE then 5 1 18 2008 PO4EG K PO4EG K E POAER K I PO4ER K I ERR K JE EPC K 1 JE EP JE endif W2 Version 3 71 Release Notes 76 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 or Enhance or GUI ment Type ment Added 30 w2 Bottom Continuation of bug fix 24 in such places 1 18 2008 processes as New code IF K KB I THEN xdum BI K 1 BH2 K 1 SW 1 18 08 ELSE xdum BI K 1 BH2 K I 1 0 BI K 1 I BI K I ENDIF SEDAS K I SEDAS K I MAX AS 0 0 ALG K I JA xdum SW 1 18 08 SEDOMS K I poms JW LPOM 1 RPOM K 1 xdum sw 1 18 08 cb 10 22 06 IF K KB I THEN SW 1 18 08 SEDSO S050 ELSE SEDSO sedS JW SED K I K 1 I BH2 K I 1 0 1 1 Endif DO K KT KB I IF K KB I THEN xdum BI K I BH2 K I SW 1 18 08 ELSE xdum BI K I BH2 K I 1 0 BI K 1 1 B1 K 1 ENDIF DO JA 1 NAL SEDASp K 1 SEDASp I MAX AS 0 0 ap ja ALG K I JA xdum SW 1 18 08 END DO DO JE 1 NEP IF EPIPHYTON CALC JW JE LPOMEPp K I I JE EMR K I JE EPC K I JE END DO do jd 1 nbo
2. Other items that have been explored but not implemented in the release version include smarter fetch calculation algorithm Updates to the selective withdrawal algorithm for multiple withdrawals Particle transport algorithm Fish bioenergetics model and fish volitional movement model DIFFERENCES BETWEEN VERSION 3 71 AND VERSION 3 7 There is only one change in the control file between Version 3 7 and 3 71 There is a new option for outlet structures dynamic centerline elevation In the control file there is an ON OFF option after declaring the of structures for each branch EDDY VISC AZC AZSLC FBC E ARODI STRCKLR BOUNDFR TKECAL WB 1 TKE IMP 1 00000 3 9 53500 0 43100 0 00000 0 00000 IMP N STRUC NSTR DYNELEV W2 Version 3 71 Release Notes 21 BRI lom BR2 0 OFF BR3 0 OFF STR INT STRIC STRIC STRIC STRIC STRIC STRIC STRIC STRIC STRIC If these fields are missing the model will assume that DYNELEV OFF DIFFERENCES BETWEEN VERSION 3 7 AND VERSION 3 6 Even though there are some cases where a Version 3 7 executable will run Version 3 6 and Version 3 5 files fine there are updates required to the w2 con npt file that need to be made The preprocessor will catch these errors Control file changes w2 con npt The main changes to the W2 control file are additional flags to turn ON OFF new control file options and the addition of new state variables for water quality BOD N and BOD P for each BO
3. N NNN 1 0 C UO ER C Co C CH EN CH EH EA EN EN C CO CH CH CH EN 00 8 8 8 8 3 43 43 43 43 29 43 43 43 43 43 43 23 3 3 43 43 43 43 23 43 43 3 1000 1000 1000 1000 O CX ODA O CX to C CX C ODO x FMTCD CDMULT 0 0 0 O co C X 00 00 00 00 CDMIN C CX C CO CD EA EA EEN CO Q EA cO I CDMAX 25 50 25 255 25 50 1 5 15 25 St Dis 20 5 145 60 50 5i 20 9 10 10 10 0 e C C C D O C C CDPL 23 24 25 DG QO N PF CAD In Version 3 2 the user has format control of all output variables as well as MULT control see User Manual In Version 3 1 some groups had one but not the other Also in Version 3 2 the groups HNAME CNAME CDNAME were reordered W2 Version 3 71 Release Notes 47 BUG FIXES AND ENHANCEMENTS BETWEEN VERSIONS There have been many updates and bug fixes between Version 3 6 and Version 3 7 that were part of the development of Version 3 7 These have not been documented Since the release of the non beta version of Version 3 7 we have
4. 5 48 W2 V3 6 Bug Fixes Enhancements and User Manual Changes 48 W2 V3 5 Bug Fixes Enhancements and User Manual Changes 61 W2 V3 2 Bug Fixes Enhancements and User Manual Changes 78 W2 Version 3 71 Release Notes 1 W2 V3 7 RELEASE PACKAGE The current release of the model includes the w2 model and preprocessor executables source codes user manuals and example applications The model files are in 1 zipped file v37 zip Contents of this file are organized in the following subdirectories and files Examples Model application examples include DeGray Reservoir Spokane River and Columbia Slough estuary e Executables executables for the preprocessor and the model in this directory were compiled using Intel Fortran 12 1 compiler and have both 32 and 64 bit versions W2ControlGUI The W2Control GUI was compiled using Visual Basic 6 The GUI directory also has an installation routine for W2Control There is a setup exe routine that installs the Visual Basic W2 V3 7 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 proper
5. B TKE KB 1 1 1 0 5 USTARB USTARB USTARB 5 0 H KB 1 JW AT 0 0 DO T T DLT 1 K TRIDIAG AT I VT 1 CT KT EN 0 09 EN Ce DO A A A TKE KB I 1 2 0 0 CT 0 0 2 1 2 DO K KT KB I K I D LT BH1 K I SIG J AZ K 1 1 AVH1 K 1 I CT K I I BB K I S rT VT K I DT K I END DO CALL 1 0 A TKE K KB I KMX TKE I D DO p L TKE K I 1 D DO nter at cell faces K KT KB I 1 pi 1 DZ K I IG J AZ K I AV T K I CT K I 1 2 1 9 2 0 5 AZ K I AZ K 1 I Z K I MAX AZMIN AZ 1 2 MIN AZMAX JW AZ K T DZMIN FRAZDZ AZ EN D DO W2 Version 3 71 Release Notes 70 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 18 W2 Restart Added TKE to restart variables written out 10 5 07 and read in 19 GUI ET The equilibrium temperature option in the 10 9 07 drop down menu was EQT rather than ET Fixed 20 W2 Sediment The SEDIMENT subroutine did not have any 10 15 07 computational mistakes just an error in assigning all array variables to the value at K l This resulted in excessive computational time The fix is shown below OLD sedbr sedb jw sed k i NEW
6. Inflow placement bottom elevation EBTR EBTR JT lt bottom active cell elevation for tributary JT W2 Version 3 71 Release Notes 50 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 12 Pre Additional Additional error checking was added to help debug an 6 22 09 error checking error in the bathymetry file when the problem was in the branch connectivity specifically BS and BE Also a false error was given when the temperature had an isothermal initial condition constituents were OFF and an initial concentration was set to 2 This was fixed 13 Pre Command line In the windows version of the preprocessor the user 9 12 09 processing can now supply a command line argument that sets and working the working directory of the code Hence one does directory not need to copy the preprocessor into every directory In a batch file for example one can execute displayed for the following command windows preW2 ivf exe C scott w2workshop 2009 workshop waterqual problem3 The preprocessor now uses the supplied directory in double quotes as the working directory for all the files The command line argument has one blank space between the end of the executable and the first quote Also the working directory is now displayed at the top of the window Additional checks were also added for checking the grid linkage
7. 0 0 ta aa CMULT 00000 00000 00000 00000 00000 00000 00000 00 00 00 00 0 00000 0 00000 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 0 00000 2 00000 1 00000 1 HFMT 00 21 0 1 0 01 1 0 2 2 0 2 99 0 1 0 1 0 1 0 1 Quy Qul m 0 1 0 1 0 1 0 1 DMULT 00000 1 00000 1 00000 1 00000 1 00000 1 00000 1 00000 1 CMIN 0000 0000 0000 0000 0000 0000 0000 0000 000 000 0000 2000 000 000 000 000 000 000 000 000 000 000 0100 0000 0000 0000 HMIN 0000 0100 0E 06 0000 7 000 E 08 E 08 E 08 E 08 E 08 E 08 E 08 E 08 E 08 E 08 CDMIN 0000 0000 0000 0000 0000 0000 0000 formats CMAX 200 000 200 00 1000 00 5 00000 300 00 3 0000 15 0000 50 000 300 00 5 0000 10 0000 15 0000 2 00000 3 0000 4 00000 3 00000 4 00000 0000 0000 0000 0000 0000 3 0000 5 0000 0 0000 200 000 O QC HMAX 100000 0 10000 0 01000 30 000 1005 00 0 00100 01000 000 000 000 000 000 000 01000 0 0000 O C O CY C O en 0 0 0 0 0 0 0 0 1 CDMAX 3 00000 25 0000 50 0000 25 0000 25 0000 25 0000 5 00000 CPL HPL CO OO OO OO OO OO OO OO OO OO OO OO CO OO OO OO OO OO OO multipliers
8. 1 1 1 1 2 1 and array viewer O vw JO amp GA N F FH N N o Ur GA N P 1 YN ta ta N P vw GA N P 45 Total nitrogen g m 3 Dissolved organic phosphorus mg m 3 Particulate organic phosphorus mg m 3 Total organic phosphorus mg m 3 Total phosphorus mg m 3 Algal production g m 2 day Chlorophyll a mg m 3 Total algae g m 3 Oxygen Gas Saturation Total suspended Solids g m 3 Total Inorganic Suspended Solids g m 3 Carbonaceous Ultimate BOD g m 3 pH CO2 HCO3 An example of the same graph file but for Version 3 2 is shown below 00000 000 00 000 00 000 00 000 00 00000 000 00 00000 00000 00000 00000 00000 00000 00000 00000 00000 o 0000 0000 0000 0000 0000 0000 0000 0000 5 0000 0000 0000 0000 00000 0000 0000 0000 Hydrodynamic constituent and derived constituent names controls Timestep violations NVIOL Horizontal velocity U m s Vertical velocity W m s Temperature T1 lt o gt C Density RHO g m 3 Vertical eddy viscosity AZ m 2 s Velocity shear stress SHEAR 1 s 2 Internal shear ST m 3 s Bottom shear SB m 3 s Longitudinal momentum ADMX m 3 s Longitudinal momentum DM m 3 s Horizontal density gradient HDG m 3 s Vertical momentum ADMZ m 3 s Horiz
9. M B M B W2 HALF ESAT EHS ENEQN ENPR 1 50 0000 40 0000 2 0 00200 TEMP ET1 ET2 ET3 ET4 1 2 00000 5 00000 20 0000 30 0000 STOI EP EN EC ESI 1 0 00500 0 08000 0 45000 0 00000 LDOMDK RDOMDK LRDDK 0 10000 0 00100 0 00100 LPOMDK RPOMDK LRPDK POMS 0 08000 0 00100 0 00100 0 10000 STOIC ORGP ORGN ORGC ORGSI 0 00500 0 08000 0 45000 0 18000 RATE OMT1 OMT2 OMK1 OMK2 4 00000 30 0000 0 10000 0 99000 Version 3 71 Release Notes EK1 EK2 EK3 4 0 10000 0 99000 0 99000 0 10000 ECHLA EPOM 65 0000 0 80000 CBOD KBOD TBOD BOD 1 0 04180 1 01470 BOD 2 0 13020 1 01470 BOD 3 0 04690 1 01470 BOD 4 0 08800 1 01470 BOD 5 0 05000 1 01470 CBOD STOIC BODP BODN BOD 1 0 00500 0 08000 BOD 2 0 00500 0 08000 BOD 3 0 00500 0 08000 BOD 4 0 00500 0 08000 BOD 5 0 00500 0 08000 PHOSPHOR WB 0 00100 0 00000 AMMONIUM NH4R NH4DK WB 0 00100 0 50000 RATE 4 1 2 WB 5 00000 25 0000 NITRATE NO3DK NO3S WB 0 05000 0 00000 RATE NO3T1 NO3T2 WB 5 00000 25 0000 SILICA DSIR PSIS WB 0 10000 0 00000 IRON FER FES WB 0 10000 0 00000 SED CO2 CO2R WB 0 10000 STOICH 1 O2NH4 020M WB 4 57000 1 40000 STOICH 2 O2AR 02AG ALG 1 10000 1 40000 STOICH 3 O2ER 02EG EPI 1 10000 1 40000 02 LIMIT KDO 0 10000 SEDIMENT SEDC SEDPRC WB ON ON SOD RATE 50071 SODT2 WB 4 00000 30 0000 S DEMAND SOD SOD RBOD 00000 00000 00000 00000 00000 ta BODC 450
10. output file format was made No change was made in the tecplot output format DO I CUS JB DS JB WRITE CPL JW A38 9 F10 3 2X CDNAME CDN JD JW CD K I CDN JD JW CDMULT C DN JD JW K KTWB JW KB I cb 6 28 13 end do IWRITE CPL JW A38 9 F10 3 2X CDNAME CDN JD JW CD K I CDN JD JW CDMULT CDN JD JW amp SW 8 12 06 K KTWB JW KB I I CUS JB DS JB CB 1 03 05 7 31 13 31 W2 Read input file A regression error that cropped up when there were 9 or greater than 10 waterbodies has been fixed This had to do with reading in derived and flux variables in the control file DO JD 1 NDC lt 10 READ CON 8 9A8 CDNAME2 JD CDWBC JD JW JW 1 NWB If nwb gt 10 READ CON A8 9 8 8 9 8 CDNAME2 JD CDWBC JD JW JW 1 NWB cb 9 13 12 sw 2 18 13 6 16 13 END DO READ CON DO JF 1 NFL do jf 1 73 Fix this later If nwb lt 10 READ CON A8 9A8 2 JF KFWBC JF JW JW 1 NWB If nwb 10 READ A8 948 8X 9A8 2 KFWBC JF JW JW 1 NWB Lech 9 13 12 sw2 18 13 6 16 13 8 13 13 32 W2 New compiler Upgraded to the Intel XE 13 1 3 198 compiler New W2 executables for 32 bit and 64 bit 8 13 13 33 W2 INIT WL An error was fixed in the initial water level computatio
11. KB I KT B KB I I 0 000001 DX KB I I DXI JW Additional layer subtractions ZMIN JW 1000 0 DO JB BS JW BE JW DO I CUS JB DS JB ZMIN JW 2 Z I END DO END DO SUB LAYER ZMIN JW gt 0 60 H KT JW AND KT lt KTMAX SR 10 17 05 END DO END DO Also done for the initial set up of the branch geometry Upstream active segment and single layer IF SLOPE JB 0 0 THEN DO I US JB 1 DS JB 1 IF KB 1 lt KT THEN lt AND I IZMIN JW SW 10 17 05 0 000001 84 Code W2 PREW2 or GUI Fix or Enhance ment Type Description of Bug Enhancement Date Bug Fixed or Enhance ment Added 17 W2 Shade algorithm No errors just an improvement in computational efficiency Delete this from the SHADING subroutine Set the angles for which topographic shade data are available DO II 1 IANG ANG II II 1 360 0 FLOAT IANG PI 180 0 END DO GAMMA 2 PI IANG and change the 2 occurrences of gamma to gama only in shading subroutine ANG2 I J 1 TOPO I J GAMA SW 10 17 05 TOPOANG TOPO I J ANG2 ANG1 ENDIF END DO IF 200 gt ANG IANG AND 700 lt 2 PI THEN ANG1 AZ00 ANG IANG ANG2 TOPO I 1 TOPO I IANG GAMA lt SW 10 17 05 ADD line tothe module SHADEC MODULE SHADEC PARAMETER IANG
12. SP 8 27 07 Shwet Prakash 15 GUI NZOOP When of zooplankton was set equal to zero 9 17 07 there was an array dimensioning error that caused the writing of the control file to only proceed part way Fixed W2 Version 3 71 Release Notes 67 Code W2 PREW2 or GUI Fix or Enhance ment Type Description of Bug Enhancement Date Bug Fixed or Enhance ment Added 16 W2 Open channel flow Variable passed between subroutines had inconsistent declaration between routines REAL ALLOCATABLE DIMENSION Y D B V CAREA TOPW BELEV Q VOLD YOLD REAL ALLOCATABLE DIMENSION Y B V CAREA BELEV Q VOLD YOLD cb 10 1 07 ALLOCATE Y NN V NN CAREA NN BELEV NN O NN VOLD NN YOLD NN D NN NN ALLOCATE Y NN V NN CAREA NN BELEV NN Q NN VOLD NN YOLD NN B NN cb 10 1 07 DEALLOCATE Y V CAREA TOPW BELEV Q VOLD YOLD D B YT VT VPR YPR TAREA TOPWT RT INDX AL DAA DEALLOCATE Y V CAREA BELEV Q VOLD YOLD B YT VT VPR YPR TAREA TOPWT RT INDX AL DAA cb 10 1 07 10 4 07 17 W2 TKE model The TKE algorithm had several bugs that have been fixed these included making the loop over layers go to KBMIN rather than KB the original code overwrote the boundary conditions when using the Thomas algorithm the orig
13. kept a list of code fixes and enhancements Also we have included below a series of tables with code fixes for Version 3 6 and earlier versions as a reference to earlier versions W2 V3 6 BUG FIXES ENHANCEMENTS AND USER MANUAL CHANGES H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 1 W2 TKE1 model The variable STRICK was incorrectly allocated as 10 11 2008 an INTEGER rather than REAL 2 W2 PIPE Code was streamlined in the subroutine ZBRENT 10 11 2008 where calls were made directly to CDFUNC rather than through the dummy function FUNC 3 W2 Manual ZO The User Manual had 20 in an incorrect line 10 28 2008 the control file w2 con npt The write up and example control file in the User Manual were corrected 4 wW2 Longitudinal The W2 program did not read initial constituent 12 4 2008 profile input concentrations in the longitudinal profile file when CCC was OFF This has been fixed 5 W2 TECPLOT When using TECPLOT output for multiple 1 26 2009 output waterbodies the output format did not allow loading the information into TECPLOT Fixed 6 W2 Epiphyton For entering vertical profile data for periphyton there wasan 5 21 2009 input index error OLD CODE IF VERT_EPIPHYTON JW JE 1 EPIVP K JW JE NEW CODE IF IWERT EPIPHYTON JW JE EPD I JE EPIVP JW JE W2 Versio
14. 11 2 2012 an effort for the error trapping to occur before the code bombs Fixed a couple of regression 11 5 2012 errors as a result of this fix 13 Excel macro Added an Excel macro utility to aid in writing out 11 5 2012 utility input files to CE QUAL W2 W2 Version 3 71 Release Notes H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 14 W2 Withdrawal Fixed an IF test that used the wrong variable in the 11 13 12 subroutine dynamic port allocation algorithm Also added code to allow the code to test for temperatures at the outlet levels specified Deleted line of code is underlined followed by the fix DO J 1 NUMTSPLT REODERING OUTLETS SO THAT HIGHEST ELEVATION STRUCTURE ON TOP ASSUMING 2 SPLIT OUTLETS IF TCNTR J ST THEN IF TSPLTCNTR J ST THEN cb 11 11 12 F ESTR JSTSPLTT J 1 TSPLTJB J ESTR JSTSPLTT J 2 TSPLTJB J THEN JSTSPLT J 1 JSTSPLTT 3 2 JSTSPLT J 2 JSTSPLTT J 1 END IF ELSE IF TCNTR J ELSE IF TSPLTCNTR J cb 11 11 12 F EWD JSTSPLTT J 1 EWD JSTSPLTT J 2 THEN WD THEN WD THEN IF TSPLTJB J JB ST THEN AND TSPLTCNTR J QALL 0 0 DO JJ 1 NOUTS J QALL QALL QSTR JSTSPLT J JJ TSPLTJB J SUM UP ALL THE FLOWS ELR SINA JB DLX DS JB 0 5 DO K KTWB JW KB DS JB
15. DS JJJB EXIT New code IF abs UHS JB DS JJJB EXIT W2 Version 3 71 Release Notes 78 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 4 wW2 Pumpback Pumpback logic was corrected this is legacy 1 25 05 code that will probably be removed from later versions of W2 Old code DO JB 1 NBR IF JB JBP JWBP JW END DO New code DO JW 1 NWB DO JB BS JW BE JW IF JB JBP JWBP JW END DO END DO 5 W2 CPL write Switched order of implied DO loop on CPL 1 25 05 write statement for output of constituents 6 W2 PRF write Changed output format for PRF output for 1 25 05 constituents from f10 2 to e13 6 W2 Version 3 71 Release Notes 79 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 7 W2 Heat balance Added the 1450 and Jackson long wave 1 25 05 radiation equation when air temperatures are below 5C The Swinbank model underpredicts long wave incoming radiation at low air temperatures by as much as 1096 The computation of long wave atmospheric radiation is done using the approach of Swinbank 1963 unless air temperatures are less than 5 C when the Idso and Jackson 1969 formula is used Wells et al 1982 The Swingbank formula for clear sky long wave atmospheric radiati
16. ENDDO go to 28995 ob 3 26 10 3 26 10 24 Prew2 Check on LAT or DOWN Added an enhancement to do a check in case a spillway pipe pump or gate was specified as DOWN In all cases where DOWN is specified the segment that the hydraulic structure originates must be at the end of a branch Additional logic was added to check for this in all the hydraulic structures 3 26 10 25 W2 Manual Light extinction ice Added more text to the section on computation of light extinction and inserted a missing reference Revised an equation for clarity in ICE algorithm and added more explanation on how to estimate HICE 4 13 2010 26 W2 Manual Precipitation input file The units of precipitation are in m s The example precipitation input file was changed to more realistic values 4 14 2010 W2 Version 3 71 Release Notes 57 Code W2 PREW2 or GUI Fix or Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 27 W2 ICE Added code to account for the need to compute long wave radiation in case user chose the equilibrium temperature approach Fixed subscript error in ice melt computation Also made the variable TICE double precision since it is assumed double precision in the call to Surface terms New code IF ICE I THEN TICE DEL J JW 2 0 1 if tair jw ge 5 0 then SW 4 19 10 RA
17. WRITE CWDOC CN JC F8 3 2 2 SW 9 23 13 Changed format from G8 3 to F8 3 to avoid format overflow CWDOC CN JC ADJUSTR CWDOC CN 2C IF QWDO J 0 0 CDWDO CDN 2JD 2W 2 CDWDO CDN JD JW 2 QWDO I WRITE CDWDOC CDN JD 2W F8 3 CDWDO CDN JD JW 2 SW 9 23 13 Changed format from G8 3 to F8 3 to avoid format overflow CDWDOC CDN JD JW ADJUSTR CDWDOC CDN JD JW 45 W2 and Inflow Added csv file format as a new file input format for 7 15 14 PREW2 Tributary flow and temperature files for inflows tributaries and Distributary distributed tributaries Also the shade file is now in and Shade lt gt file format This enhancement includes updates to the preprocessor and W2 codes Also several minor inputs bug fixes were made on the Preprocessor 46 W2 Resuspension A resuspension formula was corrected See the code 7 15 14 of inorganic change below solids HS 0 283 U2 G 0 283 TANH COEF1 TANH COEF2 TANH COEF1 LTS E 2 0 PI U2 G 1 2 TANH COEF3 TANH COEF4 TANH COEF3 TS 2 0 PI sqrt 02 G 1 2 TANH COEF3 TANH COEF4 TANH COEF3 cb 5 9 14 47 W2 Tecplot When the user sets CPL output for Tecplot the output 7 15 14 output format when HABITAC OFF was incorrect This has been fixed W2 Version 3 71 Release Notes 18 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhance
18. then SW 4 18 07 EDAS K I AS JA 0 0 ALG K I JA BH2 K I 1 0 BI K I BI K I BI K 1 1 else k SEDAS K I SEDAS K 1 SEDAS K I MAX AS JA 0 0 ALG K I BI K 1 BH2 K I 1 0 BI K 1 1 B1 K 1 endif b 1 then SW 4 18 07 SEDOMS K I POMS JW LPOM K I RPOM K I BI K I BH2 K I SEDSO POMS JW SED K I BI 1 I BH2 K I else SEDOMS K I POMS JW LPOM K I RPOM K I BI K I 1 1 0 BI K 1 I BI K I SEDSO POMS JW SED K I BI 1 I BH2 K I K 1 I BI K I endif if k X W2 Version 3 71 Release Notes 88 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 24 Wi Algae The logic for negative settling velocities for 8 27 07 algae had an error Old code ASR K 1 JA AS JA ALG K 1 I JA B K41 1 B K I H2 K I ALG K I JA BI K I BH2 K 1 New code ASR K I JA AS ALG K 1 I JA 1 1 BH2 I ALG K I JA BI 1 BH2 K I SP 8 27 07 Shwet Prakash W2 Version 3 71 Release Notes 89
19. RPOMSET CBODDK DOAP DOAR DOEP DOER DOPOM DODOM DOOM DONITR DOCBOD DOREAR DOSED DOSOD TICAG TICEG SEDDK SEDAS SEDLPOM SEDSET SODDK CST ICON TDS Genl Gen2 Gen3 Gen4 Gen5 551 PO4 NH4 NO3 DSI PSI W2 Version 3 71 Release Notes Co OO O O EH E E EA gt CO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO np np rj a n mj mj mj mj mj mj mj o a mj m m m mj m m m m m m a CI Cz Ci Ci Ci Cm Ji 1 C2IWB 00000 00000 00000 00000 00000 00000 00000 03000 01000 30000 00000 00000 C2IWB C2IWB C2IWB C2IWB C2IWB C2IWB C2IWB C2IWB LDOM RDOM LPOM RPOM BOD1 BOD2 BOD3 BOD4 BOD5 ALG1 TIC LK y FU a FU DS Genl Gen2 Gen3 Gen4 Gen5 1551 04 03 DSI PSI z z mg zj LDOM RDOM LPOM RPOM BOD1 BOD2 BOD3 BOD4 BOD5 ALG1 TIC LK D CIN CON H DS 00000 10000 10000 10000 10000 00000 00000 00000 00000 00000 10000 oo ao EN oo CO En ka N 0 a 2 EA aa eco B O nj 1 3 d J 1 1 1 J 1 1 J 1 1 J 1 i d J 1 1 J 1 3 y Oo 1 Ki OO OO OO OO OO OO OO OO OO OO OO OO z D On 2 CPRWBC CPRWBC CPRWBC CPRWBC CPRWB
20. SW 9 5 13 Added zooplankton to TSS computation W2 Version 3 71 Release Notes 13 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 37 W2 Spillway LAT When a spillway was defined with IDSP 0 and LAT a 9 11 2013 tributary was defined incorrectly The new code is shown below IF IDSP JS 0 cb 9 11 13 JTT JTT 1 QTR JTT QSP JS ITR JTT IDSP JS PLACE OTR JTT PDSPC JS DENSITY SPECIFY QTR JTT PDSPC JS SPECIFY IF SPECIFY QTR JTT THEN ELTRT JTT ETDSP JS ELTRB JTT EBDSP JS END IF JBTR JTT JBD end if cb 9 11 13 38 W2 32 bit exe on Recompiled with new settings from Visual Studio 2012 9 11 2013 XP to hopefully run on XP systems with 32 bit OS 39 W2 End Added new close open files in the end simulation 9 25 13 Simulation subroutine This is merely cleaning up the code to be consistent in closing all open files when a Stop is executed This should have no effect on the end user Part of this new code is shown below IF SELECTC ON then SW 9 25 13 New Section on closing files ifile 1949 do jb 1 nbr if nstr jb gt then ifile ifile41 close ifile endif enddo if nwd gt then ifile ifile 1 close ifile endif endif IF DOWNSTREAM_OUTFLOW THEN JFILE 0 DO JWD 1 NIWDO CLOSE WDO JWD 1 CLOSE WDO JWD 2 IF CONSTITUENTS THEN CLO
21. Thanks to Wenwei Xu for pointing this out New code is shown below CONVERT PHIO AND PHINX TO LESS THAN 2 PI SW 2 13 15 DO WHILE PHIO JW 2 PI PHIO JW PHIO JW 2 PI ENDDO DO WHILE PHINX JW 2 PI PHINX JW PHINX JW 2 PI ENDDO IF PHIO JW PHINX JW gt PI THEN PHI JW 1 0 RATIO PHINX JW 2 8 PI RATIO PHIO JW ELSEIF PHIO JW PHINX JW PI THEN WX 2 13 15 PHI JW 1 0 RATIO PHINX JW RATIO PHIO JW 42 0 PI WX 2 13 15 ELSE PHI JW 1 0 RATIO PHINX JW RATIO PHIO JW END IF W2 Version 3 71 Release Notes 19 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 53 WI Restart output Added code to write out a restart file rso opt at the 4 9 15 end of a run if restart output is ON W2 Version 3 71 Release Notes 20 W2 PLANNED ENHANCEMENTS The following list shows planned enhancements H Item Description 1 Sediment Complex sediment diagenesis model Diagenesis 2 Simultaneous Currently water surface is solved branch by branch water level The new technique will involve solving all water solution surfaces for the system or waterbody simultaneously 3 W3 3D version of W2 4 Hypoheric flow Groundwater surface water interaction algorithm 5 Sediment channel Dynamic heat transfer between channel bottom and bottom heating stream algorithm
22. m 3 s Longitudinal momentum ADMX m 3 s Longitudinal momentum DM m 3 s Horizontal density gradient HDG m 3 s Vertical momentum ADMZ m 3 s Horizontal pressure gradient HPG m 3 s Gravity term channel slope GRAV m 3 s Age days Tracer g m 3 Bacteria col 100ml Conductivity mhos W2 Version 3 71 Release Notes Fr 110 1 p D C amp CY Qo O O 0 O0 0o 0 1 F10 F10 F10 F10 F10 1 21 2 43 43 43 43 43 43 43 43 43 43 43 3 3 3 23 3 HMULT HMIN 1 0 1000 1E 6 10 0 997 0 1 08 1 08 08 1 08 1 08 1 08 1 08 1 08 1 08 0 0 EN EN EA C a 1 CMULT CMIN 0 1 0 0 1 0 0 20 000 0 20 000 0 20 000 and array viewer multipliers HMAX HPLTC 1 0 OFF 0 15 0 0 26 0 ON 1005 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 10 0 OFF 0 0 OFF CMAX CPLTC 200 0 OFF 200 0 ON 100 0 OFF 100 0 OFF 100 0 OFF 1 WN Po ta Oo vu n P WN HP 42 Chloride mg l F10 3 0 20 000 100 0 OFF ISS g m 3 F10 3 0 20 000 100 0 OFF Phosphate g m 3 F10 3 1000 0 1 0 500 0 OFF Ammonium g m 3 F10 3 1000 0 0 1000 300 0 OFF Nitrate Nitrite g m 3 F10 3 0 0 1000 5 0 OFF Dissolved silica g m 3
23. 6CBOD ON CINBRC ON OFF OFF OFF ON OFF ON ON ON ON OFF OFF OFF ON ON ON ON ON ON ON ON ON ON CINBRC CINBRC CINBRC CINBRC CINBRC CINBRC W2 Version 3 71 Release Notes 7CBOD ON ON 8CBOD ON ON 9CBOD ON ON 10CBOD ON ON 1CBOD ON 2CBODP 3CBODP 4CBODP ON 5CBODP ON 6CBODP 7CBODP ON 8CBODP 9CBODP 10CBODP ON ON ON 2CBODN ON 3CBODN ON 4CBODN ON 5CBODN ON 6 7CBODN ON 8CBODN ON 9CBODN ON ON ALG1 ON ON ALG2 ON ON ALG3 ON ON DO ON ON TIC ON ON ALK ON ON 2001 OFF OFF LDOM P ON ON RDOM P ON ON LPOM P ON ON RPOM P ON ON LDOM N ON ON RDOM N ON ON LPOM N ON ON RPOM N ON ON CTR CON CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC TDS ON ON AGE OFF OFF TRACER ON ON COL1 ON ON Conduct ON ON Chlorine ON ON ISS1 ON ON 4 ON ON NH4 ON ON NOx ON ON DSi OFF OFF PSi OFF OFF TFe OFF OFF LDOM ON ON RDOM ON ON LPOM ON ON RPOM ON ON 1CBOD ON ON 2CBOD ON ON 3CBOD ON ON 4CBOD ON ON 5CBOD ON ON W2 Version 3 71 Release Notes 6CBOD ON ON 7CBOD ON ON 8CBOD ON ON 9CBOD ON ON 10CBOD ON ON 1CBODP ON ON 2CBODP 3CBODP ON 4CBODP ON 5CBODP ON 6CBODP 7CBODP ON 8CBODP ON 9CBODP 10CBODP ON ON ON 2CBODN ON ON 3CBODN 4CBODN ON 5CBO
24. Aseries of minor memory issues were cleaned up This 9 27 13 memory should have no impacts on current model runs These issues were usually uninitialized memory Code changes made include READ CON KFNAME2 SW 9 27 13 INITIALIZE ENTIRE ARRAY KFWBC i SW 9 27 13 INITIALIZE ENTIRE ARRAY READ CON 8X 9I8 KBWD JW JW 1 NWD TRC SW 9 27 13 INITIALIZATION SINCE ALLOCATION IS TO NTRT READ CON 8X 9A8 TRC JT JT 1 NTR EHSN JE EHSSI JE JE 1 NEPT 1JE 1 NEP SW 9 27 13 READ CON 8X 2F8 0 I8 F8 0 ESAT JE EHS ENEQN JE ENPR JE JE 1 NEPT 1JE 1 NEP SW 9 27 13 READ CON 8X 8F8 0 ET1 JE ET2 JE ET3 JE ETA JE EK1 JE EK2 JE amp EK3 JE JE 1 NEPT 1JE 1 NEP SW 9 27 13 READ CON 8X 6F8 0 EP JE EN JE 2 ESI JE ECHLA JE EPOM JE JE 1 NEPT 2 1 SW 9 27 13 READ CON 8X A8 18 A8 RSOC NRSO RSIC RSOD 0 0 SW 9 27 13 INITIALIZE SINCE ALLOCATED AS NOD BUT ONLY NRSO USED READ CON 8X 9F8 0 RSOD J J 1 NRSO READ CON 8X 18 F8 0 a8 NDLT DLTMIN DLTINTER DLTD 0 0 SW 9 28 13 INITIALIZE ARRAY TO NOD SINCE ONLY NDLT ASSIGNED READ CON 8X 9F8 0 DLTD 2 J 1 NDLT SINKC 1 NSTR JB JB SINKCT 1 NSTR JB JB POINT SINK 1 NSTR JB JB SINKC 1 NSTR JB JB POINT SW 9 27 13 END DO POINT SINK S
25. W2 error message that comes on the screen traceback error should mention the file name that has problems The W2 preprocessor should catch this potential error 7 Raising level of The preprocessor will say there is an error if the user raises the weir spillway spillway weir above gate water level control or any other hydraulic element above the current grid 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 8 Internal weirs The internal weir algorithm does not work when all vertical layers of a segment are blocked by the weir 9 Multiple dams into one Currently the code will allow one dam inflow to a downstream branch by a downstream reach user specified outflow file The code though does allow multiple dams inflowing to a common downstream branch if the outflow is specified as a hydraulic structure 10 Problems reading file Sometimes the control file or bathymetry file cannot be read properly by the in GUI 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 Copy th
26. layers as the water level decreases will cause the model to bomb if an upstream segment dries up 2 Pipes under high head The pipes algorithm does not handle well high head high speed dynamic flow conditions in a pipe 3 Time step limitation in The time step for stability in a system model is governed by the lowest time a complex system step for numerical stability If you have a very dynamic river with several model 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 WDOUT files to provide boundary conditions for downstream waterbodies 4 Partitioning The partitioning coefficient for sorption is currently constant for all organic and inorganic compartments 5 Internal weir at a Dam Putting an internal weir at a Dam segment does not affect the outflow from segment the selective withdrawal structure One must 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 Version 3 71 Release Notes Item Description 6 W2 multiple file error If the model user accidentally enters duplicate file names for an input file the check 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
27. 0 1000 SH 0 2000 0 1000 0 1000 0 1000 0 1000 0 1000 0 0100 0 0100 0 0100 0 0100 0 0100 AANA OF 20 21 22 23 and array viewer multipliers HMAX HPLTC 1 0 OFF 0 15 ON 0 0 OFF 26 0 ON 1005 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 0 0 OFF 10 0 OFF 0 0 OFF CMAX CPLTC 200 0 200 0 ON 100 0 OFF 100 0 OFF 100 0 100 0 ORE 100 0 OFF 500 0 OFF 300 0 OFF 5 0 OFF 10 0 OFF 15 0 ORE 2 0 ORE E350 OFF 4 0 OFF 3 0 OFF 4 0 OFF 3 0 OFF 3 0 OFF 3 0 OFF 3 0 OFF 3 0 OFF 300 OF WN FE ta ta N P 300442 GA M OP WN N N 46 Algae g m 3 Dissolved oxygen g m 3 Inorganic carbon g m 3 Alkalinity g m 3 Dissolved organic carbon g m 3 Particulate organic carbon g m 3 Total organic carbon g m 3 Dissolved organic nitrogen g m 3 Particulate organic nitrogen g m 3 Total organic nitrogen g m 3 Total Kheldahl Nitrogen g m 3 Total nitrogen g m 3 Dissolved organic phosphorus mg m 3 Particulate organic phosphorus mg m 3 Total organic phosphorus mg m 3 Total phosphorus mg m 3 Algal production g m 2 day Chlorophyll a mg m 3 Total algae g m 3 Oxygen Gas Saturation Total suspended Solids g m 3 Total Inorganic Suspended Solids g m 3 Carbonaceous Ultimate BOD g m 3 pH CO2 HCO3
28. 00 45000 45000 45000 45000 oc NHAK1 NHAK2 0 10000 0 99000 NO3K1 NO3K2 0 10000 0 99000 PSIDK PARTSI 0 30000 0 20000 SEDCI SEDK 0 00000 0 10000 SODK1 SODK2 0 10000 0 99000 SOD SOD W2 Version 3 71 Release Notes FSOD FSED 0 1 1 00000 1 00000 SOD SOD SOD SOD SOD oo uc c 6 OY 9 O c O REAERATION EQN WB1 LAKE 6 COEF1 COEF2 Lines removed from the V3 2 control file These are a result of eliminating the pumpback and line printer settings Here is the part of the V3 2 control file that was deleted DST TRIB DTRC Br 1 ON Br 2 ON Br 3 OFF Br 4 OFF Br 5 OFF COEF4 O c 2 O O 9 HYD PRINT HPRWBC HPRWBC HPRWBC HPRWBC HPRWBC HPRWBC HPRWBC HPRWBC HPRWBC NVIOL OFF OFF U ON ON Graph npt file changes These changes are a result of the new state variables in W2 and are highlighted below Hydrodynamic constituent and derived constituent names formats controls Timestep violations NVIOL Horizontal velocity U m s Vertical velocity W m s Temperature T1 lt o gt C Density RHO g m 3 Vertical eddy viscosity AZ m 2 s Velocity shear stress SHEAR 1 s 2 Internal shear ST m 3 s Bottom shear SB
29. 10 3 1000 0 1 0 5 0 OFF Total phosphorus mg m 3 F10 3 1000 0 1 0 20 0 OFF Algal production g m 2 day F10 3 0 1 0 5 0 OFF Chlorophyll a mg m 3 F10 3 0 5 145 0 OFF Total algae g m 3 F10 3 0 1 0 60 0 OFF Oxygen Gas Saturation F10 3 0 1 0 50 0 OFF Total suspended Solids g m 3 F10 3 0 1 0 5 0 OFF Total Inorganic Suspended Solids g m 3 F10 3 0 2150 20 0 OFF Carbonaceous Ultimate BOD g m 3 F10 3 0 5 0 9 0 OFF pH F10 3 0 1 0 10 0 OFF CO2 F10 3 0 1 0 10 0 OFF HCO3 F10 3 0 1 0 10 0 OFF 10 3 0 0 0 0 0 0 W2 Version 3 71 Release Notes 3004 N NN GA 5 Q ui YNNN V DIFFERENCES BETWEEN VERSION 3 1 AND VERSION 3 2 There are minor differences 2 input files between the 2 versions w2 con npt and the graph npt file All other files are the same between the 2 versions w2 con npt The only section where there is a slight difference in the control file is in the section where the inorganic suspended solids group settling velocities are entered In Version 3 1 this section looks like this ALG EX EXA EXA EXA EXA EXA EXA 0 10000 GENERIC CGQ10 CGODK CG1DK CGS cG 1 0 00000 1 0000 0 00000 0 00000 CG 2 0 00000 0 00000 0 00000 0 00000 CG 3 1 04000 0 00000 0 50000 0 00000 4 0 00000 0 00000 0 00000 0 00000
30. 14 W2 H of The model user can now control the of physical 9 12 09 processors processors the model uses At this point dual processor model runs have shown an improvement of about 20 over a single processor But QUAD processors usually are slower It is recommended that NPROC be set to 2 in the control file The user can experiment on his her own system If this is not set by the user or is left blank the model still runs but sets it to 2 processors GRID NWB NBR IMX KMX NPROC CLOSEC 1 1 23 22 2 ON W2 Version 3 71 Release Notes 51 Code W2 PREW2 or GUI Fix or Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 15 W2 Command line processing for windows In the windows version of the w2 model the user can now supply a command line argument that sets the working directory of the code Hence one does not need to copy the model executable into every directory In a batch file for example one can execute the following command W2 ivf exe C scott w2workshop 2009 workshop waterqual problem3 The w2 model now uses the supplied directory in double quotes as the working directory for all the files The command line argument has one blank space between the end of the executable and the first quote The working directory is displayed in a text box in the window 9 12 09 16 W2 W2 window closed at end of successful execu
31. 18 REAL PARAMETER GAMA 3 1415926 2 REAL IANG SW 10 17 05 REAL DIMENSI ON IANG ANG lt SW 10 17 05 REAL ALLOCATABLE DIMENSION A00 DECL HH TTLB TTRB e LLB CLRB lt SW 10 17 05 REAL ALLOCATABLE DIMENSION SRLB1 5 1 SRLB2 SRRB2 SRFJD1 SRFJD2 SHADEI REAL ALLOCATABLE DIMENSION TOPO LOGICAL ALLOCATABLE DIMENSION DYNAMIC_SHADE DATA ANG 0 00000 0 34907 0 69813 04720 1 39626 1 74533 2 09440 44346 amp 79253 3 14159 3 49066 3 83972 18879 4 53786 4 88692 5 23599 58505 5 93412 lt 5010 17 05 END MODULE SHADEC Delete allocation statement for ang Oi N N ALLOCATE SRLB1 IMX SRRB1 IMX SRLB2 IMX 5 RRB2 IMX SRFJD1 IMX SHADEI IMX SRFJD2 IMX ALLOCATE TOPO IMX IANG lt SW10 17 05 ALLOCATE QSW KMX NWDT CTR NCT NTRT PRWBC NHY NWB Delete ang from the deallocate statement DEALLOCATE TTLB TTRB CLLB SRLB1 2 SRRB1 SRLB2 SRRB2 SRFJD1 SHADEI SRFJD2 TOPO QSW CTR SW 10 17 05 fan 10 17 05 W2V rsion 3 71 Release Notes 85 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 18 wW2 Epiphyton Several changes were made that corrected 5 26 06 algorithm errors in shallow systems where adding and subtracting layers did not reinitialize macrophyte layers w
32. 5 0 00000 0 00000 0 00000 0 00000 S SOLIDS SSS 5585 5585 555 555 555 555 555 555 1 50000 ALGAL RATE AG AR AE AM AS AHSP AHSN AHSSI ASAT ALG1 2 00000 0 12000 0 02000 0 05000 0 04000 0 00500 0 00500 0 00000 50 0000 In Version 3 2 there is now a sediment resuspension capability for wind driven resuspension along the shores of lakes and reservoirs The Version 3 2 control file has the following lines in this same section of the control file ALG EX EXA EXA EXA EXA EXA EXA 0 10000 GENERIC CGQ1O CGODK CGIDK CGS cG 1 0 00000 1 0000 0 00000 0 00000 2 0 00000 0 00000 0 00000 0 00000 3 1 04000 0 00000 0 50000 0 00000 4 0 00000 0 00000 0 00000 0 00000 5 0 00000 0 00000 0 00000 0 00000 S SOLIDS 555 SEDRC TAUCR 5510 1 50000 OFF 0 00 ALGAL RATE AG AR AE AM AS AHSP AHSN AHSSI ASAT ALG1 2 00000 0 12000 0 02000 0 05000 0 04000 0 00500 0 00500 0 00000 50 0000 For Version 3 2 SSS is the settling velocity for particle group 1 SEDRC is the control which turns ON or OFF sediment resuspension and TAUCR is the critical shear stress at which resuspension occurs For Version 3 2 each line represents 1 SS group while in Version 3 1 each group settling velocity is in the next 8 columns moving across the page W2 Version 3 71 Release Notes 44 graph npt The graph file controls output formatting and the graphing parameters used in Array Viewer only for the PC platform The files have been rearranged significantly A Version 3 1 g
33. 5 21 2009 DA Deficit upstream DB DA 1 0 0 38 AGASSP N BGASSP N CGASSP 1 0 0 11 CGASSP N 1 0 0 046 1 DB deficit downstream SAT DB 5 21 2009 W2 Version 3 71 Release Notes 49 H Code W2 Fix or Description of Bug Enhancement Date Bug PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 10 w2 Order of flux The order of flux parameters in the User 6 2 2009 t Manual and output were incorrect The parameters control file has them in this order RPOMSET CBODDK DOPOM DODOM DOOM whereas the code assumed they were in this order RPOMSET CBODDK DOPOM DODOM DOOM This has been corrected The User Manual and control file order is now reflected in the W2 code 11 Pre False errors The preprocessor sometimes gave false errors in the 6 18 09 for inflow pre err for tributary internal weirs pipes and other location hydraulic features saying that the pipe or tributary was below the elevation of the bottom of the segment The W2 model ran fine even with this error message given in the preprocessor This has been fixed Example of OLD CODE IF EBTR JT lt EL KB ITR JT 1 ITR JT THEN CALL ERRORS WRITE ERR Inflow placement bottom elevation EBTR EBTR JT bottom active cell elevation for tributary JT New CODE IF EBTR JT EL 1 THEN CALL ERRORS WRITE ERR
34. B Layer Sub Seg eliminate parentheses which caused a sign error LE NOT TRAPEZOIDAL JW BI KT IU 1 B KTI IU H1 KT IU 1 H KT JdW Z IU 1 BH1 KT IU 1 B KTI IU 1 IU 1 EL KT IU 1 EL KTI IU 1 1 IU 1 Z 10 1 COSA JB COSA JB 2 2 6 SR 10 17 05 IF KT gt KB IU 1 BH1 KT IU 1 IU 1 H1 KT IU 1 DO K KTI IU 1 1 KT BH1 KT 10 1 KT IU 1 BH1 K IU 1 END DO ELSE 10 17 05 W2 Version 3 71 Release Notes 83 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 16 W2 SUB layer for Layer SUB improve model running in 10 17 05 shallow shallow segments systems W2 Version 3 71 Release Notes Water surface minimum thickness DO JW 1 NWB KT KTWB JW ZMIN JW 1000 0 KTMAX 2 lt SR 10 17 05 DO JB BS JW BE JW DO I CUS JB DS JB IF KB I gt KTMAX lt SR 10 17 05 IF Z I gt ZMIN JW THEN IZMIN JW I JBIZ JB END IF ZMIN JW MAX ZMIN JW Z I END DO END DO ADD LAYER ZMIN JW 0 85 H KT 1 JW AND KT 2 SUB_LAYER ZMIN JW gt 0 60 H KT JW AND KT lt KTMAX lt cec SR 10 17 05 Upstream active segment IUT US JB SLOPE JB 0 0 THEN DO I US JB 1 DS JB 1 IF KB I KT THEN lt SE SR 10 17 05
35. C Pi 1 24 28 28 0 27 0 0 5 230 0 0 065 0 1 DOWN where DYNPIPE controls whether the pipe is controlled by time series of an ON OFF or partially open gate 6 Dynamic pump PUMPS 1 IUPU IDPU EPU STRTPU ENDPU EONPU EOFFPU QPU WTHLC DYNPUMP 111 0 440 1 00 366 441 0 435 0 1 0 DOWN ON where DYNPUMP controls the EPU EONPU EOFFPU and QPU over time by reading in a time series file 7 INIT CND TEMPI WTYPEC GRIDO WB 1 1 0000 0 00000 FRESH RECT where GRIDC controls whether the grid is interpreted as rectangular in depth or trapezoidal 8 CST ACTIVE Note that this change only appears if NBOD gt 0 CST ACTIVE CAC TDS ON Genl ON Gen2 OFF Gen3 OFF Gen4 OFF Gen5 OFF 1551 ON 4 ON NH4 ON NO3 ON DSI OFF PST OFF FE OFF LDOM ON RDOM ON LPOM ON RPOM ON CBOD ON 2CBOD ON 3CBOD ON 4CBOD ON 5CBOD ON 6CBOD ON 7CBOD ON 8CBOD ON 9CBOD ON 10CBOD ON 1CBODP ON 2CBODP ON W2 Version 3 71 Release Notes ALG1 ON ALG2 ON ALG3 ON DO ON TIC ON ALK ON 7001 LDOM P ON RDOM P ON LPOM P ON RPOM P ON LDOM N ON RDOM N ON LPOM N ON RPOM N ON 9 CSTICON CST PRIN CIN CON CTR CON CDT CON and CPR CON CST ICON C2IWB C2IWB DS 0 GE TRACER COL1 Conduct Chlorine ISS1 04 Ox DSi PSi Fe LDOM RDOM LPOM RPOM 1CBOD 2CBO CBOD 4CBOD O O 22 2 0 e H 6CBOD 7CBOD 8CBOD oo co oor PP P 0 us 1n eS
36. C CPRWBC CPRWBC CPRWBC CINBRC CINBRC CINBRC CINBRC CINBRC CINBRC CINBRC CINBRC W2 Version 3 71 Release Notes Genl OFF Gen2 ON Gen3 ON Gen4 ON Gen5 ON 1551 ON 4 ON NH4 ON NO3 ON DSI OFF PSI OFF FE OFF LDOM ON RDOM ON LPOM ON RPOM ON BOD1 ON BOD2 ON BOD3 ON BOD4 ON BOD5 ON ALG1 ON DO ON TIC ON ALK ON 701 LDOM P OFF RDOM P OFF LPOM P OFF RPOM P op LDOM N OFF RDOM N OFF LPOM N OFF RPOM N OFF CTR CON TDS ON ON Gen2 ON ON Gen3 ON ON Gen4 ON ON Gen5 ON ON ISS1 ON ON PO4 ON ON NH4 ON ON NO3 ON ON DSI OFF OFF PSI OFF OFF FE OFF OFF LDOM ON ON RDOM ON ON LPOM ON ON RPOM ON ON BOD1 ON ON BOD2 ON ON BOD3 ON ON BOD4 ON ON BOD5 ON ON ALG1 ON ON DO ON ON TIC ON ON ALK ON ON 201 OFF CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC CTRTRC W2 Version 3 71 Release Notes CDT CON CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC TDS ON Genl OFF Gen2 ON Gen3 ON Gen4 ON Gen5 ON 1551 ON 4 ON NH4 ON NO3 ON DSI OFF PSI OFF FE OFF LDOM ON RDOM ON LPOM ON RPOM ON BOD1 ON BOD2 ON BOD3 ON BOD4 ON BOD5 ON ALG1 ON DO ON TIC ON ALK ON 2001 OFF LDOM P OFF RDOM P OFF LPOM P OFF RPOM P OFF LDOM OFF RDOM N OFF LPOMN OFF RPOM N CPR CON CPRBRC CPRBRC CPRBRC CPRBRC CPRBRC CPRBRC CPRBRC CPRBRC TDS ON Genl OFF Gen2 ON Gen3 ON Gen4 ON Gen5 ON 1551 ON 4 ON NH4 O
37. D K I NO3 K I NO3SUW NO3TRM K I BI K I BH2 K else NOS3SED K I NO3 K 1 NO3S JW NO3TRM K I BI K I Bl K 1 1 BH2 K 1 endif New code added in sediment routine if k kb 1 then SW 4 18 07 SODD K I 50 2 I SODTRM I K I else SODD K I SOD I BH2 I SODTRM K 1 BI K I BI K 1 1 Endif New code added in suspended solids routine if k kb i then SSR EPSILON DLX I BI K I VOL K 1 else SSR EPSILON DLX I BI K I K 1 1 VOL K 1 Endif W2 Version 3 71 Release Notes 63 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 6 W2 see above New code added for heat flux to channel 4 18 07 bottom if kt kb i then SW 4 18 07 SROSED SROOUT TSEDF JW else SROSED SROOUT 1 0 BI KT 1 I BI KT I TSEDF JW Endif if k kb i then SW 4 18 07 TFLUX CBHE JW RHOWCP TSED JW T2 K I BI K I DLX I else TFLUX JW RHOWCP TSED JW T2 K 1 BI K I BI K 1 I DLX I endif New code added for sediment subroutine if k 1 then SW 4 18 07 EDAS K I AS JA 0 0 ALG K I JA BH2 K I 1 0 BI K I BI K I BI K 1 1 else k SEDAS K I SEDAS K 1 SEDAS K I MAX AS JA 0 0 A
38. D group Below is a list of changes in the control file with the card image header for each line changed highlighted options are new in V3 7 Descriptions of these new features are in the W2 User s Manual 1 MISCELL MISCELL NDAY SELECTC HABTATC ENVIRPC AERATEC INITUWL 100 OFF ON om OFF Five new variables SELECTC HABITATC ENVIRPC AERATEC and INITUWL are 5 new control variables that turn ON OFF the use of automatic selective withdrawal fish habitat volumes environmental performance criteria artificial aeration and the initial water surface and velocity computations respectively If using an old Version 3 6 control file all of these would default to OFF if they were left blank Also the model preprocessor would flag these are missing variables 2 DLTCON DLT CON NDT DLTMIN DLTINTR 1 1 00000 OFF where DLTINTR is a control for interpolating the the time step DLTMAX and DLTF rather than use as a step function 3 BRANCHG BRANCH G US DS UHS DHS NLMIN SLOPE SLOPEC Br 1 2 59 0 0 0 0 1 0 0 0 0 where SLOPEC is the hydraulic equivalent slope for a river channel that affects the momentum equation 4 GATEWEIR GATE WEIR 1 GTB1 GTA2 GTB2 DYNVAR GTIC Gatel 1 00000 1 50000 1 00000 1 50000 FLOW ON W2 Version 3 71 Release Notes 22 where GTIC is an interpolation control for the specified DYNVAR for the GATE WEIR 5 Dynamic pipe PIPES IUPI IDPI EUPI EDPI WPI DLXPI LATPI
39. DN ON 6 ON 7CBODN ON 8CBODN ON 9CBODN ON ON ALG2 ON ON ALG3 ON ON DO ON ON TIC ON ON ALK ON ON 2001 OFF OFF LDOM P ON ON RDOM P ON ON LPOM P ON ON RPOM P ON ON LDOM N ON ON RDOM N ON ON LPOM N ON ON RPOM N ON ON CDT CON CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC CDTBRC TDS ON ON AGE OFF OFF TRACER ON ON COL1 ON ON Conduct ON ON Chlorine ON ON ISS1 ON ON 4 ON ON NHA ON ON NOx ON ON DSi OFF OFF PSi OFF OFF TFe OFF OFF LDOM ON ON RDOM ON ON LPOM ON ON RPOM ON ON 1CBOD ON ON 2CBOD ON ON 3CBOD ON ON 4CBOD ON ON W2 Version 3 71 Release Notes 28 5CBOD ON ON 6CBOD ON ON 7CBOD ON ON 8CBOD ON ON 9CBOD ON ON 10CBOD ON ON 1CBOD ON 2CBODP ON 3CBODP ON 4CBODP 5CBODP ON 6CBODP ON ON TCBODP 8CBODP 9CBODP ON 10CBODP 1CBODN 2CBODN ON 3CBODN ON 4CBODN 5CBODN ON 6CBODN 7CBODN ON ON 8CBODN 9CBODN ON ON ALG1 ON ON ALG2 ON ON ALG3 ON ON DO ON ON TIC ON ON ALK ON ON 2001 OFF OFF LDOM P ON ON RDOM P ON ON LPOM P ON ON RPOM P ON ON LDOM N ON ON RDOM N ON ON LPOM N ON ON RPOM N ON ON CPR CON CPRBRC CPRBRC CPRBRC CPRBRC CPRBRC TDS ON ON AGE OFF OFF TRACER ON ON COL1 ON ON Conduct ON ON Chlorine ON ON ISS1 ON ON 4 ON O
40. EL KB IUT IUT EL KB UHS JB UHS JB THEN IF EL KB IU IU EL KB UHS JB UHS JB THEN cb 6 12 07 KB IUT 1 E B IU 1 KB IU ELSE DO K KT KB IUT IF EL KB UHS JB UHS JB gt EL K IUT THEN KB IUT 1 K EXIT DO K KT KB IU cb 6 12 07 IF EL KB UHS JB UHS JB gt EL K IU THEN cb 6 12 07 B IU 1 EXIT cb 6 12 07 END IF W2 Version 3 71 Release Notes 66 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 11 W2 CBOD settling The CBOD settling rate earlier was not 7 23 07 converted from m d in the control file to m s in the code Added code cbods cbods day cb 7 23 07 12 W2 TSR output The surface width was not correctly being 7 26 07 output Changed BI KT to BI KTWB JW FIX BI JW I SHADE 1 1 ADJU STR C2CH JAC JAC 1 amp CB 7 26 07 13 PREW2 Pumps The pump control for DOWN or LAT was not 8 14 07 being checked properly also a check on IUPUC was incorrect Fixed 14 W2 Algae The logic for negative settling velocities for 8 27 07 algae had an error Old code ASR K I JA AS JA ALG K 1 I JA B K 1 I B K I H2 K 1 ALG K I JA I BH2 K I New code ASR K I JA AS JA ALG K 1 I JA BI K 1 1 BH2 I ALG K I JA BI K I BH2 K I
41. F10 3 0 1 0 10 0 OFF Particulate silica g m 3 F10 3 0 0 2000 15 0 OFF Total iron g m 3 F10 3 0 0 1000 2 0 OFF Labile DOM g m 3 F10 3 0 0 1000 3 0 OFF Refractory DOM g m 3 F10 3 0 0 1000 4 0 OFF Labile POM g m 3 F10 3 0 0 1000 3 0 OFF Refractory POM g m 3 F10 3 0 0 1000 4 0 OFF CBOD1 g m 3 F10 3 0 0 0100 340 OFF CBOD2 g m 3 F10 3 0 0 0100 3510 OFF CBOD3 g m 3 F10 3 0 0 0100 3 0 OFF CBOD4 g m 3 F10 3 0 0 0100 3 0 OFF CBOD5 g m 3 F10 3 0 0 0100 3 0 Algae g m 3 F10 3 0 0 0100 3 0 OFF Dissolved oxygen g m 3 F10 3 0 0 0100 1 0 OFF Inorganic carbon g m 3 F10 3 0 0 0100 340 OFF Alkalinity g m 3 F10 3 0 0 0100 3 0 OFF abu caf ioe a diu CDNAME FMTCD CDMULT CDMIN CDMAX CDPLTC Dissolved organic carbon g m 3 F10 3 0 1 0 25 0 Particulate organic carbon g m 3 F10 3 0 1 0 50 0 OFF Total organic carbon g m 3 F10 3 0 sl H 25410 OFF Dissolved organic nitrogen g m 3 F10 3 0 1 0 25 0 OFF Particulate organic nitrogen g m 3 F10 3 0 1 0 25 0 OFF Total organic nitrogen g m 3 F10 3 0 1 0 50 0 Total Kheldahl Nitrogen g m 3 F10 3 0 1 0 15 0 OFF Total nitrogen g m 3 F10 3 0 1 0 15 0 OFF Dissolved organic phosphorus mg m 3 F10 3 1000 0 1 0 2540 OFF Particulate organic phosphorus mg m 3 F10 3 1000 0 lt 1 0 1 0 Total organic phosphorus mg m 3 F
42. IF EL K DS JB ELR lt ESTR JSTSPLT J JJ TSPLTIB J ISW 10 17 01 EXIT END DO KSTR K 1 KSTRSPLT JJ MIN KSTR KB DS JB ENDDO DO JJ 1 NOUTS J cb 11 11 12 dividing total flow between outlets for temperature test OSTR JSTSPLT J JJ TSPLTJB J qall real nouts j ENDDO W2 Version 3 71 Release Notes H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 15 wW2 Reading in In case a user does not enter the units in 12 3 2012 names of WQ graph npt the code improperly parses the WQ variables variable name In this case the output name is a blank To avoid this issue extra code was added to preserve the variable name even if no units were added to the graph npt list L1 SCAN CNAME JC 2 IF L1 2 L1 43 SW 12 3 2012 Implies no comma found L2 SCAN CNAME 2C L1 43 L1 IF L2 gt 43 L2 43 SW 12 3 2012 CUNIT JC CNAME JC L1 L2 CNAME1 JC CNAME JC 1 L1 3 CNAME3 JC CNAME1 JC DO WHILE L3 lt L1 3 16 PREW2 SEDS and The variable names were switched in reading the 12 30 12 SEDK control file in the preprocessor perhaps leading to incorrect warnings errors being tagged The proper order was restored READ A8 8X 2A8 6F8 0 A8 ERR 400 AID SEDC JW PRNSC JW SEDCI JW seds jw SEDDK FSOD JW amp FSED
43. INKC POINT COLDEP ELWS I COLB MACT J KT I MACT J KT 1 1 IF MACROPHYTE_ON MACT J KT I MACT J KT 1 1 SW 9 28 13 SDKV US JB DS JB SDK JW SDKV US 2B 1 DS JB 41 SDK JW SW 9 28 13 W2 Version 3 71 Release Notes 16 H Code W2 Fix or Description of Bug Enhancement Date Bug PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 42 W2 CPL output Code was added to eliminate writing out the habitat 9 28 13 index to the CPL file for Tecplot when HABITATC is OFF IF I DS JB 1 THEN IF HABTATC ON THEN WRITE CPL JW 9999 X1 I ELWS I U K I W K I T1 K I RHO K I HAB K I C2 K I CN 2C JC 1 NAC ELSE WRITE CPL JW 9999 X1 I ELWS I U K I W K I T1 K I RHO K I C2 K I CN 2C 1 NAC ENDIF ELSE XDUM 99 0 WRITE CPL JW 9999 X1 1 ELWS I XDUM XDUM XDUM XDUM 22 1 ENDIF DO K KTWB JW KMX 1 IF I DS JB 1 AND K lt KB I THEN IF HABTATC ON THEN WRITE CPL JW 9999 X1 I ELWS I DEPTHM K I U K I W K I T1 K I RHO K I I C2 K I CN 2C JC 1 NAC ELSE WRITE CPL 2W 9999 X1 I ELWS I DEPTHM K I U K I W K I T1 K I RHO K I C2 K I CN 2C 2C 1 NAC ENDIF IF K KB I THEN IF HABTATC ON THEN WRITE CPL JW 9999 X1 I ELWS I DEPTHB K I U K I W K I T1 K I
44. J BI K I 2 0 H1 KE I DLX I else J EPD K I J BI K I BI K 1 1 2 0 H1 K I DLX I endif W2 Version 3 71 Release Notes 73 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 25 W2 CBODS If the user defined particulate CBOD that 1 18 08 settles to the bottom and had SED turned ON the conversion from oxygen to organic matter was missing in the accumulation on the channel bottom or sides OLD do jd 1 nbod SEDcb K I E SEDcb K I cbods 3d 0 0 cbod K I Jd BI K I BH2 K I 1 0 BI K 1 1 BI K I end do NEW do jd 1 nbod SEDcb K I E SEDcb K I cbods 3d 0 0 d K I Jd O20M JW 1 BH2 I 1 0 BI K 1 I BI K I 1 16 08 end do 26 W2 SEDBR Eliminated a redundant definition of SEDBR in 1 18 08 the Sediment routine since it is already defined in the Kinetic rates subroutine 27 W2 SEDDK The first order sediment decay rate is an 1 18 08 average of the decay rates of all the influxes of organic matter and their respective decay rates There was an error in computing this average decay rate for CBOD treated as particulate Code fix is shown below OLD do jd 1 nbod sedsum sedsum MAX cbods 3 0 0 cbod K I Jd K I BH2 K I 1 0 BI K 1 1 BI K 1 end do NEW d
45. LG K I BI K 1 BH2 K I 1 0 BI K 1 1 B1 K 1 endif b 1 then SW 4 18 07 SEDOMS K I POMS JW LPOM K I RPOM K I BI K I BH2 K I SEDSO POMS JW SED K I BI 1 I BH2 K I else SEDOMS K I POMS JW LPOM K I RPOM K I BI K I 1 1 0 BI K 1 I BI K I SEDSO POMS JW SED K I BI 1 I BH2 K I K 1 I BI K I endif if k X W2 Version 3 71 Release Notes 64 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 or Enhance or GUI ment Type ment Added 7 W2 Zoo plankton Several fixes in the zooplankton routine were 5 21 07 fixes made Many thanks to Dr Kellie Vache Institute for Landscape Ecology and Resources Management ILR Justus Liebig University Giessen Heinrich Buff Ring 26 35392 Giessen Germany for finding these which are documented below DO K KT KB I do jz 1 nzp zgztot 0 0 kv 5 9 2007 do jjz 1 nzp 1 zooss k i jz zmu k i jz zeff jz zrt k i jz zmt k i jz zoo k i jz Zzgz k i jz jjz zoo k i jz omnivorous zooplankton zgztot zgztot zgz k i jz jjz zoo k i jz lkv 5 9 2007 end do zooss k i jz zmu k i jz zeff jz zrt k i jz zmt k i jz zoo k i jz zgztot kv 5 9 2007 end do do jjz 1 nzp tgraze k i jz tgraze k i jz prefz jz jjz zoo k i jjz tgraze k i jz tgraze
46. M RDOM LPOM RPOM BOD1 BOD2 BOD3 BOD4 BOD5 ALG1 Hj nj rj mj nu np nj nup o np nj Uu Uu m m m 1 W2 Version 3 71 Release Notes 33 ul TIC x nj mj d Wei H ei H CDWB 1 1 y o 1 1 o 1 a a J 1 o z J a o z 1 71 l ON 1 q E 2 3 z bd hd i TI 1 1 FU Di tg 1 z 1 J ATOT 1 1 i a J 1 J Tj CBOD T 1 z i Q 02 1 1 rj nj mj mj n nj nj np np m m Uu a m m m m m 1 1 1 CST FLUX CFWB J i 1 1 1 1 J 1 1 1 1 i d J 1 1 1 1 1 i ZS PU oc COE GU p 74 9 J 1 H4POM H4DOM H4OM i 1 Tj a nj nj nj mj nu np mj n np nj m n Cup oq onu nj m m m m 22222222224 W2 Version 3 71 Release Notes CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC CDWBC CFWBC NHASED NH4SOD NO3DEN NO3AG NO3EG NO3SED DSIAG DSIEG DSIPIS DSISED DSISOD DSISET PSIAI PSINET PSID FESET FESED LDOMDK LRDOM RDOMDK LDOMAP LDOMEP LPOMDK LRPOM RPOMDK LPOMAP LPOMEP LPOMSET
47. N NHA ON ON NOx ON ON DSi OFF OFF PSi OFF OFF TFe OFF OFF LDOM ON ON RDOM ON ON LPOM ON ON RPOM ON ON 1CBOD ON ON 2CBOD ON ON 3CBOD ON ON W2 Version 3 71 Release Notes 4CBOD ON ON 5CBOD ON ON 6CBOD ON ON 7CBOD ON ON 8CBOD ON ON 9CBOD ON ON 10CBOD ON ON lCBODP 2CBODP 3CBODP ON ACBODP ON 5CBODP ON 6CBODP ON 7CBODP ON ON 8CBODP ON ON 9CBODP ON 1068000 ON lCBODN ON 2CBODN ON 4CBODN ON 5 ON 6CBODN TCBODN ON 8CBODN ON ON 9CBODN ON l0CBODN ON ALG1 ON ON ALG2 ON ON ALG3 ON ON DO ON ON TIC ON ON ALK ON ON 7001 LDOM P ON ON RDOM P ON ON LPOM P ON ON RPOM P ON ON LDOM N ON ON RDOM N ON ON LPOM N ON ON RPOM N ON ON New control files Based on the options the user turns ON or OFF new control files are required These new control files are named 1 wi selective npt new variables controlling the selective withdrawal algorithm to select temperature targets w2 habitat npt new variables controlling fish habitat limits for temperature and dissolved oxygen and surface and segment volume weighted eutrophication state variables 3 2 envirpf npt new variables controlling setting environmental performance criteria w2 aerate npt variables describing use of dissolved oxygen addition to enhance dissolved oxygen levels through diffusers Details of these new cont
48. N NO3 ON DSI OFF PSI OFF FE OFF LDOM ON RDOM ON LPOM ON W2 Version 3 71 Release Notes RPOM ON 1 ON BOD2 ON BOD3 ON BOD4 ON BOD5 ON ALG1 ON DO ON TIC ON ALK ON 2001 LDOM P RDOM P OFF LPOM P RPOM P OFF LDOM N OFF RDOM N OFF LPOM N RPOM N EX COEF EXH20 EXSS EXOM BETA WB 1 0 45000 0 01000 0 40000 0 45000 ALG EX EXA EXA EXA EXA GENERIC CGO10 CGODK CG1DK CGS 1 0 00000 1 0000 0 00000 0 00000 CG 2 0 00000 0 00000 0 00000 0 00000 CG 3 1 04000 0 00000 0 50000 0 00000 CG 4 0 00000 0 00000 0 00000 0 00000 CG 5 0 00000 0 00000 0 00000 0 00000 S SOLIDS 555 SEDRC TAUCR SS1 1 50000 OFF 0 00 ALGAL RATE AG AR AE AM AS AHSP AHSN AHSSI ALG 2 00000 0 12000 0 02000 0 05000 0 04000 0 00500 0 00500 0 00000 ALGAL TEMP 1 AT2 AT3 1 2 ALG 5 00000 12 0000 20 0000 30 0000 0 10000 0 99000 0 99000 0 10000 ALG STOI ALGP ALGN ALGC ALGSI ACHLA ALPOM ANEON ANPR ALG 0 00500 0 08000 0 45000 0 00000 65 0000 0 80000 1 0 00100 EPIPHYTE EPIC EPIC EPIC EPIC EPIC EPIC EPIC EPIC EP OFF EPI PRIN EPRC EPRC EPRC EPRC EPRC EPRC EPRC EPRC EP OFF EPI INIT EPICI EPICI EPICI EPICI EPICI EPICI EPICI EPICI EP 10 0000 EPI RATE EG ER EE EM EB EHSP EHSN EHSSI EP 2 00000 0 05000 0 02000 0 05000 0 01000 0 00200 0 00200 0 00000 ASAT 50 0000 EPIC EPRC EPICI 39 pH IO Di Ed IO FU FU U lt lt
49. NLW JW 5 31E 13 273 154 TAIR JW 6 1 040 0017 CLOUD JW 2 0 97 else RANLW JW 5 62E 8 273 154 TAIR JW 4 1 0 261 exp 7 77 4 JW 2 1 0 0 0017 CLOUD JW 2 0 97 endif RN1 SRON JW REFL RHOWCP SHADE I 1 0 ALBEDO JW BETAI JW RANLW JW SW 4 19 10 DO WHILE DEL 1 0 AND J 500 CALL SURFACE TERMS TICE RN I RE I RC I 4 19 10 RN I SRON JW REFL RHOWCP SHADE I 1 0 ALBEDO JW BETAI JW RANLW JW RB I RE JW RC I OLD CODE DEL RN I RK1 RIMT TICE I IF ABS DEL gt 1 0 TICE TICE DEL 500 0 J J41 END DO 4 19 10 28 W2 Evaporation Units for EV in the SNP file were given in m s but were actually m 3 s 4 21 10 W2 Version 3 71 Release Notes 58 H Code W2 Fix or Description of Bug Enhancement Date Bug PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 29 W2 Ice In the ice melt algorithm SRON should not have been 4 21 2010 divided by RHOCP in computing RN1 and DEL in the DO WHILE loop should have been ABS DEL rather than DEL RN1 SRON JW REFL SHADE I 1 0 ALBEDO JW BETAI JW RANLW JW SW 4 19 10 eliminate spurious divsion of SRO by RHOCP DO WHILE ABS DEL 1 0 AND J 500 SW 4 21 10 Should have been ABS of DEL CALL SURFACE TERMS TICE 30 PRE Constituent The output f
50. Oo E EA O O on OO EA aa lt 2 o 5 2 o ct D C2IWB C2IWB C2IWB C2IWB C2IWB C2IWB C2IWB 24 9CBOD 10CBOD ALGI 0 1 ALG2 0 1 ALG3 0 1 DO 12 0 TIC 5 0 ALK 19 8 2001 0 0 LDOM P 0 0005 P 0 0005 LPOM P 0 0005 RPOM P 0 0005 N 0 0080 N 0 0080 N 0 0080 RPOM N 0 0080 CST PRIN CPRWBC CPRWBC CPRWBC CPRWBC CPRWBC CPRWBC CPRWBC CPRWBC CPRWBC TDS ON AGE ON TRACER ON COL1 ON Conduct ON Chlorine ON ISS1 ON 4 ON NH4 ON NOx ON DSi OFF PSi OFF TFe OFF LDOM ON RDOM ON LPOM ON RPOM ON 1CBOD ON 2CBOD ON 3CBOD ON 4CBOD ON 5CBOD ON 6CBOD ON 7CBOD ON W2 Version 3 71 Release Notes 25 8CBOD ON 9CBOD ON 10CBOD ON 1CBODP ON 2CBODP 3CBODP 4CBODP 5CBODP GCBODP 7CBODP 8CBODP 9CBODP 10CBODP ON 2CBODN 3CBODN 4CBODN ON SCBODN 6CBODN TCBODN 8CBODN ON 9CBODN 10CBODN ALG1 ON ALG2 ON ALG3 ON DO ON TIG ON ALK ON 7001 LDOM P ON RDOM P ON LPOM P ON RPOM P ON LDOM N ON RDOM N ON LPOM N ON RPOM N ON CIN CON CINBRC TDS ON AGE OFF TRACER OFF COL1 OFF Conduct ON Chlorine OFF ISS1 ON 4 NHA ON NOx ON DSi OFF PSi OFF TFe OFF LDOM ON RDOM ON LPOM ON RPOM ON 1CBOD ON 2CBOD ON 3CBOD ON 4CBOD ON 5CBOD ON
51. RHO K I I C2 K I CN 2C JC 1 NAC ELSE WRITE CPL JW 9999 X1 I ELWS I DEPTHB K I U K I W K I T1 K I RHO K I C2 K I CN 2C 2C 1 ENDIF WRITE CPL JW TITLE CE QUAL W2 IF HABTATC ON THEN WRITE CPL JW 19233 CNAME2 CN JN JN 1 NAC ELSE WRITE CPL JW 19234 CNAME2 CN JN IN 1 ENDIF sw 9 28 13 19233 FORMAT VARIABLES Distance m Elevation me TUE TO REO HABITAT lt NAC gt A8 19234 FORMAT VARIABLES Distance m Elevation m U W T RHO lt gt 8 sw 9 28 13 W2 Version 3 71 Release Notes 17 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 43 W2 SPECIFY TRIB In specifying the elevation between top and bottom 10 3 2013 for an inflow tributary the code put the inflow 1 layer below it should have been in many cases This has been fixed by the additional code shown below IF SPECIFY_QTR JT THEN KTTR IT 2 DO WHILE EL KTTR JT I gt ELTRT JT DO WHILE EL KTTR JT I gt ELTRT JT and EL KTTR JT 1 I gt ELTRT JT SW 10 3 13 KTTR JT KTTR IT 1 END DO 44 W2 CWO or Fixed a format overflow in writing out concentrations 10 4 2013 CWDO output in a withdrawal output file IF QWDO J 0 0 CWDO CN 2C 2 CWDO CN 2C 2 QWDO 3
52. SE WDO JWD 3 END IF IF DERIVED_CALC THEN CLOSE WDO JWD 4 END IF W2 Version 3 71 Release Notes 14 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 40 W2 Pumps Fixed several sections of code in the PUMP algorithm 9 25 13 Lateral in the hydroinout f90 routine Under some conditions such as specifying Lateral the PUMP algorithm may not have moved the water from the upstream to the downstream segment correctly This has been fixed and tested Part of the code changes are shown below IF LATERAL_PUMP JP THEN ELW EL KTWB JWU IUPU JP Z IUPU JP COSA JBU JWW JWW 1 SW 9 25 13 JBWD JWW JBU IWD JWW IUPU JP ELSE ELW EL KTWB JWU IUPU JP Z IUPU JP COSA JBU SINA JBU DLX IUPU JP 0 5 JSS JBU JSS JBU 1 SW 9 25 13 END IF IF PUMPON JP THEN IF LATERAL PUMP JP THEN JLAT 1 JWW 1 SW 9 25 13 CALL LATERAL WITHDRAWAL DO K KTW JWW KBW IWW QSS K I QSS K I QSW K JWW END DO IF IDPU JP 0 THEN MOVED CODE SW 9 25 13 1 JTT JTT 1 ELSE JSS JBU JSS JBU 1 SW 9 25 13 KTSW JSS JBU JBU KTPU JP W2 Version 3 71 Release Notes 15 Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 41 W2 Clean up
53. SI units of m s to English btu ft2 d F and then back to SI W m2 C ICETH 1 0 0 ICETH DLT RIMT I ICETH RK1 1 0 HIA T2 KT 1 RIMT RHOIRL1 OLD CODE ICETH I 0 0 ICETH I DLT RIMT ICETH 1 RK1 1 0 HIA HWI JW T2 KT 1 RIMT RHOIRL1 SW 10 20 09 Revised missing HWI JW ICE I ICETH I gt 0 0 ICESW I 1 0 IF ICE I THEN TFLUX 2 392 7 RIMT T2 KT I BI KT I DLX I OLD CODE TFLUX 2 392E 7 HWI JW RIMT T2 KT I KT I DLX I SW 10 20 09 Revised missing HWI JW TSS KT I TSS KT I TFLUX TSSICE JB TSSICE JB TFLUX DLT ICESW I 0 0 END IF END IF END DO END IF END IF W2 Version 3 71 Release Notes 55 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 22 WI Gates output The following bug was found in defining which branch 3 24 10 in QWD file a gate was located This affected the output for the withdrawals at a location where there were gates that were not tied to other branches Old code JWUGT JG JW IF IDGT JG 0 THEN DO JB 1 NBR IF IDGT JG gt US JB AND IDGT JG lt DS JB EXIT END DO JBDGT JG JB DO JW 1 NWB IF JB BS JW AND BE JW EXIT END DO JWDGT JG JW else BUG FIX 9 27 07 jbdgt jp 1 jwdgt 1 l
54. STARB 5 0 H1 KT I TKE KT I 2 BH2 K T I BH1 KT I DO K KT 1 KB I 1 BOUK AZ K I G RHO K 1 I RHO K I H K JW RHOW 0 0 PRDK AZ K I 0 5 U K I U K I 1 0 1 1 0 1 1 1 JW 2 0 PRHE 10 0 GC2 1 25 ABS 0 5 U K I 0 I 1 4 0 0 5 B K I 2 0 IF MANNINGS N JW THEN wS 5 start HRAD BHR K BR K I BR K 1 1I 2 0 H K JW GC2 G FRIC I FRIC 1 HRAD 0 333 if macrophyte on and mannings n jw th en call macrophyte friction hrad fric i effri Cy kyi gc2 g effric effric hrad 0 33333333 else AT not macrophyte on and mannings n jw then gc2 g fric i fric 1 hrad 0 33333333 W2 Version 3 71 Release Notes 69 H Code W2 PREW2 or GUI Enhance ment Type or Description of Bug Enhancement Date Bug Fixed or Enhance ment Added end if 5 end END IF PRHK GC2 0 5 B K I ABS 0 5 U K I 0 K I 1 1 44 3 0 UNST UNSE TKE K I 2 TKE PRDK TKE K I 2 I 1 PRDK 1 92 TKE K I 2 TKE K I 1 TKE K I 2 TKE TKE EN US SORT 1 0 5 TKE K I 1 K 1 1 DLT UNST P TKE K I 2 K 1 2 DLT UNSE P D DO TARB GC2 ABS 0 5 U K E KB I I 1 3 33 USTARB USTAR 1 1 2 RHK BOUK RHE B I I U KB I
55. SUB 0 THEN DLEL ELIU ESP JS IF ELID ESP JS DLEL ELIU ELID SW 6 7 10 IF DLEL 0 0 THEN DLEL DLEL For gates IF A2GT JG 0 0 AND G2GT JG 0 0 DLEL ELIU G2GT JG IF ELID EGT JG DLEL ELIU ELID SW 6 7 10 IF DLEL 0 0 THEN W2 Version 3 71 Release Notes 60 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 32 W2 Branch In cases where there are branch intersections 10 30 2010 intersections between waterbodies it was possible that the variable with multiple KBI and KB were incorrectly set Here is the fix Move e the statement defining in the subroutine init waterbodies geom f90 to the place shown below delete the earlier reference IF 10 1 0 0 B K ID 1 B K 1 1D 1 IF IEXIT 1 EXIT END IF END IF END IF END DO END DO SW 1 23 06 END DO SW 1 23 06 bnew b SW 1 23 06 KBI SW 10 30 2010 Upstream active segment and single layer 1 23 06 entire section moved SW DO JW 1 NWB KT JW DO JB BS JW BE JW 33 W2 SS The code index was incorrect in the loop for 2 3 2012 resuspension computing resuspension This led in some compilers to an infinite loop The corrected code is shown below SSSS KT I J SSS J SS KT I J BI KT I BH2 KT I 558 IF SEDIMENT RESUSPENSION J THEN Thanks to James Pasle
56. SW INTEGER KTTOP 4 5 13 SW 25 W2 Output for If the LAT option was chosen the output files index for 5 17 2013 pumps JWD was incorrect This may have affected output spillways temperatures and concentrations gates 26 PRE W2 Mass loading There were cases where the preprocessor bombed 6 21 2013 calculation while calculating the mass loading for output to the pre opt file This error has been fixed 27 W2 Assorted code Minor format errors that were ignored by compiler 6 21 2013 updates update to code comments and faster code initializations to speed up model performance were performed in several subroutines input PAR f90 temperature PAR f90 transport PAR f90 update f90 and w2 37 win f90 An example of an initialization code speed up from temperature PAR f90 New code DO K KT KB I AT K I 0 000 CT K I 0 000 VT K I 0 000 SW CODE SPEEDUP 6 15 13 ENDDO Old code AT 1 0 000 1 0 000 VT I 0 000 W2 Version 3 71 Release Notes 11 Code W2 PREW2 or GUI Fix or Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 28 W2 tools Excel macro Update More robust tools release 6 21 2013 29 PRE W2 Label error A label error for one spillway error was fixed It mistakenly used gate 7 2 2013 30 W2 CPL output A slight change in output format for the raw
57. W2 VERSION 3 71 RELEASE NOTES April 9 2015 The code updates and further information on the W2 model are available from the following web page subject to change http www ce pdx edu w2 Please address questions about the code to Scott Wells Department of Civil and Environmental Engineering Portland State University P O Box 751 Portland OR 97207 0751 503 725 4276 FAX 503 725 5950 e mail wellss pdx edu Table of Contents 2 V rsiom 3 71 Release Notes au 1 W2 V3 7 Release Package ores Renee 2 WZ KNOWN ISSUES 3 W2 V3 7 Bug Fixes Enhancements and User Manual Changes 00 0 0 020 0 0 000 eene nennen nennen nnne entr nnns 5 W2 Planned 2 2 teh nnne Ra ede rhe eue Ee Er v 21 Differences between Version 3 71 and Version 3 7 21 Differences between Version 3 7 and Version 3 6 22 Differences between Version 3 6 and Version 3 5 31 Differences between Version 3 2 and Version 3 5 33 Differences between Version 3 1 and Version 3 2 44 BUG FIXES AND ENHANCEMENTS BETWEEN
58. burial rate SEDB day can be specified in the SEDIMENT card section of the control file A different burial rate can be manual specified for each water body OLD NEW line example SED K I MAX SED K 1 LPOMEP K I SEDAS I S EDOMS K I SEDNS K I SEDD K I DLT 0 0 SED K I SED K I sedem SEDAS I k i SEDOMS I SEDNS 1 SEDD K I sedbr k 1 DLT 0 0 cb 11 30 06 11 30 06 5 Control File Add rate sediment model burial for This is the change 84 above implemented in the control file The new variable SEDBR is added in f8 format after the FSED variable SEDBR sediment burial rate in units of per day SEDIMENT SEDC SEDK SEDS FSOD WB 1 ON 0 10000 0 1 1 0 SEDPRC FSED SEDCI SEDBR ON 0 00000 1 00000 1 00000 W2 Version 3 71 Release Notes 62 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 6 W2 Sediment If a model added and subtracted layers that 4 18 07 heating sediment processes resulted in segment addition and subtraction there was the possibility that sediment fluxes were incorrectly computed In the NO3 subroutine Old code NO3SED K I NO3 NO3S JW NO3TRM K I BI K 1 BI K 1 1 BH2 K I New code if k kb i then NO3SE
59. d This code is repeated similarly in many of the sediment routines W2 Version 3 71 Release Notes 77 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 31 w2 Add segment The DEPTHM and DEPTHB were not 1 27 08 initial ization initialized correctly when a segment was added this does not affect internal computations just output for SPR and SNP files OLD CODE BKT I BH1 KT 1 H1 KT 1 DEPTHB K I HILIKT DEPTHM K I H1 KT 1 0 5 NEW CODE BH1 KT 1 H1 KT 1 DEPTHB KT I H1 KT I 5 1 27 08 DEPTHM KT SW 1 27 08 H1 KT I 0 5 W2 V3 2 BUG FIXES ENHANCEMENTS AND USER MANUAL CHANGES Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 1 W2 Waterbody When there was negative velocities at a 8 31 04 waterbody waterbody waterbody connection there was connection a possibility dependent on the bathymetry of the connection at the waterbody waterbody intersection that there could be temperature or concentration anomalies 2 W2 Lateral Added limit to the DLRHOMAX function 1 25 05 withdrawal Old code DLRHOMAX MAX DLRHOT DLRHOB New code DLRHOMAX MAX DLRHOT DLRHOB 1 OE 10 3 W2 Branch Logic in branch connectivity set up was fixed 1 25 05 connectivity Old code IF UHS JB
60. dies were present as a result of an earlier bug fix The new read statements occur in 2 places READ CON A8 9A8 8X 9A8 CDNAME2 20 CDWBC JD JW JW 1 NWB Icb 9 13 12 sw 2 18 13 READ CON A8 9A8 8X 9A8 2 KFWBC JF JW JW 1 NWB Icb 9 13 12 sw2 18 13 21 PREW2 More checks Additional checks were added to warn users of gaps in 2 20 2013 added meteorological data when interpolation may be inappropriate 22 W2 User Updated Updated User Manual many small additions and 2 20 2013 Manual edits REV3 23 PREW2 Improved Updated an error check for choosing inactive 3 21 2013 error check segments for ISNP output 24 PREW2 More checks Added checks for inflow temperature and tributary 3 28 2013 added temperatures W2 Version 3 71 Release Notes 10 H Code W2 Fix or Description of Bug Enhancement Date Bug PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 25 W2 Initial WL Changed SLOPE to SLOPEC in init u elws f90 routine 4 5 2013 Calculation since the normal depth should be based on SLOPEC END IF FUNCVALUE FLOM XAREA HRAD 0 6667 SLOPEC JB 0 5 FMANN SW 4 5 2013 RETURN END SUBROUTINE MANNINGS EQN Also changed KB I 1 to 1 1 for ELWS IF ABS DX LT XACC ELWS I RTBIS EL KB I 1 I SW 4 5 13 RETURN Also changed KTTOP from REAL to an INTEGER REAL XAREA WSURF 4 5 13
61. e file to a UNIX workstation and copy it back 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 Convert all tabs to spaces W2 Version 3 71 Release Notes W2 V3 7 BUG FIXES ENHANCEMENTS AND USER MANUAL CHANGES H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 1 W2 Fish habitat Changed temperature and DO criteria from 8 7 2012 limits t2 k i fishtemph ii and t2 k i fish templ ii and o2 k i fishdo ii to t2 k i fishtemph ii and t2 k i fis htempl ii and o2 k i fishdo ii This update is reflected in the manual Hence the high temperature limit and the dissolved oxygen minimum is less than or equal to given value rather than less than 2 2 Structure Added code to ensure that if flow is 0 in an 8 13 2012 gate pump outlet structure that the corresponding pipe temperature and concentration in the outlet file withdrawal is written as 99 0 Previously this was not fully output files implemented in the code Code such as this was inserted in several places in the subroutine outputa2 f90 IF QGT JS 0 0 THEN TAVGW JWD 99 0 CAVGW JWD 99 0 CDAVGW JWD 99 0 ENDIF 3 PREW2 Format Several output updates were made for warnings 8 16 2012 updates and errors 4 Resource Compiling files Updated so
62. fluxes were incorrectly computed In the NO3 subroutine Old code NO3SED K I K 1 NO3S JW NO3TRM K I BI I BI K 1 1 BH2 K 1 New code if k kb i then NO3SED K I NO3 K 1 NO3S JIW NO3TRM K Bl K 1 BH2 K 1 else NO3SED K I NO3 K I NO3S W NO3TRM K I BI K I BI K 1 1 BH2 K endif New code added in sediment routine if k kb i then SW 4 18 07 SODD K I 50 2 SODTRM K I BI K I else SODD K I SOD I BH2 I SODTRM K 1 BI K I BI K 1 1 Endif New code added in suspended solids routine if k kb i then SSR EPSILON DLX I BI I VOL K I else SSR EPSILON DLX BI BI 1 VOL K 1 Endif W2 Version 3 71 Release Notes 87 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 23 W2 see above New code added for heat flux to channel 4 18 07 bottom if kt kb i then SW 4 18 07 SROSED SROOUT TSEDF JW else SROSED SROOUT 1 0 BI KT 1 I BI KT I TSEDF JW Endif if k kb i then SW 4 18 07 TFLUX CBHE JW RHOWCP TSED JW T2 K I BI K I DLX I else TFLUX JW RHOWCP TSED JW T2 K 1 BI K I BI K 1 I DLX I endif New code added for sediment subroutine if k 1
63. hen the current KT was below KB the epiphyton burial rate was greater than specified in the control file epiphyton that are buried become part of the 1 order organic sediment as before epiphyton mortality now becomes part of the LPOM pool based on the EPOM fraction and is settled and transported downstream rather than going into the organic 1 order sediment model directly Currently this is non photosynthesizing but we will change in the next version 19 W2 ADD SUB There was a bug in addition and subtraction 5 26 06 layers of layers that led to water quality variables not being initialized correctly during riverine shallow flow 20 User Typos The manual had a few typos that were 6 11 2006 Manual corrected corrected 21 W2 Waterbody The subroutine Upstream_velocity under 6 29 2006 waterbody specific conditions did not maintain flwo connection continuity across a waterbody waterbody connection 22 W2 SNP output The algal limiting nutrient SNP output had a 6 30 2006 bug under specific conditions in writing out the information W2 Version 3 71 Release Notes 86 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 23 W2 Sediment If a model added and subtracted layers that 4 18 07 heating sediment processes resulted in segment addition and subtraction there was the possibility that sediment
64. inal code overwrote vertical eddy viscosity at the bed during the averaging process Az changed to Az 5 TKE array was initialized to zero TKE was implemented in add sub layers like AZ Many of these fixes are a result of the work of Sam Gould Gould 2006 who wrote an MS project report at PSU entitled k e Turbulence Model Further recommendations by Gould 2006 will be incorporated into the next version of CE QUAL W2 The old code is shown below as a reference to the new code in the release version OLD CODE 10 4 07 W2 Version 3 71 Release Notes 68 H Code W2 PREW2 or GUI Fix Enhance ment Type or Description of Bug Enhancement Date Bug Fixed or Enhance ment Added ENTRY CALCULATE TKE USTAR SORT 1 25 CZ I WIND10 I 2 RHO KT I IF MANNINGS N JW THEN HRAD BHR1 KT I BR KTI I I BR KT 1 1 2 AVH1 KT I if macrophyte on and mannings n jw th en call macrophyte friction hrad fric i effri c kt i gc2 g effric effric hrad 0 33333333 else if not macrophyte on and mannings n j w then gc2 g fric i fric 1 hrad 0 33333333 end if ELSE GC2 0 0 IF FRIC 1 0 0 GC2 G FRIC 1 FRIC 1 END IF USTARB SORT GC2 ABS 0 5 U KT I YU KT I 1 TKE KT 1 1 E 0 5 3 33 USTAR USTAR USTARB USTARB TKE KT I 1 BH2 KT I BH1 KT I TKE KT 1 2 0 5 USTAR USTAR USTAR USTARB USTARB U
65. ine in a DOS box you can execute the GUI as follows C scott research corps of engineersNtomcoleNw2codeNGUI36Nw2controlN w2control36 exe C scott w2workshop 2009 workshop waterqual probleml w2_con npt The first string in quotes executes the GUI The command line argument is NOT in quotes This program was developed in VB6 and does not take quotes around the command line Note that this is different than the FORTRAN command line argument So the above command will open the GUI and load the control file automatically 3 A text box now shows the file path and name of the file that you are working on 4 In file open earlier all npt files were shown Since only w2_con npt files are loaded into the GUI only the w2_con npt file was shown for opening W2 Version 3 71 Release Notes 53 Code W2 PREW2 or GUI Fix or Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 19 W2 Gates spillways pipes Whenever DOWN was specified for a gate spillway or pump the model estimated the water level at the end of the segment rather than using the branch center water level This is important in sloping river systems where a long segment may have a water surface elevation drop between the segment center and the edge In the past this was computed assuming the slope of the channel This was updated to estimate the water surface elevation using linear inter
66. k i jz prefz jjz jz zoo k i jjz cb 5 17 2007 end do do 312 omnivorous zooplankton ZGZ k 1 332 32 Zmu K I jz ZOO K I jz prefZ jz jjz tgraze K I jz ZGZ k i 33z jz Zzmu K I jz ZOO K I jz prefZ2 jjz jz tgraze K I jz ku 5 9 2007 end do 8 PRE More checks Added checks for Sediment burial rate and 6 2 2007 some further checks on grid geometry added output on SEDS and SEDBR to the pre opt file fixed condition where NZP had to equal 1 to work W2 Version 3 71 Release Notes 65 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 9 W2 Array The deallocate command on line 7557 was 6 4 2007 deallocation commented out to avoid a deallocation error when the STOP button is pushed during execution on a PC deallocate sedbr sedbrp sedbrn sedbrc SW 6 4 07 No need to deallocate pointers 10 W2 Initialization For code setting up an external head BC the 6 17 2007 of IUT variable IUT was not initialized before it was used This was fixed below Boundary bottom layers IF UH EXTERNAL JB B IUT 1 KB IUT IF UH EXTERNAL JB B IU 1 B IU cb 6 12 07 IF UH INTERNAL JB THEN IF JBUH JB gt BS JW AND JBUH JB lt BE JW THEN 1 KB IUT 1 MIN KB UHS JB KB IUT IU MIN UHS JB KB IU cb 6 12 07 ELSE IF
67. les for DS 1 The following new code was added AVHR KT DS JB 1 1 KT DS JB 1 SW 03 08 05 W2 Version 3 71 Release Notes 81 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 9 W2 Interpola tion Possible index error if there are multiple 5 10 05 multipliers waterbodies Old code RATZ K JW AVH2 K 1 1 2 K I CURZ1 K JW 2 0 H JW 2 AVH2 K 1 1 2 K I AVH2 K 1 I CURZ2 K JW 2 0 H JW 2 AVH2 K 1 I AVH2 K 1 CURZ3 K JW 2 0 H JW 2 AVH2 K 1 1 2 AVH2 K I END DO New code RATZ K JW AVH2 K 1 05 JW AVH2 K DS BE JW CURZ1 K JW 2 0 H JW 2 AVH2 K 1 DS BE JW AVH2 K DS BE JW AVH2 K 1 DS BE JW CURZ2 K JW 2 0 H JW 2 AVH2 K 1 DS JW AVH2 K DS BE JW CURZ3 K JW 2 0 H JW 2 AVH2 K 1 DS BE JW AVH2 K DS BE JW AVH2 K DS BE JW 10 W2 Spillway and Older code in order to check if it was 5 10 05 Gates submerged or not used the elevation difference relative to the channel bed on either side of the weir rather than the weir crest Also removed code line ELDN gt ESP JS DH ELUP ELDN 11 W2 Reaeration Corrected formula errors i
68. lue of zo for example if he she leaves the spaces blank for zo using a V3 5 control file the code uses 0 001 m HYD COEF AX DX CBHE TSED FI TSEDF FRICC 70 WB 1 1 00000 1 00000 0 30000 11 5000 0 01000 1 00000 MANN 0 001 A new option for output is in the format required for TECPLOT For TECPLOT animation there is only a flag in the CPL output line This allows for easy model animation of the variables U W T RHO and all active constituents at the frequency specified by the CPL file as a function of distance and elevation CPL PLOT CPLC NCPL TECPLOT WB 1 ON 1 ON A new variable for determining the fraction of NO3 N that is diffused into the sediments that becomes organic matter or SED N was introduced According to one study only about 3796 of NOs N that diffuses into the sediments becomes incorporated into organic matter in the sediments The rest is denitrified NITRATE NO3DK NO3S FNO3SED Wb 1 0 05 0 0 0 37 Wb 2 0 05 0 0 0 37 In V3 5 the model computed an average decay coefficient of the sediments based on what was deposited The user now has the option to dynamically compute that decay rate or to have it fixed and controlled by the model user A new variable was introduced called DYNSEDK which is either ON OFF to allow or not allow dynamic computation of the sediment decay rate SEDIMENT SEDC PRNSC SEDCI SEDK SEDS FSOD FSED SEDBR Wb 1 ON ON 0 0 0 1 0 0 1 0 1 0 0 001 Wb 2 ON ON 0 0 0 1 0 0 1 0 1 0 0 001 W2 Version 3 71 Release N
69. ly set up the model bathymetry prior to using the Bathymetry editor within W2Control A user manual in pdf format is included in this directory Also a separate executable W2Control is provided in case an earlier version has already been installed Note that this GUI is a part of the install routine for W2Tools now W2Tools This is the new W2 post processor by Dynamic Solutions International LLC www ds international biz They have provided an installation routine that includes both the post processor and the W2ControlGUI When the user selects W2L output the old VPL output the resulting post processing file is used by W2Tool for all post processing tasks that include contour plots animations profile plots and time series plots A brief user manual is included showing many of the features of this post processor as well as a directory that shows how to take field data and plot field data and model results in the post processor There is a zip file with an example from DeGray reservoir on how to include model predictions versus field data for reservoir profiles Source This directory contains the source code for the preprocessor and model written in Fortran The compiler settings and files necessary to compile using the Intel compiler are also included using the Intel Fortran compiler Generally we use the following compiler settings O2 maximum speed in Intel and default real is double precision Also for the following subroutines
70. me corrupted resource files that 9 12 2012 files for W2 were used to compile the source code Also zipped up source code and compiler settings together so that file locations are correct for using the Intel compiler W2 Version 3 71 Release Notes H Code W2 Fix or PREW2 or GUI Enhancement Type Description of Bug Enhancement Date Fixed or Bug Enhancemen t added and Read csv files PREW2 By inserting the character 5 as the first character of the first line the following files can now be read in free format or csv format met lpr vpr wsc met cin ctr cdtr cpre qot and qwd This is described in a Word document that the The preprocessor has also been updated for file accompanies download package checks This is part of the Version 3 71 update 9 12 2012 Read input file An input format bug was fixed for a system with more than 9 waterbodies DO JD 1 NDC READ CON A8 9A8 CDNAME2 JD CDWBC JD JW JW 1 NWB READ A8 9A8 8X 9A8 CDNAME2 JD CDWBC JD JW JW 1 NWB cb 9 13 12 END DO READ CON DO JF 1 NFL do 1 73 Fix this later READ CON 9A8 2 JF KFWBC JF JW JW 1 NWB READ CON A8 9A8 8X 9A8 KFNAME2 JF KFWBC JW JW 1 NWB cb 9 13 12 END DO This had the effect of turning OFF output for derived constituents for
71. ment Fixed or or GUI Type Enhancemen t added 48 PREW2 Warnings Fixed a name inconsistency for developing warnings 7 15 14 for input concentrations IF NAME Residence time AND NAME Water age THEN IF NAME Residence time AND NAME AGE THEN SW 7 15 14 CALL WARNINGS 49 W2 TSR filename The filename in w2_con npt for TSR is used for the 8 22 14 output filenames In order to account for complex paths that include than one the following change was made with the BACK TRUE command which checks from the right hand side rather than left hand side of the character string L1 SCAN TSRFN L1 BACK TRUE 50 PREW2 Hydraulic Added many new hydraulic structure warnings gates 9 10 14 structure spillways pumps pipes internal weirs for cases warnings where KBSTR was less than KB and fixed a few error messages for these structure checks 51 W2 TSR output The time series file has added the surface heat flux 1 15 15 terms net short wave solar net long wave radiation net back radiation heat flux evaporation heat flux conductive heat flux to the output The manual was also updated 52 W2 Interpolation In some cases the wind direction interpolation was 2 13 15 of wind direction incorrect Code was added to reduce the wind direction angle to less than 2 pi before the interpolation is performed and to consider another possible interpolation case
72. n 3 71 Release Notes 48 H Code or GUI 2 PREW2 Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 7 Prew2 Constituent loads An enhancement was added to the Preprocessor to compute loads in kg day for all inflow tributary and distributed tributaries Also these are summed up for the model application These are shown in the file pre opt These are approximate loads since the concentration data are used to set the frequency of loading update Flow rates at the time of the concentration input data are used to compute load 5 21 2009 Gas transfer at spillways A couple code fixes in the hydroinout f90 subroutine 1 CGAS needed to be initialized in some cases to CGAS C2 K ID CN JC prior to subroutine TOTAL_DISSOLVED_GAS for use in the Butts and Evans 1983 equation NEW CODE CGAS C2 ID CN JC MM 5 21 2009 calling the 2 Change logic in several lines from IF CAC NDO ON to IF CAC NDO and CN JC NDO NEW CODE IF CN JC NDO AND CAC NDO ON AND GASSPC JS ON AND OSP JS gt 0 0 THEN MM 5 21 2009 5 21 2009 Reaeration from dams An error was found in the formulae from Butts and Evans 1983 OLD CODE DB SAT C DA DB 1 0 0 38 AGASGT N BGASGT N CGASGT N 1 0 0 11 CGASGT 1 0 0 046 SAT DA NEW CODE DA SAT C MM
73. n Thackston and 5 10 05 Krenkel formula Old code USTAR SQRT ADEPTH SLOPE JB 32 2 0 5 REAER I 24 88 1 0 SORT 0 176 UAVG SQRT ADEPTH USTAR New code USTAR SQRT ADEPTH SLOPE JB 32 2 REAER I 24 88 1 0 SORT 0 176 UAVG SQRT ADEPTH FUSTAR ADEPTH Similar changes were made to the updated Thackston model Eqn 10 12 W2 Violations NV The variable Bl and VOL was not initialized 8 25 05 properly during a time step violation W2 Version 3 71 Release Notes 82 Code W2 PREW2 or GUI Fix or Enhance ment Type Description of Bug Enhancement Date Bug Fixed or Enhance ment Added 13 W2 ADD a layer The variable BI was not initialized properly during an ADD layer 8 25 05 14 W2 TRIDIAG subroutine Insert Deallocate Statement in Tridiag SUBROUTINE TRIDIAG A V C D S E N U USE PREC INTEGER IN TENT IN ate S E N REAL R8 DIMENSION INTENT IN HE A E V E C E D E REAL DIMENSION INTENT OUT U REAL R8 ALLOCATABLE DIMENSION ve BTA GMA ALLOCATE BTA N GMA N 5 V S GMA S D S DO I S 1 E BTA I V I A I BTA I 1 C I 1 I D I A I BTA I 1 GMA I 1 END DO U E GMA E BTA E DO 1 5 1 U I GMA I C I U I 1 BTA I END DO Deallocate BTA GMA lt SW 10 17 05 END SUBROUTINE TRIDIAG 10 17 05 15 W2 SUB layer In SU
74. n program for rivers The code below should have the subscript JB instead of J DO JJW 1 NWB DO JJB BS JJW BE JJW IF DHS JB gt US JJB DHS J lt 08 JJB THEN JBD JJB END IF END DO 8 20 13 W2 Version 3 71 Release Notes 12 H Code W2 Fix or Description of Bug Enhancement Date Bug PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 34 W2 INIT WL There was an index error with gates in the initial water 8 21 2013 level computation The old code is shown below IF ELWS ID lt WSUP THEN IF ELWS IDSP JS gt WSUP WSUP ELWS IDSP JS CHECKING TO SEE IF DOWNSTREAM WS ELEVATION ISN T ALREADY HIGH ELWS ID WSUP The new code is IF ELWS IDGT JG gt WSUP WSUP ELWS IDGT JG CHECKING TO SEE IF DOWNSTREAM WS ELEVATION ISN T ALREADY HIGH WX 8 21 13 35 w2 GATE Cleaning up some code in the gate algorithm 8 21 2013 Old code IF A2GT JG 0 0 AND IDGT JG 0 0 THEN New code IF A2GT JG 0 0 AND IDGT JG 0 THEN 36 W2 TSS Updated the computation for the derived variable TSS 9 6 2013 computation to include zooplankton and the particulate form of CBOD A formula was added to the User Manual reflecting this change New code includes IF CBODS IBOD 0 0 TOTSS K I TOTSS K I K I O20M JW SW 9 5 13 Added particulate CBOD to TSS computation TOTSS K I TOTSS 700 K I JZ
75. o jd 1 nbod sedsum sedsum MAX cbods 3 0 0 cbod K I Jd K 1 2 1 1 0 BI K 1 1 BI K 1 RBOD JD CBODD K I JD O20M JW end do W2 Version 3 71 Release Notes 74 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 28 W2 SRO There are some cases when segments 1 18 2008 were added subtracted that the value of BI was not correctly initialized This code is a fix to prevent such occurrences OLD CODE SRONET SROIN SROOUT SROSED SROOUT 1 0 BI K 1 I BI K I TSEDF JW NEW CODE SRONET SROIN SROOUT if k kb i then SW 1 18 08 SROSED SROOUT 1 0 BI K 1 1 BI K 1 TSEDE JW else SROSED SROOUT TSEDF JW endif W2 Version 3 71 Release Notes 75 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 29 W2 Water Quality Added several calls prevent 1 18 2008 computation of kinetic variables if epiphyton are defined in the control file with NEP 1 or more but is not ACTIVE or turned ON If the kinetic expressions are non zero and the initial concentration is given then this could add source sink terms to the oxygen balance This is typical of the code changes since several of this type were made OLD CODE DO JE 1 NEP POAEG K I
76. on of one layer The definition of KBMIN was not updated if the model started out in some segments with only one_layer This has been fixed Added code highlighted DO I IU ID IF KB I KT lt NL JB 1 1 1 ONE LAYER I KT KB 1 END DO CUS JB IUT Areas and bottom widths LE NOT TRAPEZOIDAL JW TH 12 17 07 W2 Version 3 71 Release Notes 72 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed or PREW2 Enhance or Enhance or GUI ment Type ment Added 24 W2 Bottom This is a couple more fixes related to bug fix 12 17 2007 processes 6 above The Denitrification rate and epiphyton burial rates could be affected based unique combinations of adding subtracting segments that left the value of Bl in an inactive layer below KB defined incorrectly In order to prevent the possibility of problems the following fixes were made Old Code SedNO3 K I I NO3S JW NO3TRM K I BI K I BI K 1 1 BH2 K 1 EPM K I J EPD K I J BI K 1 BI K 1 1 2 0 H1 K 1 DLX I New code if k kb i then SW 12 16 07 SedNO3 K I K I NO3S JW NO3TRM K I BI K I BH2 K I else SedNO3 K I K I NO3S JW NO3TRM K I BI K I BI K 1 1 BH2 K 1 endif if k kb i then SW 12 16 07 EPD K I
77. on is 5 318 13 T 273 where units are W m C at 2 m height Below 40 F 5 C the formula of Idso and Jackson is recommended above 109 both equations are almost identical d o T 273 1 0 261exp 7 77E 472 where units are W m and T is in units of C The Stefan Boltzmann constant 5 62E 8 W m K W2 Version 3 71 Release Notes 80 H Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 7 W2 Layer addition Mistyped subscript K instead of I 3 2 05 algorithm Old code IF KB I KBI I THEN B KB K I 0 0 DX KB I I 0 0 KB I KB I 1 IF I DS JB 1 KBMIN I MIN KB K KB 1 1 IF 1 US JB 1 KBMIN 1 1 MIN KB I 1 KB I New Code IF KB I gt KBI I THEN B KB I I 0 0 SW 3 2 05 DX KB I I 0 0 KB 1 KB 1 1 IF 1 DS JB 1 KBMIN 1 MIN KB 1 KB I 1 SW 3 2 05 IF I US JB 1 KBMIN I 1 MIN KB I 1 KB I 8 W2 Variable In some cases when there was a layer 3 9 05 initialize tion subtraction and a time step violation immediately afterward the variable SW was not initialized properly This caused problems in the Tomas Algorithm for the water surface computation The following line of code was added to the SUB layer algorithm SW KT 1 1U 1 1D 1 0 0 3 9 05 Also the variable AVHR was defined in the Update variab
78. ontal pressure gradient HPG m 3 s Gravity term channel slope GRAV m 3 s TDS g m 3 Age days Tracer g m 3 Bacteria col 100ml Conductivity mhos Chloride mg l ISS g m 3 Phosphate g m 3 Ammonium g m 3 Nitrate Nitrite g m 3 Dissolved silica g m 3 Particulate silica g m 3 Total iron g m 3 Labile DOM g m 3 Refractory DOM g m 3 Labile POM g m 3 Refractory POM g m 3 CBOD1 g m 3 CBOD2 g m 3 CBOD3 g m 3 CBOD4 g m 3 CBOD5 g m 3 W2 Version 3 71 Release Notes FMTH I10 8 8 8 8 8 8 8 8 8 8 8 8 8 8 O N P Pi Pi P Pi PP Pi P PD P P E D EN c EN EN C C lt lt gt gt Fal WO an al van E H 0 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 N BI Pi Pi P Pi pQ Pi Pi pQ PD pQ Pi Pi P P PP Pi P Pi P D C EN EN C CH EEN EN C CH 000 EN C C5 HMULT Co OO C CO O 0 0 CMULT 1000 1000 O O Co Co C C C O DO CO OOO O OO O Oo 50 0000 15 0000 15 0000 25 0000 1 0000 5 00000 70 000 5 00000 145 000 60 0000 50 0000 20 0000 9 00000 10 0000 10 0000 10 0000 formats HMIN 1 0 1000 10 0 997 0 1E 08 1 08 1 08 1 08 1 08 1 08 1 08 1 08 1 08 0 0 1 0 SEH 20 000 20 000 20 000 20 000 20 000 St 0 0 1000
79. otes 31 The User can now specify the of processors to use on the host computer Most users find that setting NPROC 2 gets the best results Sometimes setting this greater than 2 results in slower model performance Also the CLOSEC control closes the windows dialog box after the model completes its simulation This is useful in using the windows version of the release code in batch simulations These are specified in the control file as follows GRID NWB NBR IMX KMX CLOSEC 1 4 66 117 W2 Version 3 71 Release Notes 32 DIFFERENCES BETWEEN VERSION 3 2 AND VERSION 3 5 The differences in V3 5 and V3 2 input files are found in the control file w2 con npt and in the graph npt file other files are the same between the 2 versions w2 con npt Below is an example of parts of the control file from V3 5 where all new variables are highlighted Most of these changes have to do with the new zooplankton macrophyte and new state variables added to the model See the User Manual for a list of changes between V3 2 and V 3 5 in the version history Also there were some deletions from the V3 2 w2 con npt file These are shown below New variables added to the control file are highlighted IN OUTFL NTR NST NIW NWD NGT NSP NPI NPU CONSTITU NGC NSS NAL NEP NBOD 5 1 1 1 5 MISCELL NDAY 100 CST COMP CCC LIMC CUF ON ON 10 CST ACTIVE CAC TDS OFF Genl ON Gen2 Gen3 Gen4 Gen5 551 PO4 NH4 NO3 DSI PSI FE LDO
80. polation rather than the grid slope Below is an example of the code fix in this case for GATES OLD CODE ELIU ELWS IUGT JG SINA JBUGT JG DLX IUGT JG 0 5 NEW CODE ELIU ELWS IUGT JG ELWS JG 1 0 5 DLX IUGT JG DLX IUGT JG 1 DLX IUGT JG 0 5 ELWS JG 9 25 09 20 W2 New executable A new executable was made using a new release of Intel Version 11 compiler that corrected problems with Windows 7 applications 9 25 09 W2 Version 3 71 Release Notes 54 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 21 W2 ICE cover There were a couple logic errors in the ice cover 10 20 09 algorithm algorithm These were corrected below Ixxxxxxxxxxxx x Toe thickness ICETH I ICETH 1 ICETHU ICETH1 ICETH2 IF ICETH I lt ICE TOL ICETH I 0 0 IF WINTER AND NOT ICE IN JB THEN IF NOT ALLOW ICE I ICETH I 0 0 END IF ICE I ICETH I gt 0 0 IF ICE I THEN 3 27 08 SW ICESW I 0 0 ELSE ICESW I 1 0 ENDIF ICETHU 0 0 0 0 2 0 0 IF ICETH I lt ICE TOL AND ICETH I 0 0 ICETH I ICE TOL ELSE IF TERM BY TERM JW CALL EQUILIBRIUM TEMPERATURE SW 10 20 09 Must call this first otherwise ET and CSHE are 0 HIA 0 2367 CSHE 1 5 65E 8 JM 11 08 convert
81. raph file is shown below Constituent hydrodynamic and derived constituent names controls TDS g m 3 or Salinity kg m 3 Generic Constituent g m 3 41 Generic Constituent g m 3 42 Generic Constituent g m 3 3 Generic Constituent g m 3 44 Generic Constituent g m 3 5 Suspended solids g m 3 1 Phosphate g m 3 Ammonium g m 3 Nitrate Nitrite g m 3 Dissolved silica g m 3 Particulate silica g m 3 Total iron g m 3 Labile DOM g m 3 Refractory DOM g m 3 Labile POM g m 3 Refractory POM g m 3 CBOD g m 3 CBOD g m 3 2 g m 3 3 g m 3 4 g m 3 5 Algae g m 3 1 Dissolved oxygen g m 3 Inorganic carbon g m 3 Alkalinity g m 3 Timestep violations NVIOL Horizontal velocity U m s Vertical velocity W m s Temperature T1 lt o gt C Density RHO g m 3 Vertical eddy viscosity AZ m 2 s Velocity shear stress SHEAR 1 s 2 Internal shear ST m 3 s Bottom shear SB m 3 s Longitudinal momentum ADMX m 3 s Longitudinal momentum DM m 3 s Horizontal density gradient HDG m 3 s Vertical momentum ADMZ m 3 s Horizontal pressure gradient HPG m 3 s Gravity term channel slope GRAV m 3 s Dissolved organic carbon g m 3 Particulate organic carbon g m 3 Total organic carbon g m 3 Dissolved organic nitrogen g m 3 Particulate organic nitrogen g m 3 Total organic nitrogen g m 3 Total Kheldahl Nitrogen g m 3 W2 Version 3 71 Release Notes
82. rol files are in the CE QUAL W2 User Manual W2 Version 3 71 Release Notes 30 DIFFERENCES BETWEEN VERSION 3 6 AND VERSION 3 5 Version 3 6 can be run without changing any of the input files even though the preprocessor will identify errors in the control file because of missing variables Below is a highlighted list of locations in the file w2 con npt where additional variables have been added There are no other changes in the input files for Version 3 6 The TKE algorithm has been updated with new algorithms that match experimental tank data for kinetic energy and dissipation This is based on a Master s degree project by Sam Gould at Portland State University A new user option is the 1 algorithm add addition to the legacy algorithm TKE This results in several new input variables on the following line of the w2 con npt file that are only active if TKE1 is chosen for AZC WB 1 W2 IMP 1 00000 EDDY VISC Azc AZSLC ARODI STRCKLR BOUNDER TKECAL The roughness height of the water for correction of the vertical velocity wind profile is now a user defined input zo Prior to this the model had hardwired the value of zo 0 003 m for wind speed correction at 2m for evaporation where wind height at 2 m is typical and 20 0 01 m for wind at 10 m for shear stress calculations where wind height of 10 m is typical For consistency both conversions now use the same value of roughness height If the user does not specify the va
83. rom the preprocessor in the pre opt file 5 10 10 loading for constituent loading was in kg rather than the output header of kg day The output was updated to kg day by adding the following lines of code cdtload incdt 1 NACdt Jb Jb jb cdtload incdt 1 NACdt Jb Jb jb jday tstart CB 5 10 10 Change units to kg day ctrload trcn 1 NACtr Jt Jt jt ctrload trcn 1 NACtr Jt Jt jt JDAY TSTART 5 11 10 convert to units of kg day W2 Version 3 71 Release Notes 59 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 31 W2 Gate In the case where the user has specified that the flow 6 4 10 spillways is DOWN in the case of reverse flow the model did pipes not assign the flow correctly if the user had no other tributaries or withdrawals specified in the control file For this rare event additional code was written to account for this fact Also a logic error was discovered in reverse flow for spillways and gates This was corrected New code added to hydroinout f90 JWW NWD withdrawals jww gt 0 6 4 10 SW JTT NTR tributaries jtt 0 6 4 10 SW JSS NSTR IF SPILLWAY THEN END IF tributaries jtt 0 6 4 10 SW withdrawals jww gt 0 6 4 10 SW DO JW 1 NWB KT JW DO JB BS JW BE JW New code in gate spill pipe f90 For spillway IF I
84. sedbr jw DYNSEDK JW JW 1 NWB SW 6 1 07 READ A8 8X 2A8 6F8 0 A8 ERR 400 AID SEDC JW PRNSC JW SEDCI JW SEDDK JW seds jw FSOD JW amp FSED JW sedbr jw DYNSEDK JW JW 1 NWB cb 12 30 12 17 Excel macro Integer Long Some loose ends were corrected in the Visual Basic 1 2 2013 utility variables code built into the Excel macros w2tool W2 Version 3 71 Release Notes H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 18 W2 TDG output A series of code changes were made to fix some issues 1 23 2013 that arose for computing the impact of a structure on downstream TDG These fixes were made in subroutines Withdrawal outputa2w2tools w2modules hydroinout These affected calculation of output of dissolved gas concentration for output files for spillways or gates that had dissolved gas equation 19 W2 Reading in For temperature only studies the model did not 1 28 2013 dynamic update the dynamic light extinction coefficient extinction correctly This has been fixed by the added code coefficient below DO JW 1 NWB IF READ EXTINCTION JW GAMMA US BS JW D S BE JW EXH2O JW SW 1 28 13 KTWB JW IF NOT NO HEAT JW THEN 20 2 Input format A specific input read error occurred when 9 2 18 13 when 9 WBs waterbo
85. sedbr K l sedb jw sed k i 21 W2 TKE Turbulence model had an improper averaging 12 17 07 between layers A new temporary variable was defined to temporarily store the values for AZ prior to averaging to the bottom top of the layers and the horizontal layers This also affected the computation of DZ Fixed New code defined AZT and allocated memory for it such that AZT K I 0 09 TKE K I 1 TKE K I 1 TRI K I 2 and AZ K I 0 5 AZT K I 7 K 1 1 Similarly for the horizontal averaging and for DZ Also the values of DZ were fixed to be at the bottom of a cell and AZ was fixed to be at the bottom right hand edge of a cell as shown below CE QUAL W2 coordinate system X e t Segment B WD i it z h at bottom W2 Version 3 71 Release Notes 71 Code W2 PREW2 or GUI Fix or Enhance ment Type Description of Bug Enhancement Date Bug Fixed or Enhance ment Added 22 W2 SS settling The incorrect cell width was used for SSSO BI KT I was changed to BI K I OLD CODE 550 K 1 TOTSSO FES JW FPFE K I BI K T 1 BH2 K I DO1 K I FPSS K I FPSS K I TISS K I EW CODE SSSO K I TOTSSO FES JW 1 BE2 I DO1 I FPSS K I FPSS I TISS K I z 12 17 07 23 W2 Initial izati
86. t END IF New code JWUGT JG JW IF IDGT JG gt 0 THEN DO JB 1 NBR IF IDGT JG gt US JB AND IDGT JG lt DS JB EXIT END DO JBDGT JG JB DO JW 1 NWB IF JB gt BS JW AND BE JW EXIT END DO JWDGT JG JW else BUG FIX 9 27 07 jbdgt 39 1 SW 3 24 10 jwdgt 39 1 SW 3 24 10 END IF JB lt W2 Version 3 71 Release Notes 56 Code W2 PREW2 or GUI Fix or Enhancement Type Description of Bug Enhancement Date Bug Fixed or Enhancemen t added 23 Prew2 Reading of WSC Reading in of the WSC file was limited to only 100 dates in the preprocessor This limitation was fixed by the code shown below DO J 1 100 28995 continue cb 3 26 10 READ 10F8 0 8X 9F8 0 END 29000 SDAY WSC I I 1 IMX IF SDAY SDAYO THEN CALL ERRORS WRITE ERR 3 A F0 3 Julian date SDAY previous date of V in WSCEN END IF DO 1 IF WSC I lt 0 0 THEN CALL ERRORS WRITE ERR A F0 3 A 14 A Julian date SDAY WSC AT SEG 1 1 lt 0 0 in WSCEN ENDIF IF WSC I gt 2 0 THEN CALL WARNINGS WRITE WRN A F0 3 A 14 A Julian day SDAY WSC I SEG 1 1 gt 2 0 in WSCEN END IF IF WSC I gt 0 0 and wsc i lt 0 5 THEN CALL WARNINGS WRITE WRN A F0 3 A 14 A Julian day SDAY WSC I SEG I I lt 0 5 in WSCFN END IF ENDDO SDAYO SDAY
87. tion At the end of a windows run the windows dialog box waits for the user to press close to exit the window This allows the user to examine the final run parameters In the w2_con npt file there is now an option to close this window when the run has completed If this option is not set then the dialog box will stay until the user clicks close This allows for efficient batch processing of the model especially if user in conjunction with command line processing mentioned in 15 GRID NWB NBR IMX KMX NPROC 000580 When CLOSEC is set to ON then the dialog box will disappear once the run finishes If it is set to OFF then 1 1 23 22 the dialog box will remain until the user clicks close 9 12 09 17 User Manual Updates Updates and changes to the control file 813 816 were reflected in an updated User Manual 9 12 09 W2 Version 3 71 Release Notes 52 H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 18 GUI Updates The GUI was updated with the following 9 12 09 1 new control file parameters NPROC and CLOSEC were added see 14 and 16 There is also a SELECTC that will be used in V3 7 that has been included ignore it for now 2 The GUI also can be controlled by command line passing of the working directory and file In a batch program or from the command l
88. waterbody 10 9 13 2012 7 GUI Time series elevation The GUI read in values of ETSR as integers rather than real numbers This was fixed 10 30 12 Spillways Lateral Lateral spillways when connected to other model segments were sometimes not connecting as a tributary to the downstream segment This has been fixed 10 30 12 W2Tools output In place of the Vector Plot Output VPL a new output was added that allows use of the W2Tools post processing package This is part of the Version 3 71 update 10 30 12 10 2 User Manual The User Manual has been updated with the new model features as shown in 5 and 9 above In addition a separate user manual file shows how to use the w2tools post processor This is in the directory for W2tools This is the version 3 71 update 10 30 12 W2 Version 3 71 Release Notes H Code W2 Fix or Description of Bug Enhancement Date Bug or PREW2 Enhancement Fixed or or GUI Type Enhancemen t added 11 W2 Water quality A new calculation technique was added that 10 30 2012 and eliminates calling the Tri diagonal subroutine temperature These were built into the temperature and water quality subroutines This change results in improvements in computational speed of from less than 5 to over 20 for water quality models with lots of water quality state variables 12 PREW2 More checks Added more error trapping for input files This is
89. we had to use O1 optimization init cond f90 and init u elws f90 For the preprocessor the windows source code is compiled using a QuickWin application rather than a console application We use the debug version for the released executable The generic preprocessor code should work compiled as a console application W2 Version 3 71 Release Notes 2 Waterbalance This is the windows waterbalance utility that is described in the user manual The purpose of this code is to approximate the waterbalance for a reservoir or lake by computing flows positive and negative that will allow the model predicted water level to agree to water level data for a reservoir Excel macro utility for writing files in W2 format from Excel This directory contains an Excel macro that aids in writing our CE QUAL W2 compatible files from within Excel There is a short user manual describing how to use the macro This macro was developed by Jeffrey Gregory Civil Engineer USACE Nashville District W2V3 manual371 revX pdf User Manual in searchable pdf format where X is the rev number e W2 Version 3 71 Release Notes pdf Release notes in pdf format W2 KNOWN ISSUES The following list shows known bugs and issues with the current release of the code these being addressed in the next release Item Description 1 Water levels in If water levels decrease in a waterbody shaped like a bowl the removal of bowl model
90. y for this bug report fix W2 V3 5 BUG FIXES ENHANCEMENTS AND USER MANUAL CHANGES Code W2 Fix or Description of Bug Enhancement Date Bug Fixed PREW2 Enhance or Enhance or GUI ment Type ment Added 1 wW2 Zooplank ton Sign error in the zooplankton grazing on algae 8 23 06 algae term 2 W2 Input output Format for I O was changed to allow better 8 23 06 decimal precision of output W2 Version 3 71 Release Notes 61 H Code or GUI 2 PREW2 or Description of Bug Enhancement Date Bug Fixed or Enhance ment Added 3 W2 Sediment settling rate The sediment settling rate was accidentally used for POM settling This was fixed The old and new code lines are shown below OLD sedsum sedsum seds JW LPOM K I 1pomdk jw RPOM I rpomdk jw BI K I BH2 K I 1 0 K 1 1 BI K I NEW sedsum sedsum poms JW LPOM K I lpomdk jw I rpomdk jw K I BH2 K 1 1 0 K 1 1 BI K I cb 10 22 06 This was an issue in the SEDIMENT SEDIMENT C SEDIMENT P SEDIMENT and SEDIMENT DECAY subroutines 10 26 06 4 W2 Sediment burial An algorithm was added for sediment burial This is now a new parameter in the sediment part of the control file An updated user description is forthcoming The sediment

W2 Version 3 Release Notes - Civil & Environmental Engineering

Contents

Download Pdf Manuals

Related Search

Related Contents

W2 Version 3 Release Notes - Civil &amp; Environmental Engineering

Contents

Download Pdf Manuals

Related Search

Related Contents

W2 Version 3 Release Notes - Civil & Environmental Engineering