c      cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc          qbcalc routine                                 ccc
c      cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       subroutine qbcalc( rho , tem , tf , fityp, qb , gamma , nn )
       implicit real*8 (a-h,o-z)
       integer fityp
       dimension qb(nn) , tem(nn), rho(nn), tf(nn) , fityp(nn)
     &          ,gamma(nn)
       common /nika2/ a , z , ic 
c
       do 100 i = 1, ic
          if      ( fityp( i ) .eq. 1 ) then 
              call fitgas( a , z , tem( i ) , rho( i ) , qb( i ) )
          else if ( fityp( i ) .eq. 2 ) then
              call fitliq( a , z 
     &               ,gamma( i ) , tem( i ) , rho( i ) , qb( i ) )
          else if ( fityp( i ) .eq. 3 ) then
              call fitcry( a , z 
     &               ,gamma( i ) , tem( i ) , rho( i ) , qb( i ) )
          else
              call errcall(1)
          endif
 100   continue
        return
        end