|
<br />COMMON/RESTRICT/LIMITVUE,INSIDE,TOPIN,BOTIN,INSUM,INTOT,INFLAG, @@@
<br />INQUALIF,DURATION,FRACTION @@@
<br />
<br />(SUBROUTINE HINDZERO: lines 970-973)
<br />
<br />DO 10 L=1,230
<br />PTOT(L) = 0 !! zero the precip total at all ranges
<br />IF(L.GT.LIMITVUE)PTOT(L)=-l ! no confidence in results @@@
<br />10 CONTINUE
<br />
<br />(SUBROUTINE NOSPIKES: lines 1586-1587)
<br />
<br />COMMON/RESTRICT/LIMITVUE,INSIDE,TOPIN,BOTIN,INSUM,INTOT,INFLAG, @@@
<br />INQUALIF,DURATION,FRACTION @@@
<br />
<br />(SUBROUTINE NOSPIKES: lines 1606-1608, top of executables)
<br />
<br />C
<br />
<br />Set virga detection parameters
<br />INSUM=O
<br />INTOT=O
<br />
<br />@@@
<br />@@@
<br />@@@
<br />
<br />(SUBROUTINE NOSPIKES: lines 1616-1629)
<br />
<br />DO 80 N=1,4 !! step thru all four tilts (if present)
<br />IF (JAT(N) .EQ.O) GOTO 80 !! skip if no data for tilt N
<br />C Do a virga check
<br />DO 5 IBIN=2,INSIDE !! for nearly all bins within virga
<br />Z=BEAMHT(FLOAT(IBINl,ELT(Nl)
<br />IF(Z.GE.BOTIN.AND.Z.LE.TOPINlTHEN inside virga box
<br />DO 3 M=l,JAT(N) ! all azimuths at tilt
<br />INTOT=INTOT+1 ! sum all bins in box
<br />NDBZ = NDN5(M,IBIN+NDN(N)) ! get bin reflectivity
<br />IF (ISHFT(NDBZ,-7l.GE.NMIN) INSUM=INSUM+1 ! sum precip
<br />3 CONTINUE ! end of azimuth loop
<br />ENDIF ! altitude limit check
<br />5 CONTINUE! end of range loop
<br />C Isolation testing
<br />
<br /> @@@
<br />box @@@
<br /> @@@
<br /> @@@
<br />N @@@
<br /> @@@
<br /> @@@
<br />bins @@@
<br /> @@@
<br /> @@@
<br /> @@@
<br />
<br />(SUBROUTINE NOSPIKES: lines 1755-1765)
<br />
<br />C
<br />
<br />80 CONTINUE !! tilt loop, N
<br />Evaluate virga box contents
<br />IF(INTOT.GT.O)THEN
<br />F=FLOAT(INSUMl/FLOAT(INTOTl
<br />IF(F.GE.FRACTIONlTHEN ! found
<br />INFLAG=l
<br />INQUALIF=IENDSCAN ! end time
<br />ENDIF ! qualifying test
<br />ENDIF ! zero-divide protection
<br />RETURN
<br />END
<br />
<br />fraction of precip echoes in box
<br />precip within virga box
<br />
<br />of this volume scan
<br />
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />
<br />(SUBROUTINE PRECIP: lines 1874-1886)
<br />
<br />C /RESTRICT/ variables for virga reductions: @@@
<br />C LIMITVUE 1*4 Maximum range for precip accumulations @@@
<br />C INSIDE 1*4 Range within which significant echo must occur @@@
<br />C TOPIN R*4 Altitude (km) above radar for top of qualifying box @@@
<br />C BOTIN R*4 Altitude (km) above radar for bottom of qualif. box @@@
<br />C INSUM 1*4 Range bin count of echoes >=10 dBZ within box @@@
<br />C INTOT 1*4 Range bin count of all data within qualifying box @@@
<br />C INFLAG 1*4 0, else 1 if enough precip echoes inside box @@@
<br />C INQUALIF 1*4 Last qualifying time, hhmmss @@@
<br />C DURATION R*4 Duration, fracto hrs, allowed since last INFLAG=l @@@
<br />C FRACTION R*4 Qualifying fraction, INSUM/INTOT, to set INFLAG @@@
<br />COMMON/RESTRICT/LIMITVUE,INSIDE,TOPIN,BOTIN,INSUM,INTOT,INFLAG, @@@
<br />INQUALIF,DURATION,FRACTION @@@
<br />
<br />(SUBROUTINE PRECIP: lines 1933-1939)
<br />
<br />CALL TIMESPAN(JTIMSTRT,JDATSTRT,JTIMSTOP,JDATSTOP,SPANT) time
<br />C Test elapsed time for last non-virga in virga box
<br />Vo = IHMS2S(INQUALIF)/3600. !! last time INFLAG set
<br />VT=T2-VO ! fractional hours
<br />IF(VT.LT.OlVT=VT+24.
<br />IF(VT.GT.DURATION)INFLAG=O ! no precip in virga box; stop accum.
<br />C increment volume scan index: IVS is 0 to 39
<br />
<br />period
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />@@@
<br />
<br />(SUBROUTINE PRECIP: lines 1968-1987)
<br />
<br />C calculate precipitation contribution, IP
<br />
<br />28
<br />
|