Update the function to read the netCDF variables to do so slightly more intelligently wherein the last dimension is assumed to be time rather than hard coding the 4th dimension as time (which is not always the case). Also add support for an alternative time variable (t_1); eventually this will probably have to be extended to more time variables, but for me and for the time being, this seems to work OK.