obj = THG(caNl, ...)
  caNl is the n_l in [m2/W^l]
 
  options ae same as in SPM