BROOK90 Version 4 Updates

[Last revised - February 2, 2011]

To the BROOK90 Model page

64-bit Windows - February 2, 2011

Main page and documentation modified to describe required placement and registration of comdlg32.ocx for 64-bit versions of Windows. Thanks to Eloy Benito Reyes for bringing this up.

Version 4.4 - Known Issues

These will be corrected if I ever write a Version 5.

BROOK90 Version 4.4g - November 17, 2010

Documentation - November 17, 2010

Documentation - October 7, 2010

BROOK90 Version 4.4f - September 7, 2010

BROOK90 Version 4.4e - December 14, 2004

BROOK90 Version 4.4e becomes freeware.

Version 4.4d - June 23, 2002

Version 4.4c - May 2, 2002

Change in recommended KF value. Previously I have suggested a hydraulic conductivity of 2 mm/d for field capacity. My recent research suggests this is somewhat low (see Unpublished Results) and that a value of 5 mm/d is better. B90 Help and HTML documentation have been changed to reflect this.

Minor bugs fixed:

Version 4.4b - March 1, 2002

MESFL bug fixed. A recoding error in Version 4.4 and 4.4a produced zero MESFL when CheckData was turned off. EVAL thus failed to run.

WETINF is now forced to be > b/(b+1) and <= 0.999 as required by the Clapp-Hornberger algorithm.

Input solar radiation (SOLRAD) is now prevented from being larger than potential insolation (I0HDAY).

Cancel button. A cancel button has been added to parameter screens to ensure quitting without changes or lockup.

Help modified slightly in various places.

Version 4.4a - January 22, 2002

KPT Graph bug fixed.

Version 4.4 - October 18, 2001

WETINF message fixed. The annoyingly persistent warning message about WETINF when BEXP exceeds 11.5 has been fixed in Version 4.4. The version number has not been changed. The program now tests for WETINF greater than BEXP/(1+BEXP) and resets the value if necessary.

Version 4.4 - October 4, 2001

Interlayer K Weighting and Gradient

All BROOK90 versions through 3.3 and 4.3 used interpolated interlayer conductivities and potential gradients based on layer thickness in subroutine VERT. The original equations produced too little flow with a thick wet layer over a thin dry one. The fix attempted in Versions 3.3 and 4.2/4.3 reversed the problem and gave too little flow with a thin wet layer over a thick dry one. The solution seems to be to treat the interface as if it cannot "see" the thickness of a thick layer. So the equations now are:

KKMEANi = exp((ln(KKi) + ln(KKi+1)) / 2)

GRADi = (PSITi - PSITi+1) / RMIN(THICKi, THICKi+1)

Results do not differ among versions for simulations in which all layers have the same thickness. Differences are generally small unless both layer thickness and dryness differ greatly.

International awareness - comma as decimal. I have recoded parts of BROOK90 so it should now be internationally aware, responding properly to changing Windows-ControlPanel-RegionalSettings. The earlier incomplete test for a comma as decimal has been improved. Data files now may have either comma or period as the decimal point with space separators, or may be a .CSV file with period as decimal and comma as separator. Both saved and printed output will conform to the Windows RegionalSetting. However, parameter files must always include a comma separator after the parameter name and decimal periods, as in previous versions. (Thanks to Ute Assmann for raising this issue.)

Output-Append added. A new option "Append" has been added to the Output menu for each interval. This appends the current output to any existing file, whereas "Save" requires creation of a new file. "Append" asks whether the header line should be included or not.

Font selection. Omitted appearance of SelectPrinterFont window on startup when CourierNew 10 was not available. Initial font and size is thus the system default. The printer font is now only tested when a print is requested. (Thanks to Ute Assmann)

Minor errors fixed.

Documentation Corrections - July 21, 2001

On page B90-6 "Darcy-Richards equation" should be "Darcy equation".

On page KPT-3 under FDPSIDW, the subscripts in the first equation should all be f, not x.

On page EVL-2, N in the MDE equation should be n (two places).

On page PET-11 in equation 34, Lp/(Lp+Sp) in the first term should be replaced by 1/F and (Lp+Sp) in the second term by FLp, where F = (Lp+Sp/2)/Lp. This arises from the (unnecessary?) assumption that only half of Sp is distributed proportionally to L' (and the other half is below the leaves in a "stem space").

On page KPT-3 under SOILVAR, ψm should be ψ, and three lines higher AWATRZ should be AWAT.

On page WAT-3 under SRFLFR, the second paragaph should just say that BROOK90 effectively simulates the mean water budget and mean water content over the whole watershed.

