PRO generate_products, varID, GIANTtile, date, time, xdim, ydim, year, month, day, workdir

; 1. set up defaults
; 2. compute contours of interest
; 3. write out GIS shape files, point and polygons

;varID :  one of CREF, ETP18, MEHS, VIL, 1HGC, POSH, SHI

; restore a few things from save file

CASE varID OF  
   'CREF': BEGIN
             var      = 'compref_mosaic'
             thresh   = 50.0 ; dBZ
             clevs=[50.0,55.0,60.0,65.0,70.0,75.0]
             outfile  = 'RadarMax'
             varunit = "dBZ"
             description = 'Daily maximum reflectivity > ' + STRTRIM(thresh,2) + ' ' + varunit
             varname="RadarMax"
           END
   'MREF3D31L': BEGIN
             var      = 'mrefl_mosaic'
             thresh   = 40.0 ; dBZ
             clevs=[40.0,45.0,50.0,55.0,60.0,65.0,70.0,75.0]
             outfile  = 'RadarLow'
             varunit = "dBZ"
             description = 'Daily maximum 2500 m AGL reflectivity > ' + STRTRIM(thresh,2) + ' ' + varunit
             varname="RadarLow"
           END
    'ETP18': BEGIN
             var      = 'etp' 
             thresh   = 12.0 ; km
             outfile  = 'EchoTop'
             clevs=[12.0,16.0]
             outfile  = 'EchoTop'
             varunit = 'km'
             description = 'Daily maximum echo top > ' + STRTRIM(thresh,2) + ' (' + varunit + ')'
             varname='EchoTop'
           END 
   'MEHS': BEGIN
             var      = 'mehs'
             thresh   = 1 ; inches
             outfile  = 'HailSize'
             clevs=[1.0,2.0,3.0,4.0,5.0,6.0]
             varunit = 'inches' 
             description = 'Daily max estimated hail size > ' + STRTRIM(thresh,2) + ' ' + varunit
             varname="HailSize"
           END
'LIQPCP': BEGIN
             var      = 'liqpcp' 
             outfile  = 'LiqPcp'
             clevs=[60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Liquid Precipitation (' + varunit + ')'
             varname="LiqPcp"
           END 
'FRZPCP': BEGIN
             var      = 'frzpcp' 
             outfile  = 'FrzPcp'
             clevs=[5.0,60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Frozen Precipitation (' + varunit + ')'
             varname="FrzPcp"
           END 
'RAINTYPE': BEGIN
             var      = 'raintype' 
             outfile  = 'RainType'
             clevs=[60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Rain Precip Type (' + varunit + ')'
             varname="RainType"
           END 
'SNOWTYPE': BEGIN
             var      = 'snowtype' 
             outfile  = 'SnowType'
             clevs=[5.0,60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Snow Precip Type (' + varunit + ')'
             varname="SnowType"
           END 
'STRATYPE': BEGIN
             var      = 'stratype' 
             outfile  = 'StraType'
             clevs=[5.0,60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Stratiform Precip Type (' + varunit + ')'
             varname="StraType"
           END 
'CONVTYPE': BEGIN
             var      = 'convtype' 
             outfile  = 'ConvType'
             clevs=[5.0,60.0,180.0,360.0,720.0,1440.0]
             varunit='minutes'
             description = 'Convective Precip Type (' + varunit + ')'
             varname="ConvType"
           END 
'HAILDUR': BEGIN
             var      = 'haildur' 
             outfile  = 'HailDur'
             clevs=[5.0,10.0,15.0,20.1,30.0,60.0]
             varunit='minutes'
             description = 'Hail Duration (' + varunit + ')'
             varname="HailDur"
           END 
'PCPFLAG': BEGIN
             var      = 'pcp_flag' 
             outfile  = 'PCPflag'
             clevs=[20.0, 40.0, 60.0, 80.0]
             varunit='percent'
             description = 'Precipitation Flag (' + varunit + ')'
             varname="PCPflag"
           END 
'PCPTYPE': BEGIN
             var      = 'PCPtype'
             outfile  = 'PCPtype'
             clevs=[20.0, 40.0, 60.0, 80.0]
             varunit='catego'
             description = 'Precipitation Type (' + varunit + ')'
             varname="PCPtype"
           END 
   'POSH': BEGIN
             var      = 'posh' 
             outfile  = 'HailProb'
             clevs=[10.0, 30.0, 50.0, 70.0, 90.0]
             varunit='percent'
             description = 'Daily maximum probability of severe hail (' + varunit + ')'
             varname="HailProb"
           END 
    'SHI': BEGIN
             var      = 'shi' 
             outfile  = 'HailIndx'
             clevs=[30.0,60.0,100.0,150.0,210.0]
             varunit='unitless'
             description = 'Daily maximum severe hail index (' + varunit + ')'
             varname="HailIndx"
           END 
    'VIL': BEGIN
             var      = 'vil' 
             thresh   = 10 ; kg/m^2
             outfile  = 'LiqidH2O'
             clevs=[20.0,30.0,40.0,50.0, 60.0]
             varunit = 'kg/m2'
             description = 'Daily maximum vertically integrated liquid (' + varunit + ')'
             varname="LiqidH2O"
           END 
   'VILD': BEGIN
             var      = 'vil' 
             thresh   = 1 ; kg/m^3
             outfile  = 'VILdensi'
             clevs=[1.0,2.0,3.0,4.0,5.0,6.0]
             varunit = 'kg/m3'
             description = 'Daily maximum VIL density (' + varunit + ')'
             varname="VILdensi"
           END 
   '1HGC': BEGIN
             var      = 'q2_gc/q2rad_hsr_gc/1h_acc' 
             thresh   =  2 ; inches of precip
             outfile  = 'Rainfall'
             clevs=[3.0,4.0,6.0,9.0,12.0]
             varunit='inches'
             description = 'Daily accumulated precipitation > ' + STRTRIM(thresh,2) + ' ' + varunit
             varname="Rainfall"
           END
ENDCASE  


=====================================================================
;Set the attribute definitions for the new Shapefile
mynewshape->AddAttribute, STRTRIM(varname,2), 5, 5, PRECISION=1
mynewshape->IDLffShape::AddAttribute, 'dateTime', 7, 16, $
   PRECISION=0
mynewshape->IDLffShape::AddAttribute, 'hours_CST', 5, 5, $
   PRECISION=2
mynewshape->IDLffShape::AddAttribute, 'units', 7, 6, $
   PRECISION=0
mynewshape->IDLffShape::AddAttribute, 'descriptio', 7, 128, $
   PRECISION=0