Sample submission to the ICC
Notes on the COLA GSWP data tape:
Model: Simplified SiB (SSiB, Xue et al 1991). 0-Dimensional version modified to 2-D for this project.
Contact:
Paul A. Dirmeyer
COLA/IGES
4041 Powder Mill Road, Suite 302
Calverton, MD 20705 USA
(301) 902-1254 fax: (301) 595-9793
dirmeyer@cola.iges.org
Header(1): COLA
File extension: .cola
Wilting point (mm) is a function of soil AND vegetation, calculated as:
WP = n * D*(-exp(c1)/psi_s)**(-1/b)
n = porosity from GSWP Table 1
psi_s = soil tension at saturation (GSWP Table 1)
b = soil wetness exponent (GSWP Table 1)
D = depth of soil (mm)
c1 = from Xue et al 1991, Table 2 (function of vegetation)
Field capacity:
Calculated by saturating the soil at every point and draining for 120 hours (5 days).
Convergence:
This is what I got for the number of unconverged points out of 14637 points:
After year: |
Points unconverged: |
| at 5% or 5 mm |
10% or 10 mm |
| 4 |
3007 |
1118 |
| 5 |
2226 |
883 |
| 6 |
1637 |
795 |
| 7 |
1489 |
673 |
| 8 |
1365 |
653 |
| 9 |
1356 |
641 |
| 10 |
1289 |
560 |
Areas which did not fully spin up were almost all in arid regions
where slow drainage was occuring. In most of these areas the soil
moisture was below wilting point by year 10, so little impact on
surface moisture or energy fluxes is expected.
Known problems and irregularities:
- Month and day in the date/time line in the header (Header(4)) were not updated - but file names indicate the correct times.
- Snapshot fields for 00Z are at the end of the timestep at 23Z, so the forcing data saved are for 23Z.
- There is no differentiation of the form of water storage (liquid versus snow), so snow temperature is always reported the same as ground cover temperature.
- Soil layers 4 and 5 are always reported as -999 (only 3 layers in SSiB)
- Soil ice is -999 for all layers