On page WAT-4 in the SAFRAC paragraph, all QFF should be changed to QFFC.

On page SNO-3 in the middle, "resistances ra and rs" should be "resistances raa and ras"

Version 4.3 - July 19, 2001

DENSEF and CS parameters swapped. The parameter CS has been changed from a Fixed parameter to a Canopy parameter. This clarifies that CS, which is the ratio of stem area index to height, may need to change for different canopy types. The parameter DENSEF has been changed from a Canopy parameter to a Fixed parameter. Normally DENSEF should be set to 1 and MAXLAI, CS, MXRTLN, and MXKPL values should represent the actual canopy being simulated. DENSEF should only be used to simulate canopy thinning, for which multiplying the other four parameters by DENSEF < 1 is appropriate. Earlier Version 4 Canopy and Fixed parameter files will not work with Version 4.3. The files can be changed manually by replacing the DENSEF line in Canopy files with CS (normally = 0.035), and replacing the CS line in Fixed files with DENSEF (normally = 1.0). Help for DENSEF and the parameters to which it applies has been modified.

Daily precip interval. A Precipitation Interval file with only one line per day is now acceptable. This allows a user to provide different sequences of daily precipitation as precip interval files without having to change the Data file.

Errors fixed. Some Help messages were not properly changed for Version 4; at least some of these are now corrected. A few cosmetic errors are fixed. ROOTDEN is now only saved to two decimal places, the same as it is displayed.

HTML Documentation - July 7, 2001

Documentation of BROOK90, written in HTML, for Versions 4, 3.3, and 3.2x became available on the web. Documentation for Version 3.1 became obsolete and should not be used.

Documentation Error - April 27, 2001

On page SUN-3 about half way down the page, it should say that 30 to 50 W/m2 is equivalent to 1.1 to 1.8 mm/d of evaporated water.

Version 4.2 - April 2, 2001

Versions 4.2/4.3 changed the K weighting but badly - see Version 4.4 above. They should be replaced by Version 4.4.

32-bit BROOK90 Version 4.1 - February 2, 2001

Finally I have upgraded from Visual Basic 3 (16-bit) to Visual Basic 6 (32-bit). BROOK90 Version 4.1 is a 32-bit compilation of the Visual Basic 3 Version 4.0a code. Version 4.1 will run under Windows 95/98/NT.This did not require ANY code changes!! It seems to run about twice as fast as Version 4.0a, especially when Run-time Graph is turned off.

Version 4.0a

. Comma-separated Values (.csv) files saved from spreadsheets like Excel can be used for BROOK90 data files. The .csv file must have the correct 10 columns for daily input files and 6 columns for precipitation interval files. The .csv file must not remain open in Excel or BROOK90 won't be able to access it.

Some bugs involving printing and cancelling printing have been fixed

Version 4.0 - June 15, 2000

Version 4.0 separates the single parameter file into six files (Location, Flow, Canopy, Soil, Fixed, and Initial) and facilitates batch processing using "Shortcut for Windows" software.

Catchment parameters are separated into Location parameters and Flow parameters. Root distribution has been separated from the Soil layers and included with Canopy parameters. Initial values are now all in one place. I chose to leave the Fixed Parameter file separate rather than to include those parameters in the other files. Fixed parameters are not checked for abnormal values. The default extension for data and parameter files is .TXT.

The phenology of RELLAI and RELSAI is included with Location parameters but also depends on Canopy type; this means that a separate Location Parameter file is needed for each cover type at a location.

QLAYER and ILAYER in Version 3 have been replaced by QDEPTH and IDEPTH in Version 4. This makes these parameters independent of the specification of soil layers.

QFPAR = 0 (<0.01) now gives a bucket model with no SRFL or BYFL when QDEPTH water is below field capacity and all excess becoming SRFL or BYFL when above field capacity. QFFC is not used.

The documentation uses z0g when it really means either 1) z0g , 2) CZS * HEIGHT, or 3) z0s. Similar confusion existing in the program has been fixed.

Visual Basic code has been considerably changed, hopefully for the better.

To facilitate batch processing, menus can be operated with Alt-characters and parameter screens have a Tab order. Shortcuts have been added to the menu for those who like Function keys

Menus are somewhat changed.

Output-Save becomes checked if output has been saved.

All printer output now has top and side margins.

The EVAL output screen has been changed so it can be saved.

Version 3

Earlier updates can be found at BROOK90 Version 3 Updates.

[Compass Brook logo]

COMPASS BROOK