# Presentation: 826862

## Design Approach ReviewMDA Phase III

- Brian Klein, RSIS
- 12 October 2004

## Contents

- CCRs.....................................................................3
- Concept of Operations.................................................4
- Impacts to Top Level Architecture.......................................7
- Impacts to WSR-88D Architecture................................8
- MDA Rapid Update Data Flow..................................................9
- Design Approach...............................................................11
- Products.........................................................................15
- Development Approach.............................................26
- Test Approach.........................................................27
- Integration Challenges/Risks......................................28
- Proposed Schedule..................................................29
- Action Items............................................................30

## CCRs

- Four New CCRs
- NA04-28004 MDA Phase III, DEA Configuration Data and DMD Header Completion.
- NA04-28002: Elevated Circulation Cores Not Always Detected by MDA.
- CCR Pending: Additional 2D Component Attributes for DMD product
- NA04-28003: MD Product Incorrectly Applies Stationary Feature Logic.

## Current Build 7 MDA:
Three DEA parameters (URC LOCA)
2D component parameters in DMD product:
Elevation Time
Height
Diameter
Rotational Velocity
Shear
Gate-to-Gate Velocity Difference
Strength Rank
DMD attributes redundant with PDB stubbed out. 
NSSL version of core search.
Stationary feature logic incorrect in MD product.

- Three DEA parameters (URC LOCA)
- 2D component parameters in DMD product:
- Elevation Time
- Height
- Diameter
- Rotational Velocity
- Shear
- Gate-to-Gate Velocity Difference
- Strength Rank
- DMD attributes redundant with PDB stubbed out.
- NSSL version of core search.
- Stationary feature logic incorrect in MD product.

**Concept of Operations****(after IEEE Std 1362-1998)**

## Concept of Operations(continued)

- Justification for Changes:
- Transfer 32 hard-coded parameters to DEA configuration data to support parameter optimization and expected modifications in anticipation of base data resolution changes with ORDA
- Addition of point location and elevation angles for 2D component centroids will facilitate building of 3D feature graphics with correctly placed 2D components.
- Some “generic” product parameters are incomplete.
- SEC testing after Build 6 delivery showed surprisingly low 3D feature strength ranks in rare cases.
- Stationary feature logic does not meet original requirement.

## Concept of Operations(concluded)

- Concepts for Build 8 MDA:
- Adjustments to the MDA required by increased resolution base data from ORDA will be possible by editing DEA data and without affecting source code.
- MDL plans to use elevation angle information for an AWIPS OB6 meso table display.
- Third-party graphics programs can construct 3D images of mesocyclones with properly positioned 2D components.
- DMD users will be able to get all needed product information from the generic product packet alone without also accessing the PDB.

## Impacts to Top-Level Architecture

- External Display Systems
- AWIPS software modifications are required to support new 2D component location and elevation angle data in the DMD product.
- Benefits
- Correct strength rank values for 3D features.
- Data available for accurate 3D graphics display.
- Correct stationary feature information.
- Complete DMD product header.

## Impacts to WSR-88D Architecture

- Product Size
- Uncompressed DMD product using F-LOAD data will increase by about 10 Kbytes to 95 -100 Kbytes total.
- MD product size unaffected.
- New API
- Update to rpgcs_latlon functions.
- Completion of generic product header similar to RPGC_prod_desc_block function.

## Slide 9

## Slide 10

## Design Approach

- Change strength rank search logic to find “hidden” elevated cores.
- New “global” fields to DMD product.
- Elevation angles
- Elevation times (was in 2D component data)
- New 2D component fields in DMD product.
- Arrays for 2D component locations
- X/Y or azimuth/range or latitude/longitude
- Array of elevation indices
- Change stationary feature thresholds in MD product to 2.5 m/sec.
- Add 32 new DEA parameters that relate to spatial resolution

## Design Approach MDA AEL strength rank logic change

