new surface-data output; renamed output variables (pt to theta, rho_air to rho, rho_ocean to rho_sea_water)