fix: replace function now also replaces unit parameter if applicable

Author: diogomsmiranda

pandas/_libs/tslibs/timestamps.pyx

@@ -59,6 +59,7 @@ from pandas._libs.tslibs.conversion cimport (
     maybe_localize_tso,
 )
 from pandas._libs.tslibs.dtypes cimport (
+    abbrev_to_npy_unit,
     npy_unit_to_abbrev,
     npy_unit_to_attrname,
     periods_per_day,
@@ -2439,10 +2440,12 @@ default 'raise'
             datetime ts_input
             tzinfo_type tzobj
             _TSObject ts
+            NPY_DATETIMEUNIT rep_reso
 
         # set to naive if needed
         tzobj = self.tzinfo
         value = self._value
+        rep_reso = self._creso
 
         # GH 37610. Preserve fold when replacing.
         if fold is None:
@@ -2465,41 +2468,52 @@ default 'raise'
             return v
 
         if year is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_Y
             dts.year = validate("year", year)
         if month is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_M
             dts.month = validate("month", month)
         if day is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_D
             dts.day = validate("day", day)
         if hour is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_h
             dts.hour = validate("hour", hour)
         if minute is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_m
             dts.min = validate("minute", minute)
         if second is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_s
             dts.sec = validate("second", second)
         if microsecond is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_ms
             dts.us = validate("microsecond", microsecond)
         if nanosecond is not None:
+            rep_reso = NPY_DATETIMEUNIT.NPY_FR_ns
             dts.ps = validate("nanosecond", nanosecond) * 1000
         if tzinfo is not object:
             tzobj = tzinfo
+        
+        if rep_reso < self._creso:
+            rep_reso = self._creso
 
         # reconstruct & check bounds
         if tzobj is None:
             # We can avoid going through pydatetime paths, which is robust
             #  to datetimes outside of pydatetime range.
             ts = _TSObject()
             try:
-                ts.value = npy_datetimestruct_to_datetime(self._creso, &dts)
+                ts.value = npy_datetimestruct_to_datetime(rep_reso, &dts)
             except OverflowError as err:
                 fmt = dts_to_iso_string(&dts)
                 raise OutOfBoundsDatetime(
                     f"Out of bounds timestamp: {fmt} with frequency '{self.unit}'"
                 ) from err
             ts.dts = dts
-            ts.creso = self._creso
+            ts.creso = rep_reso
             ts.fold = fold
             return create_timestamp_from_ts(
-                ts.value, dts, tzobj, fold, reso=self._creso
+                ts.value, dts, tzobj, fold, reso=rep_reso
             )
 
         elif tzobj is not None and treat_tz_as_pytz(tzobj):
@@ -2518,10 +2532,10 @@ default 'raise'
             ts_input = datetime(**kwargs)
 
         ts = convert_datetime_to_tsobject(
-            ts_input, tzobj, nanos=dts.ps // 1000, reso=self._creso
+            ts_input, tzobj, nanos=dts.ps // 1000, reso=rep_reso
         )
         return create_timestamp_from_ts(
-            ts.value, dts, tzobj, fold, reso=self._creso
+            ts.value, dts, tzobj, fold, reso=rep_reso
         )
 
     def to_julian_date(self) -> np.float64: