Changeset 8377 for branches/UKMO

Timestamp:

2017-07-26T15:48:43+02:00 (7 years ago)

Author:

frrh

Message:

Prevent attempt to compute meaningless mean values in trcrst when appropriate
time level data is not available. Output "N/A" instead.
I'm not convinced this is a complete solution. fse3t_a may be undefined
when it is used following a restart in which case there's more to do.

File:

• branches/UKMO/dev_r5518_fix_NaNs/NEMOGCM/NEMO/TOP_SRC/trcrst.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_fix_NaNs/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

@@ -745 +745 @@
       !!---------------------------------------------------------------------
       INTEGER  :: jk, jn
+      CHARACTER (LEN=18) :: text_zmean
       REAL(wp) :: ztraf, zmin, zmax, zmean, areasf
       REAL(wp), DIMENSION(jpi,jpj) :: zvol
@@ -761 +762 @@
          CALL mpp_max( zmax )      ! max over the global domain
       END IF
-      zmean  = ztraf / areasf
-      IF(lwp) WRITE(numout,9002) TRIM( names ), zmean, zmin, zmax
-      !
-      IF(lwp) WRITE(numout,*)
-9002  FORMAT(' tracer name :',a10,'    mean :',e18.10,'    min :',e18.10, &
+
+      text_zmean = "N/A"
+      ! Avoid divide by zero when this routine is called from parts of the code
+      ! which may not have necessary time level variables to calculate areasf.
+      IF  (areasf > 0.0) THEN 
+         zmean = ztraf / areasf
+         WRITE(text_zmean,'(e18.10)') zmean
+      ENDIF
+
+      IF(lwp) WRITE(numout,9002) TRIM( names ), text_zmean, zmin, zmax
+
+  9002  FORMAT(' tracer name :',A,'    mean :',A,'    min :',e18.10, &
       &      '    max :',e18.10 )
       !