- _COMPUTE_ 3D Feature(Core Base)(34), 3D Feature(Core Top)(35), 3D Feature(Core Depth)(36)
- _ASSIGN_ CORE BASE = Missing
- _DO FOR ALL_ 2D Features in a 3D Feature (i.e. from bottom to top)
- 	_**IF**_** 2D Feature(Strength Rank) >= CORE RANK **_**THEN**_
**	****     **_**IF**_** Next Higher 2D Feature(Strength Rank) >= CORE RANK **_**OR**_** **
**                ****BEAM WIDTH >= THRESHOLD(Meso Min Depth) **_**THEN**_
- 	          CORE BASE = 2D Feature(Height)
- 	          _BREAK_
- 	     _**END IF**_
**	**_**END IF**_
- _END DO_
- _IF_ CORE BASE != Missing _THEN_
- 	_DO FOR ALL_ 2D Features in a 3D Feature _FROM_ Second _TO_ Last (i.e. from bottom to top)
- 	          _IF_ 2D Feature(Strength Rank) < CORE RANK _THEN_
- 	               CORE TOP = previous, lower 2D Feature(Height)
- 	               _BREAK_
- 	          _ELSE IF_ 2D Feature is the last _AND_
- 	                        2D Feature 2D Feature(Strength Rank) >= CORE RANK _THEN_
- 	               CORE TOP = current 2D Feature(Height)
- 	          _END IF_
- 	END DO
- END IF

## Design Approach MDA AEL strength rank logic change

## Design Approach

**New DEA Configuration Data**

