Dynamic FORESTBGC .DAY Output File Format


Updated    Mon Feb 04 14:37:31 1991.
By         Joe Glassy, PLDS SWG, Univ. Montana NTSG Group
           (glassy@hps1.umt.edu)

(Note:  See also OTTER SIMULATION NOTES, by Steve Running, at the end
           of this format documentation.)

Overview


This documentation applies to the "Dynamic Model" of FORESTBGC, as 
implemented within the SIMLAT5 (turbo pascal) processor.  The field (column) 
formats shown below assume that the xxxxx.DAY output file is to read as ASCII 
text, fixed (SDF) format file, and that if read using the FORTRAN or C 
formats each field is separated by one implied blank.


Field    - - Column Range - -
Name     Start    End  Length   FOR/C Fmt    Field Description        BGC Var

DAY      001      003    3       I3 , %3d    Julian Day 1..365          z[2]
SNOW     005      011    7       I7 , %7d    Snowpack ppt, in mm        x[1]   
SOIL     013      018    6       I6 , %6d    Soil h20, in mm            x[2]
OUTFLOW  020      025    6       I6 , %6d    Overland flow, in mm       x[3]
EVAP     027      034    8       I8 , %8d    Evaporation, in mm         x[5]
TRANS    036      041    6       I6 , %6d    Transpiration, in mm       x[4]
B        043      048    6       F6.3 %6.3lf Predawn leaf h20 potential g[10]
PSN      050      055    6       I6 , %6d    Net Photosynthesis         x[6]
MRESP    057      063    7       I7 , %7d    Maintenance respiration    g[100]
DECOMP   065      072    8       F8.3 %8.3lf Decomposition ratio        g[8]



Sample of DAY Output


Below is a small sample of a .DAY output file for days 1..4, with a "column
ruler" superimposed after the standard heading and after the data.



 Simlat5 V3.6  6/13/88      Copyright 1988, J.C.Coughlan.
day    snow   soil  outflow   evap  trans    b<     psn   m. resp   decomp

12345678901234567890123456789012345678901234567890123456789012345678901234567890
0--------1---------2---------3---------4---------5---------6---------7---------8
       
  1       0   2350      0        0      0  0.500      0       2    0.537
  2       0   2350    113        4      0  0.500     -8       4    0.537
  3       0   2349    124       13      1  0.500    -11       6    0.551
  4       0   2348    124       16      3  0.500    -14       9    0.565

12345678901234567890123456789012345678901234567890123456789012345678901234567890
0--------1---------2---------3---------4---------5---------6---------7---------8



Sample AWK Program script to Extract header from a .DAY file


For users in computing environments that provide AWK text processor, the
following script is offered, to help reformat the data. The following is an 
AWK program used to extract the text header (1st 7 lines) from a .DAY file, 
prior to ingesting into an SQL relational database management environment.


# extrct1.awk
# updated February 4, 1991
# use this script to extract header text and output as a free field,
# blank delimited ASCII data file.  To run, execute this awk program
# and produce dayfile.dat, enter the following at the operating 
# system prompt:
#     awk -f extrct1.awk dayfile.DAY > dayfile.DAT
{
 if(NR >= 8) {
   print $1,$2,$3,$4,$5,$6,$7,$8,$9,$10
 }
}             


# extrct2.awk
# updated February 4, 1991
# use this script to extract header text and output in a fixed format,
# using the formats described in dayfile.doc.  Note that the formatting
# preserves at least one blank between fields.
{
 if(NR >=8) {
   printf("%3d %7d %6d %6d %8d %6d %6.3lf %6d %7d %8.3lf\n",
                               $1,$2,$3,$4,$5,$6,$7,$8,$9,$10)
 {
}


*****************************************************************************

	O T T E R   S I M U L A T I O N   N O T E S

                    by Steve Running

                       31 Jan 1991

*****************************************************************************

      Here are the draft simulation runs promised for the 4 sites I promised
at the Ames meeting. Santiam did not have a complete climate file so is not
included. Note that the CASCHD, WARINGWD, and SCIO runs all were initiated
with full soil water on Jan 1, and the Metolius with 5.34cm. LAI was defined
(total) of 12 for CASCHD, WARINGWD and SCIO and 4 for Metolius.
 These are lower than measured,
but the model currently over-predicts canopy ppt interception, and possibly
leaf conductance for high LAI stands. The fix for this problem should not be
done casually, as it will help define an important research question. I did
not define site specific specific leaf areas, soil water capacities, stem
and other stand carbon variables etc. There was insufficient time, and in-
complete data to build the simulations that carefully.
         So, treat these results as a first warmup to look at relative
seasonal timing of ecosystem activity, relative differences between sites,
and a practice for all of us at my producing, and you retrieving and analyzing
these simulations. More "offficial" runs will follow.

PS. To translate the .DAY output numbers to recognizable units, all water
variables, when divided by 100 become 1-dimwensional water, cm, as expressed in
hydrology. Ex. a soil water of 534 is actually 5.34cm. All carbon variables,
when divided by 1000, become Mg/ha/yr of CARBON, not biomass!! So, 12000 PSN
is 12Mg/ha/yr photosynthetic production. The B# variable is pre-dawn leaf
water potential in MPA.