private
!
  public :: initialize_density_methods, getrho, getrho1, getlnrho, getrho_s
  public :: getdlnrho_x, getdlnrho_y, getdlnrho_z
  public :: putrho, putlnrho
  public :: getderlnrho_z