c ... vfwood = volume fraction of wood = 0.34
c ... vfwat = volume fraction of water = 0.56
C------------------------------------------------------------
WRITE(*,*) 'Enter the volume fraction of wood '
READ(30,1098) ANS,VFWOOD
WRITE(*,*) 'Enter the volume fraction of water '
READ(30,1098) ANS,VFWAT
C------------------------------------------------------------
c ... mean sap wood radius in cm
C------------------------------------------------------------
WRITE(*,*) 'Enter the sap-wood radius of each stem/root [cm] '
READ(30,1098) ANS, (SWR(J), J=1,NHPV)
WRITE(*,1098) ANS, (SWR(J), J=1,NHPV)
WRITE(*,*) 'Enter the heartwood radius of each stem/root [cm] '
READ(30,1098) ANS, (HWR(J), J=1,NHPV)
WRITE(*,1098) ANS, (HWR(J), J=1,NHPV)
WRITE(*,*)
C------------------------------------------------------------
c ... read in the probe depths
C------------------------------------------------------------
WRITE(*,*) 'Enter the thermistor depths of each probe [cm] '
DO 30 J=1,NHPV
READ(30,1098) ANS, (PD(I,J),I=1,NTC(J))
WRITE(*,1098) ANS, (PD(I,J),I=1,NTC(J))
READ(30,1098) ANS, (PS(I,J),I=1,NTC(J))
WRITE(*,1098) ANS, (PS(I,J),I=1,NTC(J))
DO 30 I=1,NTC(J)
30 PD(I,J) = ABS(SWR(J)-PD(I,J))
WRITE(*,*)
C-------------------------------------------------------------
c ... read plot range
C-------------------------------------------------------------
WRITE(*,*) 'Enter DOY range: Tmin, Tmax'
READ(30,1098) ANS, TMIN,TMAX
WRITE(*,1098) ANS, TMIN,TMAX
WRITE(*,*) 'Enter SAPFLOW range: Emin, Emax'
READ(30,1098) ANS, (EMAX(I), I=1,NHPV)
WRITE(*,1098) ANS, (EMAX(I), I=1,NHPV)
EMIN = 0.0
CALL SETTEXTPOSITION(33 ,35,CURPOS)
WRITE(*,*) '***** Click <MOUSE> to END *****'
CALL MOUSECLICK(IX4,IY4)
CALL CLEARSCREEN($GCLEARSCREEN)
STATUS = SETCOLOR(8)
STATUS = RECTANGLE($GFILLINTERIOR,0,0,NXPIX,NYPIX)
STATUS = SETCOLOR(2)
c------------------------------------------------------------
c ... main Loop
c------------------------------------------------------------
IDAY=0
100 READ(30,*,END=1000) CHANID,DAY,HRMIN,
* ((TZ(I,J),I=1,NTC(J)),J=1,NHPV)
IF(BEGYR.LT.0) THEN
TMIN = DAY
BEGYR = 9999
ENDIF
IF(HRMIN.LT.5) THEN
IDAY=IDAY+1
DO 120 J=1, NHPV
DO 110 I=1,ITOD
X1(I) = T(I)
Y1(I) = SAPFLOWV(J,I)
110 CONTINUE
X1MIN = TMIN !T(1)
X1MAX = TMIN+5. !T(ITOD)
Y1MIN = EMIN
Y1MAX = EMAX(1)
CALL NEWONEPLOT(X1,Y1,ITOD,'Day of Year', ''
&,'Water_Use_[L/h]','', X1MIN,X1MAX,Y1MIN,Y1MAX,.TRUE.,.FALSE.,8+J)
120 CONTINUE
(4 pages)
Manymanuals.com
Manymanuals.de
Manymanuals.fr
Manymanuals.it
Manymanuals.pl
Manymanuals.cz
Manymanuals.es
Manymanuals-pt.com
Comments to this Manuals