; Computes Volumenanteile :
feis = 0
fvac = [0.01, 0.1, 0.2,0.3,0.4,0.5]

; rhosi=3.5, rhoc=1.6:
; c = fsi/fac = rhoc/rhsi * m(si)/m(c)

msi_c = [0.7, 1.3, 2.3]   ; Fig. S17


;Mass ratio aC/Si of large grains:
  Simass/aCmass = 565./176. = 3.21 for single cloud in ISM

;Volumne ratio aC/Si of lkarge grains:  
;  rhc/rhSi *Simass/aCmass = 1.6/3.5 * 566./176. = 1.47 !
  
msi_c = 3.21 ; Tab.2 mean of single cloud ISM in S17 (LIPS.pdf)
c     = msi_c *1.6/3.5

fvac = 0.2
fsi  = 0.48
faC  = 0.32



print, 'fsi, faC,  fvac'
for j = 0, n_elements(fvac)-1 do begin
 for i = 0, 2 do begin
 fsi = (1-fvac(j))/ (1+1./c(i))
 fac = (1-fvac(j)) - (1-fvac(j))/ (1+1./c(i))
 if(abs(fvac(j)+ fsi + faC-1.) gt 1e-3 ) then stop
 print, format='(3f7.2)', fsi, faC,  fvac(j)
endfor
endfor

end