part way through get_model_pressure_profile

Needs logic to deal with obs below lowest level
Needs static data to calculate rho
Using two different values for gravity see #404 (comment)
I'm not sure about the check for 0.0 for surface pressure. Does this happen?
For model_rho_t checking speficically for MU variable name rather than qty_pressure