checks if the density field is a valid adiabatic re-arrangement of the base state
flag = isDensityInValidRange()
flag
This is probably best re-defined as a dynamical variable.