def METHOD_NAME(self):
    self.net = op.network.Cubic(shape=[3, 3, 3])
    self.phase = op.phase.Species(network=self.net, species='h2o')
    self.phase['pore.temperature'] = 300*np.ones(self.phase.Np,)
    self.phase['pore.salinity'] = np.zeros((self.phase.Np,))