Laserfiche WebLink
<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 />