Changeset 15518

Timestamp:

2021-11-16T17:57:25+01:00 (2 years ago)

Author:

edblockley

Message:

Correction to rheology convergence timeseries output to ensure that all iterations on last time-step are output and that the time-axis only includes point from the current run cycle (i.e., starting from nit000+1)

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90 (modified) (2 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90

@@ -1006 +1006 @@
 
       ! time
-      it = ( kt - 1 ) * kitermax + kiter
+      it = ( kt - nit000 ) * kitermax + kiter
       
       ! convergence
@@ -1026 +1026 @@
          istatus = NF90_PUT_VAR( ncvgid, nvarid, (/zresm/), (/it/), (/1/) )
          ! close file
-         IF( kt == nitend - nn_fsbc + 1 )   istatus = NF90_CLOSE(ncvgid)
+         IF( kt == nitend - nn_fsbc + 1 .AND. kiter == kitermax )   istatus = NF90_CLOSE(ncvgid)
       ENDIF