| RADAR(Beam Width) | THRESHOLD (ending Range Linear Reduction) | THRESHOLD (maximum 2D Construction Radial | THRESHOLD (3D Feature Low Level Convergence Height) |
| --- | --- | --- | --- |
| THRESHOLD (maximum Shear Segment Height) | THRESHOLD (maximum Shear Segment Length) | THRESHOLD (2D Feature minimum Range) | THRESHOLD (3D Feature Mid Level Convergence Height1) |
| THRESHOLD (maximum Shear Segment Distance) | THRESHOLD (maximum Core Shear Segment Length) | THRESHOLD (minimum 2D Association Distance) | THRESHOLD (3D Feature Mid Level Convergence Height2) |
| THRESHOLD (maximum # 1D Features) | THRESHOLD (maximum # Shear Segments) | THRESHOLD (maximum 2D Association Distance) | THRESHOLD (max number of elevations in 3D Feature) |
| THRESHOLD (maximum # 2D Features) | THRESHOLD (maximum Aspect Ratio) | THRESHOLD (Convergence maximum number of Vectors) | THRESHOLD(Shallow Allowed) |
| THRESHOLD (maximum # 3D Features) | THRESHOLD (minimum Aspect Ratio) | THRESHOLD (Convergence Look Ahead) | THRESHOLD (Meso min Rank Shallow) |
| THRESHOLD (minimum # Shear Segments) | THRESHOLD (minimum 2D Radial Diameter) | THRESHOLD (Convergence minimum Length) | THRESHOLD (Meso min Depth Shallow) |
| THRESHOLD (beginning Range Linear Reduction) | THRESHOLD (maximum 2D Diameter) | THRESHOLD (Convergence Buffer Zone) | THRESHOLD (Meso max Top Shallow) |

## Products MD PSB & GAB

**Mesocyclone Detection**
**    ****(MD, 141)**

## Products (cont’d)MD TAB

## Products (cont’d)DMD

- product id 149, name: Digital Mesocyclone Detection, desc: Data array product output of Mesocyclone Detection Algorithm.
-     type 1, time 0
-     Radar name: ***KLWX***, lat 37.9753, lon -77.4778, height   ***354***
-     Vol time ***1692***, elev time ***1796***, vol number ***80***, elev number 7
-     VCP 12, Op mode ***0***, elev 0.5000
-     compress type 0, size 78750

## Products (cont’d)DMD

- # ***5*** params:
- param[0]: id: avg_dir, attrs: Name=Average Direction of Tracked Features; Type=float; Value=249.0; Unit=deg;
- param[1]: id: avg_spd, attrs: Name=Average Speed of Tracked Features; Type=float; Value= 17.4; Unit=m/s;
- param[2]: id: last_elev_flag, attrs: Name=Last Elevation Flag; Type=int; Value=1;
***param[3]: id: elev_angle, attrs: Name=Elevation Angle; Type=float; Value=0.5,1.5,2.4,3.4,4.5; Unit=deg;***
***param[4]: id: elev_time, attrs: Name=Elevation Time; Type=int; Value=1692,1730,1750,1774,1796; Unit=s;***
- number of components 14

## Products (cont’d)DMD

- AREA: (3), area type 1
- # 4***3*** params:
- param[0]: id: meso_event_id, attrs: Name=Meso ID; Type=int; Value=82;
- param[1]: id: storm_id, attrs: Name=Associated Storm ID; Type=string; Value=A5;
- param[2]: id: age, attrs: Name=Age; Type=int; Value=3058; Unit=s;
- param[3]: id: strength_rank_type, attrs: Name=Strength Rank Type; Type=string; Value=L;
- param[4]: id: detection_status, attrs: Name=Detection Status; Type=string; Value=   ;
- param[5]: id: num2D, attrs: Name=Number of 2D Features in this 3D Feature; Type=int; Value=6;

## Products (cont’d)DMD

- param[6]: id: tvs_near, attrs: Name=Associated TVS; Type=string; Value=Y;
- param[7]: id: prop_spd, attrs: Name=Speed; Type=float; Value= 14.3; Unit=m/s;
- param[8]: id: prop_dir, attrs: Name=Direction; Type=float; Value=251.0; Unit=deg;
- param[9]: id: ll_azm, attrs: Name=Base Azimuth; Type=float; Value=152.8; Unit=deg;
- param[10]: id: ll_rng, attrs: Name=Base Range; Type=float; Value= 91.6; Unit=km;
- param[11]: id: lowest_elev, attrs: Name=Base on Lowest Elevation; Type=string; Value=Y;
- param[12]: id: base, attrs: Name=Base Height; Type=float; Value= 1.3; Unit=km;

## Products (cont’d)DMD

- param[13]: id: depth, attrs: Name=Depth; Type=float; Value=11.4; Unit=km;
- param[14]: id: num_past_pos, attrs: Name=Number of Past Positions; Type=int; Value=10;
- param[15]: id: past_lat, attrs:  Name=Past Latitude; Type=float; Value= 38.2449, 38.2186, 38.1719, 38.1716, 38.1496, 38.1197, 38.0872, 38.0616, 38.0816, 38.0797; Unit=deg;
- param[16]: id: past_lon, attrs:  Name=Past Longitude; Type=float; Value= -77.0285, -77.1031, -77.1618, -77.1979, -77.2563, -77.3720, -77.3672, -77.4831, -77.5406, -77.5807; Unit=deg;
- param[17]: id: num_fcst_pos, attrs: Name=Number of Forecast Positions; Type=int; Value=6;

## Products (cont’d)DMD

- param[18]: id: fcst_lat, attrs:  Name=Forecast Latitude; Type=float; Value= 38.2546, 38.2669, 38.2791, 38.2914, 38.3036, 38.3159; Unit=deg;
- param[19]: id: fcst_lon, attrs:  Name=Forecast Longitude; Type=float; Value= -76.9532, -76.9069, -76.8605, -76.8142, -76.7678, -76.7215; Unit=deg;
- param[20]: id: strength_rank, attrs: Name=Strength Rank; Type=int; Value=7;
- param[21]: id: ll_diam, attrs: Name=Base Diameter; Type=float; Value= 9.5; Unit=km;
- param[22]: id: ll_rot_vel, attrs: Name=Base Rotational Velocity; Type=float; Value= 29.3; Unit=m/s;
- param[22]: id: ll_rot_vel, attrs: Name=Base Rotational Velocity; Type=float; Value= 29.3; Unit=m/s;

## Products (cont’d)DMD

- param[23]: id: max_rot_vel, attrs: Name=Max Rotational Velocity; Type=float; Value= 29.3; Unit=m/s;
- param[24]: id: height_max_rot_vel, attrs: Name=Height of Max Rotational Velocity; Type=float; Value= 1.3; Unit=km;
- param[25]: id: ll_shear, attrs: Name=Base Shear; Type=float; Value= 6.2; Unit=m/s/km;
- param[26]: id: max_shear, attrs: Name=Max Shear; Type=float; Value= 9.7; Unit=m/s/km;
- param[27]: id: height_max_shear, attrs: Name=Height of Max Shear; Type=float; Value= 2.9; Unit=km;
- param[28]: id: ll_gtg_vel_dif, attrs: Name=Base Gate-to-Gate Velocity Difference; Type=float; Value= 34.7; Unit=m/s;
- param[29]: id: display_filter, attrs: Name=Overlaps Lower Feature; Type=string; Value=N;

## Products (cont’d)DMD

- param[30]: id: msi, attrs: Name=MSI; Type=int; Value=4834;
- param[31]: id: storm_rel_depth, attrs: Name=Storm Relative Depth; Type=int; Value=70; Unit=%;
- param[32]: id: ll_convergence, attrs: Name=0-2 km ARL Convergence; Type=float; Value= 27.5; Unit=m/s;
-  param[33]: id: ml_convergence, attrs: Name=2-4 km ARL Convergence; Type=float; Value= 25.2; Unit=m/s;
***param[34]: id: tilt_num, attrs: Name=Elevation Index; Type=int; Value=2,3,4,5,6,7;***
- param[35]: id: height, attrs:  Name=2D Height; Type=float; Value= 1.3, 2.9, 4.3, 6.0, 7.4, 9.2; Unit=km;
- param[36]: id: diam, attrs: Name=2D Diameter; Type=int; Value=9,5,9,10,14,7; Unit=km;

## Products (cont’d)DMD

- param[37]: id: rot_vel, attrs: Name=2D Rotational Velocity; Type=int; Value=29,25,22,23,21,19; Unit=m/s;
- param[38]: id: shear, attrs: Name=2D Shear; Type=int; Value=6,10,5,5,3,6; Unit=m/s/km;
- param[39]: id: gtgmax, attrs: Name=2D Gate-to-Gate Velocity Difference; Type=int; Value=35,26,21,23,18,15; Unit=m/s;
***param[40]: id: 2d_lat, attrs: Name=2D Latitude; Type=float; Value= ******38.2546, 38.2669, 38.2791, 38.2914, 38.3036, 38.3159******; Unit=deg;***
***param[41]: id: 2d_lon, attrs: Name=2D Longitude; Type=float; Value=***
***   ******-76.9532, -76.9069, -76.8605, -76.8142, -76.7678, -76.7215******; Unit=km;***
- param[42]: id: rank, attrs: Name=2D Strength Rank; Type=int; Value=10,7,6,5,4,5;
- $ 1 points:
- point[0]: lat =    38.25, lon =   -76.95

## Development Approach

- ANSI C conforming to RPG Coding Guidelines.
- Develop and test with Build 6v3 on dev3 Blade with connectivity to AWIPS lite.
- Coordinate DMD changes with FSL and MDL.

## Test Approach

- MD Product
- Inspection of cvg display (PSB, GAB and TAB).
- DMD Product
- cvt printouts.
- AWIPS lite displays when ready.
- Prototypes to MDL.
- Test Cases
- KCCX, 14 July 2004
- Full Load, 21 April 1996
- Still looking for stationary feature case.
- Verify results based on pre-DEA version.

## Integration Challenges/Risks

- Strength rank logic and stationary feature cases are rare.
- Uncompressed DMD product is large--up to 100 KBytes in full load--but compresses to < 9 KBytes.
- Expect increase of up to 10 Kbytes (uncompressed) for EOV DMD product.

## Proposed Schedule Summary(as of 09/17/04)

| Task Name | Duration | Start | Finish |
| --- | --- | --- | --- |
| Modify Strength Rank Logic | 43 days | 8/2/04 | 9/30/04 |
| Add DEA Data and DMD Header | 54 days | 7/22/04 | 10/6/04 |
| Add 2D Locations/Angles in DMD | 28 days | 9/8/04 | 10/18/04 |
| Correct Stationary Feature Logic | 27 days | 8/19/04 | 9/27/04 |
| Design Approach Review | 8 days | 9/13/04 | 9/22/04 |
| Integration Readiness Review | 5 days | 10/26/04 | 11/1/04 |
| Deliver Software and Documentation | 1 day | 11/2/04 | 11/2/04 |

## Action Items

- 	_______________________________
- 	_______________________________
- 	_______________________________
- 	_______________________________
- 	_______________________________