gsi=(Na+Nb)/(na+nb) for open shell implemented