Instruction manual

35
ITOD = 1
IF(IDAY.GT.5) THEN
IDAY = 1
TMIN = TMIN+5
TMAX = TMAX+5
CALL SETTEXTPOSITION(33 ,35,CURPOS)
WRITE(*,*) '***** Click <MOUSE> to CONTINUE *****'
CALL MOUSECLICK(IX4,IY4)
CALL CLEARSCREEN($GCLEARSCREEN)
STATUS = SETCOLOR(8)
STATUS = RECTANGLE($GFILLINTERIOR,0,0,NXPIX,NYPIX)
STATUS = SETCOLOR(2)
ENDIF
ELSE
ITOD = ITOD+1
ENDIF
C------------------------------------------------------------
c ... convert hrmin to time of day
C------------------------------------------------------------
HR = INT(HRMIN/100)
MIN = HRMIN - 100*HR
TOD = DAY + (HR+MIN/60)/24
T(ITOD) = TOD
C------------------------------------------------------------
c ... For each set of probes find the corrected SV profile,
C ... then compute flux using three different methods
C------------------------------------------------------------
DO 300 J=1,NHPV
DO 200 I=1,NTC(J)
IF( TZ(I,J).gt.(0.0) ) THEN
HPV(I) = PS(I,J)/(2.*TZ(I,J))*3600.
ELSE
HPV(I) = 0.0
ENDIF
C------------------------------------------------------------
c ... compute the corrected SV
C------------------------------------------------------------
IJ=(J-1)*NTC(J)+I
IF(SWANCOR) THEN
SV(I) =(WC1(5,IW(j))*HPV(I)*HPV(I)*HPV(I)
* + WC1(4,IW(j))*HPV(I)*HPV(I)
* + WC1(3,IW(j))*HPV(I)
* + WC1(2,IW(j)))*(0.505*VFWOOD + VFWAT)
ELSE
SV(I) =(WC2(5,IW(j))*HPV(I)*HPV(I)*HPV(I)
* + WC2(4,IW(j))*HPV(I)*HPV(I)
* + WC2(3,IW(j))*HPV(I)
* + WC2(2,IW(j)))*(0.505*VFWOOD + VFWAT)
ENDIF
R(I) = PD(I,J)
SVIJ(IJ) = SV(I)
SF(I) = SV(I)*R(I)
200 CONTINUE
N = NTC(J)
CALL REGRESS(N, R, SV, A0, A1, A2)
CALL REGRESS(N, R, SF, F0, F1, F2)
C------------------------------------------------------------
c ... compute total flow, making sure that it remains positive !!!
C------------------------------------------------------------
FLUXV(J) = 0.0
FLUXF(J) = 0.0
ID = HWR(J)-0.01
DO WHILE(ID.LT.SWR(J))
ID = ID + 0.01
MI = ID + 0.005
FLUXVEL = TWOPI*MI*(A0 + A1*MI + A2*MI*MI)*0.01
IF(FLUXVEL.GT.(0.0)) FLUXV(J) = FLUXV(J) + FLUXVEL
FLUXFLX = TWOPI*(F0 + F1*MI + F2*MI*MI)*0.01
IF(FLUXFLX.GT.(0.0)) FLUXF(J) = FLUXF(J) + FLUXFLX
END DO
FLUXV(J) = FLUXV(J)/1000
FLUXF(J) = FLUXF(J)/1000