Instruction manual

36
SAPFLOWV(J,ITOD) = FLUXV(J)
SAPFLOWF(J,ITOD) = FLUXF(J)
c GOTO 299
C------------------------------------------------------------
C ... compute sum of sap velocity times relative sapwood area
C------------------------------------------------------------
R1P = SWR(J)
R1M = (R(1)+R(2))/2.
R2M = (R(2)+R(3))/2.
R3M = (R(3)+R(4))/2.
R4M = HWR(J)
A1=(R1P*R1P-R1M*R1M)*PI
A2=(R1M*R1M-R2M*R2M)*PI
A3=(R2M*R2M-R3M*R3M)*PI
A4=(R3M*R3M-R4M*R4M)*PI
SAPFLOWS(J,ITOD) = A1*SV(1)+A2*SV(2)+A3*SV(3)
SAPFLOWS(J,ITOD) = (SAPFLOWS(J,ITOD)+A4*SV(4))/1000.
299 CONTINUE
C------------------------------------------------------------
c ... print out the results
C------------------------------------------------------------
300 CONTINUE
301 CONTINUE
WRITE(31,995) T(ITOD), (SAPFLOWV(J,ITOD) , J=1,NHPV)
C * ,(SAPFLOWV(1,ITOD)+SAPFLOWV(2,ITOD))/2.
WRITE(32,995) T(ITOD), (SAPFLOWF(J,ITOD) , J=1,NHPV)
C * ,(SAPFLOWF(1,ITOD)+SAPFLOWF(2,ITOD))/2.
WRITE(33,995) T(ITOD), (SAPFLOWS(J,ITOD) , J=1,NHPV)
C * ,(SAPFLOWS(1,ITOD)+SAPFLOWS(2,ITOD))/2.
GOTO 100
1000 CONTINUE
999 FORMAT(2A40)
996 FORMAT(5x,'Time [d] = ',F8.3,' Sap flow [L/h] = ',8F8.3
* ,6X,1H.)
995 FORMAT(1X,F8.3,8F8.3)
C------------------------------------------------------------
CALL SETTEXTPOSITION(33 ,35,CURPOS)
WRITE(*,*) '***** Click <MOUSE> to END *****'
CALL MOUSECLICK(IX4,IY4)
CALL CLEARSCREEN($GCLEARSCREEN)
C------------------------------------------------------------
STOP
END
c============================================================
SUBROUTINE REGRESS(N, X, Y, A0, A1, A2)
C ... fit a parabola through the SFD data using N data points
C of the form ... y = A0 + A1.x + A2.x^2
c============================================================
IMPLICIT NONE
REAL X(4), Y(4), A0, A1, A2
REAL SX, SY, SX2, SXY, SX2Y, SX3, SX4
REAL SSX2, D
INTEGER N, I
IF(N.GT.2) THEN
SX = 0
SY = 0
SX2 = 0
SXY = 0
SX3 = 0
SX2Y = 0
SX4 = 0
DO 100 I=1,N
SX = SX + X(I)
SY = SY + Y(I)
SXY = SXY + X(I)*Y(I)
SX2 = SX2 + X(I)*X(I)
SX2Y= SX2Y+ X(I)*X(I)*Y(I)
SX3 = SX3 + X(I)**3
SX4 = SX4 + X(I)**4
100 CONTINUE
SSX2= SX2
SX4 = SX4 - SX2*SX2/N