f := zone.density(z)
zone.density(z) := f

Get/set the density of the zone. Note that updates to the gravitational and inertial mass of connected gridpoints may not occur until an update occurs. This normally happens at the start of cycling and at the geometry update interval in large-strain mode. The zone.force.update intrinsic can be used to force an udpate at the start of the next cycle.

Returns:f - the current zone density
Accepts:f - the current zone density
Arguments:z - a zone pointer