Laserfiche WebLink
<br />c NTOT(IZ)=NTOT(IZ)+l <br />c 30 CONTINUE !! azimuth loop, M <br />c ENDIF ! elevation limit test <br />c 35 CONTINUE !! range loop, L <br />c 40 CONTINUE !! tilt loop, N <br />c RETURN <br />C <br />C Output image <br />c ENTRY PROFDUMP <br />C Label the graph <br />c IYR=IDATE/lOOOO <br />c lMAG(l,l)=IYR/lOO <br />c IMAG(1,2)=MOD(IYR,100) <br />c IMAG(1,3)=MOD(IDATE/IOO, 100) <br />c IMAG(1,4)=MOD(IDATE,100) <br />c lMAG(1,6)=ITIME/10000 <br />c IMAG(1,7)=MOD(ITIME/lOO,100) <br />C Fill the graph <br />c DO 60 J=1,108 ! lines and heights <br />c IF(J.GE.5.AND.J.LE.I05)THEN ! in graph interior <br />c IZ=106-J ! reverse order, prints from top to bottom in altitude <br />c KSUM=O <br />c LASTK=4 <br />c DO 45 ID=4l,1,-1 <br />c KSUM=KSUM+NSUM(IZ,ID) ! cumulative sum <br />c NK=4 <br />c IF(NTOT(IZ) .GT.O) <br />c + NK=4+NINT(200.*(FLOAT(KSUM)/FLOAT(NTOT(IZ))))! cumulative percent *2 <br />c IF(NK.GT.204)NK=204 ! protection against array overflow <br />c IF (LASTK.LT.NK) THEN <br />c DO 42 I=LASTK+l,NK <br />c 42 IMAG(J,I)=ID-l! write dBZ values into image <br />c LASTK=NK <br />c ENDIF <br />c 45 CONTINUE! dBZ loop, ID <br />c ENDIF ! graph interior test <br />C Output the graph <br />c DO 55 1=1,208 ! columns <br />c 55 LINE(I)=IMAG(J,I) <br />c 60 WRITE(l)LINE <br />c <br />C Initialize sums <br />c ENTRY PROFINIT <br />c DO 62 IZ=l,lOl <br />c DO 61 ID=1,4l <br />c 61 NSUM(IZ,ID)=O <br />c 62 NTOT(IZ)=O <br />C Initialize upper and lower image edges <br />c DO 65 J=1,4 ! lines <br />c DO 63 1=1,208 ! columns <br />c IMAG(104+J,I)=-1 <br />c 63 IMAG( J,I)=-l or 255 white <br />c DO 64 1=5,205,20 <br />c IMAG(104+J,I)=50 will color black <br />c 64 IMAG( J,I)=50 cumulative percent tics at 10% intervals <br />c 65 CONTINUE ! then upper and lower edges finished <br />C Initialize side image edges <br />c DO 68 1=1,4 ! columns <br />c DO 66 J=5,104 ! lines <br />c IMAG(J,I+204)=-1 <br />c 66 IMAG(J,I )=-1 or 255 white <br />c DO 67 J=4,104,10 <br />c IMAG(J,I+204)=50 will color black <br />c 67 IMAG(J,I )=50 elevation tics at 0.5 km intervals <br />c 68 CONTINUE then right and left edges finished <br />c RETURN <br />c END $$$ end of temporary addition <br /> <br />(lines 2091-2098, in SUBROUTINE HINDSUM) <br /> <br />IF (O.LT.ITOP) CALL TOPHOUR !! if doing any top-of-hour files <br />C CALL PROFDUMP ! $$$ and PROFINIT as well to output, initialize image <br />ENDIF <br />IYMDHB4 = NDATSTRT*lOO+NEWHR <br />IF (O.LT.IENDSUM) RETURN !! special HENDSUM call at end <br />ENDIF <br />C update the storm and 3-hr and 1-hr totals <br />DO 25 M=1,360 !! azimuth loop <br /> <br />22 <br />