Structures of rapidly rotating baroclinic stars – II. An extended numerical method for realistic stellar models