|
<br />APPENDIX A. TEMPORARY CODING TO
<br />CREATE VERTICAL PROFILE GRAPHS
<br />
<br />In an effort to understand the variations in the vertical profile of reflectivity as a stonn passes through an
<br />area, changing from virga to precipitation to vilrga, RADAR11.FOR was modified with additional coding,
<br />presented below. (The line numbers are from the experimental version of the program, but enough
<br />neighboring statements are included for location. All extra lines are commented out and many are
<br />identified with a $$$ identifier.)
<br />
<br />(the following lines through 658 are from the main PROGRAM RADARll)
<br />
<br />(lines 50-52)
<br />
<br />C 98/12 Ed: Temporary vertical profile 01' reflectivity for virga study.
<br />CAll corrunented out. Can find it by looking for $$$ on single lines and
<br />C within the PROFILE subroutine.
<br />
<br />(lines 137-139)
<br />
<br />C
<br />C
<br />
<br />COMMON/PROF/NSUM(101,41),NTOT(101) ,LINE(208) ,IMAG(108,208)
<br />BYTE LINE,IMAG ! $$$
<br />COMMON/RADIAL/ VAL(6), NUM(20) , DBZ(O:459), IFLAG
<br />
<br />$$$
<br />
<br />(lines 210-212)
<br />
<br />C
<br />C
<br />
<br />OPEN(l,FILE='PROFILE.BYT',STATUS='UNF~OWN',FORM='UNFORMATTED')
<br />CALL PROFINIT ! to initialize profile sums and image $$$
<br />CALL STARTUP
<br />
<br />$$$
<br />
<br />(lines 336-339)
<br />
<br />*
<br />
<br />22 Vol-scan starts were accepted.
<br />CALL PROFDUMP ! $$$ and PROFINIT as well to output, initialize image
<br />STOP '(passing wanted end file)'
<br />ENDIF
<br />
<br />C
<br />
<br />(lines 656.-658)
<br />
<br />C
<br />
<br />CALL OCCULT
<br />CALL NOSPIKES
<br />CALL PROFILE
<br />
<br />!! apply occult corrections
<br />!! try removing lingering clutter & blemishes
<br />! $$$ to make vertical profile image; temporary
<br />
<br />(lines 1767-1864)
<br />
<br />c SUBROUTINE PROFILE ! $$$ and followi.ng to END statement
<br />C makes an image of the vertical profile of reflectivity, contoured by
<br />C frequency of dBZ values at each height to 5.0 km
<br />c COMMON/ANGLES/ AZT(4,370), ELT(4), SHCELT(4) , JAT(4), TILTLAST
<br />c COMMON/CLOCK/ IDATE, ITIME, KDATE, K1'IME, MDATE, MTIME,
<br />c IDATl, ITIM1, ITIM9, ITIMSC
<br />c COMMON/PROF /NSUM (101,41) ,NTOT (101) ,LINE (208) , IMAG (108,208) ! $$$
<br />c BYTE LINE,IMAG ! $$$
<br />c COMMON/STORE/ N005(370,230) , ND15(370,230), ND25(370,230),
<br />c ND35(370,230), NDN(5)
<br />c INTEGER*2 N005, ND15, ND25, ND35, NDN5(370,920)'
<br />c EQUIVALENCE (ND05,NDN5)
<br />C calculates altitude, km msl, of radar beam from range, R, elevation, E
<br />c BEAMHT(R,E)=R*SIND(E)+5.8869E-5*R*R*COSD(E)*COSO(E)
<br />c SIND(A)=SIN(A*.017453293)
<br />c COSD(A)=COS(A*.017453293)
<br />C Sum the reflectivity occurances
<br />c DO 40 N=1,4 !! tilt loop
<br />c DO 35 L=3,230 !! range loop
<br />c Z=BEAMHT(FLOAT(L),ELT(N))
<br />c IZ=1+NINT(20.*Z)
<br />c IF(IZ.LE.101)THEN ! 5.0 km limit for elevation consideration
<br />c DO 30 M=l,JAT(N) !! azimuth loop
<br />C Convert biased reflectivity to dBZ
<br />C DBZ = (NDBZ -2)/2-32 was expectation, but Ra changed it:
<br />C ID=1+(NDN5(M,L+NDN(N))-2)/2-32
<br />C then dBz = (LV2128-8448)/256 elsewhere
<br />c IO=1+(NON5(M,L+NON(N))-8448)/256
<br />c IF(ID.LT.1)IO=1
<br />c IF(ID.GT.41)ID=41
<br />c NSUM(IZ,ID)=NSUM(IZ,ID)+l
<br />
<br />21
<br />
|