The density is added to the node_dof records. Now this density is used in the equations, and the data item record materi_density in the data part needs not to be specified. The density is governed by the mass conservation law
where is the velocity in the direction of the space coordinate . Don't forget to initialize the density by specifying node_dof records. Condition: the -materi differential group_type should be present.