type spherical_1D_solver!! encapsulate all grid dataprivatereal(r8k),allocatable::v(:,:)!! v(:,1) = r, v(:,2) = T, shape = [nr,2]real(r8k),allocatable::rho(:),cp(:)!! density and specific heat (size = nr)real(r8k),allocatable::T_analytical(:)!! expected solution (size = nr)contains procedure::set_vprocedure::set_material_properties_sizeprocedure::set_expected_solution_sizeprocedure::set_rhoprocedure::set_cpprocedure::time_advance_heat_equationend type spherical_1D_solver