use yaeos, only: Groups, setup_unifac, UNIFAC
 use yaeos__consistency_gemodel, only: individual_ge_calls

 type(UNIFAC) :: model

 integer, parameter :: nc = 4

 type(Groups) :: molecules(nc)

 real(pr) :: n(nc), T

 T = 303.15
 n = [400.0, 100.0, 300.0, 200.0] ! always test with sum(n) > 1

 ! Hexane [CH3, CH2]
 molecules(1)%groups_ids = [1, 2]
 molecules(1)%number_of_groups = [2, 4]

 ! Ethanol [CH3, CH2, OH]
 molecules(2)%groups_ids = [1, 2, 14]
 molecules(2)%number_of_groups = [1, 1, 1]

 ! Toluene [ACH, ACCH3]
 molecules(3)%groups_ids = [9, 11]
 molecules(3)%number_of_groups = [5, 1]

 ! Cyclohexane [CH2]
 molecules(4)%groups_ids = [2]
 molecules(4)%number_of_groups = [6]

 model = setup_unifac(molecules)

 call individual_ge_calls(model, n, T, passed)

 ! If TRUE, then all individual calls are correct
 print *, passed