; docformat = 'rst'
;
; NAME: 
;    ssoddCompPoro
; PURPOSE:
;    Compute the porosity of a SSO from its density and analog meteorite
;
;+
; :Description:
;    Compute the porosity of a SSO from its density and analog meteorite
;
; :Categories:
;    Database, SSODD, Asteroid, Density
;    
; :Params:
;
; :Keywords:
;
; :Examples:
;
; :Uses:
;
; :Author:
;    B.Carry (OCA)
;
; :History:
;    Change History::
;      Written in April 2018, B. Carry (OCA)
;-
function ssoddCompPoro, density, meteorite
  COMPILE_OPT hidden, idl2
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---    I -- Initialization And Input Verification              -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;

 ;--I.1-- Check Output Location  ---------------------------------------------------------------;


;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;
;---  TAG ---   II -- Database Step 1: Census of Objects                 -----------------------;
;-----------------------------------------------------------------------------------------------;
;-----------------------------------------------------------------------------------------------;



  poro = 100. *( 1. - density.val/meteorite.grain.val )
  dPoro = sqrt( (meteorite.grain.unc/meteorite.grain.val)^2. + (density.unc/density.val)^2. )
   dPoro *= sqrt(2.) * poro

  macro = 100. *( 1. - density.val/meteorite.bulk.val )                                        
  dMacro = sqrt( (meteorite.bulk.unc/meteorite.bulk.val)^2. + (density.unc/density.val)^2. )
    dMacro *= sqrt(2.) * macro


  return, dens

end