Changeset 12350

Timestamp:

2020-01-29T12:08:08+01:00 (3 years ago)

Author:

smueller

Message:

Rectification of a defect that halved the maximum length of the interval between restarts in model runs with enabled temporal interpolation of forcing fields (field-information flag ln_tint set to .true.)

File:

• NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/fldread.F90 (modified) (1 diff)

NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/fldread.F90

@@ -906 +906 @@
       ENDIF
       !
-      IF( sdjf%ln_tint ) THEN   ! record time defined in the middle of the record
-         sdjf%nrecsec(1:sdjf%nreclast) = 0.5 * ( sdjf%nrecsec(0:sdjf%nreclast-1) + sdjf%nrecsec(1:sdjf%nreclast) )
+      IF( sdjf%ln_tint ) THEN   ! record time defined in the middle of the record, computed using an implementation
+                                ! of the rounded average that is valid over the full integer range
+         sdjf%nrecsec(1:sdjf%nreclast) = sdjf%nrecsec(0:sdjf%nreclast-1) / 2 + sdjf%nrecsec(1:sdjf%nreclast) / 2 + &
+            & MAX( MOD( sdjf%nrecsec(0:sdjf%nreclast-1), 2 ), MOD( sdjf%nrecsec(1:sdjf%nreclast), 2 ) )
       END IF
       !