Portland State University

Program Overview

 

2009 CE-QUAL-W2 Workshop in Portland, Oregon

Noteworthy features of the new version of CE-QUAL-W2 Version 3.6:

  1. The code has been rewritten into smaller subroutines to allow better code compilation and optimization.
  2. The code has been revised with the goal of improving the computational speed. This new compiled code, using Intel Visual Fortran 10.1, should be faster on a single processor than the V3.5 code compiled on a PC with CVF 6.6.
  3. The code now has OPENMP commands embedded to allow for limited parallelization of some of the routines. Current tests show that going from 1 processor to 2 can result in up to 20-40% speed improvement.
  4. The TKE algorithm has been updated with new algorithms that match experimental tank data for kinetic energy and dissipation.
  5. The roughness height of the water for correction of the vertical velocity wind profile is now a user-defined input, z0.
  6. The Windows user interface no longer uses Array Viewer. The dialog box and PC executable no longer require installation of Array Viewer (which is now obsolete) nor do they need the Array Viewer DLL. The Dialog box has some minor improvements: model run directory displayed and a progress bar.
  7. Fixed error with Algae/chlorophyll a ratio in user manual and fixed pre-processor. The earlier language in the user manual discussed an Algae/Chlorophyll a ratio but presented information that was the ratio of chlorophyll a/algae – this has been revised and fixed in User Manual and in preprocessor.
  8. Spreadsheet output: in earlier versions put in an alphanumeric character as a space for the spreadsheet to preserve the formatting. This was changed to a default value of -99 to facilitate numeric data processing. Also, the “–Depth” output value was changed to just “Depth” since modern plotting programs can reverse an axis.
  9. Preprocessor improvements. Added variable checks for new parameters, fixed bugs, new check for wsc.npt file (not checked in earlier versions).
  10. For the generic constituent, added temperature dependence on 0th order decay and fixed errors in User Manual for units of zero order decay coefficient.
  11. Added the kinetic flux rates to the TSR file output for easier analysis using a spreadsheet of the flux terms for specific locations in the modeled system.
  12. Revised the computation of the drag coefficient for low wind speeds so that the model now agrees better with theory in this wind speed range.
  13. The light extinction coefficient (in m-1) is now included as an output variable in the TSR opt file. Because the model internally computes the light extinction coefficient based on water, SS, POM, algae, zooplankton, and macrophytes, this is an important parameter understanding the internal light transmission predicted by the model. This variable replaces the equilibrium temperature as an output variable.
  14. 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.
  15. 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 37% of NO3-N that diffuses into the sediments becomes incorporated into organic matter in the sediments. The rest is denitrified.
  16. 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.
  17. Added Kinetic flux output that sums up fluxes for all cells of a waterbody at the output frequency specified in the kinetic flux output. The output filenames are called “kflux_jw#.opt” where # is the waterbody number. All active fluxes are summed for the waterbody. This is an important overall diagnostic tool to evaluate the important fluxes in the waterbody over time. Instantaneous fluxes are output in the TSR file for individual cells and a series of fluxes at given segments are shown in the Flux output file which is similar in format as the SNP file. This new file is easy to import into a spreadsheet for analysis.
  18. Gary Hauser sent in a bug fix for the W2 code so that the selective withdrawal algorithm more closely follows the Corps’ code SELECT. The code fix is described below.
    “The problem was that DLRHOMAX (max of DLRHOT and DLRHOB) is used to calculate the relative velocity profile.  We have adjusted it to use DLRHOT and DLRHOB above and below the intake centerline, respectively.  ….  We also replaced the ABS() with a better restriction on the range of VNORM(K) (i.e., 0 to 1).”
  19. Fixed error in ice algorithm that allowed wind to influence the water surface even when ice was present.