How can I reproject a NetCDF file from its native curvilinear grid to a regular lat-lon grid while preserving all time steps and data variables using Python? Could you provide a Python code example or suggest libraries and methods for achieving this interpolation while maintaining data accuracy?