Instruction manual

32
For well-behaved velocity profiles (i.e. small curvature at large radii) all three methods
yield similar results, but for profiles where the curvature near the cambium method 3 is
recommended (Hatton et al, 1990).
6.4. Listing of the ANALYSIS program
The following is a listing of the FORTRAN program HPV-2000.FOR used to analyze the
tz-data collected by the HPV data logger. The analysis procedure follows that outlined in
the Background Section of this document. This is a fairly big program – only because it
includes routines for graphical output to the screen.
c============================================================
c ... Program HPV2000.FOR to convert tz data to correct sap
c ... velocities using numerical simulation results for a
c ... given wound width and probe spacings of (-0.5,0,1.0)
C ...........................................................
C
C ... ------ LAST modified AUG, 2000, S.R. Green ------
C
c============================================================
USE DFLIB
IMPLICIT NONE
c include files for MS Fortran graphics
CHARACTER*50 CTITLE, ans
REAL*8 NXPIX, NYPIX
REAL*8 X1(101), Y1(101), X2(101), Y2(101)
REAL*8 X1MIN, X1MAX, Y1MIN, Y1MAX
INTEGER IX4, IY4, CWIN/0/
LOGICAL STATUS, END_CALC
TYPE (RCCOORD) CURPOS
TYPE (WINDOWCONFIG) WINC
TYPE (QWINFO) QWIN
REAL TZ(4,4), HPV(4), PD(4,4), PS(4,4)
REAL TMIN, TMAX, EMIN, F0, F1, F2, BEGYR
REAL R1P, R1M, R2M, R3M, R4M, A3, A4
INTEGER IDAY, ITOD, IJ
REAL WC1(5,7), WC2(5,7)
REAL R(4), SV(4), SF(4)
REAL SVIJ(16),T(100), Y(100)
REAL SAPFLOWV(4,100), SAPFLOWF(4,100), SAPFLOWS(4,100)
REAL SWR(4), HWR(4), WW(4), EMAX(4)
REAL VFWOOD, VFWAT, TOD
REAL A0, A1, A2, FLUXVEL, FLUXFLX, FLUXV(4), FLUXF(4)
REAL ID, MI, PI, TWOPI
REAL CHANID,DAY,HRMIN, HR, MIN
INTEGER NHPV, N, I, ICORR, IW(4), J, NCOL(4), NTC(4), NC
CHARACTER*40 FILNAM
CHARACTER*4 STR
LOGICAL SWANCOR
INTEGER*2 DUMMY2, FGD
INTEGER*4 DUMMY4, BGD
COMMON WINC
COMMON QWIN
C-------------------------------------------------------------
c ... Set the wound corrections of Swanson and Whitfield (1983)
C-------------------------------------------------------------
DATA WC1/ 0, 0, 1, 0, 0.0,
* 16, 0.393, 1.356, 0.036, 0.0,
* 20, 0.807, 1.203, 0.058, 0.0,
* 24, 1.184, 1.072, 0.087, 0.0,
* 28, 1.524, 0.964, 0.124, 0.0,