<br />29
<br />C apply
<br />
<br />DO 30 L=4,230
<br />DO 30 L=4,LIMITVUE
<br />IP = 0 !! 97/9
<br />N = HYBRD(L)
<br />IF (JAN(N) .EQ.-99.) GOTO 29 !! skip vacant azimuths
<br />line would be .GT.1 but Ra scaled *128 for fine adjustment
<br />IF (NDN5(JAN(N),L+NDN(N)) .GT.128 .AND. ~.LE.NDN(N+1)-NDN(N))
<br />IP = RATE_TABLE(ISHFT(NDN5(JAN(N),L+NDN(N)),-7))
<br />FP = FLOAT (IP) *SUMT !! rnm * 100
<br />range correction: default RC1=1, RC2=RC3=O for no range cor.
<br />IF (FP.GT.O. .AND. RC4.LE.FLOAT(I")) THEN! cor. for RC4 <= range
<br />IF (INFLAG. EQ .1.0R. L. LE. INSIDE) ~rHEN ! minimal virga hazard
<br />FP=FP*(RC1+RC2*FLOAT(L)+RC3*FLOAT(L*L))
<br />ELSE
<br />FP=O.
<br />ENDIF
<br />ENDIF
<br />PPI2(L)
<br />CONTINUE
<br />
<br />! radius loop; avoid first 3 (was 2) km from radar
<br />! radius loop; avoid first 3 (was 2) km from radar @@@
<br />-Ra: Dr. Super said for nids change from L=3,230
<br />
<br />C
<br />
<br />C next
<br />
<br />@@@
<br />
<br />@@@
<br />! zero it beyond virga box edge if little precip inside @@@
<br />! @@@
<br />
<br />30
<br />
<br />NINT(FP) !! now save nearest integer result
<br />!! end of rad:lus loop, range L
<br />
<br />(SUBROUTINE HINDSUM: lines 2009-2010)
<br />
<br />COMMON/RESTRICT/LIMITVUE,INSIDE,TOPIIJ,BOTIN,INSUM,INTOT,INFLAG, @@@
<br />INQUALIF,DURATION,FRACTION @@@
<br />
<br />(SUBROUTINE HINDSUM: lines 2097-2119, range restriction)
<br />
<br />C update the storm and 3-hr and 1-hr totals
<br />DO 25 M=1,360 !! azimuth loop
<br />Call 230 km, precip in interval, mm*100
<br />READ(9, REC=NVS+M) PPI2
<br />READ(10, REC=M) PTEMP !! l:ead storm total precip
<br />C DO L=l, 230 !! radius loop, calculate depths
<br />DO L=l,LIMITVUE !! radius loop, calculate depths @@@
<br />PTEMP(L) = PTEMP(L)+PPI2(L) !! add to storm total precip
<br />DPI2 (L) = NINT (FLOAT (PPI2 (L)) *FLUFFNIlO.)
<br />ENDDO !! depth in interval, mrn*10
<br />WRITE (10, REC=M) PTEMP !! write storm total precip
<br />IF (SNSTD.NE.' ') THEN !! :If want snow depth files
<br />READ (11, REC=M) DTEMP !! storm total depth
<br />C DO L=1,230 !! add to storm total depth
<br />DO L=l,LIMITVUE !! add to storm total depth @@@
<br />DTEMP(L) = DTEMP(L)+DPI2(L)
<br />ENDDO
<br />WRITE(ll, REC=M) DTEMP !! ~~rite storm total depth
<br />ENDIF
<br />READ(9, REC=K3P+M) PTEMP !! l:ead 3-hr precip
<br />READ(9, REC=K3D+M) DTEMP !! l:ead 3-hr depth
<br />C DO 22 L=1,230 !!
<br />DO 22 L=l, LIMITVUE !! G!@@
<br />(SUBROUTINE HINDSUM: lines 2026-2127, another ran\1e limit)
<br />
<br />C
<br />
<br />DO 23 L=1,230
<br />DO 23 L=l,LIMITVUE
<br />
<br />! !
<br />!! G!@@
<br />
<br />(SUBROUTINE HINDSUM! lines 2222-2223, another range limit)
<br />
<br />C
<br />
<br />DO 30 L=1,230 ! radius loop
<br />DO 30 L=l,LIMITVUE ! radius loop @@@
<br />
<br />(SUBROUTINE HINDSUM: lines 2262-2263, another range limit)
<br />
<br />C
<br />
<br />DO 40 L=1,230
<br />DO 40 L=l,LIMITVUE
<br />
<br />!! radius loop
<br />!! radius loop @@@
<br />
<br />(SUBROUTINE FULLHOUR, lines 2443-2444)
<br />
<br />COMMON/RESTRICT/LIMITVUE,INSIDE,TOPI1J,BOTIN,INSUM,INTOT,INFLAG, @@@
<br />INQUALIF,DURATION,FRACTION @@@
<br />
<br />(SUBROUTINE FULLHOUR, lines 2487-2508, two more range restrictions)
<br />
<br />IF (1.LT.ITOP) THEN !! if want 3-hr files
<br />DO M=1,360 !! all azimuths
<br />READ(9, REC=K3P+M) PPI2
<br />DO L=1,230
<br />DO L=l,LIMITVUE
<br />C extrap. to full 3-hour period, precip
<br />PPI2(L) = NINT(F*FLOAT(PPI2(L):)
<br />ENDDO
<br />WRITE(16, REC=M) PPI2
<br />
<br />!! (ill 230 km
<br />!! all ranges
<br />!! (ill ranges @@@
<br />
<br />C
<br />
<br />!! 1111 230 km
<br />
<br />29
<br />
|