function get_chandra_rainrate,dbz,zdr,kdp,rhv,hid,method,MISSING=missing
;
; *** Function that accepts a dbz, zdr and kdp value and then
; *** returns a rain rate value and method.  This is the algorithm
; *** provided by Chandra to Dave Wolff.
;
    if(NOT KEYWORD_SET(MISSING)) then missing=-32767.0
    if(dbz eq missing) then begin
        r = missing
	method = 0
        return, r
    endif
    
    dbz_lin = 10^(0.1D * dbz)
    zdr_lin = 10^(0.1D * zdr)
;
; *** Yes -- Proceed with algorithm calculations
;        
    if(kdp gt 0.3 and dbz gt 38) then begin
        if(zdr gt 0.5) then begin
            r = 90.8 * kdp^(0.93) * zdr^(-1.69)
            method = 1          ; Kdp/Zdr
        endif else begin
            r = 50.7 * kdp^(0.85)
            method = 2          ; Kdp only
        endelse
    endif else begin
        if(zdr gt 0.5) then begin
            r = 0.0067 * dbz_lin^(0.927) * zdr^(-3.433)
            method = 3          ; Zh(linear)/Zdr
        endif else begin
            r = 0.017 * dbz_lin^(0.7143)
            method = 4          ; Zh linear only
        endelse
    endelse

    return, r
    end